From a humorous post to a detailed quantum-chemical study: isocyanate synthesis revisited.

Isocyanates play an essential role in modern manufacturing processes, especially in polyurethane production. There are numerous synthesis strategies for isocyanates both under industrial and laboratory conditions, which do not prevent searching for alternative highly efficient synthetic protocols. Here, we report a detailed theoretical investigation of the mechanism of sulfur dioxide-catalyzed rearrangement of phenylnitrile oxide into phenyl isocyanate, which was first reported in 1977. The DLPNO-CCSD(T) method and up-to-date DFT protocols were used to perform a highly accurate quantum-chemical study of the rearrangement mechanism. An overview of various organic and inorganic catalysts has revealed other potential catalysts, such as sulfur trioxide and selenium dioxide. Furthermore, the present study elucidated how substituents in phenylnitrile oxide influence reaction kinetics. This study was performed by a self-organized collaboration of scientists initiated by a humorous post on the VK social network.

Ultrafast dynamics of fluorene initiated by highly intense laser fields.

We present an investigation of the ultrafast dynamics of the polycyclic aromatic hydrocarbon fluorene initiated by an intense femtosecond near-infrared laser pulse (810 nm) and probed by a weak visible pulse (405 nm). Using a multichannel detection scheme (mass spectra, electron and ion velocity-map imaging), we provide a full disentanglement of the complex dynamics of the vibronically excited parent molecule, its excited ionic states, and fragments. We observed various channels resulting from the strong-field ionization regime. In particular, we observed the formation of the unstable tetracation of fluorene, above-threshold ionization features in the photoelectron spectra, and evidence of ubiquitous secondary fragmentation. We produced a global fit of all observed time-dependent photoelectron and photoion channels. This global fit includes four parent ions extracted from the mass spectra, 15 kinetic-energy-resolved ionic fragments extracted from ion velocity map imaging, and five photoelectron channels obtained from electron velocity map imaging. The fit allowed for the extraction of 60 lifetimes of various metastable photoinduced intermediates.

Resting-state functional connectivity correlates of brain structural aging in schizophrenia.

A large body of research has shown that schizophrenia patients demonstrate increased brain structural aging. Although this process may be coupled with aberrant changes in intrinsic functional architecture of the brain, they remain understudied. We hypothesized that there are brain regions whose whole-brain functional connectivity at rest is differently associated with brain structural aging in schizophrenia patients compared to healthy controls. Eighty-four male schizophrenia patients and eighty-six male healthy controls underwent structural MRI and resting-state fMRI. The brain-predicted age difference (b-PAD) was a measure of brain structural aging. Resting-state fMRI was applied to obtain global correlation (GCOR) maps comprising voxelwise values of the strength and sign of functional connectivity of a given voxel with the rest of the brain. Schizophrenia patients had higher b-PAD compared to controls (mean between-group difference + 2.9 years). Greater b-PAD in schizophrenia patients, compared to controls, was associated with lower whole-brain functional connectivity of a region in frontal orbital cortex, inferior frontal gyrus, Heschl's Gyrus, plana temporale and polare, insula, and opercular cortices of the right hemisphere (rFTI). According to post hoc seed-based correlation analysis, decrease of functional connectivity with the posterior cingulate gyrus, left superior temporal cortices, as well as right angular gyrus/superior lateral occipital cortex has mainly driven the results. Lower functional connectivity of the rFTI was related to worse verbal working memory and language production. Our findings demonstrate that well-established frontotemporal functional abnormalities in schizophrenia are related to increased brain structural aging.

On the thermodynamic stability of polycations.

We present a simple approximation to estimate the largest charge that a given molecule can hold until fragmentation into smaller charged species becomes more energetically favorable. This approximation solely relies on the ionization potentials, electron affinities of the parent and fragment species, and also on the neutral parent's dissociation energy. By parameterizing these quantities, it is possible to obtain analytical phase diagrams of polycationic stability. We demonstrate the applicability of this approach by discussing the maximal charge dependence on the size of the molecular system. A numerical demonstration for linear polyenes, monocyclic annulenes, and helium clusters is provided.

