QM Subsystem Research Articles

Assessing the Partial Hessian Approximation in QM/MM-Based Vibrational Analysis.

The partial Hessian approximation is often used in vibrational analysis of quantum mechanics/molecular mechanics (QM/MM) systems because calculating the full Hessian matrix is computationally impractical. This approach aligns with the core concept of QM/MM, which focuses on the QM subsystem. Thus, using the partial Hessian approximation implies that the main interest is in the local vibrational modes of the QM subsystem. Here, we investigate the accuracy and applicability of the partial Hessian vibrational analysis (PHVA) approach as it is typically used within QM/MM, i.e., only the Hessian belonging to the QM subsystem is computed. We focus on solute-solvent systems with small, rigid solutes. To separate two of the major sources of errors, we perform two separate analyses. First, we study the effects of the partial Hessian approximation on local normal modes, harmonic frequencies, and harmonic IR and Raman intensities by comparing them to those obtained using full Hessians, where both partial and full Hessians are calculated at the QM level. Then, we quantify the errors introduced by QM/MM used with the PHVA by comparing normal modes, frequencies, and intensities obtained using partial Hessians calculated using a QM/MM-type embedding approach to those obtained using partial Hessians calculated at the QM level. Another aspect of the PHVA is the appearance of normal modes resembling the translation and rotation of the QM subsystem. These pseudotranslational and pseudorotational modes should be removed as they are collective vibrations of the atoms in the QM subsystem relative to a frozen MM subsystem and, thus, not well-described. We show that projecting out translation and rotation, usually done for systems in isolation, can adversely affect other normal modes. Instead, the pseudotranslational and pseudorotational modes can be identified and removed.

QM/CG-MM: Systematic Embedding of Quantum Mechanical Systems in a Coarse-Grained Environment with Accurate Electrostatics.

Quantum Mechanics/Molecular Mechanics (QM/MM) can describe chemical reactions in molecular dynamics (MD) simulations at a much lower cost than ab initio MD. Still, it is prohibitively expensive for many systems of interest because such systems usually require long simulations for sufficient statistical sampling. Additional MM degrees of freedom are often slow and numerous but secondary in interest. Coarse-graining (CG) is well-known to be able to speed up sampling through both reduction in simulation cost and the ability to accelerate the dynamics. Therefore, embedding a QM system in a CG environment can be a promising way of expediting sampling without compromising the information about the QM subsystem. Sinitskiy and Voth first proposed the theory of Quantum Mechanics/Coarse-grained Molecular Mechanics (QM/CG-MM) with a bottom-up CG mapping. Mironenko and Voth subsequently introduced the DFT-QM/CG-MM formalism to couple a Density Functional Theory (DFT) treated QM system and to an apolar environment. Here, we present a more complete theory that addresses MM environments with significant polarity by explicitly accounting for the electrostatic coupling. We demonstrate our QM/CG-MM method with a chloride-methyl chloride SN2 reaction system in acetone, which is sensitive to solvent polarity. The method accurately recapitulates the potential of mean force for the substitution reaction, and the reaction barrier from the best model agrees with the atomistic simulations within sampling error. These models also have generalizability. In two other reactive systems that they have not been trained on, the QM/CG-MM model still achieves the same level of agreement with the atomistic QM/MM models. Finally, we show that in these examples the speed-up in the sampling is proportional to the acceleration of the rotational dynamics of the solvent in the CG system.

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. The adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems. Program SummaryProgram Title: QMMM 2023.CPC Library link to program files:https://doi.org/10.17632/zmkc9db4zj.1.Licensing provisions: Apache-2.0 for the program and CC-BY-4.0 for the manual.Programming languages: Fortran90, Perl, and shell scripting.Nature of problem: Multiscale modeling and simulations by the combined quantum mechanical and molecular mechanical (QM/MM) method.Solution method: The program performs single-point calculations, geometry optimizations, transition state optimizations, and molecular dynamics (MD) by the QM/MM method by interfacing calls to binary executables in a QM package (GAMESS-US, Gaussian, or ORCA) and an MM package (TINKER) for the required calculations. QMMM 2023 supports various types of embedding schemes and multiple ways of treating a QM-MM boundary passing through a covalent bond.Additional comments including restrictions and unusual features: (1) QMMM 2023 supports polarizable-embedding and flexible-embedding schemes based on the electronegativity equalization principle and allowing charge transfer between the QM and MM regions. (2) For a QM-MM boundary passing through covalent bonds, all link-atom-based treatments preserve the charge of the MM region (and MM bond dipoles if desired) while effectively avoiding the over-polarization of the QM region. Further tuning is available to account for the charge penetration effects at the QM-MM boundary. (3) For dynamics simulations, Hamiltonian-based adaptive-partitioning algorithms allow on-the-fly reclassification of atoms or groups into the QM or MM regions without abrupt changes in the energy of the entire system or in the forces on all atoms.

