Metabolism of 32-oxo-24,25-dihydrolanosterols by partially purified cytochrome P-450(14DM) from rat liver microsomes

Metabolism of 32-oxo-24,25-dihydrolanosterols (3 beta-hydroxylanost-8-en-32- al (4,delta 8-CHO) and 3 beta-hydroxylanost-7-en-32-al (5,delta 7-CHO)) was studied in a reconstituted system consisting of rat liver partially purified cytochrome P-450, which catalyzes lanosterol 14-demethylation (P-450(14DM)), and reduced nicotineamide adenine dinucleotide phosphate (NADPH)-cytochrome P-450 reductase. The reconstituted system converted delta 8-CHO to 4,4-dimethyl-5 alpha-cholesta-8,14-dien-3 beta-ol (2, 8, 14-Diene), which corresponds to the 14-deformylated product. delta 7-CHO, the isomer of delta 8-CHO, was not converted to the corresponding 14-deformylated product. The apparent Km value of cytochrome P-450(14DM) for delta 8-CHO was about 1/20 of that for 24,25-dihydrolanosterol (1, DHL). The metabolism of delta 8-CHO was inhibited by 7-oxo-24,25-dihydrolanosterol (6, 7-oxo-DHL), which is a potent inhibitor of cholesterol biosynthesis from lanosterol or DHL. However, the metabolism of delta 8-CHO was less inhibited by 7-oxo-DHL than that of DHL.     

Production of meiosis-activating sterols from metabolically engineered yeast

Meiosis-activating sterols (MAS), a class of potent regulators of reproductive processes, are difficult to obtain by chemical synthesis or isolation from natural sources. We demonstrate the development of metabolically engineered strains of Saccharomyces cerevisiae that accumulate MAS as the predominant sterol product. Homologous recombination was used to construct an erg24Delta erg25Delta hem1Delta mutant RXY4.3, which lacked sterol Delta14 reductase, C-4 oxidase, and delta-aminolevulinate synthase. The HEM1 deletion allowed sterol import and rendered RXY4.3 viable under aerobic conditions. This mutant accumulated 4,4-dimethyl-5alpha-cholesta-8,14,24-trien-3beta-ol (FF-MAS), and a similar erg25Delta hem1Delta mutant produced 4,4-dimethyl-5alpha-cholesta-8,24-dien-3beta-ol (T-MAS). Based on consistent yields of approximately 5 mug of FF-MAS per mL of culture, fermentation of genetically modified yeast compares favorably with other approaches to produce MAS.     

Conformational comparison of five follicular fluid meiosis-activating sterol-related active and inactive compounds

The crystal structures of five follicular fluid meiosis-activating sterol-related Δ^8,14-sterol compounds are presented. These are 4,4-di­methyl-23-phenyl-24-nor-5α-chola-8,14-dien-3β-ol, C₃₁H₄₄O, 4,4-di­methyl-22-phenyl-23,24-dinor-5α-chola-8,14-dien-3β-ol, C₃₀H₄₂O, (20R)-4,4-di­methyl-22-oxa-5α,20-chol­es­ta-8,14,24-trien-3β-ol, C₂₈H₄₄O₂, 4,4-di­methyl-23-phenyl-22-oxa-24-nor-5α-chola-8,14-dien-3β-ol–water (4/1), 4C₃₀H₄₂O₂·H₂O, and 4,4-di­methyl-5α-cholesta-8,14-dien-3-one, C₂₉H₄₆O. Two of the derivatives are inactive and three are active as agonists. Preliminary structure–activity relationship studies showed that the positions of the double bonds in the skeleton and the structures of the side chains are important determinants for activity. The conformations of the skeletons were compared with double-bond isomers retrieved from the Cambridge Structural Database [Allen & Kennard (1993). Chem. Des. Autom. News, 8 , 1, 31–37]; no significant differences were found. Thus, conformational changes induced by the double bonds are not discriminative with respect to the activity of the compounds. Comparisons of the side-chain conformations of active and inactive structures revealed that the crystal structures were not conclusive as far as correlation of conformation and activity of the side chains were concerned.     

Studies on taxol biosynthesis. Preparation of taxa-4(20),11(12)-dien-5 alpha-acetoxy-10 beta-ol by deoxygenation of a taxadiene tetraacetate obtained from Japanese yew

Taxa-4(20),11(12)-dien-5 alpha-acetoxy-10 beta-ol 6 has been identified as an early stage intermediate involved in the biosynthesis of taxol (Paclitaxel). This compound has been efficiently prepared by Barton deoxygenation of the C-2- and C-14-hydroxyl groups on a derivative semisynthetically prepared from taxa-4(20),11(12)-dien-2 alpha,5 alpha,10 beta-triacetoxy-14 beta-(2-methyl)butyrate (7), a major taxoid metabolite isolated from Japanese yew heart wood. The synthetic methodology is amenable for the preparation of isotopically labeled congeners that will be useful to probe further intermediate steps in the biosynthesis of taxol.     

Asporyergosterol, a new steroid from an algicolous isolate of Aspergillus oryzae

Asporyergosterol (1), a new steroid with an E double bond between C-17 and C-20, was identified from the culture extracts of Aspergillus oryzae, an endophytic fungus isolated from the marine red alga Heterosiphonia japonica. Moreover, four known steroids including (22E,24R)-ergosta-4,6,8(14),22-tetraen-3-one (2), (22E,24R)-3beta-hydroxyergosta-5,8,22-trien-7-one (3), (22E,24R)-ergosta-7,22-dien-3beta,5alpha,6beta-triol (4), and (22E,24R)-5alpha,8alpha-epidioxyergosta-6,22-dien-3beta-ol (5) were isolated. Structures of these compounds were unambiguously established by spectroscopic techniques and by comparison with literature values. All the isolates exhibited low activity to modulate acetylcholinesterase (AChE).     

Five-membered 2,3-dioxo heterocycles: XCII. Reaction of 3-aroyl-1H-pyrrolo[2,1-c][1,4]benzoxazine-1,2,4-triones with ethyl (2Z)-(3,3-dimethyl-8-oxo-2-azaspiro[4.5]deca-6,9-dien-1-ylidene)ethanoate. Synthesis of bridged analogs of pyrrolizidine alkaloids

3-Aroyl-1H-pyrrolo[2,1-c][1,4]benzoxazine-1,2,4-triones reacted with ethyl (2Z)-(3,3-dimethyl-8-oxo-2-azaspiro[4.5]deca-6,9-dien-1-ylidene)acetate to give ethyl 6′-aryl-2′-(2-hydroxyphenyl)-11′,11′-dimethyl-3′,4,4′,13′-tetraoxospiro[2,5-cyclohexadiene-1,9′-(7′-oxa-2′,12′-diazatetracyclo[6.5.1.0^1,5.0^8,12]tetradec-5′-ene)]-14′-carboxylates whose structure was confirmed by X-ray analysis. The products may be regarded as bridged analogs of pyrrolizidine alkaloids, 7′-oxa-2′,12′-diazatetracyclo[6.5.1.01,5.08,12]tetradecanes.     

