Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes ₃ B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes ₃ B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes ₃ B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes ₃ B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

Di- and polynuclear silver(I) saccharinate complexes of tertiary diphosphane ligands: synthesis, structures, in vitro DNA binding, and antibacterial and anticancer properties

A series of new silver(I) saccharinate (sac) complexes, [Ag₂(sac)₂(μ-dppm)H₂O]·H₂O (1), {[Ag₂(μ-sac)₂(μ-dppe)]·3H₂O·CH₂Cl₂} _n (2), [Ag₂(μ-sac)₂(μ-dppp)] _n (3), and [Ag(sac)(μ-dppb)] _n (4) [dppm is 1,1-bis(diphenylphosphino)methane, dppe is 1,2-bis(diphenylphosphino)ethane, dppp is 1,3-bis(diphenylphosphino)propane, and dppb is 1,4-bis(diphenylphosphino)butane], have been synthesized and characterized by C, H, N elemental analysis, IR spectroscopy, ¹H NMR, ¹³C NMR, and ³¹P NMR spectroscopy, electrospray ionization mass spectrometry, and thermogravimetry–differential thermal analysis. Single-crystal X-ray studies show that the diphosphanes act as bridging ligands to yield a dinuclear complex (1) and one-dimensional coordination polymers (2 and 4), whereas the sac ligand adopts a μ₂-N/O bridging mode in 2, and is N-coordinated in 1 and 4. The interaction of the silver(I) complexes with fish sperm DNA was investigated using UV–vis spectroscopy, fluorescence spectroscopy, and agarose gel electrophoresis. The binding studies indicate that the silver(I) complexes can interact with fish sperm DNA through intercalation, and complexes 1 and 3 have the highest binding affinity. The gel electrophoresis assay further confirms the binding of the complexes with the pBR322 plasmid DNA. The minimum inhibitory concentrations of the complexes indicate that complex 1 exhibits very high antibacterial activity against standard bacterial strains of Escherichia coli, Salmonella typhimurium, and Staphylococcus aureus, being much higher than those of AgNO₃, silver sulfadiazine, ciprofloxacin, and gentamicin. Moreover, complexes 1–3 exhibit very high cytotoxic activity against A549 and MCF-7 cancer cell lines, compared with AgNO₃ and cisplatin. The bacterial and cell growth inhibitions of the silver(I) complexes are closely related to their DNA binding affinities.     

New names,combinations, and lectotypifications in Heterotheca (Compositae: Astereae)

The following new names and combinations are proposed:Heterotheca barbata (Rydb.) Semple,H. horrida subsp.cinerascens (S. F. Blake) Semple,H. fulcrata vararizonica Semple,H. fulcrata var.senilis (Wooton & Standley) Semple,H. oregona var.compacta (Keck) Semple,H. oregona var.rudis (Greene) Semple,H. oregona var.scaberrima (A. Gray) Semple,H. pumila (Greene) Semple,H. villosa var.pedunculata (Greene) V. Harms ex Semple, andH. zionensis Semple. The following chromosome numbers are reported for the first time:H. fulcrata var.arizonica, 2n=9 _II ;H. horrida subsp.cinerascens, 2n=18 _II ;H. pumila, 2n=9 _II ,2n=18 _II ;H. zionensis, 2n=9 _II .     

Protonation,alkylation, addition and redox reactions of 16, 17 and 18 valence electron [Mo(L)(NO)('S4')] complexes [L = SPh,PMe3, NO; 'S4'2– = 1,2-Bis-(2-mercaptophenylthio)ethane(2-)]