Metadynamics simulations with Bohmian-style bias potential.

Here, we present a parametrization of the metadynamics simulations for reactions involving breaking the chemical bonds along a single collective variable coordinate. The parameterization is based on the similarity between the bias potential in metadynamics and the quantum potential in the de Broglie-Bohm formalism. We derive the method and test it on two prototypical reaction types: proton transfer and breaking of the cyclohexene cycle (reversed Diels-Alder reaction).

Modulation of Slow Desensitization (Tachyphylaxis) of Acid-Sensing Ion Channel (ASIC)1a.

Among the proton-activated channels of the ASIC family, ASIC1a exhibits a specific tachyphylaxis phenomenon in the form of a progressive decrease in the response amplitude during a series of activations. This process is well known, but its mechanism is poorly understood. Here, we demonstrated a partial reversibility of this effect using long-term whole-cell recording of CHO cells transfected with rASIC1a cDNA. Thus, tachyphylaxis represents a slow desensitization of ASIC1a. Prolonged acidifications provided the same recovery from slow desensitization as short acidifications of the same frequency. Slow desensitization and steady-state desensitization are independent processes although the latter attenuates the development of the former. We found that drugs which facilitate ASIC1a activation (e.g., amitriptyline) cause an enhancement of slow desensitization, while inhibition of ASIC1a by 9-aminoacridine attenuates this process. Overall, for a broad variety of exposures, including increased calcium concentration, different pH conditions, and modulating drugs, we found a correlation between their effects on ASIC1a response amplitude and the development of slow desensitization. Thus, our results demonstrate that slow desensitization occurs only when ASIC1a is in the open state.

Analysis of residue-residue interactions in the structures of ASIC1a suggests possible gating mechanisms.

The gating mechanism of acid-sensitive ion channels (ASICs) remains unclear, despite the availability of atomic-scale structures in various functional states. The collapse of the acidic pocket and structural changes in the low-palm region are assumed to trigger activation. For the acidic pocket, protonation of some residues can minimize repulsion in the collapsed conformation. The relationship between low-palm rearrangements and gating is unknown. In this work, we performed a Monte Carlo energy optimization of known ASIC1a structures and determined the residue-residue interactions in different functional states. For rearrangements in the acidic pocket, our results are consistent with previously proposed mechanisms, although significant complexity was revealed for the residue-residue interactions. The data support the proposal of a gating mechanism in the low-palm region, in which residues E80 and E417 share a proton to activate the channel.

Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling.

In this work we discuss the generally applicable Wigner sampling and introduce a new, simplified Wigner sampling method, for computationally effective modeling of molecular properties containing nuclear quantum effects and vibrational anharmonicity. For various molecular systems test calculations of (a) vibrationally averaged rotational constants, (b) vibrational IR spectra and (c) photoelectron spectra have been performed. The performance of Wigner sampling has been assessed by comparing with experimental data and with results of other theoretical models, including harmonic and VPT2 approximations. The developed simplified Wigner sampling method shows advantages in application to large and flexible molecules.

Regularized weighted sine least-squares spectral analysis for gas electron diffraction data.

Here, we present a new approach for obtaining radial distribution functions (RDF) from the electron diffraction data using a regularized weighted sine least-squares spectral analysis. It allows for explicitly transferring the measured experimental uncertainties in the reduced molecular scattering function to the produced RDF. We provide a numerical demonstration, discuss the uncertainties and correlations in the RDFs, and suggest a regularization parameter choice criterion. The approach is also applicable for other diffraction data, e.g., for x-ray or neutron diffraction of liquid samples.

Light-Sensitive Open Channel Block of Ionotropic Glutamate Receptors by Quaternary Ammonium Azobenzene Derivatives.

