Hydrogen-bonding behavior of various conformations of the HNO<Subscript>3</Subscript>…(CH<Subscript>3</Subscript>OH)<Subscript>2</Subscript> ternary system

Nine minima were found on the intermolecular potential energy surface for the ternary system HNO₃(CH₃OH)₂ at the MP2/aug-cc-pVDZ level of theory. The cooperative effect, which is a measure of the hydrogen-bonding strength, was probed in these nine conformations of HNO₃…(CH₃OH)₂. The results are discussed here in terms of structures, energetics, infrared vibrational frequencies, and topological parameters. The cooperative effect was observed to be an important contributor to the total interaction energies of the cyclic conformers of HNO₃…(CH₃OH)₂, meaning that it cannot be neglected in simulations in which the pair-additive potential is applied.

Open image in new window

Graphical abstract The H-bonding behavior of various conformations of the HNO₃(CH₃OH)₂ trimer was investigated

     

Different effects of SNP and GSNO on mitochondrial O<Stack><Subscript>2</Subscript><Superscript>.−</Superscript></Stack>/H<Subscript>2</Subscript>O<Subscript>2</Subscript> production

Sodium Nitroprusside (SNP) and S-Nitrosoglutathione (GSNO) differently affect mitochondrial H₂O₂ release at Complex-I. mM SNP increases while GSNO decreases the release induced by succinate alone or added on top of NAD-linked substrates. Stimulation likely depends on Nitric Oxide ( ^. NO) (released by SNP but not by GSNO) inhibiting cytochrome oxidase and mitochondrial respiration. Preincubations with SNP or high GSNO (10 mM plus DTE to increases its ^. NO release) induces an inhibition of the succinate dependent H₂O₂ production consistent with a ^. NO dependent covalent modification. However maximal inhibition of the succinate dependent H₂O₂ release is obtained in the presence of low GSNO (20–100 μM), but not with SNP. This inhibition appears independent of ^. NO release since μM GSNO does not affect mitochondrial respiration, or the H₂O₂ detection systems and its effect is very rapid. Inhibition may be partly due to an increased removal of O₂^.− since GSNO chemically competes with NBT and cytochrome C in O₂^.− detection.     

Electronic and optical properties of C<Subscript>24</Subscript>, C<Subscript>12</Subscript>X<Subscript>6</Subscript>Y<Subscript>6</Subscript>, and X<Subscript>12</Subscript>Y<Subscript>12</Subscript> (X = B,Al and Y = N,P)

Utilizing first-principles calculations, we studied the electronic and optical properties of C₂₄, C₁₂X₆Y₆, and X₁₂Y₁₂ fullerenes (X?=?B, Al; Y?=?N, P). These fullerenes are energetically stable, as demonstrated by their negative cohesive energies. The energy gap of C₂₄ may be tuned by doping, and the B₁₂N₁₂ fullerene was found to have the largest energy gap. All of the fullerenes had finite optical gaps, suggesting that they are optical semiconductors, and they strongly absorb UV radiation, so they could be used in UV light protection devices. They could also be used in solar cells and LEDs due to their low reflectivities.

Open image in new window

Graphical abstract Possible applications of doped C₂₄ fullerene

     

Presence of Ornithine in the Urate-binding α<Subscript>1</Subscript>—<Subscript>2</Subscript>Globulin

THE urate-binding α₁–α₂ globulin has been isolated from human plasma in a highly purified state¹. The protein was purified by DEAE-‘Sephadex’, ammonium sulphate precipitation and semi-preparative Polyacrylamide gel electrophoresis. The urate-binding α₁–α₂ globulin is a rod-shaped glycoprotein, containing 12.1% carbohydrate, with an isoelectric point of 4.6 and a molecular weight of 67,000 ± 4,000. Amino-acid analysis indicated an unknown basic compound which appeared as an extra peak just in front of lysine¹. To identify this compound, high voltage paper electrophoresis has been carried out on a plate electrophoresis apparatus in pyridine-acetate buffer pH 3.5. A spot separated out corresponding to ornithine. Amino-acid analysis on a BC-200 automatic analyser (Bio-Cal Instruments Co., West Germany), with a 54 cm column at 55° C and with 0.35 M sodium citrate buffer, pH 5.28, as elution buffer at a flow-rate of 150 ml./h, showed that ornithine was present. The presence of ornithine in the protein hydrolysate was also verified by gas chromatography/mass spectrometry².     

Structures and stability of N<Subscript>13</Subscript><Superscript>+</Superscript> and N<Subscript>13</Subscript><Superscript>−</Superscript> clusters