Challenges in Protein QM/MM Simulations with Intra-Backbone Link Atoms.

Hybrid quantum mechanical/molecular mechanical (QM/MM) simulations fuel discoveries in many fields of science including computational biochemistry and enzymology. Development of more convenient tools leads to an increase in the number of works in which mechanical insights into enzymes' mode of operation are obtained. Most commonly, these tools feature hydrogen-capping (link atom) approach to provide coupling between QM and MM subsystems across a covalent bond. Extensive studies were conducted to provide a solid foundation for the correctness of such an approach when a bond to a nonpolar MM atom is considered. However, not every task may be accomplished this way. Certain scenarios of using QM/MM in computational enzymology encourage or even necessitate the incorporation of backbone atoms into the QM region. Two out of three backbone atoms are polar, and in QM/MM with electrostatic embedding, a neighboring link atom will be hyperpolarized. Several schemes to mitigate this effect were previously proposed alongside a rigorous assessment of quantitative effects on model systems. However, it was not clear whether they may translate into qualitatively different results and how link atom hyperpolarization may manifest itself in a real-life enzymological scenario. Here, we show that the consequences of such an artifact may be severe and may completely overturn the conclusions drawn from the simulations. Our case advocates for the use of charge redistribution schemes whenever intra-backbone QM/MM boundaries are considered. Moreover, we addressed how different boundary types and charge redistribution schemes influence backbone dynamics. We showed that the results are heavily dependent on which boundary MM terms are retained, with charge alteration being of secondary importance. In the worst case, only three intra-backbone boundaries may be used with relative confidence in the adequacy of resulting simulations, irrespective of the hyperpolarization mitigation scheme. Thus, advances in the field are certainly needed to fuel new discoveries. As of now, we believe that issues raised in this work might encourage authors in the field to report what boundaries, boundary MM terms, and charge redistribution schemes they are using, so their results may be correctly interpreted.

Automatic Generation of Local Vibrational Mode Parameters: From Small to Large Molecules and QM/MM Systems.

LModeAGen, a new protocol for the automatic determination of a nonredundant, complete set of local vibrational modes is reported, which is based on chemical graph concepts. Whereas local mode properties can be calculated for a selection of parameters targeting specific local modes of interest, a complete set of nonredundant local mode parameters is requested for the adiabatic connection scheme (ACS), relating each local vibrational mode with a normal mode counterpart, and for the decomposition of normal modes (CNM) in terms of local mode contributions, a unique way to analyze vibrational spectra. So far, nonredundant parameter sets have been generated manually following chemical intuition or from a set of redundant parameters in a trial-and-error fashion, which has hampered the study of larger systems with hundreds of parameters. LModeAGen was successfully applied for a test set of 11 systems, ranging from small molecules to the large QM (>100 atoms) subsystem of carbomonoxy-neuroglobin protein, described with a hybrid QM/MM method. The ωB97X-D/aug-cc-pVDZ, M06L/def2-TZVP, and QM/MM ωB97X-D/6-31G(d,p)/AMBER model chemistries were adopted for the description of the molecules in the test set. Our new protocol is an important step forward for a routine ACS and CNM analysis of the vibrational spectra of complex and large systems with hundreds of atoms, providing new access to important encoded electronic structure information.

Conical Intersections in Solution with Polarizable Embedding: Integral-Exact Direct Reaction Field.

