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.

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.

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.

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.

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.

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.

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.

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.

The structure of metallo-DNA with consecutive thymine-HgII-thymine base pairs explains positive entropy for the metallo base pair formation.

We have determined the three-dimensional (3D) structure of DNA duplex that includes tandem Hg(II)-mediated T-T base pairs (thymine-Hg(II)-thymine, T-Hg(II)-T) with NMR spectroscopy in solution. This is the first 3D structure of metallo-DNA (covalently metallated DNA) composed exclusively of 'NATURAL' bases. The T-Hg(II)-T base pairs whose chemical structure was determined with the (15)N NMR spectroscopy were well accommodated in a B-form double helix, mimicking normal Watson-Crick base pairs. The Hg atoms aligned along DNA helical axis were shielded from the bulk water. The complete dehydration of Hg atoms inside DNA explained the positive reaction entropy (ΔS) for the T-Hg(II)-T base pair formation. The positive ΔS value arises owing to the Hg(II) dehydration, which was approved with the 3D structure. The 3D structure explained extraordinary affinity of thymine towards Hg(II) and revealed arrangement of T-Hg(II)-T base pairs in metallo-DNA.

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.

Formation of a thymine-Hg(II)-thymine metal-mediated DNA base pair: proposal and theoretical calculation of the reaction pathway.

A reaction mechanism that describes the substitution of two imino protons in a thymine:thymine (T:T) mismatched DNA base pair with a Hg(II) ion, which results in the formation of a (T)N3-Hg(II)-N3(T) metal-mediated base pair was proposed and calculated. The mechanism assumes two key steps: The formation of the first Hg(II)-N3(T) bond is triggered by deprotonation of the imino N3 atom in thymine with a hydroxo ligand on the Hg(II) ion. The formation of the second Hg(II)-N3(T) bond proceeds through water-assisted tautomerization of the remaining, metal-nonbonded thymine base or through thymine deprotonation with a hydroxo ligand of the Hg(II) ion already coordinated to the thymine base. The thermodynamic parameters ΔGR =-9.5 kcal mol(-1), ΔHR =-4.7 kcal mol(-1), and ΔSR =16.0 cal mol(-1) K(-1) calculated with the ONIOM (B3LYP:BP86) method for the reaction agreed well with the isothermal titration calorimetric (ITC) measurements by Torigoe et al. [H. Torigoe, A. Ono, T. Kozasa, Chem. Eur. J. 2010, 16, 13218-13225]. The peculiar positive reaction entropy measured previously was due to both dehydration of the metal and the change in chemical bonding. The mercury reactant in the theoretical model contained one hydroxo ligand in accord with the experimental pKa  value of 3.6 known for an aqua ligand of a Hg(II) center. The chemical modification of T:T mismatched to the T-Hg(II)-T metal-mediated base pair was modeled for the middle base pair within a trinucleotide B-DNA duplex, which ensured complete dehydration of the Hg(II) ion during the reaction.

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.

Charge carrier mobility in poly[methyl(phenyl)silylene] studied by time-resolved terahertz spectroscopy and molecular modelling.

Time-resolved terahertz spectroscopy and combination of quantum chemistry modeling and molecular dynamics simulations were used for the determination of charge carrier mobility in poly[methyl(phenyl)silylene]. Using time-resolved THz spectroscopy we established the on-chain charge carrier drift mobility in PMPSi as 0.02 cm(2) V(-1) s(-1). This value is low due to the formation of polarons: the hole is self-trapped in a potential formed by local chain distortion and the transient THz conductivity spectra show signatures of its oscillations within this potential well. This view is supported by the agreement between experimental and calculated values of the on-chain charge carrier mobility.

Charge transfer in porphyrin-calixarene complexes: ultrafast kinetics, cyclic voltammetry, and DFT calculations.

Transient absorption spectroscopy, cyclic voltammetry, and DFT calculations were used to describe charge transfer processes in a series of 5,10,15,20-tetrakis(N-methylpyridinium-n-yl) porphyrins (TMPyPn, n = 4,3,2) and TMPyPn/p-sulfonatocalix[m]arene (clxm, m = 4,6,8) complexes. Excitation of TMPyPn is accompanied by an increasing electron density at the methylpyridinium substituents in the order TMPyP2 < TMPyP3 < TMPyP4. The quenching of the excited singlet states of the complexes increases with the number of ionized phenolic groups of clxm and can be correlated with the partial transfer of the electron density from O(-) to the peripheral methylpyridinium substituents rather than to the porphyrin ring.

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.

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.

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.

Femtosecond fluorescence and intersystem crossing in rhenium(I) carbonyl-bipyridine complexes.

Ultrafast electronic-vibrational relaxation upon excitation of the singlet charge-transfer b (1)A' state of [Re(L)(CO) 3(bpy)] ( n ) (L = Cl, Br, I, n = 0; L = 4-Et-pyridine, n = 1+) in acetonitrile was investigated using the femtosecond fluorescence up-conversion technique with polychromatic detection. In addition, energies, characters, and molecular structures of the emitting states were calculated by TD-DFT. The luminescence is characterized by a broad fluorescence band at very short times, and evolves to the steady-state phosphorescence spectrum from the a (3)A" state at longer times. The analysis of the data allows us to identify three spectral components. The first two are characterized by decay times tau 1 = 85-150 fs and tau 2 = 340-1200 fs, depending on L, and are identified as fluorescence from the initially excited singlet state and phosphorescence from a higher triplet state (b (3)A"), respectively. The third component corresponds to the long-lived phosphorescence from the lowest a (3)A" state. In addition, it is found that the fluorescence decay time (tau 1) corresponds to the intersystem crossing (ISC) time to the two emissive triplet states. tau 2 corresponds to internal conversion among triplet states. DFT results show that ISC involves electron exchange in orthogonal, largely Re-localized, molecular orbitals, whereby the total electron momentum is conserved. Surprisingly, the measured ISC rates scale inversely with the spin-orbit coupling constant of the ligand L, but we find a clear correlation between the ISC times and the vibrational periods of the Re-L mode, suggesting that the latter may mediate the ISC in a strongly nonadiabatic regime.