Ozonolysis of alkenes and study of reactions of polyfunctional compounds: LXVIII.* synthesis of ω-carboxy derivatives of 20-hydroxyecdysone diacetonide

Ozonolysis of ω-anhydro-20-hydroxyecdysone diacetonide gave a mixture of 24- and 25-oxo derivatives, and only the first of these (23-carbaldehyde) reacted with malonic acid according to Knoevenagel to give 14α-hydroxy-2β,3β: 20,22-bis(isopropylidenedioxy)-6-oxo-27-nor-5β-cholesta-7,24-dien-26-oic acid. The oxidation of 23-carbaldehyde with ozone, followed by treatment with diazomethane, afforded 20-hydroxy-25,26,27-trinorecdysone-23-carboxylic acid methyl ester diacetonide.     

High-performance thin-layer chromatographic method for quantitative determination of 4alpha-methyl-24beta-ethyl-5alpha-cholesta-14,25-dien-3beta-ol, 24beta-ethylcholesta-5,9(11),22E-trien-3beta-ol, and betulinic acid in Clerodendrum inerme

A sensitive, selective, precise, and robust high-performance thin-layer chromatography method was developed and validated for analysis of two new recently isolated sterols, 4alpha-methyl-24beta-ethyl-5alpha-cholesta-14,25-dien-3beta-ol (1) and 24beta-ethylcholesta-5,9(11),22E-trien-3beta-ol (2), and a triterpene, betulinic acid (3), in Clerodendrum inerme extract. The method employed HPTLC plates precoated with silica gel 60F(254 )as the stationary phase. To achieve good separation, an optimised mobile phase consisting of toluene-acetone (94:06, v/v) was used (R(f )0.48, 0.34, and 0.22 for compounds 1, 2, and 3, respectively). Densitometric determination of the above compounds was carried out in reflection/absorption mode at 620 nm. Optimised chromatographic conditions provide well separated compact spots for the compounds 1, 2, and 3. The calibration curves were linear in the concentration range of 100-2500 ng/spot. The method was validated for precision, robustness, and recovery. The limits of detection and quantitation were 5, 6, and 10 microg/mL and 14, 18, and 29 microg/mL, respectively, for 1, 2, and 3. The method reported here is reproducible and convenient for quantitative analysis of these compounds in the aerial parts of C. inerme.     

Quantitation of lanosterol and its major metabolite FF-MAS in an inhibition assay of CYP51 by azoles with atmospheric pressure photoionization based LC-MS/MS

Azoles affect the steroid balance in all biological systems and may therefore be called endocrine disrupters. Lanosterol 14alpha-demethylase (CYP51) is an enzyme inhibited by azoles. Only few data have been reported showing their inhibitory potency since an assay in an in vitro system is not available so far. In the present work an inhibition assay using human recombinant CYP51, coexpressed with human P450 oxido-reductase by the baculovirus/insect cell expression system, and LC-MS/MS as analytical method is described. Atmospheric pressure photoionization (APPI) and atmospheric pressure chemical ionization (APCI) sources were used with a triple quadrupole mass spectrometer to compare quantitation of lanosterol (substrate) and 4,4-dimethyl-5alpha-cholesta-8,14,24-triene-3beta-ol (FF-MAS) (product of CYP51) with d(6)-2,2,3,4,4,6-cholesterol (d(6)-cholesterol) as internal standard. Optimization of analytical parameters resulted in a LC-APPI-MS/MS method with a LOQ of 10 pg on column for FF-MAS. The sensitivity of the method (LOD 0.5 ng/ml) makes it possible to analyze supernatants of inhibition experiments after precipitation of proteins by isopropanol without any sample enrichment. The coefficient of variation of the analytical method was <20% (n = 5) for FF-MAS, lanosterol and d(6)-cholesterol. The external calibration curve was linear from 1 to 10,000 ng/ml with R(2) >/= 0.999 and an accuracy of 94-115%. Compared with APCI, APPI provides a ten- to 500-fold increase in sensitivity for the analytes in this study. IC(50) values of epoxiconazole and miconazole-two widely used azole fungicides used in agriculture and in human medicine, respectively-were 1.95 microM and 0.057 microM.     

The total synthesis and anti-fertility activity of 6-silasteroids

4,4-Dimethyl-6-methoxy-4-sila-1-tetralone (2) was prepared by a modified literature procedure and converted to 3-methoxy-6,6-dimethyl-6-silaestra-1,3,5(10),8,14-pentaen-17β-yl acetate (5c). Catalytic hydrogenation of 5c gave 3-methoxy-6,6-dimethyl-6-silaestra-1,3,5(10),8-tetraen-17β-yl acetate (6b), and its 14-iso- and Δ^{1,3,5(10),8(14)} isomers, the proportions varying with the catalyst and solvent. Reduction of 6b with lithium-liquid ammonia, and O-demethylation, gave 6,6-dimethyl-6-silaestradiol (8b). Reduction of the 3-methyl ether of 8b with lithium-liquid ammonia-t-butanol and hydrolysis afforded 3-keto-6,6-dimethyl-6-silaestr-1(10)-en-17β-ol (15), which was catalytically reduced to its 1,10α-dihydro derivative 17. The 5,6 SiC bond of 8b, 15 and their derivatives was cleaved by boron tribromide, aq. ethanolic hydrogen fluoride, and other reagents, providing a series of 5,6-seco-6,6-dimethyl-6-silasteroids. X-ray crystallographic analysis of 17 and the 17α-ethynyl derivative of 15 confirmed the stereochemical assignments. None of the compounds which were subjected to uterotropic, anti-uterotropic, or post-coital assays, showed significant activity. A partially completed synthesis of 6-silaestradiol (21a) is described.     

The mass spectra of organic compounds. 7th Communication. 4,4-Dimethyl-1-pentene

