Stack bonding in polyaromatic hydrocarbons.

Parallel displacement of π-stacked component molecules enhances the efficiency of organic semiconductors by maximizing interpenetration of the π-densities. Dimers of symmetric polyaromatic hydrocarbons coronene, hexabenzo[bc,de,gh,kl,no,qr]coronene, circumcoronene, kekulene, and circumcircumcoronene are examined using density functional theory from the stack bonding perspective which considers π-stacking interactions in terms of contributions of monomer π-orbital overlap to the character of dimer orbitals. Energetically favored parallel displaced and/or twisted dimer conformations are consistent with patterns of mixing of the monomer molecular orbitals (MOs) that maximize interpenetration of the π densities. The multiple minima found along parallel displacement (PD) coordinates coincide with the formation of dimer MOs formally antibonding between the monomers at the sandwich conformation to bonding at the PD minima. Minima identified with favorable stack bonding are consistent with polymorphs found in large polyaromatic hydrocarbons.

Molecular Dynamics Simulations of Reduced and Oxidized TFIIIA Zinc Fingers Free and Interacting with 5S RNA.

Interactions of zinc finger (ZF) proteins with nucleic acids and proteins play an important role in DNA transcription and repair, biochemical recognition, and protein regulation. The release of Zn2+ through oxidation of cysteine thiolates is associated with disruption of gene expression and DNA repair, preventing tumor growth. Multi-microsecond molecular dynamics (MD) simulations were carried out to examine the effect of Cys oxidation on the ZF456 fragment of transcription factor III A (TFIIIA) and its complex with 5S RNA. In the absence of 5S RNA, the reduced ZF456 peptide undergoes conformational changes in the secondary structure due to the reorientation of the intact ZF domains. Upon oxidation, the individual ZF domains unfold to various degrees, yielding a globular ZF456 peptide with ZF4 and ZF6, responsible for base-specific hydrogen bonds with 5S RNA, losing their ßßα-folds. ZF5, on the other hand, participates in nonspecific interactions through its α-helix that conditionally unravels early in the simulation. In the presence of RNA, oxidation of the ZF456 peptide disrupts the key hydrogen bonding interactions between ZF5/ZF6 and 5S RNA. However, interactions with ZF4 are dependent on the protonation state of His119.

Halogen Bonding Interactions of Polychlorinated Biphenyls and the Potential for Thyroid Disruption.

Polychlorinated biphenyl (PCB) flame retardants are persistent pollutants and inhibit neurodevelopment, particularly in the early stages of life. Halogen bonding (XB) to the iodothyronine deiodinases (Dio) that modulate thyroid hormones (THs) is a potential mechanism for endocrine disruption. Cl⋅⋅⋅Se XB interactions of PCBs with SeMe- , a small model of the Dio active site selenocysteine, are compared with previous results on polybrominated diphenylethers (PBDEs) and THs using density functional theory. PCBs generally display weaker XB interactions compared to PBDEs and THs, consistent with the dependence of XB strength on the size of the halogen (I>Br>Cl). PCBs also do not meet a proposed energy threshold for substrates to undergo dehalogenation, suggesting they may behave as competitive inhibitors of Dio in addition to other mechanisms of endocrine disruption. XB interactions in PCBs are position-dependent, with ortho interactions slightly more favorable than meta and para interactions, suggesting that PCBs may have a greater effect on certain classes of Dio. Flexibility of PCBs around the biphenyl C-C bond is limited by ortho substitutions relative to the biphenyl linkage, which may contribute to the ability to inhibit Dio and other TH-related proteins.

A Halogen Bonding Perspective on Iodothyronine Deiodinase Activity.

