Ca+ Ions Solvated in Helium Clusters

We present a combined experimental and theoretical investigation on Ca+ ions in helium droplets, HeNCa+. The clusters have been formed in the laboratory by means of electron-impact ionization of Ca-doped helium nanodroplets. Energies and structures of such complexes have been computed using various approaches such as path integral Monte Carlo, diffusion Monte Carlo and basin-hopping methods. The potential energy functions employed in these calculations consist of analytical expressions following an improved Lennard-Jones formula whose parameters are fine-tuned by exploiting ab initio estimations. Ion yields of HeNCa+ -obtained via high-resolution mass spectrometry- generally decrease with N with a more pronounced drop between N=17 and N=25, the computed quantum HeNCa+ evaporation energies resembling this behavior. The analysis of the energies and structures reveals that covering Ca+ with 17 He atoms leads to a cluster with one of the smallest energies per atom. As new atoms are added, they continue to fill the first shell at the expense of reducing its stability, until N=25, which corresponds to the maximum number of atoms in that shell. Behavior of the evaporation energies and radial densities suggests liquid-like cluster structures.     

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.     

Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative

Human immunodeficiency virus (HIV) is a life life-threatening and serious infection caused by a virus that attacks CD4+ T-cells, which fight against infections and make a person susceptible to other diseases. It is a global public health problem with no cure; therefore, it is highly important to study and understand the intricate phenomena of HIV. In this article, we focus on the numerical study of the path-tracking damped oscillatory behavior of a model for the HIV infection of CD4+ T-cells. We formulate fractional dynamics of HIV with a source term for the supply of new CD4⁺ T-cells depending on the viral load via the Caputo–Fabrizio derivative. In the formulation of fractional HIV dynamics, we replaced the constant source term for the supply of new CD4+ T-cells from the thymus with a variable source term depending on the concentration of the viral load, and introduced a term that describes the incidence of the HIV infection of CD4+ T-cells. We present a novel numerical scheme for fractional view analysis of the proposed model to highlight the solution pathway of HIV. We inspect the periodic and chaotic behavior of HIV for the given values of input factors using numerical simulations.     

TiCl4 Dissolved in Ionic Liquid Mixtures from Аb Initio Molecular Dynamics Simulations

To gain a deeper understanding of the TiCl4 solvation effects in multi-component ionic liquids, we performed ab initio molecular dynamics simulations of 1-butyl-3-methylimidazolium [C4C1Im]+, tetrafluoroborate [BF4]−, chloride [Cl]− both with and without water and titanium tetrachloride TiCl4. Complex interactions between cations and anions are observed in all investigated systems. By further addition of water and TiCl4 this complex interaction network is extended. Observations of the radial distribution functions and number integrals show that water and TiCl4 not only compete with each other to interact mainly with [Cl]−, which strongly influences the cation-[BF4]− interaction, but also interact with each other, which leads to the fact that in certain systems the cation-anion interaction is enhanced. Further investigations of the Voronoi polyhedra analysis have demonstrated that water has a greater impact on the nanosegregated system than TiCl4 which is also due to the fact of the shear amount of water relative to all other components and its higher mobility compared to TiCl4. Overall, the polar network of the IL mixture collapses by including water and TiCl4. In the case of [Cl]− chloride enters the water continuum, while [BF4]− remains largely unaffected, which deeply affects the interaction of the ionic liquid (IL) network.     

Understanding Topological Insulators in Real Space

A real space understanding of the Su–Schrieffer–Heeger model of polyacetylene is introduced thanks to delocalization indices defined within the quantum theory of atoms in molecules. This approach enables to go beyond the analysis of electron localization usually enabled by topological insulator indices—such as IPR—enabling to differentiate between trivial and topological insulator phases. The approach is based on analyzing the electron delocalization between second neighbors, thus highlighting the relevance of the sublattices induced by chiral symmetry. Moreover, the second neighbor delocalization index, δi,i+2, also enables to identify the presence of chirality and when it is broken by doping or by eliminating atom pairs (as in the case of odd number of atoms chains). Hints to identify bulk behavior thanks to δ1,3 are also provided. Overall, we present a very simple, orbital invariant visualization tool that should help the analysis of chirality (independently of the crystallinity of the system) as well as spreading the concepts of topological behavior thanks to its relationship with well-known chemical concepts.     