Loss of CH, CH₄, C₂H₄, C₃H, C₃H₆ and C₃H₇ from the molecular ions of a number of ¹³C-labeled analogs of 4,4-dimethyl-1-pentene was studied both in normal (source) 70-eV electron impact (EI) spectra dn in metastable spectra. For loss of CH in the source, 96% of the methyl comes frm positions of 5, 5′ and 5″, while the remainder comes from position 1. In the metastable spectra, loss of C-1 (16%) and C-3 (9%) is increasing in importance. The loss of ethylene is a particular case: either C-1 or C-3 are lost with any other C-atom from positions 2,5,5′, and 5″ (8 × 10%) in the metastable spectra, the probability for simultaneous loss of C-1 and C-3 being 6%. If C-1 seems to these two positions become completely equivalent in the metastable time range. The T-values (kinetic energy release) for the different positions show small, but statisticaly different values and a small isotope effect. Loss of C₃H₅ (allylic cleavage) is 100% C-1, C-2 and C-3, i.e., no evidence for skeletal rearrangement is seen. This is also true for loss of C₃C₆ (McLafferty rearrangement) within the source, but in metastable decay the other positions gain in importance. The neutral fragment C₃H appears to be the the result of consecutive loss of CH and C₃H₄, rather than a one-step loss of propyl radical or the inverse reactions sequence. No metastable reaction can be seen for this reaction. Decomposition of labeled C₆H and C₅H secondary ions occurs in an essentially random fashion.     

Synthesis and stereochemistry of highly crowded N-benzylpiperidones

A series of N-benzylated 3,5-diakyl-2,6-diarylpiperidin-4-ones 4–8 were conveniently synthesized in significant yields of 68–88% by N-benzylation of the corresponding 2,6-diaryl-3,5-dimethylpiperidin-4-ones 1–3 using different benzyl bromides. Initially, the new piperidone 2,6-bis(4-ethoxyphenyl)-3,5-dimethylpiperidin-4-one 3 was synthesized by the condensation of 1:1:2 M ratio of 3-pentanone, ammonium acetate and para-ethoxybenzaldehyde in ethanolic medium. All the synthesized new compounds 3–8 were characterized by their analytical and spectral (IR, ¹H and ¹³C NMR) interpretations. The stereochemistry of the new piperidone 3 was elucidated as chair conformation with an equatorial orientation of all substituents, suggested by its vicinal couplings from ¹H NMR spectrum. To investigate the impact on piperidone stereochemistry as well as NMR chemical shifts, all the N-benzylated products 4–8 were compared with their corresponding precursors, and as a result, it is clearly established that all the synthesized N-benzyl piperidones exist in the chair conformation with an equatorial orientation of all the substituents at C-2, C-3, C-5, C-6 and N. Contrary to the probability all N-benzylated compounds retain the same conformation and configuration as their precursors, however, a remarkable change on the chemical shifts are observed. For the further unambiguous confirmation of stereochemistry, the 1-benzyl-3,5-dimethyl-2,6-diphenylpiperidin-4-one 4 was examined by single-crystal X-ray diffraction. The compound 4, C₂₆H₂₇NO, crystallized in a P-1 space group under triclinic system with unit cell dimensions a, b, c (Å) and α, β, γ (°) of 10.156(2), 11.002(2), 11.348(4) and 116.74(4), 100.81(3), 100.17(3), respectively.     

Crown ethers with converging neutral binding sites: Synthesis and complexation with t-butylammonium hexafluorophosphate

The influence of substituents in close proximity to crown ether cavities, on the stability of complexes of the crown ethers with t-butylammonium salts, has been investigated. Crown ethers with intra-annular donor substituents (2–4) were prepared by the reaction of 2-acetylresorcinol (1) with polyethylene glycol ditosylates and subsequent modification of the acetyl group. Crown ethers with substituents above and below the plane of the crown ether 0 atoms were synthesized by the reaction of 2,2'-dihydroxy-1,1'-biphenyls with polyethylene glycol ditosylates. Chloromethylation of 5,5'-dimethyl-1,1'-biphenyl crown ethers (6) yielded 4,4'-bis(chloromethyl)-1,1'-biphenyl crown ethers (10). 3,3'-Disubstituted-1,1'-biphenyl crown ethers (13–24) were synthesized by the reaction of 3,3'-diallyl-2,2'-dihydroxy-1,1'-biphenyl (12) with polyethylene glycol ditosylates. The allyl groups of 13 were isomerized with sodium hydride to propen- 1-yl groups. Ozonolysis of 13 and 14 gave the corresponding dialdehydes (15 and 18) which were converted into other 3,3'-disubstituted biphenyl-20-crown-6 derivatives (RCH₂COOMe, CH₂COOH, CH₂OH, CH₂Cl, CH₂OMe, OH and Me) by standard operations. The thermodynamic stability of the complexes of these functionalized crown ethers with t-butylammonium hexafluorophosphate has been studied in deuterochloroform in competition experiments with m-xyleno-18-crown-5 and benzo-15-crown-5 as the reference compounds. The nature of the 2-substituents in the crown ethers 2 and 3 has little effect on the stability of the complexes. The stability of the complexes of 3,3'-disubstituted biphenyl crown ethers depends of ringsize and the size and nature of the substituents. The most stable complexes are those of 24 (R = Me) and 14 (R=CH=CHMe).The Me groups in 24 represent the optimum between relief of O-O repulsion in the polyether ring and steric hindrance of complexation. The propen-1-yl substituents of 14 stabilize the complex because they provide extended π-electron donor stabilization. Substitution at the 4- and 4'-positions of the aryl groups has little effect on the stability of the complexes.     

Combined biotransformations of 4(20),11-taxadienes

Taxuyunnanine C (1) and its analogs (2 and 3), the C-14 oxygenated 4(20), 11-taxadienes from callus cultures of Taxus sp., were regio- and stereo-selectively hydroxylated at the 7β position by a fungus, Abisidia coerulea IFO 4011, and it was interesting that the longer the alkyl chain of the acyloxyl group at C-14 became, the higher the yield of 7β-hydroxylated product was. Besides the three 7β-hydroxylated products (5, 9, 17), other nine new products (7, 11, 12, 14, 15, 16, 18, 20 and 21) and six known products (4, 6, 8, 10, 13 and 19) were obtained. Subsequently, the acetylated derivatives (24 and 27) of 7β-and 9α-hydroxylated products of 1 were regio- and stereo-specifically hydroxylated at the 9α position by Ginkgo cells and 7β position by A. coerulea, respectively. Thus, the two specific oxidations have been combined. These bioconversions would provide not only valuable intermediates for the semi-synthesis of paclitaxel or other bioactive taxoids from 1 and its analogs, but also some useful hints for the biosynthetic pathway of taxoid in the natural Taxus plant.     

Reactions of Sulfenamides with Compounds of Trivalent Phosphorus