Glutamate ionotropic receptors mediate fast excitation processes in the central nervous system of vertebrates and play an important role in synaptic plasticity, learning, and memory. Here, we describe the action of two azobenene-containing compounds, AAQ (acrylamide-azobenzene-quaternary ammonium) and QAQ (quaternary ammonium-azobenzene-quaternary ammonium), which produced rapid and fully reversible light-dependent inhibition of glutamate ionotropic receptors. The compounds demonstrated voltage-dependent inhibition with only minor voltage-independent allosteric action. Calcium-impermeable AMPA receptors had weaker sensitivity compared to NMDA and calcium-permeable AMPA receptors. We further revealed that the compounds bound to NMDA and calcium-permeable AMPA receptors in different modes. They were able to enter the wide selectivity filter of AMPA receptors, and strong negative voltages caused permeation into the cytoplasm. The narrow selectivity filter of the NMDA receptors did not allow the molecules to bypass them; therefore, QAQ and AAQ bound to the shallow channel site and prevented channel closure by a foot-in-the-door mechanism. Computer simulations employing available AMPA and NMDA receptor structures readily reproduced the experimental findings, allowing for the structure-based design of more potent and selective drugs in the future. Thus, our work creates a framework for the development of light-sensitive blockers of calcium-permeable AMPA receptors, which are desirable tools for neuroscience.

Possible Compensatory Role of ASICs in Glutamatergic Synapses.

Proton-gated channels of the ASIC family are widely distributed in central neurons, suggesting their role in common neurophysiological functions. They are involved in glutamatergic neurotransmission and synaptic plasticity; however, the exact function of these channels remains unclear. One problem is that acidification of the synaptic cleft due to the acidic content of synaptic vesicles has opposite effects on ionotropic glutamate receptors and ASICs. Thus, the pH values required to activate ASICs strongly inhibit AMPA receptors and almost completely inhibit NMDA receptors. This, in turn, suggests that ASICs can provide compensation for post-synaptic responses in the case of significant acidifications. We tested this hypothesis by patch-clamp recordings of rat brain neuron responses to acidifications and glutamate receptor agonists at different pH values. Hippocampal pyramidal neurons have much lower ASICs than glutamate receptor responses, whereas striatal interneurons show the opposite ratio. Cortical pyramidal neurons and hippocampal interneurons show similar amplitudes in their responses to acidification and glutamate. Consequently, the total response to glutamate agonists at different pH levels remains rather stable up to pH 6.2. Besides these pH effects, the relationship between the responses mediated by glutamate receptors and ASICs depends on the presence of Mg2+ and the membrane voltage. Together, these factors create a complex picture that provides a framework for understanding the role of ASICs in synaptic transmission and synaptic plasticity.

Glutamate potentiates heterologously expressed homomeric acid-sensing ion channel 1a.

Acid-sensing ion channels (ASICs) participate in synaptic transmission due to the acidic content of synaptic vesicles, but their contribution to postsynaptic currents is small. This has stimulated attempts to find endogenous ASIC potentiators that could enhance ASIC-mediated currents to physiologically relevant values. Here we demonstrate that glutamate, which serves as a neurotransmitter, potentiates recombinant ASIC1a in the submillimolar concentration range. The effect of glutamate is especially interesting as ASIC's presence has been shown in glutamatergic synapses. At pH=6.5 glutamate had maximum potentiation of 87% with an EC50 value of 0.65 mM. The mechanism of potentiation is due to a shift of pH-dependent activation to less acidic values, with 0.5 mM glutamate increasing pH50 from 6.04 to 6.43. Due to this mechanism, ASIC1a in glutamatergic synapses might be intrinsically potentiated. Furthermore, this effect could compensate for the inhibition of ionotropic glutamate receptors by extracellular acidification during synaptic transmission.

The kinetic energy of PAH dication and trication dissociation determined by recoil-frame covariance map imaging.

We investigated the dissociation of dications and trications of three polycyclic aromatic hydrocarbons (PAHs), fluorene, phenanthrene, and pyrene. PAHs are a family of molecules ubiquitous in space and involved in much of the chemistry of the interstellar medium. In our experiments, ions are formed by interaction with 30.3 nm extreme ultraviolet (XUV) photons, and their velocity map images are recorded using a PImMS2 multi-mass imaging sensor. Application of recoil-frame covariance analysis allows the total kinetic energy release (TKER) associated with multiple fragmentation channels to be determined to high precision, ranging 1.94-2.60 eV and 2.95-5.29 eV for the dications and trications, respectively. Experimental measurements are supported by Born-Oppenheimer molecular dynamics (BOMD) simulations.