Magnetic-Field-Orientation Dependent Thermal Entanglement of a Spin-1 Heisenberg Dimer: The Case Study of Dinuclear Nickel Complex with an Uniaxial Single-Ion Anisotropy

The bipartite entanglement in pure and mixed states of a quantum spin-1 Heisenberg dimer with exchange and uniaxial single-ion anisotropies is quantified through the negativity in a presence of the external magnetic field. At zero temperature the negativity shows a marked stepwise dependence on a magnetic field with two abrupt jumps and plateaus, which can be attributed to the quantum antiferromagnetic and quantum ferrimagnetic ground states. The magnetic-field-driven phase transition between the quantum antiferromagnetic and quantum ferrimagnetic ground states manifests itself at nonzero temperatures by a local minimum of the negativity, which is followed by a peculiar field-induced rise of the negativity observable in a range of moderately strong magnetic fields. The rising temperature generally smears out abrupt jumps and plateaus of the negativity, which cannot be distinguished in the relevant dependencies above a certain temperature. It is shown that the thermal entanglement is most persistent against rising temperature at the magnetic field, for which an energy gap between a ground state and a first excited state is highest. Besides, temperature variations of the negativity of the spin-1 Heisenberg dimer with an easy-axis single-ion anisotropy may exhibit a singular point-kink, at which the negativity has discontinuity in its first derivative. The homodinuclear nickel complex [Ni2(Medpt)2(μ-ox)(H2O)2](ClO4)2·2H2O provides a suitable experimental platform of the antiferromagnetic spin-1 Heisenberg dimer, which allowed us to estimate a strength of the bipartite entanglement between two exchange-coupled Ni2+ magnetic ions on the grounds of the interaction constants reported previously from the fitting procedure of the magnetization data. It is verified that the negativity of this dinuclear compound is highly magnetic-field-orientation dependent due to presence of a relatively strong uniaxial single-ion anisotropy.     

MHD Stagnation Point on Nanofluid Flow and Heat Transfer of Carbon Nanotube over a Shrinking Surface with Heat Sink Effect

This study is to investigate the magnetohydrodynamic (MHD) stagnation point flow and heat transfer characteristic nanofluid of carbon nanotube (CNTs) over the shrinking surface with heat sink effects. Similarity equations deduced from momentum and energy equation of partial differential equations are solved numerically. This study looks at the different parameters of the flow and heat transfer using first phase model which is Tiwari-Das. The parameter discussed were volume fraction nanoparticle, magnetic parameter, heat sink/source parameters, and a different type of nanofluid and based fluids. Present results revealed that the rate of nanofluid (SWCNT/kerosene) in terms of flow and heat transfer is better than (MWCNT/kerosene) and (CNT/water) and regular fluid (water). Graphically, the variation results of dual solution exist for shrinking parameter in range λc<λ≤−1 for different values of volume fraction nanoparticle, magnetic, heat sink parameters, and a different type of nanofluid. However, a unique solution exists at −1<λ<1, and no solutions exist at λ<λc which is a critical value. In addition, the local Nusselt number decreases with increasing volume fraction nanoparticle when there exists a heat sink effect. The values of the skin friction coefficient and local Nusselt number increase for both solutions with the increase in magnetic parameter. In this study, the investigation on the flow and heat transfer of MHD stagnation point nanofluid through a shrinking surface with heat sink effect shows how important the application to industrial applications.     

Zhang–Zhang Polynomials of Multiple Zigzag Chains Revisited: A Connection with the John–Sachs Theorem