Chloro-dimethylamino-phenyl-p-tolylthio-phosphonium chloride 2 , dimethylamino-diphenyl-p-tolylthio-phosphonium chloride 3 , bis(diethylamino)-dimethylamino-p-tolylthiophosphonium chloride 4 and tert-butyl-dimethylamino-phenyl-p-tolylthio-phosphonium chloride 5 were prepared by the reaction of N,N-dimethylamino-p-tolylsulfenamide 1 with PhPCl₂, Ph₂PCl, (Et₂N)₂PCl and ^tBu(Ph)PCl, respectively. The reaction of N,N′-dimethyl-N,N′-bis(trimethylsilyl)urea 9 and N-methyl-N′-phenyl-N,N′-bis(trimethylsilyl)urea 10 with phenylsulfenyl chloride 6 or p-nitrophenylsulfenyl chloride 8 furnished the N-arylthio-N,N′-diorgano-N′-(trimethylsilyl)-ureas 11 – 14 . The reaction of 11 – 14 and of the previously known compounds 15 and 16 with MePCl₂, ClCH₂PCl₂, ^tBuPCl₂ and PhPCl₂ resulted in the formation of the 2-arylthio-2-chloro-1,2,3-triorgano-1,3,2λ⁵-diazaphosphetidin-4-ones 17 – 26 . 1,3-Dimethyl-2-(1,1,1,3,3,3-hexafluoro-2-propoxy)-2-phenyl-2-phenylthio-1,3,2λ⁵-diazaphosphetidin-4-one 29 and the 2-arylthio-1,3-dimethyl-2-(p-nitrophenoxy)-2-organyl-1,3,2λ⁵-diazaphosphetidin-4-ones 30 – 32 were obtained in the reactions of compounds 17, 24 and 27 with 1,1,1,3,3,3-hexafluoro-2-propanol or p-nitrophenol in the presence of triethylamine. The reaction of compound 21 with thiophenol in the presence of triethylamine resulted in a mixture of products, from which 1,3,4,5,7-pentamethyl-1,3,5,7-tetraaza-4λ⁵-phosphaspiro[3.3]heptan-2,6-dione 33 was isolated. The identity and structure of all the new compounds were established by ¹H-, ¹³C- and ³¹P-NMR spectroscopy and by elemental analysis. A possible mechanism of reaction of sulfenamides with compounds of trivalent phosphorus is discussed. For the compounds 5a, 32 and 33 X-ray structure analyses were conducted. The cation of compound 5a involves four-coordinate phosphorus (essentially tetrahedral geometry) and is a rare example of a P–S single bond in such a system (P–S 207.37(9) pm). In 32 the geometry at phosphorus is distorted trigonal bipyramidal, with axial positions occupied by oxygen and nitrogen atoms. In the spirophosphorane 33 the geometry at phosphorus is intermediate between trigonal bipyramidal and square pyramidal, with essentially planar four-membered rings.     

Linear homobimetallic palladium complexes

The oxidative addition of C₆H₄-1,4-I₂ (1) to Pd(PPh₃)₄ (2) gives mononuclear trans-(Ph₃P)₂Pd(C₆H₄-4-I)(I) (3), which can be converted to trans-(Ph₃P)₂Pd(C₆H₄-4-I)(OTf) (5) by its reaction with [AgOTf] (4). Complex 5 can be used in the high-yield preparation of a series of unique cationic mono- and dinuclear palladium complexes of structural type [trans-(Ph₃P)₂Pd(C₆H₄-4-I)(L)]⁺ (7, L = C₄H₄N₂; 9a, L = C₅H₄N-4-CN; 9b, L = NC-4-C₅H₄N) and [trans-(C₆H₄-4-I)(Ph₃P)₂Pd ← N^∩N → Pd(PPh₃)₂(C₆H₄-4-I)]²⁺ (14a, N^∩N = C₆H₄-1,4-(CN)₂; 14b, N^∩N = (C₆H₄-4-CN)₂; 14c, N^∩N = 4,4′-bipyridine (=bipy)). Complexes 7, 9 and 14 rearrange in solution to give [trans-(Ph₃P)₂Pd(C₆H₄-4-PPh₃)(L)]²⁺ (10, L = C₄H₄N₂; 12a, L = C₅H₄N-4-CN; 12b, L = NC-4-C₅H₄N) and [trans-(C₆H₄-4-PPh₃)(Ph₃P)₂Pd ← N^∩N → Pd(PPh₃)₂(C₆H₄-4-PPh₃)]⁴⁺ (15a, N^∩N = C₆H₄-1,4-(CN)₂; 15b, N^∩N = (C₆H₄-4-CN)₂) along with {[(Ph₃P)₂(Ph₃P-4-C₆H₄)Pd(μ-I)]₂}²⁺ (11).The solid state structures of 3, 9a, 10, 11 and 15b are reported. Most characteristic for all complexes is the square-planar coordination geometry of palladium with trans-positioned PPh₃ ligands. In 3 the iodide and the 4-iodo-benzene are linear oriented laying with the palladium atom on a crystallographic C₂ axes. In 9a this symmetry is broken by steric interactions of the PPh₃ ligands with the 4-cyanopyridine and 4-iodobenzene groups. Compound 11 contains two μ-bridging iodides with different Pd-I separations showing that the ligand induces a stronger trans-influence than PPh₃. In 15b, the Ph₃PC₆H₄Pd ← NCC₆H₄C₆H₄CN → PdC₆H₄PPh₃ building block is rigid-rod structured with the C₆H₄ units perpendicular oriented to the Pd coordination plane, while the biphenylene connecting moiety is in-plane bound.     

Synthesis, structural characterization, and initial electroluminescent properties of bis-cycloiridiated complexes of 2-(3,5-bis(trifluoromethyl)phenyl)-4-methylpyridine