Iodothyronine deiodinases (Dios) are involved in the regioselective removal of iodine from thyroid hormones (THs). Deiodination is essential to maintain TH homeostasis, and disruption can have detrimental effects. Halogen bonding (XB) to the selenium of the selenocysteine (Sec) residue in the Dio active site has been proposed to contribute to the mechanism for iodine removal. Polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls (PCBs) are known disruptors of various pathways of the endocrine system. Experimental evidence shows PBDEs and their hydroxylated metabolites (OH-BDEs) can inhibit Dio, while data regarding PCB inhibition are limited. These xenobiotics could inhibit Dio activity by competitively binding to the active site Sec through XB to prevent deiodination. XB interactions calculated using density functional theory (DFT) of THs, PBDEs, and PCBs to a methyl selenolate (MeSe-) arrange XB strengths in the order THs > PBDEs > PCBs in agreement with known XB trends. THs have the lowest energy C-X*-type unoccupied orbitals and overlap with the Se lp donor leads to high donor-acceptor energies and the greatest activation of the C-X bond. The higher energy C-Br* and C-Cl* orbitals similarly result in weaker donor-acceptor complexes and less activation of the C-X bond. Comparison of the I···Se interactions for the TH group suggest that a threshold XB strength may be required for dehalogenation. Only highly brominated PBDEs have binding energies in the same range as THs, suggesting that these compounds may inhibit Dio and undergo debromination. While these small models provide insight on the I···Se XB interaction itself, interactions with other active site residues are governed by regioselective preferences observed in Dios.

Conformation dynamics of cyclic disulfides and selenosulfides in CXXC(U) (X = Gly, Ala) tetrapeptide redox motifs.

Thioredoxin fold proteins often contain a Cys-(Xxx)n -Cys(Sec) or CXn C(U) motif, where the active cysteine (C) or selenocysteine (U) is bridged by X residues, which vary with protein function. The effect of the X residues on the conformation space of the oxidized disulfide and selenosulfide forms of the CXXC(U) motif has been investigated using molecular dynamics (MD) and density functional theory. Multi-microsecond-length MD simulations of the CGGC, CGAC, and CAGC cyclic peptides show that CGGC rings readily exchange between several conformations over the course of the simulation, but steric interactions with the methyl group of Ala limit the conformation space available to the cyclic peptide, especially for CGAC. The potential for the motif to be reduced, as measured by the energy of the lowest unoccupied molecular orbitals, is dependent upon the ring conformation. These results suggest that control of available conformations by the bridging residues and the protein tertiary structure may be important for defining the function of the CXXC motif. Theoretical 77 Se chemical shifts of the selenosulfide moiety are dependent upon the conformation and/or intramolecular Se···O interactions with the backbone carbonyl group of the C-terminal U residue.

Trigger bond analysis of nitroaromatic energetic materials using wiberg bond indices.

The identification of trigger bonds, bonds that break to initiate explosive decomposition, using computational methods could help direct the development of novel, "green" and efficient high energy density materials (HEDMs). Comparing bond densities in energetic materials to reference molecules using Wiberg bond indices (WBIs) provides a relative scale for bond activation (%ΔWBIs) to assign trigger bonds in a set of 63 nitroaromatic conventional energetic molecules. Intramolecular hydrogen bonding interactions enhance contributions of resonance structures that strengthen, or deactivate, the CNO2 trigger bonds and reduce the sensitivity of nitroaniline-based HEDMs. In contrast, unidirectional hydrogen bonding in nitrophenols strengthens the bond to the hydrogen bond acceptor, but the phenol lone pairs repel and activate an adjacent nitro group. Steric effects, electron withdrawing groups and greater nitro dihedral angles also activate the CNO2 trigger bonds. %ΔWBIs indicate that nitro groups within an energetic molecule are not all necessarily equally activated to contribute to initiation. %ΔWBIs generally correlate well with impact sensitivity, especially for HEDMs with intramolecular hydrogen bonding, and are a better measure of trigger bond strength than bond dissociation energies (BDEs). However, the method is less effective for HEDMs with significant secondary effects in the solid state. Assignment of trigger bonds using %ΔWBIs could contribute to understanding the effect of intramolecular interactions on energetic properties. © 2018 Wiley Periodicals, Inc.

Halogen-Bonding Interactions of Polybrominated Diphenyl Ethers and Thyroid Hormone Derivatives: A Potential Mechanism for the Inhibition of Iodothyronine Deiodinase.