The structures and stabilities of eleven N₁₃ ⁺ and N₁₃ ^– isomers have been investigated with second-order Møller–Plesset (MP2) and density functional theory (DFT) methods. Five N₁₃ ⁺ isomers and six N₁₃ ^– isomers are all reasonable local minima on their potential energy hypersurfaces. The most stable N₁₃ ⁺ cation is structure C-2 with C_2v symmetry, which contains a pentazole ring and two N₄ open chains. It is different from those of the N₇ ⁺ and N₉ ⁺ clusters, but similar to the N₁₁ ⁺ cluster. Meanwhile, the most stable N₁₃ ^– structure A-2 is composed of a pentazole ring and a six-membered ring connected by two nitrogen atoms. It is not only different from those of the N₇ ^– and N₉ ^– clusters, but also from the N₁₁ ^– cluster. The decomposition pathways of structures C-2 and A-2 were investigated at the B3LYP/(aug)-cc-pVDZ level. From the barrier heights of the structures C-2 and A-2 decomposition processes, it is suggested that C-2 is difficult to observe experimentally and A-2 may be observed as a short-lived species. Figure Optimized geometrical parameters of N₁₃ ⁺ isomer C-2      

The neural γ<Subscript>2</Subscript>α<Subscript>1</Subscript>β<Subscript>2</Subscript>α<Subscript>1</Subscript>β<Subscript>2</Subscript> gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models—one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2–M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.     

Solid phase and solution synthesis of NvocLys(CO(CH<Subscript>2</Subscript>)<Subscript>5</Subscript>NH–NBD)OCH<Subscript>2</Subscript>CN,a trifunctional fluorescent lysine derivative

Herein, we describe a general strategy for the facile synthesis of a multifunctional amino acid derivative bearing both fluorescent and photolabile groups such as the lysine derivative NvocLys(CO(CH₂)₅NH–NBD)OCH₂CN (1) that can be used as a biophysical tool for studying protein structure. The synthetic strategy involves functionalization of the amine groups while the amino acid is attached to a solid support, followed by esterification of the carboxylic acid in solution. The solid support protects the caboxylic acid, preventing a side reaction associated with the synthesis in solution and obviating the need for chromatographic purification of several intermediates. This synthetic strategy can be used for the preparation of a variety of amino acid derivatives with unusual α-amine and side chain functionalities.     

Synthesis,characterization, and studies on DNA binding of the complex Fe(Sal<Subscript>2</Subscript>dienNO<Subscript>3</Subscript>·H<Subscript>2</Subscript>O)

A new metal complex, Fe(Sal₂dienNO₃·H₂O) (where Sal is salicylaldehyde and dien is diethylenetriamine), has been synthesized and characterized. The interactions between the Fe(III) complex and calf thymus DNA has been investigated using UV and fluorescence spectra, viscosity, thermal denaturation, and molecular modeling. The cleavage reaction on plasmid DNA has been monitored by agarose gel electrophoresis. The experimental results show that the mode of binding of the complex to DNA is classical intercalation and the complex can cleave pBR322 DNA.     

Insight into π-hole interactions containing the inorganic heterocyclic compounds S<Subscript>2</Subscript>N<Subscript>2</Subscript>/SN<Subscript>2</Subscript>P<Subscript>2</Subscript>

Similar to σ-hole interactions, the π-hole interaction has attracted much attention in recent years. According to the positive electrostatic potentials above and below the surface of inorganic heterocyclic compounds S₂N₂ and three SN₂P₂ isomers (heterocyclic compounds 1–4), and the negative electrostatic potential outside the X atom of XH₃ (X = N, P, As), S₂N₂/SN₂P₂?XH₃ (X = N, P, As) complexes were constructed and optimized at the MP2/aug-cc-pVTZ level. The X atom of XH₃ (X = N, P, As) is almost perpendicular to the ring of the heterocyclic compounds. The π-hole interaction energy becomes greater as the trend goes from 1?XH₃ to 4?XH₃. These π-hole interactions are weak and belong to “closed-shell” noncovalent interactions. According to the energy decomposition analysis, of the three attractive terms, the dispersion energy contributes more than the electrostatic energy. The polarization effect also plays an important role in the formation of π-hole complexes, with the contrasting phenomena of decreasing electronic density in the π-hole region and increasing electric density outside the X atom of XH₃ (X = N, P, As).

Open image in new window

Graphical abstract Computed density difference plots for the complexes 3?NH ₃ (a ₁), 3?PH ₃ (b ₁), 3?AsH ₃ (c ₁) and electron density shifts for the complexes 3?NH ₃ (a ₂), 3?PH ₃ (b ₂),3?AsH ₃ (c ₂) on the 0.001 a.u. contour

     

