Enzymes exist in ensembles of states that encode the energetics underlying their catalysis. Conformational ensembles built from 1231 structures of 17 serine proteases revealed atomic-level changes across their reaction states. By comparing the enzymatic and solution reaction, we identified molecular features that provide catalysis and quantified their energetic contributions to catalysis. Serine proteases precisely position their reactants in destabilized conformers, creating a downhill energetic gradient that selectively favors the motions required for reaction while limiting off-pathway conformational states. The same catalytic features have repeatedly evolved in proteases and additional enzymes across multiple distinct structural folds. Our ensemble-function analyses revealed previously unknown catalytic features, provided quantitative models based on simple physical and chemical principles, and identified motifs recurrent in nature that may inspire enzyme design.

The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics, using a combination of constant-pH molecular dynamics and semi-classical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena Sensory Rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties such as the isomerization quantum yield or decay rates. We decompose our results in single titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.

Red fluorescent protein (RFP) variants are highly sought after for in-vivo imaging since longer wavelengths improve depth and contrast in fluorescence imaging. However, the lower energy emission wavelength usually correlates with a lower fluorescent quantum yield than their green emitting counterparts. To guide the rational design of bright variants, we have theoretically assessed two variants (mScarlet and mRouge) which are reported to have very different brightness. Using an α-CASSCF QM/MM framework (chromophore and all protein residues within 6 Å of it in the QM region, for a total of more than 450 QM atoms), we identify key points on the ground and first excited state potential energy surfaces. The brighter variant mScarlet has a rigid scaffold, and the chromophore stays largely planar on the ground state. The dimmer variant mRouge shows more flexibility and can accommodate a pre-twisted chromophore conformation which provides easier access to conical intersections. The main difference between the variants lies in the intersection seam regions, which appear largely inaccessible in mScarlet but partially accessible in mRouge. This observation is mainly related with changes in the cavity charge distribution, the hydrogen-bonding network involving the chromophore and a key ARG/THR mutation (which changes both charge and steric hindrance).

The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. We describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google's protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ∼250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI-)-a model chromophore in green fluorescent proteins-on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.

The ab initio nanoreactor has previously been introduced to automate reaction discovery for ground state chemistry. In this work, we present the nonadiabatic nanoreactor, an analogous framework for excited state reaction discovery. We automate the study of nonadiabatic decay mechanisms of molecules by probing the intersection seam between adiabatic electronic states with hyper-real metadynamics, sampling the branching plane for relevant conical intersections, and performing seam-constrained path searches. We illustrate the effectiveness of the nonadiabatic nanoreactor by applying it to benzene, a molecule with rich photochemistry and a wide array of photochemical products. Our study confirms the existence of several types of S0/S1 and S1/S2 conical intersections which mediate access to a variety of ground state stationary points. We elucidate the connections between conical intersection energy/topography and the resulting photoproduct distribution, which changes smoothly along seam space segments. The exploration is performed with minimal user input, and the protocol requires no previous knowledge of the photochemical behavior of a target molecule. We demonstrate that the nonadiabatic nanoreactor is a valuable tool for the automated exploration of photochemical reactions and their mechanisms.

Over the last decade, artificial intelligence has been propelled forward by advances in machinelearning algorithms and computational hardware, opening up myriad new avenues for scientificresearch. Nevertheless, virtual assistants and voice control have yet to be widely utilized in thenatural sciences. Here, we present ChemVox, an interactive Amazon Alexa skill that uses speechrecognition to perform quantum chemistry calculations. This new application interfaces Alexawith cloud computing and returns the results through a capable device. ChemVox paves the wayto making computational chemistry routinely accessible to the wider community

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Anabaena Sensory Rhodopsin (ASR), a microbial photoactive protein featuring the retinal chromophore in two different conformations, exhibits a pH-dependent electronic absorption spectrum. Using the recently developed CpHMD-then-QM/MM multiscale protocol applied to ASR embedded in a membrane model, the pH-induced changes in its maximum absorption wavelength have been reproduced and analyzed. While the acidic tiny red-shift is essentially correlated with the deprotonation of an aspartic acid located on the ASR extracellular side, the larger blue-shift experimentally reported at pH values larger than 5 involves a cluster of titrating residues sitting on the cytoplasmic side. The ASR pH-dependent spectrum is the consequence of the competitive stabilization of retinal ground and excited states by the protein electrostatic potential.

The relaxation dynamicsof thymine and its derivatives thymidine and thymidine monophosphate werestudied using time-resolved photoelectron spectroscopy applied to a watermicrojet. Two absorption bands were studied, the first is a bright ππ* statewhich was populated using tunable-ultraviolet light in the range of 4.74 – 5.17eV and probed using a 6.20 eV probe pulse. By reversing the order of thesepulses, a band containing multiple ππ* states was populated by the 6.20 eVpulse and the lower energy pulse served as the probe. The lower lying ππ* statewas found to decay in ~400 fs in both thymine and thymidine independent of pumpphoton energy while thymidine monophosphate decays varied from 670-840 fs withsome pump energy dependence. The application of a computational QM/MM scheme atthe XMS-CASPT2//CASSCF/AMBER level of theory suggeststhat conformational differences existing between thymidine and thymidinemonophosphate in solution accounts for this difference. Thehigher lying ππ* band was found to decay in ~600 fs in all three cases, but wasonly able to be characterized when using the 5.17 eV probe pulse. Notably, nolong-lived signal from an np* state could be identified in eitherexperiment on any of the three molecules.