In the quest for complexes modelling functional characteristics of metal sulfur oxidoreductases, a series of molybdenum nitrosyl complexes with sulfur-dominated coordination sphere was synthesized. Treatment of the 16, 17 and 18 valence electron (VE) complexes [Mo(L)(NO)('S₄')] (1–3) [L?=?SPh (1), PMe₃ (2), NO (3), 'S₄'^2–?=?1,2-bis-(2-mercaptophenylthio) ethane(2-)] with the Brönsted acid HBF₄ resulted in formation of different types of products. 1 and 3 were reversibly protonated at one thiolate atom of the 'S₄'^2– ligand;2, however, yielded the phosphonium salt [HPMe₃]BF₄ and the dinuclear [Mo(NO)('S₄')]₂. Alkylation of 1, 2 and 3 by Me₃OBF₄ or Et₃OBF₄ uniformly resulted in high yields of [Mo(L)(NO)(R-'S₄')]BF₄ complexes [L?=?SPh: R?=?Me (5), Et (6); L?=?PMe₃: R?=?Me (7); L?=?NO: R?=?Me (8), Et (9)] in which one thiolate atom of the 'S₄'^2– ligand had become alkylated; the NMR spectra of 5, 6, 8 and 9 indicated that only one out of four theoretically possible diastereoisomers had formed. 5 and 6 were characterized also by single-crystal X-ray structure analyses. A comparison of ν(NO) bands and redox potentials (cyclic voltammetry) of parent complexes and alkylated derivatives showed that alkylation leads to a decrease in electron density at the molybdenum center and to a positive shift in redox potentials. The 16 VE complex 1 could be reduced, also chemically, to give the corresponding 17 VE anion [1]^–, and inserted elemental sulfur into the Mo-SPh bond, forming the 18 VE phenylperthio complex [Mo(η²–SSPh)(NO)('S₄')] (11) which, upon reaction with PPh₃, gave SPPh₃ and regenerated the parent complex 1. These results are discussed with regard to the sequence of proton and electron transfer steps occurring in substrate conversions catalyzed by metal sulfur oxidoreductases.     

Theoretical description of halogen bonding – an insight based on the natural orbitals for chemical valence combined with the extended-transition-state method (ETS-NOCV)

In the present study we have characterized the halogen bonding in selected molecules H₃N–ICF₃ (1-NH ₃ ), (PH₃)₂C–ICF₃ (1-CPH ₃ ), C₃H₇Br–(IN₂H₂C₃)₂C₆H₄ (2-Br), H₂–(IN₂H₂C₃)₂C₆H₄ (2-H ₂ ) and Cl–(IC₆F₅)₂C₇H₁₀N₂O₅ (3-Cl), containing from one halogen bond (1-NH ₃ , 1-CPH ₃ ) up to four connections in 3-Cl (the two Cl–HN and two Cl–I), based on recently proposed ETS-NOCV analysis. It was found based on the NOCV-deformation density components that the halogen bonding C–X ^… B (X-halogen atom, B-Lewis base), contains a large degree of covalent contribution (the charge transfer to X ^… B inter-atomic region) supported further by the electron donation from base atom B to the empty σ*(C–X) orbital. Such charge transfers can be of similar importance compared to the electrostatic stabilization. Further, the covalent part of halogen bonding is due to the presence of σ-hole at outer part of halogen atom (X). ETS-NOCV approach allowed to visualize formation of the σ-hole at iodine atom of CF₃I molecule. It has also been demonstrated that strongly electrophilic halogen bond donor, [C₆H₄(C₃H₂N₂I)₂][OTf]₂, can activate chemically inert isopropyl bromide (2-Br) moiety via formation of Br–I bonding and bind the hydrogen molecule (2-H ₂ ). Finally, ETS-NOCV analysis performed for 3-Cl leads to the conclusion that, in terms of the orbital-interaction component, the strength of halogen (Cl–I) bond is roughly three times more important than the hydrogen bonding (Cl–HN).

Theoretical study on aluminum carbide endohedral fullerene-Al4C@C80

The possibility of a new endohedral fullerene with a trapped aluminum carbide cluster, Al₄C @C₈₀-I _h , was theoretical investigated. The geometries and electronic properties of it were investigated using density functional theory methods. The Al₄C unit formally transfers six electrons to the C₈₀ cage which induces stabilization of Al₄C@C₈₀. A favorable binding energy, relatively large HOMO-LUMO gap, electron affinities and ionization potentials suggested the Al₄C@C₈₀ is rather stable. The analysis of vertical ionization potential and vertical electron affinity indicate Al₄C@C₈₀ is a good electron acceptor.