A common strategy to exploring the properties and reactivity of complex systems is to use quantum mechanics/molecular mechanics (QM/MM) embedding, wherein a QM region is defined and treated with electronic structure theory, and the remainder of the system is treated with a force field. Important to the description of electronic excited states, especially those of charge-transfer character, is the treatment of the coupling between the QM and MM subsystems. The state of the art is to use a polarizable force field for the MM region and mutually couple the QM wavefunction and MM induced dipoles, in addition to the usual electrostatic embedding, yielding a polarizable embedding (QM/MM-Pol) approach. However, we showed previously that current popular QM/MM-Pol approaches exhibit issues of root flipping and/or incorrect descriptions of electronic crossings in multistate calculations [J. Chem. Theory Comput. 14, 2137 (2018)]. Here, we demonstrate a solution to these problems with an integral-exact reformulation of the direct reaction field approach of Thole and Van Duijnen (QM/MM-IEDRF). The resulting embedding potential includes one- and two-electron operators and many-body dipole-induced dipole interactions and thus includes a natural description of the screening of electron-electron interactions by the MM induced dipoles. Pauli repulsion from the environment is mimicked by effective core potentials on the MM atoms. Inherent to the DRF approach is the assumption that MM dipoles respond instantaneously to the positions of the QM electrons; therefore, dispersion interactions are captured approximately. All electronic states are eigenfunctions of the same Hamiltonian, while the polarization induced in the environment and the associated energetic stabilization are unique to each state. This allows for a consistent definition of transition properties and state crossings. We demonstrate QM/MM-IEDRF by exploring the influence of a (polarizable) inert xenon matrix environment on the conical intersection underlying the photoisomerization of ethylene.

Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase.

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the “oxyanion hole”, and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.

Three-Layer Multiscale Approach Based on Extremely Localized Molecular Orbitals to Investigate Enzyme Reactions.

Quantum mechanics/molecular mechanics (QM/MM) calculations are widely used embedding techniques to computationally investigate enzyme reactions. In most QM/MM computations, the quantum mechanical region is treated through density functional theory (DFT), which offers the best compromise between chemical accuracy and computational cost. Nevertheless, to obtain more accurate results, one should resort to wave function-based methods, which however lead to a much larger computational cost already for relatively small QM subsystems. To overcome this drawback, we propose the coupling of our QM/ELMO (quantum mechanics/extremely localized molecular orbital) approach with molecular mechanics, thus introducing the three-layer QM/ELMO/MM technique. The QM/ELMO strategy is an embedding method in which the chemically relevant part of the system is treated at the quantum mechanical level, while the rest is described through frozen ELMOs. Since the QM/ELMO method reproduces results of fully QM computations within chemical accuracy and with a much lower computational effort, it can be considered a suitable strategy to extend the range of applicability and accuracy of the QM/MM scheme. In this paper, other than briefly presenting the theoretical bases of the QM/ELMO/MM technique, we will also discuss its validation on the well-tested deprotonation of acetyl coenzyme A by aspartate in citrate synthase.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials.

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme–substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Atomistic Polarizable Embeddings: Energy, Dynamics, Spectroscopy, and Reactivity.

ConspectusThe computational modeling of realistic extended systems, relevant in, e.g., Chemistry and Biophysics, is a fundamental problem of paramount importance in contemporary research. Enzymatic catalysis and photoinduced processes in pigment–protein complexes are typical problems targeted by computer-aided approaches, to complement experiments as interpretative tools at a molecular scale. The daunting complexity of this task lies in between the opposite stringent requirements of results’ reliability for structural/dynamical properties and related intermolecular interactions, and a mandatory principle of realism in the modeling strategy. Therefore, in practice, a truly realistic computational model of a biologically relevant system can easily fail to meet the accuracy requirement, in order to balance the excessive computational cost necessary to reach the desired precision.To address such an “accuracy vs reality” dualistic requirement, mixed quantum mechanics/classical mechanics approaches within Atomistic (i.e., preserving the discrete particle configuration) Polarizable Embeddings (QM/APEs) methods have been proposed over the years. In this Account, we review recent developments in the design and application of general QM/APE methods, targeting situations where a local intrinsically quantum behavior is coupled to a large molecular system (i.e., an environment), often involving processes with different dynamical time scales, in order to avoid brute-force, unpractical quantum chemistry calculations on the complete system.In the first place, our interest is devoted to the available APEs models presently implemented in computational software, highlighting the quantum chemistry methods that can be used to treat the QM subsystem. We review the coupling strategy between the QM subsystem and the APE, which requires to examine the way the QM/MM mutual interactions are accounted for and how the polarization of the classical environment is considered with respect to (wrt) the quantum variables. Because of the need of reliable molecular and macromolecular structures, a pivotal aspect to address here is the handling of the system dynamics (i.e., gradients wrt nuclear positions are required), especially for large molecular assemblies composed by an overwhelming number of atoms, exploring many conformations on a complex energy landscape.Alongside, we highlight our views on the necessary steps to take toward more accurate general-purposes and transferable explicit embeddings. The main objective to achieve here is to design a more physically grounded multiscale approach. To do so, one should apply advanced new generation classical models to account for refined induction effects that are able to (i) improve the quality of QM/MM interaction energies; (ii) enhance transferability by avoiding the compulsory partial (or total) reparameterization of the classical model. Moreover, the extension of recent developments originating from the field of advanced classical molecular dynamics (MD) to the realm of QM/APE methods is a key direction to improve both speed and efficiency for the phase space exploration of systems of growing size and complexity.Lastly, we point out specific research topics where an advanced QM/APE dynamics can certainly shed some light. For example, we discuss chemical reactions in “harsh” environments and the case of spectroscopic theoretical modeling where the inclusion of refined environment effects is often mandatory.

