Theoretical Design, Synthesis, and In Vitro Neurobiological Applications of a Highly Efficient Two-Photon Caged GABA Validated on an Epileptic Case.

In this paper, we present an additional, new cage-GABA compound, called 4-amino-1-(4'-dimethylaminoisopropoxy-5',7'-dinitro-2',3'-dihydro-indol-1-yl)-1-oxobutane-γ-aminobutyric acid (iDMPO-DNI-GABA), and currently, this compound is the only photoreagent, which can be applied for GABA uncaging without experimental compromises. By a systematic theoretical design and successful synthesis of several compounds, the best reagent exhibits a high two-photon efficiency within the 700-760 nm range with excellent pharmacological behavior, which proved to be suitable for a complex epileptic study. Quantum chemical design showed that the optimal length of the cationic side chain enhances the two-photon absorption by 1 order of magnitude due to the cooperating internal hydrogen bonding to the extra nitro group on the core. This feature increased solubility while suppressing membrane permeability. The efficiency was demonstrated in a systematic, wide range of in vitro single-cell neurophysiological experiments by electrophysiological as well as calcium imaging techniques. Scalable inhibitory ion currents were elicited by iDMPO-DNI-GABA with appropriate spatial-temporal precision, blocking both spontaneous and evoked cell activity with excellent efficiency. Additionally, to demonstrate its applicability in a real neurobiological study, we could smoothly and selectively modulate neuronal activities during artificial epileptic rhythms first time in a neural network of GCaMP6f transgenic mouse brain slices.

Reply to the 'Comment on "Penicillin's catalytic mechanism revealed by inelastic neutrons and quantum chemical theory"' by S. A. Glover, Phys. Chem. Chem. Phys., 2019, 21, 18012.

Herein we comparatively comment on the molecular metric 'amidicity', a descriptor of amide reactivity, and differing methods to determining it; with focus on lactam-rings. Specifically, our established amidicity percentage (AM%) approach is quantitatively contrasted with the transamidation (TA) method. This comment is organised into two sections, firstly addressing the differing methods in context of the computational bases of amidicity. This is followed by the quantitative demonstration that although both the AM% and HRS methods provide estimates of resonance enthalpy (ΔHRE), the former is more reliable across a wider set of systems. The robustness of the AM% approach is affirmed by quantitative matching of experimental NMR and kinetics measurements tracking changes in amide reactivities, including in Penicillin arising from modulation of its amide group and environmental effects.

Oxidatively-mediated in silico epimerization of a highly amyloidogenic segment in the human calcitonin hormone (hCT15-19).

In order to study the effects of peptide exposure to oxidative attack, we chose a model reaction in which the hydroxyl radical discretely abstracts a hydrogen atom from the α-carbon of each residue of a highly amyloidogenic region in the human calcitonin hormone, hCT15-19. Based on a combined Molecular Mechanics / Quantum Mechanics approach, the extended and folded L- and D-configuration and radical intermediate hCT15-19 peptides were optimized to obtain their compactness, secondary structure and relative thermodynamic data. The results suggest that the epimerization of residues is generally an exergonic process that can explain the cumulative nature of molecular aging. Moreover, the configurational inversion induced conformational changes can cause protein dysfunction. The epimerization of the central residue to the D-configuration induced a hairpin structure in hCT15-19, concomitant with a possible oligomerization of human calcitonin into Aß(1-42)-like amyloid fibrils present in patients suffering from Alzheimer's disease.

Synthesis and spectroscopic characterization of novel GFP chromophore analogues based on aminoimidazolone derivatives.

In order to improve the fluorescence properties of the green fluorescent protein chromophore, p­HOBDI ((5­(4­hydroxybenzylidene)­2,3­dimethyl­3,5­dihydro­4H­imidazol­4­one), sixteen dihydroimidazolone derivates were synthesized from thiohydantoin and arylaldehydes. The synthesis developed is an efficient, novel, one-pot procedure. The study provides a detailed description of the spectroscopic characteristics of the newly synthesized compounds, using p­HOBDI as a reference. The new compounds all exhibited significantly stronger fluorescence than p­HOBDI, up to 28 times higher quantum yields. An experimental and theoretical investigation of the relationship of the fluorescence properties with the molecular structure was also carried out. A good correlation was found between the emission wavenumber and the Hammett constant of the functional group, which suggests the intermolecular charge transfer (ICT) mechanism between the aromatic groups.

Formation Mechanism of Benzo(a)pyrene: One of the Most Carcinogenic Polycyclic Aromatic Hydrocarbons (PAH).

The formation of polycyclic aromatic hydrocarbons (PAHs) is a strong global concern due to their harmful effects. To help the reduction of their emissions, a crucial understanding of their formation and a deep exploration of their growth mechanism is required. In the present work, the formation of benzo(a)pyrene was investigated computationally employing chrysene and benz(a)anthracene as starting materials. It was assumed a type of methyl addition/cyclization (MAC) was the valid growth mechanism in this case. Consequently, the reactions implied addition reactions, ring closures, hydrogen abstractions and intramolecular hydrogen shifts. These steps of the mechanism were computed to explore benzo(a)pyene formation. The corresponding energies of the chemical species were determined via hybrid density funcional theory (DFT), B3LYP/6-31+G(d,p) and M06-2X/6-311++G(d,p). Results showed that the two reaction routes had very similar trends energetically, the difference between the energy levels of the corresponding molecules was just 6.13 kJ/mol on average. The most stable structure was obtained in the benzo(a)anthracene pathway.

Amide Activation in Ground and Excited States.

Not all amide bonds are created equally. The purpose of the present paper is the reinterpretation of the amide group by means of two concepts: amidicity and carbonylicity. These concepts are meant to provide a new viewpoint in defining the stability and reactivity of amides. With the help of simple quantum-chemical calculations, practicing chemists can easily predict the outcome of a desired process. The main benefit of the concepts is their simplicity. They provide intuitive, but quasi-thermodynamic data, making them a practical rule of thumb for routine use. In the current paper we demonstrate the performance of our methods to describe the chemical character of an amide bond strength and the way of its activation methods. Examples include transamidation, acyl transfer and amide reductions. Also, the method is highly capable for simple interpretation of mechanisms for biological processes, such as protein splicing and drug mechanisms. Finally, we demonstrate how these methods can provide information about photo-activation of amides, through the examples of two caged neurotransmitter derivatives.

High efficiency two-photon uncaging coupled by the correction of spontaneous hydrolysis.

Two-photon (TP) uncaging of neurotransmitter molecules is the method of choice to mimic and study the subtleties of neuronal communication either in the intact brain or in slice preparations. However, the currently available caged materials are just at the limit of their usability and have several drawbacks. The local and focal nature of their use may for example be jeopardized by a high spontaneous hydrolysis rate of the commercially available compounds with increased photochemical release rate. Here, using quantum chemical modelling we show the mechanisms of hydrolysis and two-photon activation, and synthesized more effective caged compounds. Furthermore, we have developed a new enzymatic elimination method removing neurotransmitters inadvertently escaping from their compound during experiment. This method, usable both in one and two-photon experiments, allows for the use of materials with an increased rate of photochemical release. The efficiency of the new compound and the enzymatic method and of the new compound are demonstrated in neurophysiological experiments.

Protein Stability and Unfolding Following Glycine Radical Formation.

Glycine (Gly) residues are particularly susceptible to hydrogen abstraction; which results in the formation of the capto-dative stabilized Cα-centered Gly radical (GLR) on the protein backbone. We examined the effect of GLR formation on the structure of the Trp cage; tryptophan zipper; and the villin headpiece; three fast-folding and stable miniproteins; using all-atom (OPLS-AA) molecular dynamics simulations. Radicalization changes the conformation of the GLR residue and affects both neighboring residues but did not affect the stability of the Trp zipper. The stability of helices away from the radical center in villin were also affected by radicalization; and GLR in place of Gly15 caused the Trp cage to unfold within 1 µs. These results provide new evidence on the destabilizing effects of protein oxidation by reactive oxygen species.

Radical Formation Initiates Solvent-Dependent Unfolding and ß-sheet Formation in a Model Helical Peptide.

We examined the effects of Cα-centered radical formation on the stability of a model helical peptide, N-Ac-KK(AL)10KK-NH2. Three, 100 ns molecular dynamics simulations using the OPLS-AA force field were carried out on each α-helical peptide in six distinct binary TIP4P water/2,2,2-trifluoroethanol (TFE) mixtures. The α-helicity was at a maximum in 20% TFE, which was inversely proportional to the number of H-bonds between water molecules and the peptide backbone. The radial distribution of TFE around the peptide backbone was highest in 20% TFE, which enhanced helix stability. The Cα-centered radical initiated the formation of a turn within 5 ns, which was a smaller kink at high TFE concentrations, and a loop at lower TFE concentrations. The highest helicity of the peptide radical was measured in 100% TFE. The formation of hydrogen bonds between the peptide backbone and water destabilized the helix, whereas the clustering of TFE molecules around the radical center stabilized the helix. Following radical termination, the once helical structure converted to a ß-sheet rich state in 100% water only, and this transition did not occur in the nonradical control peptide. This study gives evidence on how the formation of peptide radicals can initiate α-helical to ß-sheet transitions under oxidative stress conditions.

Glutathione as a prebiotic answer to α-peptide based life.

The energetics of peptide bond formation is an important factor not only in the design of chemical peptide synthesis, but it also has a role in protein biosynthesis. In this work, quantum chemical calculations at 10 different levels of theory including G3MP2B3 were performed on the energetics of glutathione formation. The strength of the peptide bond is found to be closely related to the acid strength of the to-be N-terminal and the basicity of the to-be C-terminal amino acid. It is shown that the formation of the first peptide activates the amino acid for the next condensation step, manifested in bacterial protein synthesis where the first step is the formation of an N-formylmethionine dipeptide. The possible role of glutathione in prebiotic molecular evolution is also analyzed. The implications of the thermodynamics of peptide bond formation in prebiotic peptide formation as well as in the preference of α- instead of ß- or γ-amino acids are discussed. An empirical correction is proposed for the compensation of the error due to the incapability of continuum solvation models in describing the change of the first solvation shell when a peptide bond is formed from two zwitterions accompanied by the disappearance of one ion pair.

Molecular ageing: free radical initiated epimerization of thymopentin--a case study.

The epimerization of amino acid residues increases with age in living organisms. In the present study, the structural consequences and thermodynamic functions of the epimerization of thymopentin (TP-5), the active site of the thymic hormone thymopoietin, were studied using molecular dynamics and density functional theory methods. The results show that free radical-initiated D-amino acid formation is energetically favoured (-130 kJmol(-1)) for each residue and induces significant changes to the peptide structure. In comparison to the wild-type (each residue in the L-configuration), the radius of gyration of the D-Asp(3) epimer of the peptide decreased by 0.5 Å, and disrupted the intramolecular hydrogen bonding of the native peptide. Beyond establishing important structural, energetic and thermodynamic benchmarks and reference data for the structure of TP-5, these results disseminate the understanding of molecular ageing, the epimerization of amino acid residues.

Atropisomerism of the Asn α radicals revealed by Ramachandran surface topology.

C radicals are typically trigonal planar and thus achiral, regardless of whether they originate from a chiral or an achiral C-atom (e.g., C-H + (•)OH → C• + H2O). Oxidative stress could initiate radical formation in proteins when, for example, the H-atom is abstracted from the Cα-carbon of an amino acid residue. Electronic structure calculations show that such a radical remains achiral when formed from the achiral Gly, or the chiral but small Ala residues. However, when longer side-chain containing proteogenic amino acid residues are studied (e.g., Asn), they provide radicals of axis chirality, which in turn leads to atropisomerism observed for the first time for peptides. The two enantiomeric extended backbone structures, •ßL and •ßD, interconvert via a pair of enantiotopic reaction paths, monitored on a 4D Ramachandran surface, with two distinct transition states of very different Gibbs-free energies: 37.4 and 67.7 kJ/mol, respectively. This discovery requires the reassessment of our understanding on radical formation and their conformational and stereochemical behavior. Furthermore, the atropisomerism of proteogenic amino acid residues should affect our understanding on radicals in biological systems and, thus, reframes the role of the D-residues as markers of molecular aging.

Glutathione--hydroxyl radical interaction: a theoretical study on radical recognition process.

Non-reactive, comparative (2 × 1.2 µs) molecular dynamics simulations were carried out to characterize the interactions between glutathione (GSH, host molecule) and hydroxyl radical (OH(•), guest molecule). From this analysis, two distinct steps were identified in the recognition process of hydroxyl radical by glutathione: catching and steering, based on the interactions between the host-guest molecules. Over 78% of all interactions are related to the catching mechanism via complex formation between anionic carboxyl groups and the OH radical, hence both terminal residues of GSH serve as recognition sites. The glycine residue has an additional role in the recognition of OH radical, namely the steering. The flexibility of the Gly residue enables the formation of further interactions of other parts of glutathione (e.g. thiol, α- and ß-carbons) with the lone electron pair of the hydroxyl radical. Moreover, quantum chemical calculations were carried out on selected GSH/OH(•) complexes and on appropriate GSH conformers to describe the energy profile of the recognition process. The relative enthalpy and the free energy changes of the radical recognition of the strongest complexes varied from -42.4 to -27.8 kJ/mol and from -21.3 to 9.8 kJ/mol, respectively. These complexes, containing two or more intermolecular interactions, would be the starting configurations for the hydrogen atom migration to quench the hydroxyl radical via different reaction channels.

Penicillin's catalytic mechanism revealed by inelastic neutrons and quantum chemical theory.

Penicillin, travels through bodily fluids, targeting and acylatively inactivating enzymes responsible for cell-wall synthesis in gram-positive bacteria. Somehow, it avoids metabolic degradation remaining inactive en route. To resolve this ability to switch from a non-active, to a highly reactive form, we investigated the dynamic structure-activity relationship of penicillin by inelastic neutron spectroscopy, reaction kinetics, NMR and multi-scale theoretical modelling (QM/MM and post-HF ab initio). Results show that by a self-activating physiological pH-dependent two-step proton-mediated process, penicillin changes geometry to activate its irreversibly reactive acylation, facilitated by systemic intramolecular energy management and cooperative vibrations. This dynamic mechanism is confirmed by the first ever reported characterisation of an antibiotic by neutrons, achieved on the TOSCA instrument (ISIS facility, RAL, UK).

Application of the systems chemistry approach on the ammonolysis of 1-ethoxycarbonyl- and 1-phenoxycarbonyl-3-(2-thienyl)oxindoles. A method to predict reactivity.

The routine prediction of the reactivity of a complex, multifunctional molecule is a challenging and time-consuming procedure. In the last step of the synthesis of the well-known drug substance tenidap, a nonexpected difference was observed between the reactivities of two closely related carbamate moieties, the N-ethoxycarbonyl and the N-phenoxycarbonyl group. A detailed kinetic study, necessitating a significant computational effort, is described in the present paper for this reaction step. On the other hand, the systems chemistry concept, by analyzing the details of the electronic structure and the connections between functional groups in a fast and simple way, is also able to answer this question using various "-icity" parameters (aromaticity, carbonylicity, olefinicity). The complete systems chemistry approach involves all these conjugativicity parameters, while its further simplified version is based on only one key parameter, which is carbonylicity in the present case. The above methods were compared in terms of their predictive power. The results show that the systems chemistry concept, even its one-parameter version, is applicable for the characterization of this challenging reactivity issue.

Qualitative assessment of microstructure and Hertzian indentation failure in biocompatible glass ionomer cements.

Discs of biocompatible glass ionomer cements were prepared for Hertzian indentation and subsequent fracture analyses. Specifically, 2 × 10 mm samples for reproducing bottom-initiated radial fracture, complemented by 0.2 × 1 mm samples for optimal resolution with X-ray micro tomography (µCT), maintaining dimensional ratio. The latter allowed for accurate determination of volumetric-porosity of the fully cured material, fracture-branching through three Cartesian axes and incomplete bottom-initiated cracking. Nanocomputed tomography analyses supported the reliability of the µCT results. Complementary 2-dimensional fractographic investigation was carried out by optical and scanning electron microscopies on the larger samples, identifying fracture characteristics. The combined 3-D qualitative assessment of microstructure and fractures, complemented by 2-D methods, provided an increased understanding of the mechanism of mechanical failure in these cements. Specifically, cracks grew to link pores while propagating along glass-matrix interfaces. The methodological development herein is exploitable on related biomaterials and represents a new tool for the rational characterisation, optimisation and design of novel materials for clinical service.

Conformation-dependent ˙OH/H2O2 hydrogen abstraction reaction cycles of Gly and Ala residues: a comparative theoretical study.

To determine if (•)OH can initiate the unfolding of an amino acid residue, the elementary reaction coordinates of H abstraction by (•)OH different conformations (ß(L), γ(L), γ(D), α(L), and α(D)) of Gly and Ala dimethyl amides were computed using first-principles quantum computations. The MPWKCIS1K/6-311++G(3df,2p)//BHandHLYP/6-311+G(d,p) level of theory was selected after different combinations of functionals and basis sets were compared. The structures of Gly and Ala in the elementary reaction steps were compared to the conformers of the Gly, Gly(•), Ala, and Ala(•) structures in the absence of (•)OH/H(2)O, which were identified by optimizing the minima of the respective potential energy surfaces. A dramatic change in conformation is observed in the Gly and Ala conformers after conversion to Gly(•) and Ala(•), respectively, and this change can be monitored along the minimal energy pathway. The ß(L) conformer of Gly (-0.3 kJ mol(-1)) and Ala (-1.6 kJ mol(-1)) form the lowest-lying transition states in the reaction with (•)OH, whereas the side chain of Ala strongly destabilizes the α conformers compared to the γ conformers, which could cause the lower reactivity shown in Ala. This effect shown in Ala could affect the abstraction of hydrogen from Ala and the other chiral amino acid residues in the helices. The energy of subsequent hydrogen abstraction reactions between Ala(•) and Gly(•) and H(2)O(2) remains approximately 90 kJ mol(-1) below the entrance level of the (•)OH reaction, indicating that the (•)OH radical can initiate an α to ß transition in an amino acid residue if a molecule such as H(2)O(2) can provide the hydrogen atom necessary to re-form Gly and Ala. This work delineates the mechanism of the rapid (•)OH-initiated unfolding of peptides and proteins which has been proposed in Alzheimer's and other peptide misfolding diseases involving amyloidogenic peptides.

Exploiting diverse stereochemistry of ß-amino acids: toward a rational design of sheet-forming ß-peptide systems.

Due to the two methylene groups in their backbone, ß-amino acids can adopt numerous secondary structures, including helices, sheets and nanotubes. Chirality introduced by the additional side chains can significantly influence the folding preference of ß-peptides composed of chiral ß-amino acids. However, only conceptual suggestions are present in the literature about the effect of chirality on folding preferences. Summarizing both the experimental and computational results, Seebach (Chem Biodivers 1:1111-1240, 2004) has proposed the first selection rule on the effect of side chain chirality, on the folding preference of ß-peptides. In order to extend and fine-tune the aforementioned predictions of Seebach, we have investigated its validity to the novel type of apolar sheet proposed recently (Pohl et al. in J Phys Chem B 114:9338-9348, 2010). In order to facilitate the rational design of sheet-like structures, a systematic study on the effect of chirality on "apolar" sheet stability is presented on disubstituted [HCO-ß-Ala-ß(2,3)-hAla-ß-Ala-NH(2)](2) model peptides calculated at the M05-2X/6-311++G(d,p)//M05-2X/6-31G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) levels of theory both in vacuum and in polar and apolar solvents. In addition, both types of "apolar" sheets were investigated; the one with two strands of identical (AA) and enantiomeric (AB) backbone structure. Our results show that heterochirally disubstituted sheets have the greatest preference for sheet formation (ΔG ~ -11 kcal mol(-1)). However, in contrast to Seebach's predictions, "homochiral disubstitution" itself does not necessarily disrupt the sheet structure, rather it could result stable fold (ΔG ~ -5 kcal mol(-1)). Results indicate that both the methyl group orientation and the local conformational effect of substitution affects sheet stability, as point chirality was found to have influence only on the backbone torsional angles. These results enabled us to extend and generalize Seebach's predictions and to propose a more general and accurate "rule of thumb" describing the effect of chirality on sheet stability. This offers an easy-to-use summary on how to design ß-peptide sheet structures. We conclude that heterochirally disubstituted models are the best candidates for sheet formation, if the two strands are substituted in a way to create identical torsional angle sets on the two backbones for ideal hydrogen-bonding pattern. With adequately selected side chains, homochirally disubtituted derivatives may also form sheet structures, and the position of methyl groups would prevent assembly of more than two strands making it ideal to create hairpins.

The Effect of Newly Developed OPLS-AA Alanyl Radical Parameters on Peptide Secondary Structure.

Recent studies using ab initio calculations have shown that Cα-centered radical formation by H-abstraction from the backbone of peptide residues has dramatic effects on peptide structure and have suggested that this reaction may contribute to the protein misfolding observed in Alzheimer's and Parkinson's diseases. To enable the effects of Cα-centered radicals to be studied in longer peptides and proteins over longer time intervals, force-field parameters for the Cα-centered Ala radical were developed for use with the OPLS force field by minimizing the sum of squares deviation between the quantum chemical and OPLS-AA energy hypersurfaces. These parameters were used to determine the effect of the Cα-centered Ala radical on the structure of a hepta-alanyl peptide in molecular dynamics (MD) simulations. A negligible sum-of-squares energy deviation was observed in the stretching parameters, and the newly developed OPLS-AA torsional parameters showed a good agreement with the LMP2/cc-pVTZ(-f) hypersurface. The parametrization also demonstrated that derived force-field bond length and bond angle parameters can deviate from the quantum chemical equilibrium values, and that the improper torsional parameters should be developed explicitly with respect to the coupled torsional parameters. The MD simulations showed planar conformations of the Cα-containing residue (Alr) are preferred and these conformations increase the formation of γ-, α-, and π-turn structures depending on the position in the turn occupied by the Alr residue. Higher-ordered structures are destabilized by Alr except when this residue occupies position "i + 1" of the 310-helix. These results offer new insight into the protein-misfolding mechanisms initiated by H-abstraction from the Cα of peptide and protein residues.

The role of entropy in initializing the aggregation of peptides: a first principle study on oligopeptide oligomerization.

The initiation and progression of Alzheimer's disease is coupled to the oligo- and polymerization of amyloid peptides in the brain. Amyloid like aggregates of protein domains were found practically independent of their primary sequences. Thus, the driving force of the transformation from the original to a disordered amyloid fold is expected to lie in the protein backbone common to all proteins. In order to investigate the thermodynamics of oligomerization, full geometry optimizations and frequency calculations were performed both on parallel and antiparallel ß-pleated sheet model structures of [HCO-(Ala)(1-6)-NH(2)](2) and (For-Ala(1-2)-NH(2))(1-6) peptides, both at the B3LYP and M05-2X/6-311++G(d,p)//M05-2X/6-31G(d) levels of theory, both in vacuum and in water. Our results show that relative entropy and enthalpy both show a hyperbolic decrease with increasing residue number and with increasing number of strands as well. Thus, di- and oligomerization are always thermodynamically favored. Antiparallel arrangements were found to have greater stability than parallel arrangements of the polypeptide backbones. During our study the relative changes in thermodynamic functions are found to be constant for long enough peptides, indicating that stability and entropy terms are predictable. All thermodynamic functions of antiparallel di- and oligomers show a staggered nature along the increasing residue number. By identifying and analyzing the 6 newly emerging dimer vibrational modes of the 10- and 14-membered building units, the staggered nature of the entropy function can be rationalized. Thus, the vanishing rotational and translational modes with respect to single strands are converted into entropy terms "holding tight" the dimers and oligomers formed, rationalizing the intrinsic adherence of natural polypeptide backbones to aggregate.