IF<Subscript>1</Subscript> distribution in HepG2 cells in relation to ecto–F<Subscript>0</Subscript>F<Subscript>1</Subscript>ATPsynthase and calmodulin

F₀F₁ATPsynthase is now known to be expressed as a plasma membrane receptor for several extracellular ligands. On hepatocytes, ecto–F₀F₁ATPsynthase binds apoA–I and triggers HDL endocytosis concomitant with ATP hydrolysis. Considering that inhibitor protein IF₁ was shown to regulate the hydrolytic activity of ecto–F₀F₁ATPsynthase and to interact with calmodulin (CaM) in vitro, we investigated the subcellular distributions of IF₁, calmodulin (CaM), OSCP and β subunits of F₀F₁ATPsynthase in HepG2 cells. Using immunofluorescence and Western blotting, we found that around 50% of total cellular IF₁ is localized outside mitochondria, a relevant amount of which is associated to the plasma membrane where we also found Ca²⁺–CaM, OSCP and β. Confocal microscopy showed that IF₁ colocalized with Ca²⁺–CaM on plasma membrane but not in mitochondria, suggesting that Ca²⁺–CaM may modulate the cell surface availability of IF₁ and thus its ability to inhibit ATP hydrolysis by ecto–F₀F₁ATPsynthase. These observations support a hypothesis that the IF₁–Ca²⁺–CaM complex, forming on plasma membrane, functions in the cellular regulation of HDL endocytosis by hepatocytes.     

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α.     

The effect of NBD-Cl in nucleotide-binding of the major subunit α and B of the motor proteins F<Subscript>1</Subscript>F<Subscript>O</Subscript> ATP synthase and A<Subscript>1</Subscript>A<Subscript>O</Subscript> ATP synthase

Subunit α of the Escherichia coli F₁F_O ATP synthase has been produced, and its low-resolution structure has been determined. The monodispersity of α allowed the studies of nucleotide-binding and inhibitory effect of 4-Chloro-7-nitrobenzofurazan (NBD-Cl) to ATP/ADP-binding. Binding constants (K _d ) of 1.6 μM of bound MgATP-ATTO-647N and 2.9 μM of MgADP-ATTO-647N have been determined from fluorescence correlation spectroscopy data. A concentration of 51 μM and 55 μM of NBD-Cl dropped the MgATP-ATTO-647N and MgADP-ATTO-647N binding capacity to 50% (IC₅₀), respectively. In contrast, no effect was observed in the presence of N,N′-dicyclohexylcarbodiimide. As subunit α is the homologue of subunit B of the A₁A_O ATP synthase, the interaction of NBD-Cl with B of the A-ATP synthase from Methanosarcina mazei Gö1 has also been shown. The data reveal a reduction of nucleotide-binding of B due to NBD-Cl, resulting in IC₅₀ values of 41 μM and 42 μM for MgATP-ATTO-647N and MgADP-ATTO-647N, respectively.     

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².     

Thyroid hormone receptor α<Subscript>1</Subscript>–β<Subscript>1</Subscript> expression in epididymal epithelium from euthyroid and hypothyroid rats

The objectives of the present work were to assess whether epithelial cells from the different segments of epididymis express TRα₁–β₁ isoforms, to depict its subcellular immunolocalization and to evaluate changes in their expression in rats experimentally submitted to a hypothyroid state by injection of 131I. In euthyroid and hypothyroid groups, TR protein was expressed in epididymal epithelial cells, mainly in the cytoplasmic compartment while only a few one showed a staining in the nucleus as well. A similar TR immunostaining pattern was detected in the different segments of the epididymis. In hypothyroid rats, the number of TR-immunoreactive epithelial cells as well as the intensity of the cytoplasmic staining significantly increased in all sections analyzed. In consonance to the immunocytochemical analysis, the expression of TRα₁–β₁ isoforms, assessed by Western blot revealed significantly higher levels of TR in cytosol compared to the nuclear fractions. Furthermore, TR expression of both α₁ and β₁ isoforms and their mRNA levels were increased by the hypothyroid state. The immuno-electron-microscopy showed specific reaction for TR in principal cells associated with eucromatin, cytosolic matrix and mitochondria. The differences in expression levels assessed in control and thyroidectomized rats ascertain a specific function of TH on this organ.     