Keto‐enol tautomerism of the 4,5‐dimethyl‐2‐(2′‐hydroxyphenyl)imidazole in water solution: Modeling equilibrium between neutral forms and accurate assignment of the absorption bands

AbstractTautomerizarion between keto, cis‐ and trans‐enol forms of the 4,5‐dimethyl‐2‐(2′‐hydroxyphenyl)imidazole (DMHI) is studied using molecular dynamics simulations with the quantum mechanics/molecular mechanics (QM/MM) potentials. We applied the live solvent selection approach that presumes division of water molecules to either QM or MM subsystems on the fly depending on their distances from the chromophore. This allowed us to treat a water shell around the chromophore at the QM level and to describe properly all chromophore‐water interactions. According to the calculated Gibbs energy profiles the most stable is a keto form and a cis‐enol form is ~1 kcal/mol higher. The trans‐enol form is much higher in energy and should not be observed in experiment. Vertical S0,min − S1 excitation energies of all neutral tautomers were calculated at the XMCQDPT2 level: both enol forms have similar values ~4 eV and keto form is red‐shifted, with the excitation energy of 3.3 eV. Thus, we conclude that two absorption bands in experimental absorption spectrum of DMHI in water solution are due to the presence of cis‐enol and keto forms of the chromophore.

Density Functional Theory-Based Quantum Mechanics/Coarse-Grained Molecular Mechanics: Theory and Implementation.

Quantum mechanics/molecular mechanics (QM/MM) is a standard computational tool for describing chemical reactivity in systems with many degrees of freedom, including polymers, enzymes, and reacting molecules in complex solvents. However, QM/MM is less suitable for systems with complex MM dynamics due to associated long relaxation times, the high computational cost of QM energy evaluations, and expensive long-range electrostatics. Recently, a systematic coarse graining of the MM part was proposed to overcome these QM/MM limitations in the form of the quantum mechanics/coarse-grained molecular mechanics (QM/CG-MM) approach. Herein, we recast QM/CG-MM in the density functional theory formalism and, by employing the force-matching variational principle, assess the method performance for the two model systems: QM CCl4 in the MM CCl4 liquid and the reaction of tert-butyl hypochlorite with the benzyl radical in the MM CCl4 solvent. We find that density functional theory (DFT)-QM/CG-MM accurately reproduces DFT-QM/MM radial distribution functions and three-body correlations between the QM and CG-MM subsystems. The free-energy profile of the reaction is also described well, with an error <1-2 kcal/mol. DFT-QM/CG-MM is a general, systematic, and computationally efficient approach to include chemical reactivity in coarse-grained molecular models.

The Impact of Electron Correlation on Describing QM/MM Interactions in the Attendant Molecular Dynamics Simulations of CO in Myoglobin

The impact of the dispersion and electron correlation effects on describing quantum mechanics/molecular mechanics (QM/MM) interactions in QM/MM molecular dynamics (MD) simulations was explored by performing a series of up to 2 ns QM/MM MD simulations on the B states of the myoglobin–carbon monoxide (MbCO) system. The results indicate that both dispersion and electron correlations play significant roles in the simulation of the ratios of two B states (B1/B2), which suggests that the inclusion of the electron correlation effects is essential for accurately modeling the interactions between QM and MM subsystems. We found that the QM/MM interaction energies between the CO and the surroundings statistically present a linear correlation with the electric fields along the CO bond. This indicates that QM/MM interactions can be described by a simple physical model of a dipole with constant moment under the action of the electric fields. The treatment provides us with an accurate and effective approach to account for the electron correlation effects in QM/MM MD simulations.