Polybrominated diphenyl ethers (PBDEs) increase flame resistance in consumer goods, but these compounds and their hydroxylated derivatives (OH-BDEs) impair normal thyroid function. Halogen bonding (XB) of PBDEs to an active site selenocysteine may prevent iodothyronine deiodinase(Dio)-catalyzed activation/deactivation of thyroid hormone (TH) derivatives. In this study, we compare the strength of the XB interactions of TH derivatives, iodine-based contrast agents and PBDEs/OH-BDEs with a methylselenolate model of the Dio active site using density functional theory calculations. The strength of the XB interaction depends upon the acceptor halide, the position of the halide, the number of ring substituents, and the proximity of hydroxyl groups to the XB site. The weaker Se⋅⋅⋅Br interactions relative to Se⋅⋅⋅I interactions are consistent with a model of competitive inhibition that blocks binding of THs at elevated PBDE/OH-BDE concentrations. XB interactions were generally more favorable at ortho and meta positions and in substrates with more electron-withdrawing substituents. PBDEs/OH-BDEs that mimic the binding behavior of THs, that is, containing ortho and meta bromides and adjacent hydroxyl groups, may be the most effective inhibitors. Highly-brominated PBDEs/OH-BDEs have comparable interaction energies to THs and may undergo debromination. These results may also suggest that XB strength must exceed a threshold value in order for PBDEs/OH-BDEs to undergo nucleophilic attack by Dio.

Redox active motifs in selenoproteins.

Selenoproteins use the rare amino acid selenocysteine (Sec) to act as the first line of defense against oxidants, which are linked to aging, cancer, and neurodegenerative diseases. Many selenoproteins are oxidoreductases in which the reactive Sec is connected to a neighboring Cys and able to form a ring. These Sec-containing redox motifs govern much of the reactivity of selenoproteins. To study their fundamental properties, we have used (77)Se NMR spectroscopy in concert with theoretical calculations to determine the conformational preferences and mobility of representative motifs. This use of (77)Se as a probe enables the direct recording of the properties of Sec as its environment is systematically changed. We find that all motifs have several ring conformations in their oxidized state. These ring structures are most likely stabilized by weak, nonbonding interactions between the selenium and the amide carbon. To examine how the presence of selenium and ring geometric strain governs the motifs' reactivity, we measured the redox potentials of Sec-containing motifs and their corresponding Cys-only variants. The comparisons reveal that for C-terminal motifs the redox potentials increased between 20-25 mV when the selenenylsulfide bond was changed to a disulfide bond. Changes of similar magnitude arose when we varied ring size or the motifs' flanking residues. This suggests that the presence of Sec is not tied to unusually low redox potentials. The unique roles of selenoproteins in human health and their chemical reactivities may therefore not necessarily be explained by lower redox potentials, as has often been claimed.

Metal Ion Capture Mechanism of a Copper Metallochaperone.

A novel cation-π interaction between the bound Cu(+) metal ion and Trp44 in the periplasmic Cu(+)/Ag(+) metallochaperone Escherichia coli CusF protects Cu(+) from the oxidative influence of the periplasm. In a popular model of metal ion transfer, a conformational change in the metal binding loop disrupts the cation-π interaction and moves Trp44 aside to provide access to the occluded metal ion binding site in an "open" conformation. In this study, our molecular dynamics simulations support this putative mechanism of metal ion transfer. We find that the apoprotein undergoes a transition back and forth from the crystallographically observed "closed" state to the hypothesized open conformation over multiple microseconds. In agreement with nuclear magnetic resonance data, our simulations show that similar transitions are prohibited in Cu(+)·CusF, suggesting that the conformational transitions are gated by a metal ion-mediated second-shell hydrogen bond between metal binding residue His36 and Asp37 of the metal binding loop region. Ab initio quantum mechanical calculations indicate that metal ion binding strengthens this interaction significantly, much like what is found in the case of other metalloproteins. The study builds toward a common evolutionary role of metal ion-mediated second-shell hydrogen bonds in metalloprotein structure and function.

Predicting Trigger Bonds in Explosive Materials through Wiberg Bond Index Analysis.

