Conical Intersection Accessibility Dictates Brightness in Red Fluorescent Proteins.

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 compared to 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 pretwisted 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).

Modeling pH-Dependent Biomolecular Photochemistry.

The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical 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 into 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.

TeraChem protocol buffers (TCPB): Accelerating QM and QM/MM simulations with a client-server model.

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 non-adiabatic nanoreactor: towards the automated discovery of photochemistry.

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.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

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.

Voice-controlled quantum chemistry.

Over the past decade, artificial intelligence has been propelled forward by advances in machine learning algorithms and computational hardware, opening up myriads of new avenues for scientific research. Nevertheless, virtual assistants and voice control have yet to be widely used in the natural sciences. Here, we present ChemVox, an interactive Amazon Alexa skill that uses speech recognition to perform quantum chemistry calculations. This new application interfaces Alexa with cloud computing and returns the results through a capable device. ChemVox paves the way to making computational chemistry routinely accessible to the wider community.

Relaxation Dynamics of Hydrated Thymine, Thymidine, and Thymidine Monophosphate Probed by Liquid Jet Time-Resolved Photoelectron Spectroscopy.

The relaxation dynamics of thymine and its derivatives thymidine and thymidine monophosphate are studied using time-resolved photoelectron spectroscopy applied to a water microjet. Two absorption bands are studied; the first is a bright ππ* state which is populated using tunable-ultraviolet light in the range 4.74-5.17 eV and probed using a 6.20 eV probe pulse. By reversing the order of these pulses, a band containing multiple ππ* states is populated by the 6.20 eV pulse and the lower energy pulse serves as the probe. The lower lying ππ* state is found to decay in ∼400 fs in both thymine and thymidine independent of pump photon energy, while thymidine monophosphate decays vary from 670 to 840 fs with some pump energy dependence. The application of a computational quantum mechanical/molecular mechanical scheme at the XMS-CASPT2//CASSCF/AMBER level of theory suggests that conformational differences existing between thymidine and thymidine monophosphate in solution account for this difference. The higher lying ππ* band is found to decay in ∼600 fs in all three cases, but it is only able to be characterized when the 5.17 eV probe pulse is used. Notably, no long-lived signal from an nπ* state can be identified in either experiment on any of the three molecules.

CpHMD-Then-QM/MM Identification of the Amino Acids Responsible for the Anabaena Sensory Rhodopsin pH-Dependent Electronic Absorption Spectrum.

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.

Media Framing and the Threat of Global Pandemics: The Ebola Crisis in UK Media and Policy Response.

Pandemics pose new and difficult challenges. Risks associated with the spread of pandemics generate intense speculation in Western media. Taking the 2014-2015 Ebola outbreak as a case study, the article critically analyses how the risk of contagion in the US, Europe, and the UK has been constructed in UK media and policy discourse. Drawing on the importance of media framing in shaping a given problem definition, causal interpretation and treatment recommendation, the article critically assesses the impacts of the British newspaper framing of Ebola, questioning the rationale of a UK domestic political response based on containment and border screenings. The article also takes a comparative angle, engaging with constructions of previous pandemics. Underscoring the importance of a sociological analysis of these framings, the article critically reflects on the role of media communication in reproducing certain topoi, which reduce the scope for open public debate around best responses to a pandemic emergency.

Mapping the ultrafast vibrational dynamics of all-trans and 13-cis retinal isomerization in Anabaena Sensory Rhodopsin.

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.

Sampling the protonation states: the pH-dependent UV absorption spectrum of a polypeptide dyad.

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.

pH-Dependent absorption spectrum of a protein: a minimal electrostatic model of Anabaena sensory rhodopsin.

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.

Design, Synthesis, and Dynamics of a Green Fluorescent Protein Fluorophore Mimic with an Ultrafast Switching Function.

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.

Directionality of Double-Bond Photoisomerization Dynamics Induced by a Single Stereogenic Center.

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.

Risky individuals and the politics of genetic research into aggressiveness and violence.

New genetic technologies promise to generate valuable insights into the aetiology of several psychiatric conditions, as well as a wider range of human and animal behaviours. Advances in the neurosciences and the application of new brain imaging techniques offer a way of integrating DNA analysis with studies that are looking at other biological markers of behaviour. While candidate 'genes for' certain conditions, including schizophrenia and bipolar disorders, are said to be 'un-discovered' at a faster rate than they are discovered, many studies are being conducted on personality traits such as aggressiveness and anti-social traits. The clinical applicability and implications of these studies are often discussed within the scientific community. However, little attention has so far been paid to their possible policy implications in relation to criminality management and to Criminal Law itself. Similarly, the related ethical issues arising in the field of crime control, and the tensions between enhancing security for society and protecting civil liberties, are currently under-explored. This paper investigates these ethical issues by focusing on the views of those professionals - including judges, lawyers, probation officers and social workers - who work with individuals 'deemed at risk' of violent and aggressive behaviours. It also discusses and problematizes mainstream rhetoric and arguments around the notion of 'risky individuals'.