When a chromophore interacts with several titratable molecular sites, the modeling of its photophysical properties requires to take into account all their possible protonation states. We have developed a multi-scale protocol, based on constant-pH molecular dynamics simulations coupled to QM/MM excitation energy calculations, aimed at sampling both the phase space and protonation state space of a short polypeptide featuring a tyrosine-tryptophan dyad interacting with two aspartic acid residues. We show that such a protocol is accurate enough to help in the interpretation of the experimental tyrosine UV absorption spectrum at both acidic and basic pH. Moreover, it is confirmed that radical tryptophan probably contributes to the peptide spectrum, thanks to a UV-induced electron transfer from tyrosine to tryptophan, ultimately shedding light on the complex pH-dependent behavior of the peptide spectrum.

Discrepancies in the isomerization dynamics and quantum yields of the trans and cis retinal protonated Schiff base is a well-known issue in the context of retinal photochemistry. Anabaena Sensory Rhodopsin (ASR) is a microbial retinal protein that comprises a retinal chromophore in two ground state (GS) conformations: all-trans, 15-anti (AT) and 13-cis, 15-syn (13C). In this study, we applied impulsive vibrational spectroscopic techniques (DFWM, pump-DFWM and pump-IVS) to ASR to shed more light on how the structural changes take place in the excited state within the same protein environment. Our findings point to distinct features in the ground state structural conformations as well as to drastically different evolutions in the excited state manifold. The ground state vibrational spectra show stronger Raman activity of the C14-H out-of-plane wag (at about 805 cm-1) for the 13C isomer than that for the AT isomer, which hints at a pre-distortion of 13C in the ground state. Evolution of the Raman frequency after interaction with the actinic pulse shows a blue-shift for the C[double bond, length as m-dash]C stretching and CH3 rocking mode for both isomers. For AT, however, the blue-shift is not instantaneous as observed for the 13C isomer, rather it takes more than 200 fs to reach the maximum frequency shift. This frequency blue-shift is rationalized by a decrease in the effective conjugation length during the isomerization reaction, which further confirms a slower formation of the twisted state for the AT isomer and corroborates the presence of a barrier in the excited state trajectory previously predicted by quantum chemical calculations.

A minimal electrostatic model is introduced which aims at reproducing and analyzing the visible-light absorption energy shift of a protein with pH. It relies on the existence of a protein structure, the prediction of titratable amino-acid pKa values and a very limited set of parameters. Applied to the case of the photochromic Anabaena sensory rhodopsin protein, the model succeeds in reproducing qualitatively the reported experimental data, confirming the importance of aspartic acid 217 in the observed blue shift in the λmax of ASR at neutral pH. It also suggests for the first time the role of two other amino acids, glutamic acid 36 at basic pH and aspartic acid 120 at acidic pH.

While rotary molecular switches based on neutral and cationic organic π-systems have been reported, structurally homologous anionic switches providing complementary properties have not been prepared so far. Here we report the design and preparation of a molecular switch mimicking the anionic p-HBDI chromophore of the green fluorescent protein. The investigation of the mechanism and dynamics of the E/Z switching function is carried out both computationally and experimentally. The data consistently support axial rotary motion occurring on a sub-picosecond time scale. Transient spectroscopy and trajectory simulations show that the nonadiabatic decay process occurs in the vicinity of a conical intersection (CInt) between a charge transfer state and a covalent/diradical state. Comparison of our anionic p-HBDI-like switch with the previously reported cationic N-alkyl indanylidene pyrrolinium switch mimicking visual pigments reveals that these similar systems translocate, upon vertical excitation, a similar net charge in the same axial direction.

In light-driven single-molecule rotary motors, the photoisomerization of a double bond converts light energy into the rotation of a moiety (the rotor) with respect to another (the stator). However, at the level of a molecular population, an effective rotary motion can only be achieved if a large majority of the rotors rotate in the same, specific direction. Here we present a quantitative investigation of the directionality (clockwise vs counterclockwise) induced by a single stereogenic center placed in allylic position with respect to the reactive double bond of a model of the biomimetic indanylidene-pyrrolinium framework. By computing ensembles of nonadiabatic trajectories at 300 K, we predict that the photoisomerization is >70% unidirectional for the Z → E and E → Z conversions. Most importantly, we show that such directionality, resulting from the asymmetry of the excited state force field, can still be observed in the presence of a small (ca. 2°) pretwist or helicity of the reactive double bond. This questions the validity of the conjecture that a significant double-bond pretwist (e.g., >10°) in the ground state equilibrium structure of synthetic or natural rotary motors would be required for unidirectional motion.