Theoretical descriptions of novel triplet germylenes M<Subscript>1</Subscript>-Ge-M<Subscript>2</Subscript>-M<Subscript>3</Subscript> (M<Subscript>1</Subscript> = H,Li, Na,K; M<Subscript>2</Subscript> = Be,Mg, Ca; M<Subscript>3</Subscript> = H,F, Cl,Br)

In a quest to identify new ground-state triplet germylenes, the stabilities (singlet–triplet energy differences, ΔE_S–T) of 96 singlet (s) and triplet (t) M₁-Ge-M₂-M₃ species were compared and contrasted at the B3LYP/6–311++G**, QCISD(T)/6–311++G**, and CCSD(T)/6–311++G** levels of theory (M₁?=?H, Li, Na, K; M₂?=?Be, Mg, Ca; M₃?=?H, F, Cl, Br). Interestingly, F-substituent triplet germylenes (M₃?=?F) appear to be more stable and linear than the corresponding Cl- or Br-substituent triplet germylenes (M₃?=?Cl or Br). Triplets with M₁?=?K (i.e., the K-Ge-M₂-M₃ series) seem to be more stable than the corresponding triplets with M₁?=?H, Li, or Na. This can be attributed to the higher electropositivity of potassium. Triplet species with M₃?=?Cl behave similarly to those with M₃?=?Br. Conversely, triplets with M₃?=?H show similar stabilities and linearities to those with M₃?=?F. Singlet species of formulae K-Ge-Ca-Cl and K-Ge-Ca-Br form unexpected cyclic structures. Finally, the triplet germylenes M₁-Ge-M₂-M₃ become more stable as the electropositivities of the α-substituents (M₁ and M₂) and the electronegativity of the β-substituent (M₃) increase.     

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.     

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.     

Should maize doubled haploids be induced among F<Subscript>1</Subscript> or F<Subscript>2</Subscript> plants?

Maize (Zea mays L.) doubled haploid lines are typically produced from F₁ plants. Studies have suggested that the low frequency of recombinants in doubled haploids may reduce the response to selection. My objective was to determine if, for sustaining long-term response, doubled haploids should be induced in F₁ or F₂ plants during maize inbred development. In simulation experiments, I examined the response to multiple cycles of testcross selection among doubled haploid lines derived from F₁ plants (denoted by DH), doubled haploid lines derived from F₂ plants (DH_F2), and recombinant inbred (RI) lines derived by single-seed descent. For a trait controlled by 100 or more quantitative trait loci (QTL), the cumulative responses to selection were up to 4–6% larger among DH_F2 lines than among DH lines. The cumulative responses were up to 5–8% larger among RI lines than among DH lines. The QTL become unlinked as the number of QTL in a finite genome decreases, and the responses among RI, DH, and DH_F2 lines were equal or nearly equal when only 20 QTL controlled the trait. Metabolic-flux epistasis reduced the differences in the response among RI, DH, and DH_F2 lines. Overall, the results indicated that doubled haploids should be induced from F₂ plants rather than from F₁ plants. If year-round nurseries are used and new F₁ crosses for inbred development are initially created on a speculative basis, the development of doubled haploids from F₂ rather than F₁ plants should not cause a delay in inbred development.     

Control of rotation of the F<Subscript>1</Subscript>F<Subscript>O</Subscript>-ATP synthase nanomotor by an inhibitory α-helix from unfolded ε or intrinsically disordered ζ and IF<Subscript>1</Subscript> proteins

