Pencil graphite electrodes for in-situ spectroelectrochemical sensing of reaction intermediates and products in organic solvents.

BACKGROUND: Spectroelectrochemistry (SEC) is a valuable analytical tool providing insights to reaction mechanisms and the structure of species involved in charge transfer reactions. Most of commercial SEC setups are based on platinum working electrodes where the adsorption of species involved in reactions often complicates their analysis. RESULTS: In this work, we employ an array of pencil graphite rods as an optically transparent working electrode in a custom-made air-tight thin-layer cell suitable for the SEC analysis performed here in acetonitrile as a representative non-aqueous solvent. The functionality of the device was demonstrated by UV-Vis SEC sensing of charge transfer reactions of ruthenium acetylacetonate, ferrocene and ethylviologen dibromide redox probes performed employing the cyclic voltammetry. The SEC response obtained for all three probes confirmed no adsorption and the absence of oxygen in the cell. Furthermore, we have developed and utilized finite element method numerical simulations considering charge transfer reactions coupled with the diffusional mass transport to model the cyclic voltammetric response and the reaction conversion in the thin-layer SEC cell. SIGNIFICANCE: Our work paves the way for easy-to-assemble customized air-tight adsorption-free SEC devices with the manufacturing costs well below those of commercially available platforms. Developed computational approaches have the predictive power for optimizing reaction conditions and the geometry of the SEC cell.

Self-Assembled Monolayers of Molecular Conductors with Terpyridine-Metal Redox Switching Elements: A Combined AFM, STM and Electrochemical Study.

Self-assembled monolayers (SAMs) of terpyridine-based transition metal (ruthenium and osmium) complexes, anchored to gold substrate via tripodal anchoring groups, have been investigated as possible redox switching elements for molecular electronics. An electrochemical study was complemented by atomic force microscopy (AFM) and scanning tunneling microscopy (STM) methods. STM was used for determination of the SAM conductance values, and computation of the attenuation factor ß from tunneling current-distance curves. We have shown that SAMs of Os-tripod molecules contain larger adlayer structures compared with SAMs of Ru-tripod molecules, which are characterized by a large number of almost evenly distributed small islands. Furthermore, upon cyclic voltammetric experimentation, Os-tripod films rearrange to form a smaller number of even larger islands, reminiscent of the Ostwald ripening process. Os-tripod SAMs displayed a higher surface concentration of molecules and lower conductance compared with Ru-tripod SAMs. The attenuation factor of Os-tripod films changed dramatically, upon electrochemical cycling, to a higher value. These observations are in accordance with previously reported electron transfer kinetics studies.

Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA.

Interpretation of 3JP,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of 3JP,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of 3JP,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD 3JP,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of 3JP,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.

The Ad-MD method to calculate NMR shift including effects due to conformational dynamics: The 31 P NMR shift in DNA.

A method for averaging of NMR parameters by molecular dynamics (MD) has been derived from the method of statistical averaging in MD snapshots, benchmarked and applied to structurally dynamic interpretation of the 31 P NMR shift (δ31P ) in DNA phosphates. The method employs adiabatic dependence of an NMR parameter on selected geometric parameter(s) that is weighted by MD-calculated probability distribution(s) for the geometric parameter(s) (Ad-MD method). The usage of Ad-MD for polymers is computationally convenient when one pre-calculated structural dependence of an NMR parameter is employed for all chemically equivalent units differing only in dynamic behavior. The Ad-MD method is benchmarked against the statistical averaging method for δ31P in the model phosphates featuring distinctively different structures and dynamic behavior. The applicability of Ad-MD is illustrated by calculating 31 P NMR spectra in the Dickerson-Drew DNA dodecamer. δ31P was calculated with the B3LYP/IGLO-III/PCM(water) and the probability distributions for the torsion angles adjacent to the phosphorus atoms in the DNA phosphates were calculated using the OL15 force field.

Raman spectroscopy and DFT calculations of PEDOT:PSS in a dipolar field.

The conductive polymer-electrolyte interface plays an important role in many electrochemical devices. An unusual situation arises when a solvent-free ionic liquid (SF-IL) is used as the electrolyte because it behaves as a molten salt rather than an electrolyte solution. On the basis of Raman spectra, it was found that the presence of ion pairs of SF-IL in the vicinity of poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) results in a decrease in the oxidation level of the polymer and an increase in the HOMO-LUMO gap. The process of polymer "dedoping" and the modification of the electronic structure of the polymer are illustrated by quantum chemical calculations.

Environmental Control of Single-Molecule Junction Evolution and Conductance: A Case Study of Expanded Pyridinium Wiring.