Different non-host resistance responses of two rice subspecies,japonica and indica,to Puccinia striiformis f. sp. tritici

Key message

Japonica and indica have different non-host resistance (NHR) abilities to Puccinia striiformis f. sp. tritici ( Pst ), and hydrogen peroxide (H ₂ O ₂ ) has a positive function in NHR to japonica against Pst.

Abstract

Non-host interactions between Puccinia striiformis f. sp. tritici (Pst) and two rice subspecies were characterized using 23 rice varieties, including 11 japonica and 12 indica. Results showed that the infected fungal structures were easily produced in the leaves of indica, whereas only several substomatal vesicles and primary infection hyphae were observed in the leaves of japonica. This result indicated that indica is less resistant or more susceptible to Pst than japonica. Hydrogen peroxide accumulated in the initial phase of japonica–Pst interaction but not in indica–Pst interaction. A set of reactive oxygen species (ROS)-related genes was also induced in response to Pst infection, suggesting that ROS activation is one of the major mechanisms of non-host resistance of rice to Pst.     

Characterization of bio-related vanadium and zinc complexes containing tetradentate dithiolate-disulfide, -diamine and -amine-amide ligands

The reaction of [VCl₃(PMe₂Ph)₃] _{with HSS′S′SH (where the HS are thiophenolate and the S′ thioether functions, respectively), H2–1, yields [VCl(μ-SS′S′S)]2 (3) with one of the thiolate groups of each of the two ligands in the bridging mode. Reaction of Na2–1 with [VOCl2(thf)2] leads to a polymeric product of composition [VO(SS′S′S)]x (4). The products obtained from the reaction between [VOCl2(thf)2] and NaSNNSNa, Na2–2, (S is thiophenolate, N the amine function) depend on subtle changes in the diamine backbone of this ligand: If the amine functions are linked by -CH2CH2– (2a), the tetranuclear V^IV complex [V(SNNS)μ-O]4 (5) is formed alongside the V^III complex [VCl(SNNS)]. If the backbone is -CH(Me)CH(Me)- (2b), [VO(SNNS)] (7) and the dinuclear, asymmetrically oxo-bridged V^IV complex [{(SNN ^– S)(thf)V}μ-O{V(SNN ^– S)}] (8) are obtained. In 8, one amine of each of the two ligands is deprotonated to the amide group. In either case, the complexation is accompanied by oxidation of the thiolates to disulfides, leading to the generation of teraazatetrathio-cycloeicosanes (6a/b). Compounds 5 and 8·2THF have been structurally characterized by X-ray analyses. The connectivities have further been established for 3·2CH2Cl2 and for 6b, which exhibits the same conformation as formally characterized 6a. The cluster compound 5 is stabilized by an extended intramolecular N-H...O and N-H...S) hydrogen-bonding network. In 7·2THF, one of the THFs of crystallization is hydrogen-bonded to the NH of the penta-coordinated {VO(SNN ^– S)} moiety; further, there is an intramolecular hydrogen bond between one of the thiolates of this tetragonal-pyramidal half of the molecule and the NH of the octahedral {VO(SNN ^– S)thf} half. The generation of the ligand 2b from its precursor compound, the zinc complex [Zn(SNNS)] (9) leads to the structural characterization of 9·CH3OH with a large SZnS bite angle and a strong hydrogen bond between the methanolic OH and one of the thiolate sulfurs. The relevance of these compounds in biological systems is discussed.}     

Driving and retarding forces in a chemical reaction

The reaction force F(ξ) is the negative gradient of the potential energy of a chemical process along the intrinsic reaction coordinate ξ. We extend the rigorous concept of F(ξ) to the “activation strain model” of Bickelhaupt et al., to formulate the “strain” force F _str(ξ) that retards a reaction and the “interaction” force F _int(ξ) that drives it. These are investigated for a group of Diels-Alder cycloadditions. The results fully support the interpretation of the minimum of F(ξ) as defining the beginning of the transition from deformed reactants to eventual products.     