The ATP synthase is a ubiquitous nanomotor that fuels life by the synthesis of the chemical energy of ATP. In order to synthesize ATP, this enzyme is capable of rotating its central rotor in a reversible manner. In the clockwise (CW) direction, it functions as ATP synthase, while in counter clockwise (CCW) sense it functions as an proton pumping ATPase. In bacteria and mitochondria, there are two known canonical natural inhibitor proteins, namely the ε and IF₁ subunits. These proteins regulate the CCW F₁F_O-ATPase activity by blocking γ subunit rotation at the α_DP/β_DP/γ subunit interface in the F₁ domain. Recently, we discovered a unique natural F₁-ATPase inhibitor in Paracoccus denitrificans and related α-proteobacteria denoted the ζ subunit. Here, we compare the functional and structural mechanisms of ε, IF₁, and ζ, and using the current data in the field, it is evident that all three regulatory proteins interact with the α_DP/β_DP/γ interface of the F₁-ATPase. In order to exert inhibition, IF₁ and ζ contain an intrinsically disordered N-terminal protein region (IDPr) that folds into an α-helix when inserted in the α_DP/β_DP/γ interface. In this context, we revised here the mechanism and role of the ζ subunit as a unidirectional F-ATPase inhibitor blocking exclusively the CCW F₁F_O-ATPase rotation, without affecting the CW-F₁F_O-ATP synthase turnover. In summary, the ζ subunit has a mode of action similar to mitochondrial IF₁, but in α-proteobacteria. The structural and functional implications of these intrinsically disordered ζ and IF₁ inhibitors are discussed to shed light on the control mechanisms of the ATP synthase nanomotor from an evolutionary perspective.     

G protein β<Subscript>1</Subscript>γ<Subscript>2</Subscript> subunits purification and their interaction with adenylyl cyclase

A preliminary study on the interaction of G protein (guanine triphosphate binding pro- tein) β1γ2 subunits and their coupled components in cell signal transduction was conducted in vitro. The insect cell lines, Sf9 (Spodoptera frugiperda) and H5 (Trichoplusia ni) were used to express the recombinant protein Gβ1γ2. The cell membrane containing Gβ1γ2 was isolated through affinity chromatography column with Ni-NTA agarose by FPLC method, and the highly purified protein was obtained. The adenylyl cyclase 2 (AC2) activity assay showed that the purified Gβ1γ2 could significantly stimulate AC2 activity. The interaction of β1γ2 subunits of G protein with the cytoplasmic tail of various mammalian adenylyl cyclases was monitored by BIAcore technology using NTA sensor chip, which relies on the phenomenon of surface plasmon resonance (SPR). The experiments showed the direct binding of Gβ1γ2 to the cytoplasmic tail C2 domain of AC2. The specific binding domain of AC2 with Gβ1γ2 was the same as AC2 activity domain which was stimulated by β1γ2.     

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

     

Guanidino Acids Act as ρ1 GABA<Subscript>C</Subscript> Receptor Antagonists

GABA_C receptors play a role in myopia, memory-related disorders and circadian rhythms signifying a need to develop potent and selective agents for this class of receptors. Guanidino analogs related to glycine, β-alanine and taurine were evaluated at human ρ₁GABA_C receptors expressed in Xenopus oocytes using 2-electrode voltage clamp methods. Of the 12 analogs tested, 8 analogs were active as antagonists and the remaining were inactive. (S)-2-Guanidinopropionic acid (IC₅₀ = 2.2 μM) and guanidinoacetic acid (IC₅₀ = 5.4 μM; K _B = 7.75 μM [pK _B = 5.11 ± 0.06]) were the most potent being competitive antagonists at this receptor. In contrast, the β-alanine and GABA guanidino analogs showed reduced activity, indicating the distance between the carboxyl carbon and terminal nitrogen of the guanidino group is critical for activity. Substituting the C2-position of guanidinoacetic acid with various alkyl groups reduced activity indicating that steric effects may impact on activity. The results of this study contribute to the structure–activity-relationship profile required in developing novel therapeutic agents.     

Co-transformation of <Emphasis Type="Italic">Panax</Emphasis> major ginsenosides Rb<Subscript>1</Subscript> and Rg<Subscript>1</Subscript> to minor ginsenosides C–K and F<Subscript>1</Subscript> by <Emphasis Type="Italic">Cladosporium cladosporioides</Emphasis>