Efficient Analytic Second Derivative of Electrostatic Embedding QM/MM Energy: Normal Mode Analysis of Plant Cryptochrome.

Analytic second derivatives of electrostatic embedding (EE) quantum mechanics/molecular mechanics (QM/MM) energy are important for performing vibrational analysis and simulating vibrational spectra of quantum systems interacting with an environment represented as a classical electrostatic potential. The main bottleneck of EE-QM/MM second derivatives is the solution of coupled perturbed equations for each MM atom perturbation. Here, we exploit the Q-vector method [J. Chem. Phys., 2019, 151, 041102] to workaround this bottleneck. We derive the full analytic second derivative of the EE-QM/MM energy, which allows us to compute QM, MM, and QM-MM Hessian blocks in an efficient and easy to implement manner. To show the capabilities of our method, we compute the normal modes for the full Arabidopsis thaliana plant cryptochrome. We show that the flavin adenine dinucleotide vibrations (QM subsystem) strongly mix with protein modes. We compute approximate vibronic couplings for the lowest bright transition, from which we extract spectral densities and the homogeneous broadening of FAD absorption spectrum in protein using vibrationally resolved electronic spectrum simulations.

Simulating the Nonadiabatic Relaxation Dynamics of 4-(N,N-Dimethylamino)benzonitrile (DMABN) in Polar Solution.

The compound 4-(N,N-dimethylamino)benzonitrile (DMABN) represents the archetypal system for dual fluorescence, a rare photophysical phenomenon in which a given fluorophore shows two distinct emission bands. Despite extensive studies, the underlying mechanism remains the subject of debate. In the present contribution, we address this issue by simulating the excited-state relaxation process of DMABN as it occurs in polar solution. The potential energy surfaces for the system are constructed with the use of the additive quantum mechanics/molecular mechanics (QM/MM) method, and the coupled dynamics of the electronic wave function and the nuclei is propagated with the semiclassical fewest switches surface hopping method. The DMABN molecule, which comprises the QM subsystem, is treated with the use of the second-order algebraic diagrammatic construction (ADC(2)) method with the imposition of spin-opposite scaling (SOS). It is verified that this level of theory achieves a realistic description of the excited-state potential energy surfaces of DMABN. The simulation results qualitatively reproduce the main features of the experimentally observed fluorescence spectrum, thus allowing the unambiguous assignment of the two fluorescence bands: the normal band is due to the near-planar locally excited (LE) structure of DMABN, while the so-called "anomalous" second band arises from the twisted intramolecular charge transfer (TICT) structure. The transformation of the LE structure into the TICT structure takes place directly via intramolecular rotation, and is not mediated by another excited-state structure. In particular, the oft-discussed rehybridized intramolecular charge transfer (RICT) structure, which is characterized by a bent nitrile group, does not play a role in the relaxation process.

Dynamical properties of enzyme-substrate complexes disclose substrate specificity of the SARS-CoV-2 main protease as characterized by the electron density descriptors.

A dynamical approach is proposed to discriminate between reactive (rES) and nonreactive (nES) enzyme-substrate complexes taking the SARS-CoV-2 main protease (Mpro) as an important example. Molecular dynamics simulations with the quantum mechanics/molecular mechanics potentials (QM(DFT)/MM-MD) followed by the electron density analysis are employed to evaluate geometry and electronic properties of the enzyme with different substrates along MD trajectories. We demonstrate that mapping the Laplacian of the electron density and the electron localization function provides easily visible images of the substrate activation that allow one to distinguish rES and nES. The computed fractions of reactive enzyme-substrate complexes along MD trajectories well correlate with the findings of recent experimental studies on the substrate specificity of Mpro. The results of our simulations demonstrate the role of the theory level used in QM subsystems for a proper description of the nucleophilic attack of the catalytic cysteine residue in Mpro. The activation of the carbonyl group of a substrate is correctly characterized with the hybrid DFT functional PBE0, whereas the use of a GGA-type PBE functional, that lacks the admixture of the Hartree-Fock exchange fails to describe activation.

Polarizable Embedding with a Transferable H2O Potential Function I: Formulation and Tests on Dimer.