A series of bis-cyclometalated Ir(III) complexes (8-10, 12, 15, 17, 19, 21, 23, 25, 28, 29 and 33) bearing two chromophoric N^∧C cyclometalated ligands derived from 2-(3,5-bis(trifluoromethyl)phenyl)-4-methylpyridine (1) and a third nonchromophoric ligand has been synthesized. A palladium-catalyzed cross-coupling reaction between 2-chloro-4-methylpyridine (2) and 3,5-bis(trifluoromethyl)phenylboronic acid (3) was used to prepare 2-(3,5-bis(trifluoromethyl)phenyl)-4-methylpyridine (1). Cyclometalation of (1) by IrCl₃ was carried out in (MeO)₃PO, with the formation of chloro-bridged dimer [N^∧C]₂Ir(μ-Cl)₂Ir[C^∧N]₂ (8). Reaction of (8) with lithium 2,4-pentanedionate, lithium 2,2,6,6-tetramethyl-heptane-3,5-dionate (13), dipivaloyltrimethylsilylphosphine (14), 2,2-dimethyl-6,6,7,7,8,8,8-heptafluoro-3,5-octadione (16), 1,1,1,3,3,3-hexafluoro-2-pyridin-2-yl-propan-2-ol (18), 1,1,1,3,3,3-hexafluoro-2-pyrazol-1-ylmethyl-propan-2-ol (20), 2-diphenylphosphanylethanol (22), and 1-diphenylphosphanylpropan-2-ol (24), afforded octahedral iridium complexes 9, 12, 15, 17, 19, 21, 23 and 25, respectively. Complex 10, which contains three different ligands (L₁ = N^∧C of 1; L₂ = N^∧C of 4,4′-dimethyl-[2,2′]bipyridinyl 4; L₃ = O^∧O of 2,4-pentanedione), and complex 11, which contains no cyclometalated ligands (L₁ = 4; L₂ = L₃ = Cl; L₄ = O^∧O of 2,4-pentanedione) were also isolated as minor products in a one-pot reaction between a 94:5 mixture of 1 and 4, IrCl₃ and lithium 2,4-pentanedionate. Reaction of 8 with diphenylphosphanylmethanol (27) in 1,2-dichloroethane unexpectedly led to complexes 28 and 29. The reactions of 8 with benzoylformic acid resulted in the formation of hydroxyl-bridged dimer [N^∧C]₂Ir(μ-OH)₂Ir[C^∧N]₂ (33). According to X-ray analyses, Ir-to-Ir distances in the crystal cell increase from 6.86 Å for 10 to 13.31 Å for 33. The angle theta, which represents the twisting of two cyclometalated C-Ir-N planes relative to each other, varies from 97.5° for 21 to 90.76 for complex 28. OLED devices were fabricated from several Ir complexes and preliminary results are discussed.     

Construction of vicinal 4°/3°-carbons via reductive Cope rearrangement

Herein reported is a strategy for constructing vicinal 4°/3° carbons via reductive Cope rearrangement. Substrates have been designed which exhibit Cope rearrangement kinetic barriers of ∼23 kcal mol⁻¹ with isoenergetic favorability (ΔG ∼ 0). These fluxional/shape-shifting molecules can be driven forward by chemoselective reduction to useful polyfunctionalized building blocks.

Herein reported is a strategy for constructing vicinal 4°/3° carbons via reductive Cope rearrangement.