Kinetic analysis of recombinant mammalian α<Subscript>1</Subscript> and α<Subscript>1</Subscript>β glycine receptor channels

To analyze the influence of the beta-subunit on the kinetic properties of GlyR channel currents, alpha(1)-subunits and alpha(1)beta-subunits were transiently expressed in HEK 293 cells. A piezo dimorph was used for fast application of glycine to outside-out patches. The rise time of activation was dose dependent for both receptors and decreased with increasing glycine concentrations. Subunit composition had no effect on the time course of activation. Coexpression of alpha(1)- and beta-subunits resulted in a significantly lower EC(50) and a reduced slope of the dose-response curve of glycine compared with expression of alpha(1)-subunits alone. For both receptor subtypes, the time course of desensitization was concentration dependent. Desensitization was best fitted with a single time constant at 10-30 micro M, with two at 0.1 mM, and at saturating concentrations (0.3-3 mM) with three time constants. Desensitization of homomeric alpha(1)-receptor channels was significantly slower than that of alpha(1)beta-receptor channels. The time course of current decay after the end of glycine pulses was tested at different pulse durations of 1 mM glycine. It was best fitted with two time constants for both alpha(1) and alpha(1)beta GlyR channels, and increased significantly with increasing pulse duration.     

Effects of Prostaglandins F<Subscript>2α</Subscript>, E<Subscript>2</Subscript> and E<Subscript>1</Subscript> on Fertility in Mice

PROSTAGLANDIN (PG) F_2αhas antifertility effects in many species^1–3 but there are conflicting suggestions as to its mechanism of action. For example, it may cause the degeneration of the corpus luteum by decreasing blood flow in the uteroovarian vein⁴; alternatively, its action may be due to a hypersecretion of luteinizing hormone (LH) by the pituitary^3,5. I have investigated the effects of PGF_2α, E₂ and E₁ on pregnancy in mice and examined the mechanism of action of PGF_2α.     

Effects of activation of μ-, κ<Subscript>1</Subscript>-, δ<Subscript>1</Subscript>-, δ<Subscript>2</Subscript>-, and ORL<Subscript>1</Subscript>-receptors on heart resistance to the pathogenic action of delayed ischemia and reperfusion

Different subtypes of opioid receptors (OR) were activated in rats in vivo to study the activation effect on the heart’s resistance to ischemia and reperfusion. It has been established that administration of deltorphin II, a selective δ₂-OR agonist, lowered the infarct size/area at risk index (IS/AAR) by 23%. Naltrexone, naloxone methiodide (an OR inhibitor not penetrating the blood-brain barrier (BBB)), and naltriben (δ₂-antagonist) eliminated the cardioprotective effect of deltorphin II, while BNTX (a δ₁-antagonist) produced no effect on the cardioprotective action of the δ₂-agonist. The infarct-reducing effect of deltorphin II was eliminated by administration of chelerythrine (a protein kinase C (PKC) inhibitor), glibenclamide (a K_ATP-channels inhibitor), and 5-hydroxydecanoate (a mitochondrial K_ATP-channel blocker). Administration of other opioids did not reduce the IS/AAR index. It has been established that all the deltorphins manifest antiarrhythmic potency. Other opioids do not produce any effect on the incidence of arrhythmia occurrences. The antiarrhythmic effect of deltorphin II was eliminated by preliminary administration of naltrexone, naloxone methiodide, and naltriben, but BNTX did not affect the δ₂-agonist’s anti-arrhythmic effect. The preliminary administration of chelerythrine, a PKC inhibitor, eliminated the δ₂ agonist’s antiarrhythmic action. However, glibenclamide and 5-hydroxydecanoate did not alter the antiarrhythmic effect by deltorphin II. Therefore, activation of the peripheral δ₂-ORs reduces the infarct size and prevents the onset of arrhythmias. The antiarrhythmic effect of the δ₂-OR stimulation is mediated by activating PKC and opening the mitochondrial K_ATP-channels. PKC participates in the antiarrhythmic effect of the δ₂-OR activation, but this effect does not depend on the condition of K_ATP-channels.     

Succinate is the controller of O<Stack><Subscript>2</Subscript><Superscript>−</Superscript></Stack>/H<Subscript>2</Subscript>O<Subscript>2</Subscript> release at mitochondrial complex I : negative modulation by malate,positive by cyanide