Computational Analysis of the Crystal and Cryo-EM Structures of P-Loop Channels with Drugs.

The superfamily of P-loop channels includes various potassium channels, voltage-gated sodium and calcium channels, transient receptor potential channels, and ionotropic glutamate receptors. Despite huge structural and functional diversity of the channels, their pore-forming domain has a conserved folding. In the past two decades, scores of atomic-scale structures of P-loop channels with medically important drugs in the inner pore have been published. High structural diversity of these complexes complicates the comparative analysis of these structures. Here we 3D-aligned structures of drug-bound P-loop channels, compared their geometric characteristics, and analyzed the energetics of ligand-channel interactions. In the superimposed structures drugs occupy most of the sterically available space in the inner pore and subunit/repeat interfaces. Cationic groups of some drugs occupy vacant binding sites of permeant ions in the inner pore and selectivity-filter region. Various electroneutral drugs, lipids, and detergent molecules are seen in the interfaces between subunits/repeats. In many structures the drugs strongly interact with lipid and detergent molecules, but physiological relevance of such interactions is unclear. Some eukaryotic sodium and calcium channels have state-dependent or drug-induced π-bulges in the inner helices, which would be difficult to predict. The drug-induced π-bulges may represent a novel mechanism of gating modulation.

Double Proton Transfer Across a Table: The Formic Acid Dimer-Fluorobenzene Complex.

Proton transfer via tunneling is a fundamental quantum-mechanical phenomenon. We report rotational spectroscopy measurements of this process in the complex of the formic acid dimer with fluorobenzene. The assignment of the spectrum indicates that this complex exists in the form of a π-π stacked structure. Each rotational transition of the parent isotopologue exhibits splitting. Isotopic substitution experiments show that the spectral splitting results from double-proton transfer tunneling in the formic acid dimer. Presence of fluorobenzene as a neighboring molecule does not quench the double proton transfer in the formic acid dimer but decreases its tunneling splitting from 341(3) MHz to 267.608(1) MHz. Calculations suggest that the presence of the weakly bounded fluorobenzene does not influence the activation energy of the proton transfer. The fluorobenzene is reoriented with respect to the formic acid dimer during the course of the reaction, slowing down the proton transfer motion.

Optical control of L-type Ca2+ channels using a diltiazem photoswitch.

L-type Ca2+ channels (LTCCs) play a crucial role in excitation-contraction coupling and release of hormones from secretory cells. They are targets of antihypertensive and antiarrhythmic drugs such as diltiazem. Here, we present a photoswitchable diltiazem, FHU-779, which can be used to reversibly block endogenous LTCCs by light. FHU-779 is as potent as diltiazem and can be used to place pancreatic ß-cell function and cardiac activity under optical control.

Phenylalkylamines in calcium channels: computational analysis of experimental structures.

Experimental 3D structures of calcium channels with phenylalkylamines (PAAs) provide basis for further analysis of atomic mechanisms of these important cardiovascular drugs. In the crystal structure of the engineered calcium channel CavAb with Br-verapamil and in the cryo-EM structure of the Cav1.1 channel with verapamil, the ligands bind in the inner pore. However, there are significant differences between these structures. In the crystal structure the ligand ammonium group is much closer to the ion in the selectivity-filter region Site 3, which is most proximal to the inner pore, than in the cryo-EM structure. Here we used Monte Carlo energy minimizations to dock PAAs in calcium channels. Our computations suggest that in the crystal structure Site 3 is occupied by a water molecule rather than by a calcium ion. Analysis of the published electron density map does not rule out this possibility. In the cryo-EM structures the ammonium group of verapamil is shifted from the calcium ion in Site 3 either along the pore axis, towards the cytoplasm or away from the axis. Our unbiased docking reproduced these binding modes. However, in the cryo-EM structures detergent and lipid molecules interact with verapamil. When we removed these molecules, the nitrile group of verapamil bound to the calcium ion in Site 3. Models of Cav1.2 with different PAAs suggest similar binding modes and direct contacts of the ligands electronegative atoms with the calcium ion in Site 3. Such interactions explain paradoxes in structure-activity relationships of PAAs.