Multiple zigzag chains Zm,n of length n and width m constitute an important class of regular graphene flakes of rectangular shape. The physical and chemical properties of these basic pericondensed benzenoids can be related to their various topological invariants, conveniently encoded as the coefficients of a combinatorial polynomial, usually referred to as the ZZ polynomial of multiple zigzag chains Zm,n. The current study reports a novel method for determination of these ZZ polynomials based on a hypothesized extension to John–Sachs theorem, used previously to enumerate Kekulé structures of various benzenoid hydrocarbons. We show that the ZZ polynomial of the Zm,n multiple zigzag chain can be conveniently expressed as a determinant of a Toeplitz (or almost Toeplitz) matrix of size m2×m2 consisting of simple hypergeometric polynomials. The presented analysis can be extended to generalized multiple zigzag chains Zkm,n, i.e., derivatives of Zm,n with a single attached polyacene chain of length k. All presented formulas are accompanied by formal proofs. The developed theoretical machinery is applied for predicting aromaticity distribution patterns in large and infinite multiple zigzag chains Zm,n and for computing the distribution of spin densities in biradical states of finite multiple zigzag chains Zm,n.     

The Flow of Blood-Based Hybrid Nanofluids with Couple Stresses by the Convergent and Divergent Channel for the Applications of Drug Delivery

This research work aims to scrutinize the mathematical model for the hybrid nanofluid flow in a converging and diverging channel. Titanium dioxide and silver TiO2 and Ag are considered as solid nanoparticles while blood is considered a base solvent. The couple-stress fluid model is essentially use to describe the blood flow. Therefore, the couple-stress term was used in the recent study with the existence of a magnetic field and a Darcy–Forchheiner porous medium. The heat absorption/omission and radiation terms were also included in the energy equation for the sustainability of drug delivery. An endeavor was made to link the recent study with the applications of drug delivery. It has already been revealed by the available literature that the combination of TiO2 with any other metal can destroy cancer cells more effectively than TiO2 separately. Both the walls are stretchable/shrinkable, whereas flow is caused by a source or sink with α as a converging/diverging parameter. Governing equations were altered into the system of non-linear coupled equations by using the similarity variables. The homotopy analysis method (HAM) was applied to obtain the preferred solution. The influences of the modeled parameters have been calculated and displayed. The confrontation of wall shear stress and hybrid nanofluid flow increased as the couple stress parameter rose, which indicates an improvement in the stability of the base fluid (blood). The percentage (%) increase in the heat transfer rate with the variation of nanoparticle volume fraction was also calculated numerically and discussed theoretically.     

Binding Ensembles of p53-MDM2 Peptide Inhibitors by Combining Bayesian Inference and Atomistic Simulations

Designing peptide inhibitors of the p53-MDM2 interaction against cancer is of wide interest. Computational modeling and virtual screening are a well established step in the rational design of small molecules. But they face challenges for binding flexible peptide molecules that fold upon binding. We look at the ability of five different peptides, three of which are intrinsically disordered, to bind to MDM2 with a new Bayesian inference approach (MELD × MD). The method is able to capture the folding upon binding mechanism and differentiate binding preferences between the five peptides. Processing the ensembles with statistical mechanics tools depicts the most likely bound conformations and hints at differences in the binding mechanism. Finally, the study shows the importance of capturing two driving forces to binding in this system: the ability of peptides to adopt bound conformations (ΔGconformation) and the interaction between interface residues (ΔGinteraction).     

Thermal stability of hybrid nanofluid with viscous dissipation and suction/injection applications: Dual branch framework

The thermal stability of nanomaterials is quite necessary for controlling the heat and cooling phenomenon. It is worthy observed that much research has been focused scientists towards the thermal significance of nanoparticles with multidisciplinary engineering and industrial applications. On this end, this report explores the improved thermal mechanism water base material with interaction of hybrid nanofluid stretching and shrinking surface. The cooling and heat phenomenon is observed in presence of viscous dissipation. The hybrid nanofluid characteristics are inspected with combination of copper $\left(Cu\right)$ and aluminum oxide $\left(A{l}_2{O}_3\right)$ nanoparticles with stable prospective. The consideration of such hybrid nanoparticles is due to impressive thermal characteristics and stable thermal performances. Although some studies are focused by researchers on hybrid nanofluid, however the measurement of thermal stability is not claimed yet. The stretching and shrinking configuration specify the porous medium features. The problem is compiled into dimensionless structure which is further preceded via bvp4c scheme. The resultant ODEs are successfully numerically solved using the bvp4c solver technique. Under restricting conditions, numerical findings are compared to previously published results. Non-dimensional profiles of velocity and temperature are shown graphically. Furthermore, graphs and tables show the effects of the physical parameters used on the reduced skin friction and heat transfer rate. Dual branches are found in specified domain of suction factor.     