Synthesis of some quinolinyl chalcone analogues and investigation of their anticancer and synergistic anticancer effect with doxorubicin

Two derivatives of 2-(4-acetylanilino)quinolines (IIIa, b) were synthesized as scaffolds for synthesis of open chalcone analogues (Va-f) through Claisen-Schmidt condensation with a set of aromatic aldehydes (IVa-d). Derivatives (Va, b) were further manipulated into cyclic ??,??-unsaturated ketones by Michael-addition of acetylacetone and ethylacetoacetate affording derivatives (VI?CVII). Deethoxycarboxylation of derivatives (VIIa, b) afforded cyclohexenons (VIIIa, b) allowing formation of a mini library of ??,??-unsaturated ketones for screening their anticancer and synergistic anticancer effect with doxorubicin using colon cancer cell line (Caco-2). Two open enones, (Vb) and (Ve), showed significant anticancer activity with IC₅₀ of 5.0 and 2.5 ??M respectively. Only one cyclic enone, (VIa) showed synergistic anticancer activity with doxorubicin at 10 ??M.     

Insights into the structural determinants for selective inhibition of nitric oxide synthase isoforms

Selective inhibition of the nitric oxide synthase isoforms (NOS) is a promising approach for the treatment of various disorders. However, given the high active site conservation among all NOS isoforms, the design of selective inhibitors is a challenging task. Analysis of the X-ray crystal structures of the NOS isoforms complexed with known inhibitors most often gives no clues about the structural determinants behind the selective inhibition since the inhibitors share the same binding conformation. Aimed at a better understanding of the structural factors responsible for selective inhibition of NOS isoforms we have performed MD simulations for iNOS, nNOS and eNOS complexed with N^ω-NO₂-L-Arg (1), and with the aminopyridine derivatives 2 and 3. The slightly better selectivity of 1 for nNOS may be assigned to the presence of extra charge–charge interactions due to its “extended” conformation. While the high affinity of 2 for iNOS can be explained by the formation of an iNOS-specific subpocket upon binding, the lack of affinity for eNOS is associated to a conformational change in Glu363. The strong van der Waals and electrostatic interactions between 3 and the active site of nNOS are most likely responsible for its higher affinity for this isoform. Owing to the elongated and narrow binding pocket of iNOS, the correct positioning of 3 over the heme group is difficult, which may account for its lower affinity toward this isoform. Brought together, our results might help to rationalize the design of selective NOS inhibitors.

Spin-flip reactions of Zr + C2H6 researched by relativistic density functional theory

Density functional theory (DFT) with relativistic corrections of zero-order regular approximation (ZORA) has been applied to explore the reaction mechanisms of ethane dehydrogenation by Zr atom with triplet and singlet spin-states. Among the complicated minimum energy reaction path, the available states involves three transition states (TS), and four stationary states (1) to (4) and one intersystem crossing with spin-flip (marked by ?): ³ Zr + C ₂ H ₆ → ³ Zr-CH ₃ -CH ₃ ( ³ 1) → ³ TS _1/2 → ³ ZrH-CH ₂ -CH ₃ ( ³ 2) → ³ TS _2/3 ? ¹ ZrH₂-CH₂ = CH₂ ( ¹ 3) → ¹ TS _3/4 → ¹ ZrH ₃ -CH = CH ₂ ( ¹ 4). The minimum energy crossing point is determined with the help of the DFT fractional-occupation-number (FON) approach. The spin inversion leads the reaction pathway transferring from the triplet potential energy surface (PES) to the singlet’s accompanying with the activation of the second C-H bond. The overall reaction is calculated to be exothermic by about 231 kJ mol^?1. Frequency and NBO analysis are also applied to confirm with the experimental observed data.

Computational investigations on the electronic and structural properties of the unsaturated silylenoid HP=SiLiF