Mitochondrial production of H₂O₂ is low with NAD substrates (glutamate/pyruvate, 3 and 2 mM) (G/P) and increases over ten times upon further addition of succinate, with the formation of a sigmoidal curve (semimaximal value at 290 μM, maximal H₂O₂ production at 600 μM succinate). Malate counteracts rapidly the succinate induced increased H₂O₂ release and moves the succinate dependent H₂O₂ production curve to the right. Nitric oxide (NO) and carbon monoxide (CO) are cytochrome c oxidase inhibitors which increase mitochondrial ROS production. Cyanide (CN⁻) was used to mimic NO and CO. In the presence of G/P and succinate (300 μM), CN⁻ progressively increased the H₂O₂ release rate, starting at 1.5 μM. The succinate dependent H₂O₂ production curve was moved to the left by 30 μM CN⁻. The V_max was little modified. We conclude that succinate is the controller of mitochondrial H₂O₂ production, modulated by malate and CN⁻. We propose that succinate promotes an interaction between Complex II and Complex I, which activates O₂⁻ production.     

Association of β<Subscript>2</Subscript>-microglobulin with the α<Subscript>3</Subscript> domain of H-2D<Superscript>b</Superscript> heavy chain

MHC class I molecules are heterotrimeric complexes composed of heavy chain, ₂-microglobulin (₂m) and short peptide. This trimeric complex is generated in the endoplasmic reticulum (ER), where a peptide loading complex (PLC) facilitates transport from the cytosol and binding of the peptide to the preassembled ER resident heavy chain/₂m dimers. Association of mouse MHC class I heavy chain with ₂m is characterized by allelic differences in the number and/or positions of amino acid interactions. It is unclear, however, whether all alleles follow common binding patterns with minimal contributions by allele-specific contacts, or whether essential contacts with ₂m are different for each allele. While searching for the PLC binding site in the ₃ domain of the mouse MHC class I molecule H-2D^b, we unexpectedly discovered a site critical for binding mouse, but not human, ₂m. Interestingly, amino acids in the corresponding region of another MHC class I heavy chain allele do not make contacts with the mouse ₂m. Thus, there are allelic differences in the modes of binding of ₂m to the heavy chain of MHC class I.     

Negative molecular ion of deuterium D<Stack><Subscript>2</Subscript><Superscript>−</Superscript></Stack>

D₂⁻ ions produced in collisions of D⁻ ions with relative energies of 2.5–9.2 eV were detected for the first time. It is shown that the effective cross section for this reaction is no less than 1.5 × 10⁻¹⁴ cm². Along with the theoretically predicted short-lived state of negative molecular deuterium ions, a state existing for more than 1 μs was observed.     

Spectroscopic Behavior and Biological Activity of K<Subscript>2</Subscript>[VO(O<Subscript>2</Subscript>)NTA]·2H<Subscript>2</Subscript>O

The dihydrated potassium salt of the complex anion [VO(O₂)NTA]²⁻ (NTA = nitrilotriacetate anion, [N(CH₂-COO)₃]³⁻) was thoroughly characterized by electronic and vibrational (infrared and Raman) spectroscopies. The bioactivity of the complex on the cell proliferation was tested on three cell lines in culture (UMR106 rat osteosarcoma-derived cells, Caco-2 derived from a human colon adenocarcinoma, and RAW 264.7, a macrophage murine cell line).     

NMR analysis of duplex d(CGCGATCGCG)<Subscript>2</Subscript> modified by <Emphasis Type="Italic">Λ</Emphasis>- and <Emphasis Type="Italic">Δ</Emphasis>-[Ru(bpy)<Subscript>2</Subscript>(m-GHK)]Cl<Subscript>2</Subscript> and DNA photocleavage study

The interaction of the diastereomeric complexes Λ-[Ru(bpy)₂(m-GHK)]Cl₂ and Δ-[Ru(bpy)₂(m-GHK)]Cl₂ (bpy is 2,2′-bipyridine, GHK is glycine–l-histidine–l-lysine) with the deoxynucleotide duplex d(5′-CGCGATCGCG)₂ was studied by means of ¹H NMR spectroscopy. At a Δ-isomer to DNA ratio of 1:1, significant shifts for the metal complex are observed, whereas there is negligible effect on the oligonucleotide protons and only one intermolecular nuclear Overhauser effect (NOE) is present at the 2D nuclear Overhauser enhancement spectroscopy spectrum. The ¹Η NMR spectrum at ratio 2:1 is characterized by a slight shift for the Δ-isomer’s bpy aromatic protons as well as significant shifts for the decanucleotide G₄ H1′ and Η2″, A₅ H2, G₁₀ H1′, T₆ NH and G₂ NH protons. Furthermore, at ratio 2:1, 11 intermolecular NOEs are observed. The majority of the NOEs involve the sugar Η2′ and Η2″ protons sited in the major groove of the decanucleotide. Increasing the Δ-isomer to d(CGCGATCGCG)₂ ratio to 5:1 results in noteworthy spectral changes. The Δ-isomer’s proton shifts are reduced, whereas significant shifts are observed for the decanucleotide protons, especially the sugar protons, as well as for the exchangeable protons. Interaction is characterized by the presence of only one intermolecular NOE. Furthermore, there is significant broadening of the imino proton signals as the ratio of the Δ-isomer to DΝΑ increases, which is attributed to the opening of the two strands of the duplex. The Λ-isomer, on the other hand, approaches the minor groove of the oligonucleotide and interacts only weakly, possibly by electrostatic interactions. Photocleavage studies were also conducted with the plasmid pUC19 and a 158-bp restriction fragment, showing that both diastereomers cleave DNA with similar efficiency, attacking mainly the guanines of the sequence probably by generating active oxygen species. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. An erratum to this article can be found at     