Effects of Total Thermal Balance on the Thermal Energy Absorbed or Released by a High-Temperature Phase Change Material

Front tracking and enthalpy methods used to study phase change processes are based on a local thermal energy balance at the liquid–solid interface where mass accommodation methods are also used to account for the density change during the phase transition. Recently, it has been shown that a local thermal balance at the interface does not reproduce the thermodynamic equilibrium in adiabatic systems. Total thermal balance through the entire liquid–solid system can predict the correct thermodynamic equilibrium values of melted (solidified) mass, system size, and interface position. In this work, total thermal balance is applied to systems with isothermal–adiabatic boundary conditions to estimate the sensible and latent heat stored (released) by KNO3 and KNO3/NaNO3 salts which are used as high-temperature phase change materials. Relative percent differences between the solutions obtained with a local thermal balance at the interface and a total thermal balance for the thermal energy absorbed or released by high-temperature phase change materials are obtained. According to the total thermal balance proposed, a correction to the liquid–solid interface dynamics is introduced, which accounts for an extra amount of energy absorbed or released during the phase transition. It is shown that melting or solidification rates are modified by using a total thermal balance through the entire system. Finally, the numerical and semi-analytical methods illustrate that volume changes and the fraction of melted (solidified) solid (liquid) estimated through a local thermal balance at the interface are not invariant in adiabatic systems. The invariance of numerical and semi-analytical solutions in adiabatic systems is significantly improved through the proposed model.     

Kinetic Study on Alpha-Form Crystallization of Mixed-Acid Triacylglycerols POP,PPO, and Their Mixture

The crystallization behavior of the metastable α form of triacylglycerols (TAGs) plays a critical role as a precursor for the crystallization of more stable β′ and β forms for various applications in food and pharmaceutical products. However, precise analysis of the crystallization kinetics of α has not been performed, likely due to its rapid and complex behavior. This paper presents the observation results of the initial stages of the isothermal crystallization kinetics of α forms of 1,3-dipalmitoyl-2-oleoyl-glycerol (POP), 1,2-dipalmitoyl-3-oleoyl-rac-glycerol (rac-PPO), and molecular compound (MC) crystals of a POP/rac-PPO (1/1) mixture (MCPOP/PPO) using synchrotron radiation time-resolved X-ray diffraction and polarized optical microscopy. In all the TAGs, α crystals with a worm-like morphology started to grow rapidly in the first stage. Then, the α crystals slowly transformed into more stable forms in different manners for different TAG samples. In POP, the conversion was simple, as the α-2 form transformed into γ-3, whereas in rac-PPO, the lamellar distance values of the α-2 form continuously decreased with time and changed into the α-3 form. In the MCPOP/PPO crystals, in contrast, separate crystallization of α-2 of a rac-PPO fraction initially occurred, followed by the crystallization of α-2 of POP, and the two α forms merged into α-2 of MCPOP/PPO. This separate crystallization was caused by large differences in the crystallization kinetics of the α forms of POP and rac-PPO.     

Catechin Photolysis Suppression by Aluminum Chloride under Alkaline Conditions and Assessment with Liquid Chromatography–Mass Spectrometry