Constructing sterically congested vicinal quaternary–tertiary carbons (4°/3° carbons) via Cope rearrangement is currently quite limited with only a handful of papers on the subject published over the past 40 years. This stands in stark contrast to the plethora of other methods for establishing sterically congested vicinal carbons.^1–5 Central to the challenge are kinetic and thermodynamic issues associated with the transformation. In the simplest sense, Cope rearrangements proceed in the direction that results in highest alkene substitution (Fig. 1).^6,7 To forge 4°/3° motifs by Cope rearrangement, additional driving forces must be introduced to reverse the [3,3] directionality and compensate for the energetic penalty associated with the steric and torsional strain of the targeted vicinal 4°/3° motif. With limited reports in all cases, oxy-Cope substrates (Scheme 1, eqn (1)),^8–14 divinylcyclopropanes (Scheme 1, eqn (2)),^15–20 and vinylidenecyclopropane-based 1,5-dienes²¹ (Scheme 1, eqn (3)) have demonstrated favourability for constructing vicinal 4°/3° carbons. Malachowski et al. put forth a series of studies on the construction of quaternary centers via Cope rearrangement driven forward by a conjugation event (Scheme 1, eqn (4)).^22–25 In their work, a single example related to the construction of vicinal 4°/3° centers was disclosed, though kinetic (180 °C) and thermodynamic (equilibrium mixtures) challenges are also observed.²³ And of particular relevance to this work, Wigfield et al. demonstrated that 3,3-dicyano-1,5-dienes with the potential to generate vicinal 4°/3° carbons instead react via an ionic mechanism yielding the less congested products (Scheme 1, eqn (5)).²⁶Open in a separate windowFig. 1Cope equilibrium of 1,1,6-trisubstituted 1,5-dienes.Open in a separate windowScheme 1(A) Cope rearrangements for constructing vicinal 4°/3°-centers (B) this report.Our group has been examining strategies to decrease kinetic barriers and increase the thermodynamic favourability of 3,3-dicyano-1,5-diene-based Cope substrates.^27–31 Beyond the simplest, unsubstituted variants, this class of 1,5-diene is not particularly reactive in both a kinetic and thermodynamic sense (e.g.Scheme 1, eqn (5)).^26,32 Reactivity issues aside, these substrates are attractive building blocks for two main reasons: (1) they have straightforward accessibility from alkylidenemalononitriles and allylic electrophiles by deconjugative allylic alkylation.³³ (2) The 1,5-diene termini are substantially different (malononitrile vs. simple alkene) thus allowing for orthogonal functional group interconversion facilitating target and analogue synthesis.³⁴ Herein we report that a combination of 1,5-diene structural engineering^28,31 and reductive conditions (the reductive Cope rearrangement^29,30) can result in the synthesis of building blocks containing vicinal gem-dimethyl 4°/3° carbons along with orthogonal malononitrile and styrene functional groups for interconversion (Scheme 1B). On this line, malononitrile can be directly converted to amides³⁴ yielding functionally dense β-gem-dimethylamides, important pharmaceutical scaffolds.³⁵This project began during the Covid-19 pandemic lockdown (ca. March–May 2020). As such, we were not permitted to use our laboratory out of an abundance of caution. We took this opportunity to first computationally investigate a Cope rearrangement that could result in vicinal 4°/3° carbons (Scheme 2). Then, when permitted to safely return to the lab, we would experimentally validate our findings (vide infra). From our previous work, it is known that by adding either a 4-aromatic group²⁸ or a 4-methyl group³¹ to a 3,3-dicyano-1,5-diene, low barrier (rt – 80 °C) diastereoselective Cope rearrangements can occur. Notably, the 4-substituent was found to destabilize the starting material (weaken the C3–C4 bond, conformationally bias the substrate for [3,3]), and stabilize the product side of the equilibrium via resonance (phenyl group) or hyperconjugation (methyl group). In this study, we modelled substrates 1, 3, and 5 that have variable 4-substitution and would result in vicinal gem-dimethyl- and phenyl-containing 4°/3° carbons upon Cope rearrangement to 2, 4, or 6, respectively. We chose to target this motif due to likely synthetic accessibility from simple starting materials but also because of the important and profound impact that gem-dimethyl groups impart on pharmaceuticals.³⁵ Substrate 1 lacking 4-substitution had an extremely unfavourable kinetic and thermodynamic profile (ΔG^‡ = 31.6; ΔG = +5.3 kcal mol⁻¹). When a 4-methyl group was added, the kinetic barrier (ΔG^‡) dropped appreciably to 28.2 kcal mol; however, the thermodynamics were still quite endergonic (ΔG = +4.4 kcal mol⁻¹). Most excitingly, it was uncovered that the 4-phenyl group dramatically impacted the kinetics and thermodynamics: the [3,3] has a barrier of 22.9 kcal mol⁻¹ (ΔG^‡) and is ∼isoenergetic (ΔG = +0.17 kcal mol⁻¹). Thus, the reaction appears to be fluxional/shape-shifting at room temperature.^36–40 For this substrate, we also modelled the dissociative pathway (Scheme 2D). It was found that bond breakage to two allylic radical intermediates is a higher energy process than the concerted transition state (Scheme 2Cvs.Scheme 2D). Specifically, the dissociative pathway was found to be kinetically less favourable (ΔG^‡ ∼ 27.6 kcal mol; ΔG = 26.2 kcal mol⁻¹) than the concerted process (ΔG^‡ = 22.9 kcal mol⁻¹). While the dissociative pathway is less favourable than the concerted transformation, we surmised that the two-step process becomes accessible at elevated temperature (vide infra). Finally, the ionic pathway was calculated to be significantly higher for this substrate (see the ESI^†).Open in a separate windowScheme 2Computational analysis of 3,3-dicyano-1,5-diene that in theory could result in vicinal 4°/3° carbons. (A) 4-Unsubstituted 3,3-dicyano-1,5-diene. (B) 4-Methyl 3,3-dicyano-1,5-diene. (C) 4-Phenyl 3,3-dicyano-1,5-diene. (D) The dissociative mechanism for substrate 5 is higher than the closed transition state. (E) visualization of the kinetic- and thermodynamic differences of transformations (A–D).The class of substrate uncovered from our computational investigation could be accessed from γ,γ-dimethyl-alkylidenemalononitrile (7a) and 1,3-diarylallyl electrophiles (such as 8a) by Pd-catalyzed deconjugative allylic alkylation (Scheme 3A).³³ As such, model 1,5-diene 5a was prepared to verify the computational results. It was found that upon synthesis of 5a, an inseparable 21 : 79 mixture of 1,5-diene 5a and the 1,5-diene 6a was observed. The predicted ratio of 5a to 6a was 57 : 43 (Scheme 2C). These two results are within the error of the calculations (predicted; slightly endergonic, observed; slightly exergonic). To determine whether the transformation was progressing through the predicted concerted pathway (Scheme 2C) over the dissociative pathway (Scheme 2D), substrate 5b was prepared by an analogous deconjugative allylic alkylation reaction. Similarly, two Cope equilibrium isomers 5b and 6b are observed at room temperature in a 12 : 88 ratio. Upon heating at 100 °C for 3 h, the 1,5-dienes “scramble” (e.g. iso-6b is observed; 0.2 : 1.0 : 1.5 ratio of 5b : 6b : iso-6b) indicating that the dissociative pathway is only accessible at elevated temperature. This is all in good agreement with the calculated kinetics and thermodynamics of this system (Scheme 2).Open in a separate windowScheme 3(A) Observation of fluxional [3,3] and confirmation of calculated predictions. (B) Optimization of a reductive Cope rearrangement protocol for constructing vicinal 4°/3° centers. (C) The Pd-catalyzed deconjugative allylic alkylation must be regioselective.With respect to the synthetic methodology, we aimed to increase the overall efficiency and applicability of the sequence (Scheme 3B). Specifically, we wanted to avoid [3,3] equilibrium mixtures and sensitive/unstable substates and intermediates. It was found that the direct coupling of 7a with diphenylallyl alcohol 9a could take place in the presence of DMAP, Ac₂O, and Pd(PPh₃)₄. When the coupling was complete, methanol and NaBH₄ were added to drive the Cope equilibrium forward, yielding the reduced Cope rearrangement product 10a in 76% isolated yield. In terms of practicality and efficiency, this method utilizes diphenylallyl alcohols, which are more stable and synthetically accessible than their respective acetates, and the [3,3] equilibrium mixture can be directly converted dynamically to a single reduced product.With an efficient protocol in hand for constructing malononitrile–styrene-tethered building blocks featuring central vicinal 4°/3° carbons, we next examined the scope of the transformation (Scheme 4). We chose diarylallyl alcohols with the propensity to react regioselectively via an electronic bias (Scheme 3C).^41,42 The combination of p-nitrophenyl and phenyl (10b) or p-methoxyphenyl (10c) yielded regioselective outcomes with the electron-deficient arene at the allylic position. This is consistent with the expected regiochemical outcome where the nucleophile reacts preferentially at the α-position and the electrophile reacts at the allylic position bearing the donor-arene (Scheme 3C).^41,42 Then, reductive Cope rearrangement occurs to position the electron-deficient arene adjacent to the gem-dimethyl quaternary center. This is an exciting outcome as many pharmaceutically relevant (hetero)arenes are electron deficient. Thus, fluorinated arenes were installed at the allylic position of products 10d–10k. While the phenyl group resulted in poor regioselectivity (1 : 1–3 : 1), the p-methoxyphenyl group enhanced the regiomeric ratios in all cases (3 : 1–15 : 1). The degree of selectivity is correlated with the number and position of fluorine atoms. N-Heterocycles could be incorporated with excellent regioselectivity, generally speaking (10l–10q). For example, 3-chloro-4-pyridyl (10l/10m) groups were installed at the allylic position with >20 : 1 rr. 4-Chloro-3-pyridyl was poorly regioselective (10n), but the combination of 4-trifluomethyl-3-pyridyl/p-methoxyphenyl (10o) gave good regioselectivity of 11 : 1. 2-Pyridyl/p-methoxyphenyl (10q) was also a regioselective combination. We also examined a few other heterocycles including quinoline (10s) and thiazole (10t and 10u) with excellent and modest regioselectivity observed, respectively. As a general trend, when the arenes on the allylic electrophile become less polarized, poor regioselectivity is observed in the Pd-catalyzed allylic alkylation. For example, the combination of p-chlorophenyl and p-methoxyphenyl (10v) or phenyl (10w) yields regioisomeric mixtures of products. This can be circumvented by utilizing symmetric electrophiles (to 10x).Open in a separate windowScheme 4Scope of the 4°/3°-center-generating reductive Cope rearrangement.The phenyl or the p-methoxyphenyl group is necessary to achieve the 4°/3° carbon-generating Cope rearrangement: it functions as an “activator” by lowering the kinetic barrier and increasing thermodynamic favourability. These activating groups can be removed through alkene C C cleavage reactions (e.g. metathesis (Scheme 5) and ozonolysis (Scheme 6B)). In this regard, highly substituted cycloheptenes 11 were prepared by allylation and metathesis (Scheme 4).^28,43 The yields were modest to excellent over this two-step sequence. In many cases, where 10 exists as a mixture of regioisomers, the major allylation/RCM products 11 could be chromatographically separated from their minor constituents. As shown in Scheme 6A, the malononitrile can be transformed via oxidative amidation³⁴ to products 12 containing a dense array of pharmaceutically relevant functionalities (amides, gem-dimethyl, fluoroaromatics, and heteroaromatics). Following this transformation, ozonolysis terminated with a NaBH₄ quench installs an alcohol moiety on small molecule 13a.Open in a separate windowScheme 5Removal of the “activating group” by ring-closing metathesis.Open in a separate windowScheme 6(A) oxidative amidation of malononitrile. (B) Removal of “activating group” by ozonolysis.These first computational and experimental studies utilizing 3,3-dicyano-1,5-dienes as substrates for constructing vicinal 4°/3° centers sets the stage for much further examination and application. For example, while we focused our efforts on gem-dimethyl-based quaternary carbons, it is likely that other functionality can be installed at this position. For example, while unoptimized, it appears the protocol is reasonably effective at incorporating a piperidine moiety in addition to heteroarenes from the allylic electrophile (7b + 9f → 14a; Scheme 7A). Similar functional group interconversion chemistry as described in Schemes 5 and ​and66 can thus yield functionally dense building blocks 15 and 16 in good yields.Open in a separate windowScheme 7(A) The construction of 4/3° centres on piperidines. (B) Promoting endergonic [3,3] rearrangements is possible, assuming the [3,3] kinetic barrier is sufficiently low.While the 4,6-diaryl-3,3-dicyano-1,5-dienes offered the most attractive energetic profile (low kinetic barrier, isoenergetic [3,3] equillibrium; Scheme 2C), the 4-methyl analogue is also intriguing to consider as a viable substrate class for reductive Cope rearrangement (Scheme 2B). The challenge here is that the kinetics and thermodynamics are quite unfavourable (not observable by NMR), but potentially not prohibitively so. It is extremely exciting to find that Cope equilibria that are significantly endergonic in the desired, forward direction (e.g.3a to 4a) can be promoted by a related reductive protocol (Scheme 7B). While unoptimized, we were able to isolate product 17 in xx% yield by heating at 90 °C in the presence of Hantzsch ester in DMF.     