Rb₁ and Rg₁ are the major ginsenosides in protopanaxadiol and protopanaxatriol. Their content in ginsenosides was 23.8 and 17.6%, respectively. A total of 22 isolates of β-glucosidase producing microorganisms were isolated from the soil of a ginseng field using Esculin-R2A agar. Among these isolates, the strain GH21 showed the strongest activities to convert ginsenoside Rb₁ and Rg₁ to minor ginsenosides compound-K and F₁, respectively. Ginsenosides Rb₁ and Rg₁ bioconversion rates were 74.2 and 89.3%, respectively. Meanwhile, the results demonstrated that the ginsenoside Rg₁ could change the biotransformation pathway of ginsenoside Rb₁ by inhibiting the formation of the intermediate metabolite gypenoside-XVII. GH21 was identified as a Cladosporium cladosporioides species based on the internal transcribed spacers (ITS) ITS1-5.8S-ITS2 rRNA gene sequences constructed phylogenetic trees.     

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

     

Deciphering the binding mode of Zolpidem to GABA<Subscript>A</Subscript> α<Subscript>1</Subscript> receptor – insights from molecular dynamics simulation

To investigate the binding mode of Zolpidem to GABA(A) and to delineate the conformational changes induced upon agonist binding, we carried out atomistic molecular dynamics simulation using the ligand binding domain of GABA(A) α(1) receptor. Comparative molecular dynamics simulation of the apo and the holo form of GABA(A) receptor revealed that γ(2)/α(1) interface housing the benzodiazepine binding site undergoes distinct conformational changes upon Zolpidem binding. We notice that C loop of the α(1) subunit experiences an inward motion toward the vestibule and the F loop of γ(2) sways away from the vestibule, an observation that rationalizes Zolpidem as an alpha1 selective agonist. Energy decomposition analysis carried out was able to highlight the important residues implicated in Zolpidem binding, which were largely in congruence with the experimental data. The simulation study disclosed herein provides a meaningful insight into Zolpidem-GABA(A)R interactions and helps to arrive at a binding mode hypothesis with implications for drug design.     

NMR solution structure of the N-terminal domain of subunit E (E<Subscript>1–52</Subscript>) of A<Subscript>1</Subscript>A<Subscript>O</Subscript> ATP synthase from <Emphasis Type="Italic">Methanocaldococcus jannaschii</Emphasis>

The N-termini of E and H of A₁A_O ATP synthase have been shown to interact and an NMR structure of N-terminal H_1–47 has been solved recently. In order to understand the E-H assembly and the N-terminal structure of E, the truncated construct E_1–52 of Methanocaldococcus jannaschii A₁A_O ATP synthase was produced, purified and the solution structure of E_1–52 was determined by NMR spectroscopy. The protein is 60.5 Å in length and forms an α helix between the residues 8–48. The molecule is amphipathic with a strip of hydrophobic residues, discussed as a possible helix-helix interaction with neighboring subunit H.     

A theoretical study on the reaction mechanism of O<Subscript>2</Subscript> with C<Subscript>4</Subscript>H<Subscript>9</Subscript>• radical

Ab initio calculations have been performed using the complete basis set model (CBS-QB3) to study the reaction mechanism of butane radical (C₄H₉•) with oxygen (O₂). On the calculated potential energy surface, the addition of O₂ to C₄H₉• forms three intermediates barrierlessly, which can undergo subsequent isomerization or decomposition reaction leading to various products: HOO• + C₄H₈, C₂H₅• + CH₂CHOOH, OH• + C₃H₇CHO, OH• + cycle-C₄H₈O, CH₃• + CH₃CHCHOOH, CH₂OOH• + C₃H₆. Five pathways are supposed in this study. After taking into account the reaction barrier and enthalpy, the most possible reaction pathway is C₄H₉• + O₂ → IM1 → TS5 → IM3 → TS6 → IM4 → TS7 → OH• + cycle-C₄H₈O.