Icosahedral Carbaboranes with Peripheral Hydrogen-Chalcogenide Groups: Structures from Gas Electron Diffraction and Chemical Shielding in Solution.

Carbaboranes 1,2-(EH)2 -closo-1,2-C2 B10 H10 (E=S, Se) were prepared, in the case of E=Se for the first time. Their semi-experimental equilibrium molecular structures were established by the concerted use of quantum-chemical calculations and gas electron diffraction. A method was developed and implemented to quantify the contribution of experimental data to each refined structural parameter. The accuracy of the experimental structures and those calculated at the MP2 level of theory were gauged by comparison of experimental 11 B NMR chemical shifts with quantum-chemically computed values; the inclusion of electron correlation (GIAO-MP2) provided superior results. For the purpose of geometrical prediction, the remaining group 16 elements were considered, and the icosahedral structures for E=O and Te were also computed; for E=O the same theoretical approach was used as for E=S, and for E=Te a description similar to that for E=Se was employed.

Extremely Potent Block of Bacterial Voltage-Gated Sodium Channels by µ-Conotoxin PIIIA.

µ-Conotoxin PIIIA, in the sub-picomolar, range inhibits the archetypal bacterial sodium channel NaChBac (NavBh) in a voltage- and use-dependent manner. Peptide µ-conotoxins were first recognized as potent components of the venoms of fish-hunting cone snails that selectively inhibit voltage-gated skeletal muscle sodium channels, thus preventing muscle contraction. Intriguingly, computer simulations predicted that PIIIA binds to prokaryotic channel NavAb with much higher affinity than to fish (and other vertebrates) skeletal muscle sodium channel (Nav 1.4). Here, using whole-cell voltage clamp, we demonstrate that PIIIA inhibits NavBac mediated currents even more potently than predicted. From concentration-response data, with [PIIIA] varying more than 6 orders of magnitude (10-12 to 10-5 M), we estimated an IC50 = ~5 pM, maximal block of 0.95 and a Hill coefficient of 0.81 for the inhibition of peak currents. Inhibition was stronger at depolarized holding potentials and was modulated by the frequency and duration of the stimulation pulses. An important feature of the PIIIA action was acceleration of macroscopic inactivation. Docking of PIIIA in a NaChBac (NavBh) model revealed two interconvertible binding modes. In one mode, PIIIA sterically and electrostatically blocks the permeation pathway. In a second mode, apparent stabilization of the inactivated state was achieved by PIIIA binding between P2 helices and trans-membrane S5s from adjacent channel subunits, partially occluding the outer pore. Together, our experimental and computational results suggest that, besides blocking the channel-mediated currents by directly occluding the conducting pathway, PIIIA may also change the relative populations of conducting (activated) and non-conducting (inactivated) states.

Hydrophobic Amines and Their Guanidine Analogues Modulate Activation and Desensitization of ASIC3.

Acid-sensing ion channel 3 (ASIC3) is an important member of the acid-sensing ion channels family, which is widely expressed in the peripheral nervous system and contributes to pain sensation. ASICs are targeted by various drugs and toxins. However, mechanisms and structural determinants of ligands' action on ASIC3 are not completely understood. In the present work we studied ASIC3 modulation by a series of "hydrophobic monoamines" and their guanidine analogs, which were previously characterized to affect other ASIC channels via multiple mechanisms. Electrophysiological analysis of action via whole-cell patch clamp method was performed using rat ASIC3 expressed in Chinese hamster ovary (CHO) cells. We found that the compounds studied inhibited ASIC3 activation by inducing acidic shift of proton sensitivity and slowed channel desensitization, which was accompanied by a decrease of the equilibrium desensitization level. The total effect of the drugs on the sustained ASIC3-mediated currents was the sum of these opposite effects. It is demonstrated that drugs' action on activation and desensitization differed in their structural requirements, kinetics of action, and concentration and state dependencies. Taken together, these findings suggest that effects on activation and desensitization are independent and are likely mediated by drugs binding to distinct sites in ASIC3.