Understanding the explosive decomposition pathways of high-energy-density materials (HEDMs) is important for developing compounds with improved properties. Rapid reaction rates make the detonation mechanisms of HEDMs difficult to understand, so computational tools are used to predict trigger bonds-weak bonds that break, leading to detonation. Wiberg bond indices (WBIs) have been used to compare bond densities in HEDMs to reference molecules to provide a relative scale for the bond strength to predict the activated bonds most likely to break to trigger an explosion. This analysis confirms that X-NO2 (X=N,C,O) bonds are trigger linkages in common HEDMs such as TNT, RDX and PETN, consistent with previous experimental and theoretical studies. Calculations on a small test set of substituted tetrazoles show that the assignment of the trigger bond depends upon the functionality of the material and that the relative weakening of the bond correlates with experimental impact sensitivities.

Effect of Methoxy Substituents on the Activation Barriers of the Glutathione Peroxidase-Like Mechanism of an Aromatic Cyclic Seleninate.

Density functional theory (DFT) models including explicit water molecules have been used to model the redox scavenging mechanism of aromatic cyclic seleninates. Experimental studies have shown that methoxy substitutions affect the rate of scavenging of reactive oxygen species differently depending upon the position. Activities are enhanced in the para position, unaffected in the meta, and decreased in the ortho. DFT calculations show that the activation barrier for the oxidation of the selenenyl sulfide, a proposed key intermediate, is higher for the ortho methoxy derivative than for other positions, consistent with the low experimental conversion rate.

Mechanistic insights into Mg2+-independent prenylation by CloQ from classical molecular mechanics and hybrid quantum mechanics/molecular mechanics molecular dynamics simulations.

Understanding the mechanism of prenyltransferases is important to the design of engineered proteins capable of synthesizing derivatives of naturally occurring therapeutic agents. CloQ is a Mg(2+)-independent aromatic prenyltransferase (APTase) that transfers a dimethylallyl group to 4-hydroxyphenylpyruvate in the biosynthetic pathway for clorobiocin. APTases consist of a common ABBA fold that defines a ß-barrel containing the reaction cavity. Positively charged basic residues line the inside of the ß-barrel of CloQ to activate the pyrophosphate leaving group to replace the function of the Mg(2+) cofactor in other APTases. Classical molecular dynamics simulations of CloQ, its E281G and F68S mutants, and the related NovQ were used to explore the binding of the 4-hydroxyphenylpyruvate (4HPP) and dimethylallyl diphosphate substrates in the reactive cavity and the role of various conserved residues. Hybrid quantum mechanics/molecular mechanics potential of mean force (PMF) calculations show that the effect of the replacement of the Mg(2+) cofactor with basic residues yields a similar activation barrier for prenylation to Mg(2+)-dependent APTases like NphB. The topology of the binding pocket for 4HPP is important for selective prenylation at the ortho position of the ring. Methylation at this position alters the conformation of the substrate for O-prenylation at the phenol group. Further, a two-dimensional PMF scan shows that a "reverse" prenylation product may be a possible target for protein engineering.

Thiol reduction of arsenite and selenite: DFT modeling of the pathways to an as-se bond.

The reactivity of arsenite and selenite with biological thiols plays an important role in the toxicity of these elements. However, toxic effects are eliminated when the species are coadministered, due to the antagonistic relationship between selenium and arsenic. The reduction of arsenous acid and selenious acid by thiol and the formation of an As-Se species have been modeled using density functional theory (DFT) and solvent-assisted proton exchange (SAPE), a microsolvation technique that uses a network of water molecules to mimic the participation of bulk solvent in proton transfer processes. Activation barriers and relative energies were calculated for the stepwise thiol reduction of arsenite to form As(SR)3 and selenious acid to first form a selenotrisulfide (Se(SR)2) and then H2Se. Several pathways were explored for the formation of an As-Se bond: the nucleophilic attack of selenide or selenopersulfide on As(OH)3, (RS)As(OH)2, and (RS)2AsOH to form (RS)2AsSeH. On the basis of the lower activation barrier and bioavailability of (RS)2AsOH, the reaction of H2Se with (RS)2AsOH is deemed the most favorable, consistent with previous experimental studies.

Density functional theory study of the attack of ebselen on a zinc-finger model.

Density functional theory and solvent-assisted proton exchange are used to model the attack of ebselen 1 on a zinc-finger model, an important step in the regulation of zinc signaling by reducible selenium compounds. These calculations show that the formation of a selenosulfide bond from an Se···S intermediate complex between 1 and a Cys2His2 zinc-finger model can occur through a moderate activation barrier that is consistent with experimental observations of the relative rates of Zn(2+) release from zinc-finger transcription factors and metallothionein.