Environmental control of single-molecule junction evolution and conductance was demonstrated for expanded pyridinium molecules by scanning tunneling microscopy break junction method and interpreted by quantum transport calculations including solvent molecules explicitly. Fully extended and highly conducting molecular junctions prevail in water environment as opposed to short and less conducting junctions formed in non-solvating mesitylene. A theoretical approach correctly models single-molecule conductance values considering the experimental junction length. Most pronounced difference in the molecular junction formation and conductance was identified for a molecule with the highest stabilization energy on the gold substrate confirming the importance of molecule-electrode interactions. Presented concept of tuning conductance through molecule-electrode interactions in the solvent-driven junctions can be used in the development of new molecular electronic devices.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids.

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

QM and QM/MM umbrella sampling MD study of the formation of Hg(II)-thymine bond: Model for evaluation of the reaction energy profiles in solutions with constant pH.

The formation of the Hg-N3(T) bond between the 1-methylthymine (T) molecule and the hydrated Hg2+ cation was explored with the combined quantum mechanics/molecular mechanics (QM/MM) method including Free Energy Perturbation corrections. The thermodynamic properties were determined in the whole pH range, when these systems were explicitly investigated and considered as the QM part: (1) T + [Hg(H2 O)6 ]2+ , (2) T + [Hg(H2 O)5 (OH)]+ , (3) T + Hg(H2 O)4 (OH)2 , and (4) N3-deprotonated T + Hg(H2 O)4 (OH)2 . The MM part contained only solvent molecules and counterions. As a result, the dependence of Gibbs-Alberty reaction free energy on pH was obtained along the reaction coordinate. We found that an endoergic reaction in acidic condition up to pH < 4-5 becomes exoergic for a higher pH corresponding to neutral and basic solutions. The migration of the Hg2+ cation between N3 and O4/2 positions in dependence on pH is discussed as well. For the verification, DFT calculations of stationary points were performed confirming the qualitative trends of QM/MM MD simulations and NMR parameters were determined for them.

Tuning the contact conductance of anchoring groups in single molecule junctions by molecular design.

A tetraphenylmethane tripod functionalized with three thiol moieties in the para position can serve as a supporting platform for functional molecular electronic elements. A combined experimental scanning tunneling microscopy break junction technique with theoretical approaches based on density functional theory and non-equilibrium Green's function formalism was used for detailed charge transport analysis to find configurations, geometries and charge transport pathways in the molecular junctions of single molecule oligo-1,4-phenylene conductors containing this tripodal anchoring group. The effect of molecular length (n = 1 to 4 repeating phenylene units) on the charge transport properties and junction configurations is addressed. The number of covalent attachments between the electrode and the tripodal platform changes with n affecting the contact conductance of the junction. The longest homologue n = 4 adopts an upright configuration with all three para thiolate moieties of the tripod attached to the gold electrode. The contact conductance of the tetraphenylmethane tripod substituted by thiols in the para position is higher than that substituted in the meta position. Such molecular arrangement is highly conducting and allows well-defined directional positioning of a variety of functional groups.

Probabilistic mapping of single molecule junction configurations as a tool to achieve the desired geometry of asymmetric tripodal molecules.

Four molecules containing identical tripodal anchors and p-oligophenylene molecular wires of increasing length were used to demonstrate tuning of the asymmetric molecular junction to the desired geometry by probabilistic mapping of single molecule junction configurations in a scanning tunnelling microscopy break junction experiment.

The structural model of Zika virus RNA-dependent RNA polymerase in complex with RNA for rational design of novel nucleotide inhibitors.

Zika virus is a global health threat due to significantly elevated risk of fetus malformations in infected pregnant women. Currently, neither an effective therapy nor a prophylactic vaccination is available for clinical use, desperately necessitating novel therapeutics and approaches to obtain them. Here, we present a structural model of the Zika virus RNA-dependent RNA polymerase (ZIKV RdRp) in complex with template and nascent RNAs, Mg2+ ions and accessing nucleoside triphosphate. The model allowed for docking studies aimed at effective pre-screening of potential inhibitors of ZIKV RdRp. Applicability of the structural model for docking studies was illustrated with the NITD008 artificial nucleotide that is known to effectively inhibit the function of the ZIKV RdRp. The ZIKV RdRp - RNA structural model is provided for all possible variations of the nascent RNA bases pairs to enhance its general utility in docking and modelling experiments. The developed model makes the rational design of novel nucleosides and nucleotide analogues feasible and thus provides a solid platform for the development of advanced antiviral therapy.

The benchmark of 31P NMR parameters in phosphate: a case study on structurally constrained and flexible phosphate.