The incorporation of mutual polarization in multiscale simulations where different regions of the system are treated at different level of theory is important in studies of, for example, electronic excitations and charge transfer processes. We present here an energy functional for describing a quantum mechanics/molecular mechanics (QM/MM) scheme that includes reciprocal polarization between the two subsystems. The inclusion of polarization alleviates shortcomings inherent in electrostatic embedding QM/MM models based on point-charge force fields. A density functional theory (DFT) description of the QM subsystem is coupled to a single center multipole expansion (SCME) description of H2O molecules in the MM subsystem that includes anisotropic dipole and quadrupole polarizability as well as static multipoles up to and including the hexadecapole. The energy functional and the coupling scheme is general and can be extended to arbitrary order in terms of both the static and induced moments. Tests of the energy surface for the H2O dimer show that the QM/MM results lie in between the pure DFT and pure SCME values. The consistency of the many-body contributions to the energy and analytical forces is demonstrated for an H2O pentamer.

Tracking Proton Transfer through Titratable Amino Acid Side Chains in Adaptive QM/MM Simulations.

Ubiquitous throughout biological processes, proton transport usually occurs through the Grotthuss shuttling mechanism, where the hydrated proton exists as a charge or structural defect and is propagated quickly through a network of hydrogen bonds with minimal perturbation to the positions of the involved heavy atoms. The rapid reorganization of the bonding network and changing identity of the migrating proton can cause difficulties in molecular dynamics (MD) simulations. Previously, we formulated a proton indicator that tracks the proton as a structural defect for proton exchange between water molecules. In this work, we extend the proton indicator to treat proton transfer between water and amino acid side chains with titratable functional groups. Of particular interest are histidine, glutamate, and arginine, all of which have titratable groups featuring multiple protonation sites. Comparison with the modified center of excess charge (mCEC) suggests that the proton indicator and mCEC are both comparable in approximating the location of the proton. The location of the proton indicator was then used as the center of the QM subsystem in adaptive quantum-mechanical/molecular-mechanical (QM/MM) simulations of proton transport through a model channel along a path consisting of water and titratable amino acid side chains. In the adaptive QM/MM simulations, atoms were reclassified on the fly in a continuous and smooth manner as QM or MM depending on their distances from the proton indicator. Employing a small, mobile QM subsystem, the adaptive QM/MM simulations were found to be much more efficient than the conventional QM/MM simulations with a large QM subsystem that covered the entire pathway for proton relay.

Analytic QM/MM atomic charge derivatives avoiding the scaling of coupled perturbed equations with the MM subsystem size.

In electrostatic embedding mixed quantum and molecular mechanics (QM/MM) approaches, the QM charge distribution is polarized by the electrostatic interaction with the MM environment. Analytic derivatives of expectation values of operators are required to extract properties such as vibrational spectra. These derivatives usually require solving a set of coupled perturbed equations for each nucleus/atom in the system, thus becoming prohibitive when the MM subsystem contains thousands of atoms. In the context of Electrostatic Potential Fitting (ESPF) QM/MM, we can easily overcome this bottleneck by defining a set of auxiliary coupled perturbed equations called the Q-vector equations. The Q-vector method scales only with the size of the QM subsystem, producing an effective charge tensor that leads to the atomic charge derivative after contraction with the MM electrostatic potential gradient. As an example, we use the charge derivatives as an analysis tool to identify the most important chromophore-polarizing amino-acids in plant cryptochrome. This finding opens up the route of defining polarizable force fields and simulating vibrational spectroscopy using ESPF QM/MM electrostatic embedding at an affordable computational cost.

MiMiC: A Novel Framework for Multiscale Modeling in Computational Chemistry.

We present a flexible and efficient framework for multiscale modeling in computational chemistry (MiMiC). It is based on a multiple-program multiple-data (MPMD) model with loosely coupled programs. Fast data exchange between programs is achieved through the use of MPI intercommunicators. This allows exploiting the existing parallelization strategies used by the coupled programs while maintaining a high degree of flexibility. MiMiC has been used in a new electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) implementation coupling the highly efficient CPMD and GROMACS programs, but it can also be extended to use other programs. The framework can also be utilized to extend the partitioning of the system into several domains that can be treated using different models, such as models based on wave function or density functional theory as well as coarse-graining and continuum models. The new QM/MM implementation treats long-range electrostatic QM-MM interactions through the multipoles of the QM subsystem which substantially reduces the computational cost without loss of accuracy compared to an exact treatment. This enables QM/MM molecular dynamics (MD) simulations of very large systems.