Tea is rich in catechins and aluminum. In this study, the process of catechin photolysis was applied as a model for examining the effects of aluminum chloride (AlCl₃) on the structural changes of catechin and the alteration of aluminum complexes under blue light irradiation (BLI) at pH 8 using liquid chromatography and mass spectrometry techniques. Additionally, the effects of anions on catechin upon the addition of AlCl₃ and treatment with BLI were also studied. In this study, when 1 mM catechin was treated with BLI, a superoxide anion radical (O2•−) was generated in an air-saturated aqueous solution, in addition to forming a dimeric catechin (proanthocyanidin) via a photon-induced redox reaction. The relative percentage of catechin was found to be 59.0 and 95.7 for catechin treated with BLI and catechin upon the addition of 1 mM AlCl₃ treated with BLI, respectively. It suggested that catechin treated with BLI could be suppressed by AlCl₃, while AlCl₃ did not form a complex with catechin in the photolytic system. However, under the same conditions, it was also found that the addition of AlCl₃ inhibited the photolytic formation of O2•−, and reduced the generation of proanthocyanidin, suggesting that the disconnection of proanthocyanidin was achieved by AlCl₃ acting as a catalyst under treatment with BLI. The influence of 1 mM fluoride (F−) and 1 mM oxalate (C2O42−) ions on the photolysis of 1 mM catechin upon the addition of 1 mM AlCl₃ and treatment with BLI was found to be insignificant, implying that, during the photolysis of catechin, the Al species were either neutral or negatively charged and the aluminum species did not form a complex with anions in the photolytic system. Therefore, aluminum, which is an amphoteric species, has an inherent potential to stabilize the photolysis of catechin in an alkaline conditions, while suppressing the O2•− and proanthocyanidin generation via aluminum ion catalysis in the catechin/Al system under treatment with BLI.     

Two Intercalation Mechanisms of Oxazole Yellow Dimer (YOYO-1) into DNA

The oxazole yellow dye, YOYO-1 (a symmetric homodimer), is a commonly used molecule for staining DNA. We applied the brightness analysis to study the intercalation of YOYO-1 into the DNA. We distinguished two binding modes of the dye to dsDNA: mono-intercalation and bis-intercalation. Bis-intercalation consists of two consecutive mono-intercalation steps, characterised by two distinct equilibrium constants (with the average number of base pair per binding site equals 3.5): K1=3.36±0.43×107M−1 and K2=1.90±0.61×105M−1, respectively. Mono-intercalation dominates at high concentrations of YOYO-1. Bis-intercalation occurs at low concentrations.     

Electrochemical Resistive DNA Biosensor for the Detection of HPV Type 16

In this work, a low-cost and rapid electrochemical resistive DNA biosensor based on the current relaxation method is described. A DNA probe, complementary to the specific human papillomavirus type 16 (HPV-16) sequence, was immobilized onto a screen-printed gold electrode. DNA hybridization was detected by applying a potential step of 30 mV to the system, composed of an external capacitor and the modified electrode DNA/gold, for 750 µs and then relaxed back to the OCP, at which point the voltage and current discharging curves are registered for 25 ms. From the discharging curves, the potential and current relaxation were evaluated, and by using Ohm’s law, the charge transfer resistance through the DNA-modified electrode was calculated. The presence of a complementary sequence was detected by the change in resistance when the ssDNA is transformed in dsDNA due to the hybridization event. The target DNA concentration was detected in the range of 5 to 20 nM. The results showed a good fit to the regression equation ΔRtotal(Ω)=2.99 × [DNA]+81.55, and a detection limit of 2.39 nM was obtained. As the sensing approach uses a direct current, the electronic architecture of the biosensor is simple and allows for the separation of faradic and nonfaradaic contributions. The simple electrochemical resistive biosensor reported here is a good candidate for the point-of-care diagnosis of HPV at a low cost and in a short detection time.     

Carbon Atoms Speaking Out: How the Geometric Sensitivity of 13C Chemical Shifts Leads to Understanding the Colour Tuning of Phycocyanobilin in Cph1 and AnPixJ