A benchmark for structural interpretation of the 31P NMR shift and the 2JP,C NMR spin-spin coupling in the phosphate group was obtained by means of theoretical calculations and NMR measurements in diethylphosphate (DEP) and 5,5-dimethyl-2-hydroxy-1,3,2-dioxaphosphinane 2-oxide (cDEP). The NMR parameters were calculated employing the B3LYP, BP86, BPW91, M06-2X, PBE0, KT2, KT3, MP2, and HF methods, and the 6-31+G(d), Iglo-n (n = II, III), cc-pVnZ (n = D, T, Q, 5), aug-cc-pVnZ (n = D, T and Q), and pcS-n and pcJ-n (n = 1, 2, 3, 4) bases, including the solvent effects described with explicit water molecules and/or the implicit Polarizable Continuum Model (PCM). The effect of molecular dynamics (MD) on NMR parameters was MD-calculated using the GAFF force field inclusive of explicit hydration with TIP3P water molecules. Both the optimal geometries and the dynamic behaviors of the DEP and cDEP phosphates differed notably, which allowed a reliable theoretical benchmark of the 31P NMR parameters for highly flexible and structurally constrained phosphate in a one-to-one relationship with the corresponding experiment. The calculated 31P NMR shifts were referenced employing three different NMR reference schemes to highlight the effect of the 31P NMR reference on the accuracy of the calculated 31P NMR shift. The relative Δδ(31P) NMR shift calculated employing the MD/B3LYP/Iglo-III/PCM method differed from the experiment by 0.16 ppm while the NMR shifts referenced to H3PO4 and/or PH3 deviated from the experiment notably more, which illustrated the superior applicability of the relative NMR reference scheme. The 2JP,C coupling in DEP and cDEP calculated employing the MD/B3LYP/Iglo-III(DSO,PSO,SD)/cc-PV5Z(FC)/PCM method inclusive of correction due to explicit hydration differed from the experiment by 0.32 Hz and 0.15 Hz, respectively. The NMR calculations demonstrated that reliable structural interpretation of the 31P NMR parameters in phosphate must involve both the structural and the dynamical components.

Theoretical and experimental study of the antifreeze protein AFP752, trehalose and dimethyl sulfoxide cryoprotection mechanism: correlation with cryopreserved cell viability.

In this work the physico-chemical properties of selected cryoprotectants (antifreeze protein TrxA-AFP752, trehalose and dimethyl sulfoxide) were correlated with their impact on the constitution of ice and influence on frozen/thawed cell viability. The freezing processes and states of investigated materials solutions were described and explained from a fundamental point of view using ab-initio modelling (molecular dynamics, DFT), Raman spectroscopy, Differential Scanning Calorimetry and X-Ray Diffraction. For the first time, in this work we correlated the microscopic view (modelling) with the description of the frozen solution states and put these results in the context of human skin fibroblast viability after freezing and thawing. DMSO and AFP had different impacts on their solution's freezing process but in both cases the ice crystallinity size was considerably reduced. DMSO and AFP treatment in different ways improved the viability of frozen/thawed cells.

The mechanism of the glycosylase reaction with hOGG1 base-excision repair enzyme: concerted effect of Lys249 and Asp268 during excision of 8-oxoguanine.

The excision of 8-oxoguanine (oxoG) by the human 8-oxoguanine DNA glycosylase 1 (hOGG1) base-excision repair enzyme was studied by using the QM/MM (M06-2X/6-31G(d,p):OPLS2005) calculation method and nuclear magnetic resonance (NMR) spectroscopy. The calculated glycosylase reaction included excision of the oxoG base, formation of Lys249-ribose enzyme-substrate covalent adduct and formation of a Schiff base. The formation of a Schiff base with ΔG# = 17.7 kcal/mol was the rate-limiting step of the reaction. The excision of the oxoG base with ΔG# = 16.1 kcal/mol proceeded via substitution of the C1΄-N9 N-glycosidic bond with an H-N9 bond where the negative charge on the oxoG base and the positive charge on the ribose were compensated in a concerted manner by NH3+(Lys249) and CO2-(Asp268), respectively. The effect of Asp268 on the oxoG excision was demonstrated with 1H NMR for WT hOGG1 and the hOGG1(D268N) mutant: the excision of oxoG was notably suppressed when Asp268 was mutated to Asn. The loss of the base-excision function was rationalized with QM/MM calculations and Asp268 was confirmed as the electrostatic stabilizer of ribose oxocarbenium through the initial base-excision step of DNA repair. The NMR experiments and QM/MM calculations consistently illustrated the base-excision reaction operated by hOGG1.

