Influence of the Greater Protein Environment on the Electrostatic Potential in Metalloenzyme Active Sites: The Case of Formate Dehydrogenase.

The Mo/W-containing metalloenzyme formate dehydrogenase (FDH) is an efficient and selective natural catalyst that reversibly converts CO2 to formate under ambient conditions. In this study, we investigate the impact of the greater protein environment on the electrostatic potential (ESP) of the active site. To model the enzyme environment, we used a combination of classical molecular dynamics and multiscale quantum-mechanical (QM)/molecular-mechanical (MM) simulations. We leverage charge shift analysis to systematically construct QM regions and analyze the electronic environment of the active site by evaluating the degree of charge transfer between the core active site and the protein environment. The contribution of the terminal chalcogen ligand to the ESP of the metal center is substantial and dependent on the chalcogen identity, with similar, less negative ESPs for Se and S terminal chalcogens in comparison to O regardless of whether the metal is Mo or W. The orientation of the side chains and conformations of the cofactor also affect the ESP, highlighting the importance of sampling dynamic fluctuations in the protein. Overall, our observations suggest that the terminal chalcogen ligand identity plays an important role in the enzymatic activity of FDH, suggesting opportunities for a rational bioinspired catalyst design.

Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins.

Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.

Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning.

Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.

Revealing quantum mechanical effects in enzyme catalysis with large-scale electronic structure simulation.

Enzymes have evolved to facilitate challenging reactions at ambient conditions with specificity seldom matched by other catalysts. Computational modeling provides valuable insight into catalytic mechanism, and the large size of enzymes mandates multi-scale, quantum mechanical-molecular mechanical (QM/MM) simulations. Although QM/MM plays an essential role in balancing simulation cost to enable sampling with full QM treatment needed to understand electronic structure in enzyme active sites, the relative importance of these two strategies for understanding enzyme mechanism is not well known. We explore challenges in QM/MM for studying the reactivity and stability of three diverse enzymes: i) Mg2+-dependent catechol O-methyltransferase (COMT), ii) radical enzyme choline trimethylamine lyase (CutC), and iii) DNA methyltransferase (DNMT1), which has structural Zn2+ binding sites. In COMT, strong non-covalent interactions lead to long range coupling of electronic structure properties across the active site, but the more isolated nature of the metallocofactor in DNMT1 leads to faster convergence of some properties. We quantify these effects in COMT by computing covariance matrices of by-residue electronic structure properties during dynamics and along the reaction coordinate. In CutC, we observe spontaneous bond cleavage following initiation events, highlighting the importance of sampling and dynamics. We use electronic structure analysis to quantify the relative importance of CHO and OHO non-covalent interactions in imparting reactivity. These three diverse cases enable us to provide some general recommendations regarding QM/MM simulation of enzymes.

Quantum Mechanical Description of Electrostatics Provides a Unified Picture of Catalytic Action Across Methyltransferases.

Methyl transferases (MTases) are a well-studied class of enzymes for which competing enzymatic enhancement mechanisms have been suggested, ranging from structural methyl group CH···X hydrogen bonds (HBs) to electrostatic- and charge-transfer-driven stabilization of the transition state (TS). We identified all Class I MTases for which reasonable resolution (<2.0 Å) crystal structures could be used to form catalytically competent ternary complexes for multiscale (i.e., quantum-mechanical/molecular-mechanical or QM/MM) simulation of the SN2 methyl transfer reaction coordinate. The four Class I MTases studied have both distinct functions (e.g., protein repair or biosynthesis) and substrate nucleophiles (i.e., C, N, or O). While CH···X HBs stabilize all reactant complexes, no universal TS stabilization role is found for these interactions in MTases. A consistent picture is instead obtained through analysis of charge transfer and electrostatics, wherein much of cofactor-substrate charge separation is maintained in the TS region, and electrostatic potential is correlated with substrate nucleophilicity (i.e., intrinsic reactivity).

Communication: Observation of local-bender eigenstates in acetylene.

We report the observation of eigenstates that embody large-amplitude, local-bending vibrational motion in acetylene by stimulated emission pumping spectroscopy via vibrational levels of the S1 state involving excitation in the non-totally symmetric bending modes. The N(b) = 14 level, lying at 8971.69 cm(-1) (J = 0), is assigned on the basis of degeneracy due to dynamical symmetry breaking in the local-mode limit. The level pattern for the N(b) = 16 level, lying at 10 218.9 cm(-1), is consistent with expectations for increased separation of ℓ = 0 and 2 vibrational angular momentum components. Increasingly poor agreement between our observations and the predicted positions of these levels highlights the failure of currently available normal mode effective Hamiltonian models to extrapolate to regions of the potential energy surface involving large-amplitude displacement along the acetylene ⇌ vinylidene isomerization coordinate.