We present a combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics–statistical approach for the interpretation of nuclear magnetic resonance (NMR) chemical shift patterns in phycocyanobilin (PCB). These were originally associated with colour tuning upon photoproduct formation in red/green-absorbing cyanobacteriochrome AnPixJg2 and red/far-red-absorbing phytochrome Cph1Δ2. We pursue an indirect approach without computation of the absorption frequencies since the molecular geometry of cofactor and protein are not accurately known. Instead, we resort to a heuristic determination of the conjugation length in PCB through the experimental NMR chemical shift patterns, supported by quantum chemical calculations. We have found a characteristic correlation pattern of 13C chemical shifts to specific bond orders within the π-conjugated system, which rests on the relative position of carbon atoms with respect to electron-withdrawing groups and the polarisation of covalent bonds. We propose the inversion of this regioselective relationship using multivariate statistics and to apply it to the known experimental NMR chemical shifts in order to predict changes in the bond alternation pattern. Therefrom the extent of electronic conjugation, and eventually the change in absorption frequency, can be derived. In the process, the consultation of explicit mesomeric formulae plays an important role to qualitatively account for possible conjugation scenarios of the chromophore. While we are able to consistently associate the NMR chemical shifts with hypsochromic and bathochromic shifts in the Pg and Pfr, our approach represents an alternative method to increase the explanatory power of NMR spectroscopic data in proteins.     

An Effective MM/GBSA Protocol for Absolute Binding Free Energy Calculations: A Case Study on SARS-CoV-2 Spike Protein and the Human ACE2 Receptor

The binding free energy calculation of protein–ligand complexes is necessary for research into virus–host interactions and the relevant applications in drug discovery. However, many current computational methods of such calculations are either inefficient or inaccurate in practice. Utilizing implicit solvent models in the molecular mechanics generalized Born surface area (MM/GBSA) framework allows for efficient calculations without significant loss of accuracy. Here, GBNSR6, a new flavor of the generalized Born model, is employed in the MM/GBSA framework for measuring the binding affinity between SARS-CoV-2 spike protein and the human ACE2 receptor. A computational protocol is developed based on the widely studied Ras–Raf complex, which has similar binding free energy to SARS-CoV-2/ACE2. Two options for representing the dielectric boundary of the complexes are evaluated: one based on the standard Bondi radii and the other based on a newly developed set of atomic radii (OPT1), optimized specifically for protein–ligand binding. Predictions based on the two radii sets provide upper and lower bounds on the experimental references: −14.7(ΔGbindBondi)<−10.6(ΔGbindExp.)<−4.1(ΔGbindOPT1) kcal/mol. The consensus estimates of the two bounds show quantitative agreement with the experiment values. This work also presents a novel truncation method and computational strategies for efficient entropy calculations with normal mode analysis. Interestingly, it is observed that a significant decrease in the number of snapshots does not affect the accuracy of entropy calculation, while it does lower computation time appreciably. The proposed MM/GBSA protocol can be used to study the binding mechanism of new variants of SARS-CoV-2, as well as other relevant structures.     

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.     

A Computational Study of the S2 State in the Oxygen-Evolving Complex of Photosystem II by Electron Paramagnetic Resonance Spectroscopy

The S₂ state produces two basic electron paramagnetic resonance signal types due to the manganese cluster in oxygen-evolving complex, which are influenced by the solvents, and cryoprotectant added to the photosystem II samples. It is presumed that a single manganese center oxidation occurs on S₁ → S₂ state transition. The S₂ state has readily visible multiline and g4.1 electron paramagnetic resonance signals and hence it has been the most studied of all the Kok cycle intermediates due to the ease of experimental preparation and stability. The S₂ state was studied using electron paramagnetic resonance spectroscopy at X-band frequencies. The aim of this study was to determine the spin states of the g4.1 signal. The multiline signal was observed to arise from a ground state spin ½ centre while the g4.1 signal generated at ≈140 K NIR illumination was proposed to arise from a spin 52 center with rhombic distortion. The ‘ground’ state g4.1 signal was generated solely or by conversion from the multiline. The data analysis methods used involved numerical simulations of the experimental spectra on relevant models of the oxygen-evolving complex cluster. A strong focus in this paper was on the ‘ground’ state g4.1 signal, whether it is a rhombic 52 spin state signal or an axial 32 spin state signal. The data supported an X-band CW-EPR-generated g4.1 signal as originating from a near rhombic spin 5/2 of the S₂ state of the PSII manganese cluster.