Structure Determination of an Ag(I) -Mediated Cytosine-Cytosine Base Pair within DNA Duplex in Solution with (1) H/(15) N/(109) Ag NMR Spectroscopy.

The structure of an Ag(I) -mediated cytosine-cytosine base pair, C-Ag(I) -C, was determined with NMR spectroscopy in solution. The observation of 1-bond (15) N-(109) Ag J-coupling ((1) J((15) N,(109) Ag): 83 and 84 Hz) recorded within the C-Ag(I) -C base pair evidenced the N3-Ag(I) -N3 linkage in C-Ag(I) -C. The triplet resonances of the N4 atoms in C-Ag(I) -C demonstrated that each exocyclic N4 atom exists as an amino group (-NH2 ), and any isomerization and/or N4-Ag(I) bonding can be excluded. The 3D structure of Ag(I) -DNA complex determined with NOEs was classified as a B-form conformation with a notable propeller twist of C-Ag(I) -C (-18.3±3.0°). The (109) Ag NMR chemical shift of C-Ag(I) -C was recorded for cytidine/Ag(I) complex (δ((109) Ag): 442 ppm) to completed full NMR characterization of the metal linkage. The structural interpretation of NMR data with quantum mechanical calculations corroborated the structure of the C-Ag(I) -C base pair.

Benchmark Theoretical and Experimental Study on (15)N NMR Shifts of Oxidatively Damaged Guanine.

The (15)N NMR shifts of 9-ethyl-8-oxoguanine (OG) were calculated and measured in liquid DMSO and in crystal. The OG molecule is a model for oxidatively damaged 2'-deoxyguanosine that occurs owing to oxidative stress in cell. The DNA lesion is repaired with human 8-oxoguanine glycosylase 1 (hOGG1) base-excision repair enzyme, however, the exact mechanism of excision of damaged nucleobase with hOGG1 is currently unknown. This benchmark study on (15)N NMR shifts of OG aims their accurate structural interpretation and calibration of the calculation protocol utilizable in future studies on mechanism of hOGG1 enzyme. The effects of NMR reference, DFT functional, basis set, solvent, structure, and dynamics on calculated (15)N NMR shifts were first evaluated for OG in crystal to calibrate the best performing calculation method. The effect of large-amplitude motions on (15)N NMR shifts of OG in liquid was calculated employing molecular dynamics. The B3LYP method with Iglo-III basis used for B3LYP optimized geometry with 6-311++G(d,p) basis and including effects of solvent and molecular dynamic was the calculation protocol used for calculation of (15)N NMR shifts of OG. The NMR shift of N9 nitrogen of OG was particularly studied because the atom is involved in an N-glycosidic bond that is cleaved with hOGG1. The change of N9 NMR shift owing to oxidation of 9-ethylguanine (G) measured in liquid was -27.1 ppm. The calculated N9 NMR shift of OG deviated from experiment in crystal and in liquid by 0.45 and 0.65 ppm, respectively. The calculated change of N9 NMR shift owing to notable N9-pyramidalization of OG in one previously found polymorph was 20.53 ppm. We therefore assume that the pyramidal geometry of N9 nitrogen that could occur for damaged DNA within hOGG1 catalytic site might be detectable with (15)N NMR spectroscopy. The calculation protocol can be used for accurate structural interpretation of (15)N NMR shifts of oxidatively damaged guanine DNA residue.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions.

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

Structures, physicochemical properties, and applications of T-Hg(II)-T, C-Ag(I)-C, and other metallo-base-pairs.

Recently, metal-mediated base-pairs (metallo-base-pairs) have been studied extensively with the aim of exploring novel base-pairs; their structures, physicochemical properties, and applications have been studied. This trend has become more evident after the discovery of Hg(II)-mediated thymine-thymine (T-Hg(II)-T) and Ag(I)-mediated cytosine-cytosine (C-Ag(I)-C) base-pairs. In this article, we focus on the basic science and applications of these metallo-base-pairs, which are composed of natural bases.

Direct detection of the mercury-nitrogen bond in the thymine-Hg(II)-thymine base-pair with (199)Hg NMR spectroscopy.

We have observed the 1-bond (199)Hg-(15)N J-coupling ((1)J((199)Hg,(15)N) = 1050 Hz) within the Hg(II)-mediated thymine-thymine base pair (T-Hg(II)-T). This strikingly large (1)J((199)Hg,(15)N) is the first one for canonical sp(2)-nitrogen atoms, which can be a sensitive structure-probe of N-mercurated compounds and a direct evidence for N-mercuration.