Cooperation of adenosine and prostaglandin E<Subscript>2</Subscript> (PGE<Subscript>2</Subscript>) in amplification of cAMP–PKA signaling and immunosuppression

INTRODUCTION: We hypothesize that adenosine and PGE(2) could have a complementary immunosuppressive effect that is mediated via common cAMP-PKA signaling. MATERIALS AND METHODS: To test this hypothesis, the effect of adenosine and PGE(2) on the cytotoxic activity and cytokine production of lymphokine activated killer (LAK) cells was investigated. RESULTS: PGE(2) and adenosine inhibited LAK cells cytotoxic activity and production of INF-gamma, GM-CSF and TNF-alpha. In combination they showed substantially higher inhibition than each modality used alone. Using agonists and antagonists specific for PGE(2) and adenosine receptors we found that cooperation of PGE(2) and adenosine in their inhibitory effects are mediated via EP(2) and A(2A) receptors, respectively. LAK cells have 35-fold higher expression of EP(2) than A(2A). Combined PGE(2) and adenosine treatment resulted in augmentation of cAMP production, PKA activity, CREB phosphorylation and inhibition of Akt phosphorylation. Wortmannin and LY294002 enhanced the suppressive effects of adenosine and PGE(2). In contrast, Rp-8-Br-cAMPS, an inhibitor of PKA type I blocked their immunosuppressive effects, suggesting that the inhibitory effects of PGE(2) and adenosine are mediated via common pathway with activation of cAMP-PKA and inhibition of Akt. CONCLUSION: In comparison to other immunosuppressive molecules (TGF-beta and IL-10), adenosine and PGE(2) are unique in their ability to inhibit the executive function of highly cytotoxic cells. High intratumor levels of adenosine and PGE(2) could protect tumor from immune-mediated destruction by inactivation of the tumor infiltrating functionally active immune cells.     

Metal–metal bonding and aromaticity in [M<Subscript>2</Subscript>(NHCHNH)<Subscript>3</Subscript>]<Subscript>2</Subscript> (μ-E)<Subscript>2</Subscript> (E = O,S; M = Nb,Mo, Tc,Ru, Rh)

The nature of M–M bonding and aromaticity of [M₂(NHCHNH)₃]₂(μ-E)₂ (E?=?O, S; M?=?Nb, Mo, Tc, Ru, Rh) was investigated using atoms in molecules (AIM) theory, electron localization function (ELF), natural bond orbital (NBO) and molecular orbital analysis. These analyses led to the following main conclusions: in [M₂(NHCHNH)₃]₂(μ-E)₂ (E?=?O, S; M?=?Nb, Mo, Tc, Ru, Rh), the Nb–Nb, Ru–Ru, and Rh–Rh bonds belong to “metallic” bonds, whereas Mo–Mo and Tc–Tc drifted toward the “dative” side; all these bonds are partially covalent in character. The Nb–Nb, Mo–Mo, and Tc–Tc bonds are stronger than Ru–Ru and Rh–Rh bonds. The M–M bonds in [M₂(NHCHNH)₃]₂(μ-S)₂ are stronger than those in [M₂(NHCHNH)₃]₂(μ-O)₂ for M?=?Nb, Mo, Tc, and Ru. The NICS(1)_ZZ values show that all of the studied molecules, except [Ru₂(NHCHNH)₃]₂(μ-O)₂, are aromaticity molecules. O-bridged compounds have more aromaticity than S-bridged compounds.

Open image in new window

Graphical Abstract Left Molecular graph, and right electron localization function (ELF) isosurface of [M₂(NHCHNH)₃]₂(μ-E)₂(E?=?O, S; M?=?Nb, Mo, Tc, Ru, Rh)