Oxidation of biologically relevant chalcogenones and their Cu(I) complexes: insight into selenium and sulfur antioxidant activity.

Hydroxyl radical damage to DNA causes disease, and sulfur and selenium antioxidant coordination to hydroxyl-radical-generating Cu(+) is one mechanism for their observed DNA damage prevention. To determine how copper binding results in antioxidant activity, biologically relevant selone and thione ligands and Cu(+) complexes of the formula [Tpm*Cu(L)](+) [Tpm* = tris(3,5-dimethylpyrazolyl)methane; L = N,N'-dimethylimidazole selone or thione] were treated with H2O2 and the products analyzed by (1)H, (13)C{(1)H}, and (77)Se{(1)H} NMR spectroscopy, mass spectrometry, and X-ray crystallography. Upon H2O2 treatment, selone and thione binding to Cu(+) prevents oxidation to Cu(2+); instead, the chalcogenone ligand is oxidized. Thus, copper coordination by sulfur and selenium compounds can provide targeted sacrificial antioxidant activity.

Ruthenium(II) dichloride complexes of chiral, tetradentate aminosulfoxide ligands: stereoisomerism and redox-induced linkage isomerism.

Ruthenium(II) dichloride complexes of two chiral tetradentate aminosulfoxide ligands, varying only in the N-N linker, were synthesized. With each ligand, two major isomers formed, and these were structurally assigned and characterized through a combination of NMR and UV-vis spectroscopies, X-ray crystallography, and density functional theory calculations. The cis-ß geometric isomer was formed by each ligand, whereas the trans and cis-α geometric isomers were significant components for one ligand only. Cyclic voltammetry studies show that only the cis-ß isomers undergo linkage isomerism upon oxidation to ruthenium(III), whereas the trans and cis-α isomers show simple reversible redox couples.

Orbital-based insights into parallel-displaced and twisted conformations in π-π interactions.

Dispersion and electrostatics are known to stabilize π-π interactions, but the preference for parallel-displaced (PD) and/or twisted (TW) over sandwiched (S) conformations is not well understood. Orbital interactions are generally believed to play little to no role in π-stacking. However, orbital analysis of the dimers of benzene, pyridine, cytosine and several polyaromatic hydrocarbons demonstrates that PD and/or TW structures convert one or more π-type dimer MOs with out-of-phase or antibonding inter-ring character at the S stack to in-phase or bonding in the PD/TW stack. This change in dimer MO character can be described in terms of a qualitative stack bond order (SBO) defined as the difference between the number of occupied in-phase/bonding and out-of-phase/antibonding inter-ring π-type MOs. The concept of an SBO is introduced here in analogy to the bond order in molecular orbital theory. Thus, whereas the SBO of the S structure is zero, parallel displacement or twisting the stack results in a non-zero SBO and overall bonding character. The shift in bonding/antibonding character found at optimal PD/TW structures maximizes the inter-ring density, as measured by intermolecular Wiberg bond indices (WBIs). Values of WBIs calculated as a function of the parallel-displacement are found to correlate with the dispersion and other contributions to the π-π interaction energy determined by the highly accurate density-fitting DFT symmetry adapted perturbation theory (DF-DFT-SAPT) method. These DF-DFT-SAPT calculations also suggest that the dispersion and other contributions are maximized at the PD conformation rather than the S when conducted on a potential energy curve where the inter-ring distance is optimized at fixed slip distances. From these results of this study, we conclude that descriptions of the qualitative manner in which orbitals interact within π-stacking interactions can supplement high-level calculations of the interaction energy and provide an intuitive tool for applications to crystal design, molecular recognition and other fields where non-covalent interactions are important.

The effect of metalation on antimicrobial piscidins imbedded in normal and oxidized lipid bilayers.

Metalation of the N-terminal Amino Terminal Cu(ii)- and Ni(ii)-binding (ATCUN) motif may enhance the antimicrobial properties of piscidins. Molecular dynamics simulations of free and nickelated piscidins 1 and 3 (P1 and P3) were performed in 3 : 1 POPC/POPG and 2.6 : 1 : 0.4 POPC/POPG/aldo-PC bilayers (POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPG, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol; aldo-PC, 1-palmitoyl-2-(9'-oxo-nonanoyl)-sn-glycero-3-phosphocholine) bilayer models. Nickel(ii) binding decreases the conformation dynamics of the ATCUN motif and lowers the charge of the N-terminus to allow it to embed deeper in the bilayer without significantly changing the overall depth due to interactions of the charged half-helix of the peptide with the headgroups. Phe1⋯Ni2+ cation-π and Phe2-Phe1 CH-π interactions contribute to a small fraction of structures within the nickelated P1 simulations and may partially protect a bound metal from metal-centered chemical activity. The substitution of Phe2 for Ile2 in P3 sterically blocks conformations with cation-π interactions offering less protection to the metal. This difference between metalated P1 and P3 may indicate a mechanism by which peptide sequence can influence antimicrobial properties. Any loss of bilayer integrity due to chain reversal of the oxidized phospholipid chains of aldo-PC may be enhanced in the presence of metalated piscidins.

In silico insights into the dimer structure and deiodinase activity of type III iodothyronine deiodinase from bioinformatics, molecular dynamics simulations, and QM/MM calculations.

The homodimeric family of iodothyronine deiodinases (Dios) regioselectively remove iodine from thyroid hormones. Currently, structural data has only been reported for the monomer of the mus type III thioredoxin (Trx) fold catalytic domain (Dio3Trx), but the mode of dimerization has not yet been determined. Various groups have proposed dimer structures that are similar to the A-type and B-type dimerization modes of peroxiredoxins. Computational methods are used to compare the sequence of Dio3Trx to related proteins known to form A-type and B-type dimers. Sequence analysis and in silico protein-protein docking methods suggest that Dio3Trx is more consistent with proteins that adopt B-type dimerization. Molecular dynamics (MD) simulations of the refined Dio3Trx dimer constructed using the SymmDock and GalaxyRefineComplex databases indicate stable dimer formation along the ß4α3 interface consistent with other Trx fold B-type dimers. Free energy calculations show that the dimer is stabilized by interdimer interactions between the ß-sheets and α-helices. A comparison of MD simulations of the apo and thyroxine-bound dimers suggests that the active site binding pocket is not affected by dimerization. Determination of the transition state for deiodination of thyroxine from the monomer structure using QM/MM methods provides an activation barrier consistent with previous small model DFT studies.Communicated by Ramaswamy H. Sarma.

Modeling the mechanism of the glutathione peroxidase mimic ebselen.

Ebselen (1), the quintessential mimic of the antioxidant selenoenzyme glutathione peroxidase (GPx), is a potential chemopreventative for various diseases associated with oxidative stress. Density-functional theory (DFT) and solvent-assisted proton exchange (SAPE) are used to model the complex mechanism for scavenging of reactive oxygen species by 1. SAPE is a microsolvation method designed to approximate the role of bulk solvent in chemical processes involving proton transfer. Consistent with experimental studies, SAPE studies predict the reaction of 1 with thiol (RSH) to form a selenenyl sulfide 2 to be preferred under most conditions, with an alternate pathway through a selenoxide 3 possible at high reactive oxygen species (ROS) concentrations ([ROS] ≫ [RSH]). The reduction of 2 to the selenol 4, known to be rate-determining in the protein, has a high SAPE activation barrier due to a strong Se···O interaction which reduces the electrophilicity of the sulfur center of the -SeS- bond of 2. Thiols, such as dithiols and peptide-based thiols, are expected to overcome this barrier through structural features that increase the probability of attack at this sulfur. Thus, in vivo, the GPx-like pathway is the most likely mechanism for 1 under most circumstances, except, perhaps, under extreme oxidative stress where initial oxidation to 3 could compete with formation of 2. Simple thiols, used in various in vitro studies, are predicted by SAPE modeling to proceed through oxidation of 2 to a seleninyl sulfide intermediate. Overall, SAPE modeling provides a realistic interpretation of the redox mechanism of 1 and holds promise for further exploration of complex aqueous-phase reaction mechanisms.