The structures of unsaturated silylenoid HP=SiLiF were studied by density functional theory at the B3LYP/6-311+G(d,P) level. Four equilibrium structures, the three-membered ring (1), the four-membered ring (2), the “classical” silane (3), and the linear (4) structures, were located. Their energies are in the order of 4?>?3?>?1?>?2. To exploit the stability of HP=SiLiF, the insertions reaction of 2 and HP=Si into C-Cl have been investigated, respectively. The results show that the insertion of HP=Si is more favorable. To compare with the saturated silylenoid, the insertion reaction of H₂SiLiF was also investigated. The calculations indicate that the insertion of HP=SiLiF (2) is more favorable. The unsaturated siylenoid HP=SiLiF has similar reaction characters to saturated silylenoid H₂SiLiF and silylene HP=Si.     

QTL mapping for salt tolerance and domestication-related traits in Vigna marina subsp. oblonga,a halophytic species

Key message

QTL mapping in F ₂ population [ V. luteola × V. marina subsp. oblonga ] revealed that the salt tolerance in V. marina subsp. oblonga is controlled by a single major QTL.

Abstract

The habitats of beach cowpea (Vigna marina) are sandy beaches in tropical and subtropical regions. As a species that grows closest to the sea, it has potential to be a gene source for breeding salt-tolerant crops. We reported here for the first time, quantitative trait loci (QTLs) mapping for salt tolerance in V. marina. A genetic linkage map was constructed from an F₂ population of 120 plants derived from an interspecific cross between V. luteola and V. marina subsp. oblonga. The map comprised 150 SSR markers. The markers were clustered into 11 linkage groups spanning 777.6 cM in length with a mean distance between the adjacent markers of 5.59 cM. The F_2:3 population was evaluated for salt tolerance under hydroponic conditions at the seedling and developmental stages. Segregation analysis indicated that salt tolerance in V. marina is controlled by a few genes. Multiple interval mapping consistently identified one major QTL which can explain about 50 % of phenotypic variance. The flanking markers may facilitate transfer of the salt tolerance allele from V. marina subsp. oblonga into related Vigna crops. The QTL for domestication-related traits from V. marina are also discussed.     

Estimating the absorptive root area in Norway spruce by using the common direct and indirect earth impedance methods

Aims

Estimates of root absorption magnitude are needed for the balanced management of forest ecosystems, but no methods able to work on the whole tree and stand level were available. Modified earth impedance method was developed recently and here it was tested, by comparing the results with those obtained by combination of several classical methods.

Methods

We used direct (soil cores, scanning and microscopy) and indirect (sap flow patterns and modified earth impedance) methods in an attempt to estimate the absorptive root area indexes (RAI) at two sites of about 25 and 40-years-old Norway spruce. We considered the geometric surfaces of all scanned fine roots to be equal to the fine root absorptive area (RAI _scan ). To estimate the potentially physically permeable area of fine roots, we microscopically evaluated the point of secondary xylem appearance and calculated the geometric area of root portions with primary structure (RAI _micro ). We termed the area of electrically conductive root surface as the active (ion) absorptive area (RAI _mei ) and measured its extent by the modified earth impedance (MEI) method.

Results

The highest values for absorptive root areas at the two experimental sites we obtained with the scanning method (RAI _scan was considered to be 100%), followed by the RAI _micro (51%) and RAI _mei (32%). RAImei reached about 2/3 of RAImicro. The surface area of the ectomycorrhizal hyphae was an order of magnitude larger than that of all fine roots, but the MEI did not measure such increase.

Conclusions

We showed that the absorptive root area, indirectly estimated by the MEI, provides consistent results that approach the values obtained for fine roots with a primary structure estimated by traditional direct methods. The similar range of the values for the absorptive root surface area obtained by microscopy and by the MEI method indicates that this method is feasible and that it could be used to determine the extent of active absorptive root surface areas in forests.     

New combinations and one new variety among the genera Indigofera,Robinia, and Tephrosia (Leguminosae)

This report includes: a new variety,Indigofera miniata var.florida, and five new combinations.I. mucronata var.keyensis, Robinia hispida var.kelseyi, Tephrosia angustissima var.corallicola, T. angustissima var.curtissii, andT. × floridana.