Bare and ligand protected planar hexacoordinate silicon in SiSb3M3+ (M = Ca,Sr, Ba) clusters

The occurrence of planar hexacoordination is very rare in main group elements. We report here a class of clusters containing a planar hexacoordinate silicon (phSi) atom with the formula SiSb₃M₃⁺ (M = Ca, Sr, Ba), which have D_3h (¹A₁′) symmetry in their global minimum structure. The unique ability of heavier alkaline-earth atoms to use their vacant d atomic orbitals in bonding effectively stabilizes the peripheral ring and is responsible for covalent interaction with the Si center. Although the interaction between Si and Sb is significantly stronger than the Si–M one, sizable stabilization energies (−27.4 to −35.4 kcal mol⁻¹) also originated from the combined electrostatic and covalent attraction between Si and M centers. The lighter homologues, SiE₃M₃⁺ (E = N, P, As; M = Ca, Sr, Ba) clusters, also possess similar D_3h symmetric structures as the global minima. However, the repulsive electrostatic interaction between Si and M dominates over covalent attraction making the Si–M contacts repulsive in nature. Most interestingly, the planarity of the phSi core and the attractive nature of all the six contacts of phSi are maintained in N-heterocyclic carbene (NHC) and benzene (Bz) bound SiSb₃M₃(NHC)₆⁺ and SiSb₃M₃(Bz)₆⁺ (M = Ca, Sr, Ba) complexes. Therefore, bare and ligand-protected SiSb₃M₃⁺ clusters are suitable candidates for gas-phase detection and large-scale synthesis, respectively.

The global minimum of SiSb₃M₃⁺ (M = Ca, Sr, Ba) is a D_3h symmetric structure containing an elusive planar hexacoordinate silicon (phSi) atom. Most importantly, the phSi core remains intact in ligand protected environment as well.