Simplified Cartesian basis model for intrapolyad emission intensities in the bent-to-linear electronic transition of acetylene.

The acetylene emission spectrum from the trans-bent electronically excited à state to the linear ground electronic X̃ state has attracted considerable attention because it grants Franck­Condon access to local bending vibrational levels of the X̃ state with large-amplitude motion along the acetylene ⇌ vinylidene isomerization coordinate. For emission from the ground vibrational level of the à state, there is a simplifying set of Franck­Condon propensity rules that gives rise to only one zero-order bright state per conserved vibrational polyad of the X̃ state. Unfortunately, when the upper level involves excitation in the highly admixed ungerade bending modes, ν4' and ν6', the simplifying Franck­Condon propensity rule breaks down--as long as the usual polar basis (with v and l quantum numbers) is used to describe the degenerate bending vibrations of the X̃ state--and the intrapolyad intensities result from complicated interference patterns between many zero-order bright states. In this article, we show that, when the degenerate bending levels are instead treated in the Cartesian two-dimensional harmonic oscillator basis (with vx and vy quantum numbers), the propensity for only one zero-order bright state (in the Cartesian basis) is restored, and the intrapolyad intensities are simple to model, as long as corrections are made for anharmonic interactions. As a result of trans ⇌ cis isomerization in the à state, intrapolyad emission patterns from overtones of ν4' and ν6' evolve as quanta of trans bend (ν3') are added, so the emission intensities are not only relevant to the ground-state acetylene ⇌ vinylidene isomerization, they are also a direct reporter of isomerization in the electronically excited state.

Cis-trans isomerization in the S1 state of acetylene: identification of cis-well vibrational levels.

A systematic analysis of the S(1)-trans (Ã(1)A(u)) state of acetylene, using IR-UV double resonance along with one-photon fluorescence excitation spectra, has allowed assignment of at least part of every single vibrational state or polyad up to a vibrational energy of 4200 cm(-1). Four observed vibrational levels remain unassigned, for which no place can be found in the level structure of the trans-well. The most prominent of these lies at 46 175 cm(-1). Its (13)C isotope shift, exceptionally long radiative lifetime, unexpected rotational selection rules, and lack of significant Zeeman effect, combined with the fact that no other singlet electronic states are expected at this energy, indicate that it is a vibrational level of the S(1)-cis isomer (Ã(1)A(2)). Guided by ab initio calculations [J. H. Baraban, A. R. Beck, A. H. Steeves, J. F. Stanton, and R. W. Field, J. Chem. Phys. 134, 244311 (2011)] of the cis-well vibrational frequencies, the vibrational assignments of these four levels can be established from their vibrational symmetries together with the (13)C isotope shift of the 46 175 cm(-1) level (assigned here as cis-3(1)6(1)). The S(1)-cis zero-point level is deduced to lie near 44 900 cm(-1), and the ν(6) vibrational frequency of the S(1)-cis well is found to be roughly 565 cm(-1); these values are in remarkably good agreement with the results of recent ab initio calculations. The 46 175 cm(-1) vibrational level is found to have a 3.9 cm(-1) staggering of its K-rotational structure as a result of quantum mechanical tunneling through the isomerization barrier. Such tunneling does not give rise to ammonia-type inversion doubling, because the cis and trans isomers are not equivalent; instead the odd-K rotational levels of a given vibrational level are systematically shifted relative to the even-K rotational levels, leading to a staggering of the K-structure. These various observations represent the first definite assignment of an isomer of acetylene that was previously thought to be unobservable, as well as the first high resolution spectroscopic results describing cis-trans isomerization.

Design and evaluation of a pulsed-jet chirped-pulse millimeter-wave spectrometer for the 70-102 GHz region.

Chirped-pulse millimeter-wave (CPmmW) spectroscopy is the first broadband (multi-GHz in each shot) Fourier-transform technique for high-resolution survey spectroscopy in the millimeter-wave region. The design is based on chirped-pulse Fourier-transform microwave (CP-FTMW) spectroscopy [G. G. Brown, B. C. Dian, K. O. Douglass, S. M. Geyer, S. T. Shipman, and B. H. Pate, Rev. Sci. Instrum. 79, 053103 (2008)], which is described for frequencies up to 20 GHz. We have built an instrument that covers the 70-102 GHz frequency region and can acquire up to 12 GHz of spectrum in a single shot. Challenges to using chirped-pulse Fourier-transform spectroscopy in the millimeter-wave region include lower achievable sample polarization, shorter Doppler dephasing times, and problems with signal phase stability. However, these challenges have been partially overcome and preliminary tests indicate a significant advantage over existing millimeter-wave spectrometers in the time required to record survey spectra. Further improvement to the sensitivity is expected as more powerful broadband millimeter-wave amplifiers become affordable. The ability to acquire broadband Fourier-transform millimeter-wave spectra enables rapid measurement of survey spectra at sufficiently high resolution to measure diagnostically important electronic properties such as electric and magnetic dipole moments and hyperfine coupling constants. It should also yield accurate relative line strengths across a broadband region. Several example spectra are presented to demonstrate initial applications of the spectrometer.

Darling-Dennison resonance and Coriolis coupling in the bending overtones of the A 1A(u) state of acetylene, C2H2.

Rotational analyses have been carried out for the overtones of the nu(4) (torsion) and nu(6) (in-plane cis-bend) vibrations of the A (1)A(u) state of C(2)H(2). The v(4)+v(6)=2 vibrational polyad was observed in high-sensitivity one-photon laser-induced fluorescence spectra and the v(4)+v(6)=3 polyad was observed in IR-UV double resonance spectra via the ground state nu(3) (Sigma(+) (u)) and nu(3)+nu(4) (Pi(u)) vibrational levels. The structures of these polyads are dominated by the effects of vibrational angular momentum: Vibrational levels of different symmetry interact via strong a-and b-axis Coriolis coupling, while levels of the same symmetry interact via Darling-Dennison resonance, where the interaction parameter has the exceptionally large value K(4466)=-51.68 cm(-1). The K-structures of the polyads bear almost no resemblance to the normal asymmetric top patterns, and many local avoided crossings occur between close-lying levels with nominal K-values differing by one or more units. Least squares analysis shows that the coupling parameters change only slightly with vibrational excitation, which has allowed successful predictions of the structures of the higher polyads: A number of weak bands from the v(4)+v(6)=4 and 5 polyads have been identified unambiguously. The state discovered by Scherer et al. [J. Chem. Phys. 85, 6315 (1986)], which appears to interact with the K=1 levels of the 3(3) vibrational state at low J, is identified as the second highest of the five K=1 members of the v(4)+v(6)=4 polyad. After allowing for the Darling-Dennison resonance, the zero-order bending structure can be represented by omega(4)=764.71, omega(6)=772.50, x(44)=0.19, x(66)=-4.23, and x(46)=11.39 cm(-1). The parameters x(46) and K(4466) are both sums of contributions from the vibrational angular momentum and from the anharmonic force field. For x(46) these contributions are 14.12 and -2.73 cm(-1), respectively, while the corresponding values for K(4466) are -28.24 and -23.44 cm(-1). It is remarkable how severely the coupling of nu(4) and nu(6) distorts the overtone polyads, and also how in this case the effects of vibrational angular momentum outweigh those of anharmonicity in causing the distortion.

Contrasting singlet-triplet dynamical behavior of two vibrational levels of the acetylene S1 2(1)3(1)B2 polyad.

Surface electron ejection by laser-excited metastables (SEELEM) and LIF spectra of acetylene were simultaneously recorded in the regions of the A1Au-X1Sigmag+ nominal 2(1)3(1)4(2) Ka=1<--00 and 2(1)3(1)6(2) Ka=1<--00 bands near 46,140 cm(-1). The upper states of these two bands are separated by only approximately 100 cm(-1), and the two S1 vibrational levels are known to be strongly mixed by anharmonic and Coriolis interactions. Strikingly different patterns were observed in the SEELEM spectra in the regions of the 2(1)3(1)4(2) and 2(1)3(1)6(2) vibrational levels. Because the equilibrium structure of the T3 electronic state is known to be nonplanar, excitation of nu4 (torsion) and nu6 (antisymmetric in-plane bend) are expected respectively to promote and suppress vibrational overlap between low-lying S1 and T3 vibrational levels. The nearly 50:50 mixed 2(1)3(1)4(2)-2(1)3(1)6(2) character of the S1 vibrational levels rules out this simple Franck-Condon explanation for the different appearance of the SEELEM spectra. A simple model is applied to the SEELEM/LIF spectra to explain the differences between spectral patterns in terms of a T3 doorway-mediated singlet-triplet coupling model.

Observation of the A1A" state of isocyanogen.

The A1A" state of isocyanogen, CNCN, is observed using photofragment fluorescence excitation spectroscopy in a room temperature cell and in a molecular beam. The spectra are highly congested, but progressions that correspond to the Franck-Condon active C-N-C bending vibration in the excited state are evident. Linewidth measurements indicate that the excited state lifetime is <10 ps. These measurements are consistent with previous ab initio calculations, which predicted a bent excited state with a short lifetime due to predissociation. Although we do not believe that we have observed the origin band of the electronic transition, we place an upper limit of 42,523 cm(-1) on the energy of the excited state zero point level.

Beam action spectroscopy via inelastic scattering.

In this article, a new technique we call Beam Action Spectroscopy via Inelastic Scattering (BASIS) is demonstrated. BASIS takes advantage of the sensitivity of rotational state distributions in a supersonic molecular beam to inelastic scattering within the beam. We exploit BASIS to achieve increased sensitivity in two very different types of experiments. In the first, the UV photodissociation spectrum of OClO is recovered by monitoring intensity changes in the pure rotational transition of a spectator molecule (OCS) downstream from the nozzle, revealing a new vibrational structure in the region between 30,000 and 36,000 cm(-1). In the second, the mid-IR vibrational spectrum of acetylene is recorded simply by monitoring a single pure rotational transition of OCS co-expanded with acetylene. The technique may prove particularly fruitful when an excitation process produces product dark states that are not easily probed by conventional spectroscopy.

Electronic signatures of large amplitude motions: dipole moments of vibrationally excited local-bend and local-stretch states of S0 acetylene.

A one-dimensional local bend model is used to describe the variation of electronic properties of acetylene in vibrational levels that embody large amplitude local motions on the S0 potential energy surface. Calculations performed at the CCSD(T) and MR-AQCC levels of theory predict an approximately linear dependence of the dipole moment on the number of quanta in either the local bending or local stretching excitation. In the local mode limit, one quantum of stretching excitation in one CH bond leads to an increase of 0.025 D in the dipole moment, and one quantum of bending vibration in the CCH angle leads to an increase of 0.068 D. The use of a one-dimensional model for the local bend is justified by comparison to the well-established polyad model which reveals a decoupling of the large amplitude bending from other degrees of freedom in the range of Nbend = 14-22. We find that the same one-dimensional large amplitude bending motion emerges from two profoundly different representations, a one-dimensional cut through an ab initio, seven-dimensional Hamiltonian and the three-dimensional (l = 0) pure-bending experimentally parametrized spectroscopic Hamiltonian.

Millimeter-wave-detected, millimeter-wave optical polarization spectroscopy.

We report a new form of microwave optical double-resonance spectroscopy called millimeter-wave-detected, millimeter-wave optical polarization spectroscopy (mmOPS). In contrast to other forms of polarization spectroscopy, in which the polarization rotation of optical beams is detected, the mmOPS technique is based on the polarization rotation of millimeter waves induced by the anisotropy from optical pumping out of the lower or upper levels of the millimeter wave transition. By monitoring ground-state rotational transitions with the millimeter waves, the mmOPS technique is capable of identifying weak or otherwise difficult-to-observe optical transitions in complex chemical environments, where multiple molecular species or vibrational states can lead to spectral congestion. Once a transition is identified, mmOPS can then be used to record pure rotational transitions in vibrationally and electronically excited states, with the resolution limited only by the radiative decay rate. Here, the sensitivity of this nearly-background-free technique is demonstrated by optically pumping the weak, nominally spin-forbidden CS e (3)Sigma(-)-X (1)Sigma(+) (2-0) and d (3)Delta-X (1)Sigma(+) (6-0) electronic transitions while probing the CS X (1)Sigma(+) (v(")=0,J(")=2-1) rotational transition with millimeter waves. The J(')=2,N(')=2<--J(')=1,N(')=1 pure rotational transition of the CS e (3)Sigma(-) (v(')=2) state is then recorded by optically preparing the J(')=1,N(')=1 level of the e (3)Sigma(-) (v(')=2) state via the J(')=1,N(')=1<--J(")=1 transition of the e (3)Sigma(-)-X (1)Sigma(+) (2-0) band.

CH-stretching overtone spectroscopy of 1,1,1,2-tetrafluoroethane.

We have recorded the vibrational absorption spectrum of 1,1,1,2-tetrafluoroethane (HFC-134a) in the fundamental and first five CH-stretching overtone regions with the use of Fourier transform infrared, dispersive long-path, intracavity laser photoacoustic, and cavity ringdown spectroscopies. We compare our measured total oscillator strengths in each region with intensities calculated using an anharmonic oscillator local mode model. We calculate intensities with 1D, 2D, and 3D Hamiltonians, including one or two CH stretches and two CH stretches with the HCH bending mode, respectively. The dipole moment function is calculated ab initio with self-consistent-field Hartree-Fock and density functional theories combined with double- and triple-zeta-quality basis sets. We find that the basis set choice affects the total intensity more than the choice of the Hamiltonian. We achieve agreement between the calculated and measured total intensities of approximately a factor of 2 or better for the fundamental and first five overtones.