Exploring the bonding capacity of main-group elements (such as carbon or silicon) beyond the traditional tetrahedral concept has been a fascinating subject in chemistry for five decades. The 1970 pioneering work of Hoffmann and coworkers¹ initiated the field of planar tetracoordinate carbons (ptCs), or more generally, planar hypercoordinate carbons. The past 50 years have witnessed the design and characterization of an array of ptC and planar pentacoordinate carbon (ppC) species.^2–14 However, it turned out to be rather challenging to go beyond ptC and ppC systems. The celebrated CB₆²⁻ cluster and relevant species^15,16 were merely model systems because C avoids planar hypercoordination in such systems.^17,18 In 2012, the first genuine global minimum D_3h CO₃Li₃⁺ cluster was reported to have six interactions with carbon in planar form, although electrostatic repulsion between positively charged phC and Li centers and the absence of any significant orbital interaction between them make this hexacoordinate assignment questionable.¹⁹ It was only very recently that a series of planar hexacoordinate carbon (phC) species, CE₃M₃⁺ (E = S–Te; M = Li–Cs), were designed computationally by the groups of Tiznado and Merino (Fig. 1; left panel),²⁰ in which there exist pure electrostatic interactions between the negative C^δ− center and positive M^δ+ ligands. These phC clusters were achieved following the so-called “proper polarization of ligand” strategy.Open in a separate windowFig. 1The pictorial depiction of previously reported phC CE₃M₃⁺ (E = S–Te; M = Li–Cs) clusters and the present SiE₃M₃⁺ (E = S–Te and N–Sb; M = Li–Cs and Ca–Ba) clusters. Herein the solid and dashed lines represent covalent and ionic bonding, respectively. The opposite double arrows illustrate electrostatic repulsion.The concept of planar hypercoordinate carbons has been naturally extended to their next heavier congener, silicon-based systems. Although the steric repulsion between ligands decreases due to the larger size, the strength of π- and σ-bonding between the central atom and peripheral ligands dramatically decreases, which is crucial for stability. Planar tetracoordinate silicon (ptSi) was first experimentally observed in a pentaatomic C_2v SiAl₄⁻ cluster by Wang and coworkers in 2000.²¹ Very recently, this topic got a huge boost by the room-temperature, large-scale syntheses of complexes containing a ptSi unit.²² A recent computational study also predicted the global minimum of SiMg₄Y⁻ (Y = In, Tl) and SiMg₃In₂ to have unprecendented planar pentacoordinate Si (ppSi) units.²³ Planar hexacoordinate Si (phSi) systems seem to be even more difficult to stabilize. Previously, a C_2v symmetric Cu₆H₆Si cluster was predicted as the true minimum,²⁴ albeit its potential energy surface was not fully explored. A kinetically viable phSi SiAl₃Mg₃H₂⁺ cluster cation was also predicted.²⁵ However, these phSi systems^24,25 are only local minima and not likely to be observed experimentally. In 2018, the group of Chen identified the Ca₄Si₂²⁻ building block containing a ppSi center and constructed an infinite CaSi monolayer, which is essentially a two-dimensional lattice of the Ca₄Si₂ motif.²⁶ Thus, it is still an open question to achieve a phSi atom to date.Herein we have tried to find the correct combination towards a phSi system as the most stable isomer. Gratifyingly, we found a series of clusters, SiE₃M₃⁺ (E = N, P, As, Sb; M = Ca, Sr, Ba), having planar D_3h symmetry with Si at the center of the six membered ring, as true global minimum forms. Si–E bonds are very strong in all the clusters, and alkaline-earth metals interact with the Si center by employing their d orbitals. However, electrostatic repulsion originated from the positively charged Si and M centers for E = N, P, and As dominates over attractive covalent interaction, making individual Si–M contacts repulsive in nature. This makes the assignment of SiE₃M₃⁺ (E = N, P, As; M = Ca, Sr, Ba) as genuine phSi somewhat skeptical. SiSb₃M₃⁺ (M = Ca, Sr, Ba) clusters are the sole candidates which possess genuine phSi centers as both electrostatic and covalent interactions in Si–M bonds are attractive. The d orbitals of M ligands play a crucial role in stabilizing the ligand framework and forming covalent bonds with phSi. Such planar hypercoordinate atoms are, in general, susceptible to external perturbations. However, the present title clusters maintain the planarity and the attractive nature of the bonds even after multiple ligand binding at M centers in SiSb₃M₃(NHC)₆⁺ and SiSb₃M₃(Bz)₆⁺. This would open the door for large-scale synthesis of phSi as well.Two major computational efforts were made before reaching our title phSi clusters. The first one is to examine SiE₃M₃⁺ (E = S–Po; M = Li–Cs) clusters, which adopt D_3h or C_3v structures as true minima (see Table S1 in ESI^†), being isoelectronic to the previous phC CE₃M₃⁺ (E = S–Po; M = Li–Cs) clusters. In the SiE₃M₃⁺ (E = S–Po; M = Li–Cs) clusters, the Si center always carries a positive charge ranging from 0.01 to +1.03|e|, in contrast to the corresponding phC species (see Fig. 1). Thus, electrostatic interactions between the Si^δ+ and M^δ+ centers would be repulsive (Fig. 1). Given that the possibility of covalent interaction with an alkali metal is minimal, it would be a matter of debate whether they could be called true coordination. A second effort is to tune the electronegativity difference between Si and M centers so that the covalent contribution in Si–M bonding becomes substantial. Along this line, we consider the combinations of SiE₃M₃⁺ (E = N, P, As, Sb; M = Be, Mg, Ca, Sr, Ba). The results in Fig. S1^† show that for E = Be and Mg, the phSi geometry has a large out-of-plane imaginary frequency mode, which indicates a size mismatch between the Si center and peripheral E₃M₃ (E = N–Bi; M = Be, Mg) ring. On the other hand, the use of larger M = Ca, Sr, Ba atoms effectively expands the size of the cavity and eventually leads to perfect planar geometry with Si atoms at the center as minima. In the case of SiBi₃M₃⁺, the planar isomer possesses a small imaginary frequency for M = Ca. Although planar SiBi₃Sr₃⁺ and SiBi₃Ba₃⁺ are true minima, they are 2.2 and 2.5 kcal mol⁻¹ higher in energy than the lowest energy isomer, respectively (Fig. S2^†). Fig. 2 displays some selected low-lying isomers of SiE₃M₃⁺ (E = N, P, As, Sb; M = Ca, Sr, Ba) clusters (see Fig. S3–S6^† for additional isomers). The global minimum structure is a D_3h symmetric phSi with an ¹A₁′ electronic state for all the twelve cases. The second lowest energy isomer, a ppSi, is located more than 49 kcal mol⁻¹ above phSi for E = N. This relative energy between the most stable and nearest energy isomer gradually decreases upon moving from N to Sb. In the case of SiSb₃M₃⁺ clusters, the second-lowest energy isomer is 4.6–6.1 kcal mol⁻¹ higher in energy than phSi. The nearest triplet state isomer is very high in energy (by 36–53 kcal mol⁻¹, Fig. S3–S6^†) with respect to the global minimum.Open in a separate windowFig. 2The structures of low-lying isomers of SiE₃M₃⁺ (E = N, P, As, Sb; M = Ca, Sr, Ba) clusters. Relative energies (in kcal mol⁻¹) are shown at the single-point CCSD(T)/def2-TZVP//PBE0/def2-TZVP level, followed by a zero-energy correction at PBE0. The values from left to right refer to Ca, Sr, and Ba in sequence. The group symmetries and electronic states are also given.Born–Oppenheimer molecular dynamics (BOMD) simulations at room temperature (298 K), taking SiE₃Ca₃⁺ clusters as case studies, were also performed. The results are displayed in Fig. S7.^† All trajectories show no isomerization or other structural alterations during the simulation time, as indicated by the small root mean square deviation (RMSD) values. The BOMD data suggest that the global minimum also has reasonable kinetic stability against isomerization and decomposition.The bond distances, natural atomic charges, and bond indices for SiE₃Ca₃⁺ clusters are given in † for M = Sr, Ba). The Si–E bond distances are shorter than the typical Si–E single bond distance computed using the self-consistent covalent radii proposed by Pyykkö.²⁷ In contrast, the Si–M bond distance is almost equal to the single bond distance. This gives the first hint of the presence of covalent bonding therein. However, the Wiberg bond indices (WBIs) for the Si–M links are surprisingly low (0.02–0.04). We then checked the Mayer bond order (MBO), which can be seen as a generalization of WBIs and is more acceptable since the approach of WBI calculations assumes orthonormal conditions of basis functions while the MBO considers an overlap matrix. The MBO values for the Si–M links are now sizable (0.13–0.18). These values are reasonable considering the large difference in electronegativity between Si and M, and, therefore, only a very polar bond is expected between them. In fact, the calculations of WBIs after orthogonalization of basis functions by the Löwdin method gives significantly large bond orders (0.48–0.55), which is known to overestimate the bond orders somewhat. The above results indicate that the presence of covalent bonding cannot be ruled out only by looking at WBI values.Bond distances (r, in Å), different bond orders (WBIs) {MBOs} [WBI in orthogonalized basis], and natural atomic charges (q, in |e|) of SiE₃Ca₃⁺ (E = N, P, As, Sb) clusters at the PBE0/def2-TZVP level