Exchange pathways of plastoquinone and plastoquinol in the photosystem II complex

6PPD-quinone affects the photosynthetic carbon fixation in cyanobacteria by extracting photosynthetic electrons

Photosynthetic carbon fixation by cyanobacteria plays a pivotal role in the global carbon cycle but is threatened by environmental pollutants. To date, the impact of quinones, with electron shuttling properties, on cyanobacterial photosynthesis is unknown. Here, we present the first study investigating the effects of an emerging quinone pollutant, i.e., 6PPD-Q (N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine-quinone), on the cyanobacterium Synechocystis sp. over a 400-generation exposure period. Synechocystis sp. exhibited distinct sequential phases, including hormesis, toxicity, and eventual recovery, throughout this exposure. Extensive evidence, including results of thylakoid membrane morphological and photosynthetic responses, carbon fixation rate, and key gene/protein analyses, strongly indicates that 6PPD-Q is a potent disruptor of photosynthesis. 6PPD-Q accepts photosynthetic electrons at the plastoquinone QB site in photosystem II (PSII) and the phylloquinone A1 site in PSI, leading to a sustained decrease in the carbon fixation of cyanobacteria after an ephemeral increase. This work revealed the specific mechanism by which 6PPD-Q interferes with photosynthetic carbon fixation in cyanobacteria, which is highly important for the global carbon cycle.

Tuning the Thermochemistry and Reactivity of a Series of Cu-Based 4H+/4e- Electron-Coupled-Proton Buffers

Electron-coupled-proton buffers (ECPBs) store and deliver protons and electrons in a reversible fashion. We have recently reported an ECPB based on Cu and a redox-active ligand that promoted 4H+/4e- reversible transformations (J. Am. Chem. Soc. 2022, 144, 16905). Herein, we report a series of Cu-based ECPBs in which the ability of these to accept and/or donate H• equivalents can be tuned via ligand modification. The thermochemistry of the 4H+/4e- ECPB equilibrium was determined using open-circuit potential measurements. The reactivity of the ECPBs against proton-coupled electron transfer (PCET) reagents was also analyzed, and the results obtained were rationalized based on the thermochemical parameters. Experimental and computational analysis of the thermochemistry of the H+/e- transfers involved in the 4H+/4e- ECPB transformations found substantial differences between the stepwise (namely, BDFE1, BDFE2, BDFE3, and BDFE4) and average bond dissociation free energy values (BDFEavg.). Our analysis suggests that this "redox unleveling" is critical to promoting the disproportionation and ligand-exchange reactions involved in the 4H+/4e- ECPB equilibria. The difference in BDFEavg. within the series of Cu-based ECPBs was found to arise from a substantial change in the redox potential (E1/2) upon modification of the ligand scaffold, which is not fully compensated for by a change in the acidity/basicity (pKa), suggesting "thermochemical decompensation".

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Modulation of the wheat transcriptome by TaZFP13D under well-watered and drought conditions

Maintaining global food security in the context of climate changes will be an important challenge in the next century. Improving abiotic stress tolerance of major crops such as wheat can contribute to this goal. This can be achieved by the identification of the genes involved and their use to develop tools for breeding programs aiming to generate better adapted cultivars. Recently, we identified the wheat TaZFP13D gene encoding Zinc Finger Protein 13D as a new gene improving water-stress tolerance. The current work analyzes the TaZFP13D-dependent transcriptome modifications that occur in well-watered and dehydration conditions to better understand its function during normal growth and during drought. Plants that overexpress TaZFP13D have a higher biomass under well-watered conditions, indicating a positive effect of the protein on growth. Survival rate and stress recovery after a severe drought stress are improved compared to wild-type plants. The latter is likely due the higher activity of key antioxidant enzymes and concomitant reduction of drought-induced oxidative damage. Conversely, down-regulation of TaZFP13D decreases drought tolerance and protection against drought-induced oxidative damage. RNA-Seq transcriptome analysis identified many genes regulated by TaZFP13D that are known to improve drought tolerance. The analysis also revealed several genes involved in the photosynthetic electron transfer chain known to improve photosynthetic efficiency and chloroplast protection against drought-induced ROS damage. This study highlights the important role of TaZFP13D in wheat drought tolerance, contributes to unravel the complex regulation governed by TaZFPs, and suggests that it could be a promising marker to select wheat cultivars with higher drought tolerance.

Lys264 of the D2 Protein Performs a Dual Role in Photosystem II Modifying Assembly and Electron Transfer through the Quinone-Iron Acceptor Complex

Bicarbonate (HCO3-) binding regulates electron flow between the primary (QA) and secondary (QB) plastoquinone electron acceptors of Photosystem II (PS II). Lys264 of the D2 subunit of PS II contributes to a hydrogen-bond network that stabilizes HCO3- ligation to the non-heme iron in the QA-Fe-QB complex. Using the cyanobacterium Synechocystis sp. PCC 6803, alanine and glutamate were introduced to create the K264A and K264E mutants. Photoautotrophic growth was slowed in K264E cells but not in the K264A strain. Both mutants accumulated an unassembled CP43 precomplex as well as the CP43-lacking RC47 assembly intermediate, indicating weakened binding of the CP43 precomplex to RC47. Assembly was impeded more in K264E cells than in the K264A strain, but K264A cells were more susceptible to high-light-induced photodamage when assayed using PS II-specific electron acceptors. Furthermore, an impaired repair mechanism was observed in the K264A mutant in protein labeling experiments. Unexpectedly, unlike the K264A strain, the K264E mutant displayed inhibited oxygen evolution following high-light exposure when HCO3- was added to support whole chain electron transport. In both mutants, the decay of chlorophyll fluorescence was slowed, indicating impaired electron transfer between QA and QB. Furthermore, the fluorescence decay kinetics in the K264E strain were insensitive to addition of either formate or HCO3-, whereas HCO3--reversible formate-induced inhibition in the K264A mutant was observed. Exchange of plastoquinol with the membrane plastoquinone pool at the QB-binding site was also retarded in both mutants. Hence, D2-Lys264 possesses key roles in both assembly and activity of PS II.

Proteome of Agave angustifolia Haw.: Uncovering metabolic alterations, over-accumulation of amino acids, and compensatory pathways in chloroplast-deficient albino plantlets

Amino acids (AA) are essential molecules for plant physiology, acting as precursor molecules for proteins and other organic compounds. Chloroplasts play a vital role in AA metabolism, yet little is known about the impact on AA metabolism of albino plants' lack of chloroplasts. In this study, we conducted a quantitative proteome analysis on albino and variegated somaclonal variants of Agave angustifolia Haw. to investigate metabolic alterations in chloroplast-deficient plants, with a focus on AA metabolic pathways. We identified 82 enzymes involved in AA metabolism, with 32 showing differential accumulation between the somaclonal variants. AaCM, AaALS, AaBCAT, AaIPMS1, AaSHMT, AaAST, AaCGS, and AaMS enzymes were particularly relevant in chloroplast-deficient Agave plantlets. Both variegated and albino phenotypes exhibited excessive synthesis of AA typically associated with chloroplasts (aromatic AAs, BCAAs, Asp, Lys, Pro and Met). Consistent trends were observed for AaBCAT and AaCM at mRNA and protein levels in albino plantlets. These findings highlight the critical activation and reprogramming of AA metabolic pathways in plants lacking chloroplasts. This study contributes to unraveling the intricate relationship between AA metabolism and chloroplast absence, offering insights into survival mechanisms of albino plants.

Structure of Chlorella ohadii Photosystem II Reveals Protective Mechanisms against Environmental Stress

Green alga Chlorella ohadii is known for its ability to carry out photosynthesis under harsh conditions. Using cryogenic electron microscopy (cryoEM), we obtained a high-resolution structure of PSII at 2.72 Å. This structure revealed 64 subunits, which encompassed 386 chlorophylls, 86 carotenoids, four plastoquinones, and several structural lipids. At the luminal side of PSII, a unique subunit arrangement was observed to protect the oxygen-evolving complex. This arrangement involved PsbO (OEE1), PsbP (OEE2), PsbB, and PsbU (a homolog of plant OEE3). PsbU interacted with PsbO, PsbC, and PsbP, thereby stabilizing the shield of the oxygen-evolving complex. Significant changes were also observed at the stromal electron acceptor side. PsbY, identified as a transmembrane helix, was situated alongside PsbF and PsbE, which enclosed cytochrome b559. Supported by the adjacent C-terminal helix of Psb10, these four transmembrane helices formed a bundle that shielded cytochrome b559 from the surrounding solvent. Moreover, the bulk of Psb10 formed a protective cap, which safeguarded the quinone site and likely contributed to the stacking of PSII complexes. Based on our findings, we propose a protective mechanism that prevents QB (plastoquinone B) from becoming fully reduced. This mechanism offers insights into the regulation of electron transfer within PSII.

Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

Solar energy conversion by photosystem II: principles and structures

Photosynthetic water oxidation by Photosystem II (PSII) is a fascinating process because it sustains life on Earth and serves as a blue print for scalable synthetic catalysts required for renewable energy applications. The biophysical, computational, and structural description of this process, which started more than 50 years ago, has made tremendous progress over the past two decades, with its high-resolution crystal structures being available not only of the dark-stable state of PSII, but of all the semi-stable reaction intermediates and even some transient states. Here, we summarize the current knowledge on PSII with emphasis on the basic principles that govern the conversion of light energy to chemical energy in PSII, as well as on the illustration of the molecular structures that enable these reactions. The important remaining questions regarding the mechanism of biological water oxidation are highlighted, and one possible pathway for this fundamental reaction is described at a molecular level.

Structural insights into the action mechanisms of artificial electron acceptors in photosystem II

Photosystem II (PSII) utilizes light energy to split water, and the electrons extracted from water are transferred to Q_B, a plastoquinone (PQ) molecule bound to the D1 subunit of PSII. Many artificial electron acceptors (AEAs) with similar molecular structures to PQ can accept electrons from PSII. However, the molecular mechanism by which AEAs act on PSII is unclear. Here, we solved the crystal structure of PSII treated with three different AEAs, 2,5-dibromo-1,4-benzoquinone, 2,6-dichloro-1,4-benzoquinone, and 2-phenyl-1,4-benzoquinone, at 1.95-2.10 Å resolution. Our results show that all AEAs substitute for Q_B and are bound to the Q_B-binding site (Q_B site) to receive electrons, but their binding strengths are different, resulting in differences in their efficiencies to accept electrons. The acceptor 2-phenyl-1,4-benzoquinone binds most weakly to the Q_B site, and showed the highest oxygen-evolving activity, implying a reverse relationship between the binding strength and oxygen-evolving activity. In addition, a novel quinone binding site, designated the Q_D site, was discovered, which is located in the vicinity of Q_B site and close to Q_C site, a binding site reported previously. This Q_D site is expected to play a role as a channel or a storage site for quinones to be transported to the Q_B site. These results provide the structural basis for elucidating the actions of AEAs and exchange mechanism of Q_B in PSII, and also provide information for the design of more efficient electron acceptors.

CryoEM PSII structure reveals adaptation mechanisms to environmental stress in Chlorella ohadii

Performing photosynthesis in the desert is a challenging task since it requires a fast adaptation to extreme illumination and temperature changes. To understand adaptive mechanisms, we purified Photosystem II (PSII) from Chlorella ohadii , a green alga from the desert soil surface, and identified structural elements that might enable the photosystem functioning under harsh conditions. The 2.72 Å cryogenic electron-microscopy (cryoEM) structure of PSII exhibited 64 subunits, encompassing 386 chlorophylls, 86 carotenoids, four plastoquinones, and several structural lipids. At the luminal side of PSII, the oxygen evolving complex was protected by a unique subunit arrangement - PsbO (OEE1), PsbP (OEE2), CP47, and PsbU (plant OEE3 homolog). PsbU interacted with PsbO, CP43, and PsbP, thus stabilising the oxygen evolving shield. Substantial changes were observed on the stromal electron acceptor side - PsbY was identified as a transmembrane helix situated alongside PsbF and PsbE enclosing cytochrome b559, supported by the adjacent C-terminal helix of Psb10. These four transmembrane helices bundled jointly, shielding cytochrome b559 from the solvent. The bulk of Psb10 formed a cap protecting the quinone site and probably contributed to the PSII stacking. So far, the C. ohadii PSII structure is the most complete description of the complex, suggesting numerous future experiments. A protective mechanism that prevented Q B from rendering itself fully reduced is proposed.

In situ structure of the red algal phycobilisome-PSII-PSI-LHC megacomplex

In oxygenic photosynthetic organisms, light energy is captured by antenna systems and transferred to photosystem II (PSII) and photosystem I (PSI) to drive photosynthesis1,2. The antenna systems of red algae consist of soluble phycobilisomes (PBSs) and transmembrane light-harvesting complexes (LHCs)3. Excitation energy transfer pathways from PBS to photosystems remain unclear owing to the lack of structural information. Here we present in situ structures of PBS-PSII-PSI-LHC megacomplexes from the red alga Porphyridium purpureum at near-atomic resolution using cryogenic electron tomography and in situ single-particle analysis4, providing interaction details between PBS, PSII and PSI. The structures reveal several unidentified and incomplete proteins and their roles in the assembly of the megacomplex, as well as a huge and sophisticated pigment network. This work provides a solid structural basis for unravelling the mechanisms of PBS-PSII-PSI-LHC megacomplex assembly, efficient energy transfer from PBS to the two photosystems, and regulation of energy distribution between PSII and PSI.

Isolation of a novel heterodimeric PSII complex via strep-tagged PsbO

The multi-subunit membrane protein complex photosystem II (PSII) catalyzes the light-driven oxidation of water and with this the initial step of photosynthetic electron transport in plants, algae, and cyanobacteria. Its biogenesis is coordinated by a network of auxiliary proteins that facilitate the stepwise assembly of individual subunits and cofactors, forming various intermediate complexes until fully functional mature PSII is present at the end of the process. In the current study, we purified PSII complexes from a mutant line of the thermophilic cyanobacterium Thermosynechococcus vestitus BP-1 in which the extrinsic subunit PsbO, characteristic for active PSII, was fused with an N-terminal Twin-Strep-tag. Three distinct PSII complexes were separated by ion-exchange chromatography after the initial affinity purification. Two complexes differ in their oligomeric state (monomeric and dimeric) but share the typical subunit composition of mature PSII. They are characterized by the very high oxygen evolving activity of approx. 6000 μmol O₂·(mg Chl·h)^-1. Analysis of the third (heterodimeric) PSII complex revealed lower oxygen evolving activity of approx. 3000 μmol O₂·(mg Chl·h)^-1 and a manganese content of 2.7 (±0.2) per reaction center compared to 3.7 (±0.2) of fully active PSII. Mass spectrometry and time-resolved fluorescence spectroscopy further indicated that PsbO is partially replaced by Psb27 in this PSII fraction, thus implying a role of this complex in PSII repair.

Using supramolecular machinery to engineer directional charge propagation in photoelectrochemical devices

Molecular photoelectrochemical devices are hampered by electron-hole recombination after photoinduced electron transfer, causing losses in power conversion efficiency. Inspired by natural photosynthesis, we demonstrate the use of supramolecular machinery as a strategy to inhibit recombination through an organization of molecular components that enables unbinding of the final electron acceptor upon reduction. We show that preorganization of a macrocyclic electron acceptor to a dye yields a pseudorotaxane that undergoes a fast (completed within ~50 ps) 'ring-launching' event upon electron transfer from the dye to the macrocycle, releasing the anionic macrocycle and thus reducing charge recombination. Implementing this system into p-type dye-sensitized solar cells yielded a 16-fold and 5-fold increase in power conversion efficiency compared to devices based on the two control dyes that are unable to facilitate pseudorotaxane formation. The active repulsion of the anionic macrocycle with concomitant reformation of a neutral pseudorotaxane complex circumvents recombination at both the semiconductor-electrolyte and semiconductor-dye interfaces, enabling a threefold enhancement in hole lifetime.

Approaches and Perspective of Coarse-Grained Modeling and Simulation for Polymer-Nanoparticle Hybrid Systems

Molecular modeling and simulations have emerged as effective and indispensable tools to characterize polymeric systems. They provide fundamental and essential insights to design a product of the required properties and to improve the understanding of a phenomenon at the molecular level for a particular system. The polymer-nanoparticle hybrids are materials with outstanding properties and correspondingly large applications whose study has benefited from this new paradigm. However, despite the significant expansion of modern day computational powers, investigation of the long time and large length scale phenomenon in polymeric and polymer-nanoparticle systems is still a challenging task to complete through all-atom molecular dynamics (AA-MD) simulations. To circumvent this problem, a variety of coarse-grained (CG) models have been proposed, ranging from the generic CG models for qualitative properties predictions to more realistic chemically specific CG models for quantitative properties predictions. These CG models have already delivered some success stories in the study of several spatial and temporal evolutions of many processes. Some of these studies were beyond the feasibility of traditional atomistic resolution models due to either the size or the time constraints. This review captures the different types of popular CG approaches that are utilized in the investigation of the microscopic behavior of polymer-nanoparticle hybrid systems. The rationale of this article is to furnish an overview of the popular CG approaches and their applications, to review several important and most recent developments, and to delineate the perspectives on future directions in the field.

Improved Coarse-Grained Model for Nanoparticles Based on the Martini Force Field and Its Application in Molecular Dynamics Simulation on Gel Ink

Research on the dispersion and stability of nanoparticles in liquid media is one of the key subjects for nanomaterial utilization. In consideration of the preponderance of coarse-grained (CG) molecular dynamics (MD) simulation in following and understanding the structure and dynamics on the nanoscale, an improved CG model for nanoparticles based on the Martini force field is established to facilitate the more extensive applications of this simulation method and further studies on complex nanoparticle liquid systems. Gel ink is selected as the liquid system for nanoparticles to validate the improved CG model on the one hand and introduce the CGMD simulation method into the studies of this system on the other. The calculation shows that the improved model can provide relatively precise results and has good computational stability. The effect mechanisms of the thickener and disperser on the carbon black nanoparticle are similar, namely the result of a delicate balance between the interaction of the thickener/disperser with the carbon black nanoparticle and the interaction of the thickener and disperser with each other. Furthermore, the phase assimilating effect of disperser molecules is key for separating the agglomerated carbon black nanoparticles; thereafter, the space steric hindrance effect and the electrostatic hindrance effect play main roles in maintaining the dispersion of carbon black nanoparticles.

Gut bacteriome and metabolome of Ascaris lumbricoides in patients

The most frequent intestinal helminth infections in humans are attributed to Ascaris lumbricoides, and there are concerns over the anthelminthic resistance of this species. The gut microbiota has essential roles in host physiology. Therefore, discovering host-parasite-microbiota interactions could help develop alternative helminthiasis treatments. Additionally, these interactions are modulated by functional metabolites that can reveal the mechanisms of infection and disease progression. Thus, we aimed to investigate bacteriomes in the gut of helminths and fecal samples of patients via next-generation sequencing. Our results showed that infection intensity was associated with the bacterial composition of helminth guts but not with the intestinal bacteriome of human hosts. Moreover, the metabolomes of A. lumbricoides in the heavy and light ascariasis cases were characterized using ultra-high performance liquid chromatography/time-of-flight mass spectrometry. Increased levels of essential biomolecules, such as amino acids, lipids, and nucleotide precursors, were found in the guts of helminths isolated from heavily infected patients, implying that these metabolites are related to egg production and ascariasis pathogenicity. These findings are the first step towards a more complete understanding of the mechanisms by which the bacteriome of helminth guts affect their colonization and may reveal novel and more effective approaches to parasitic disease therapy.

Chemistry of Lipoquinones: Properties, Synthesis, and Membrane Location of Ubiquinones, Plastoquinones, and Menaquinones

Lipoquinones are the topic of this review and are a class of hydrophobic lipid molecules with key biological functions that are linked to their structure, properties, and location within a biological membrane. Ubiquinones, plastoquinones, and menaquinones vary regarding their quinone headgroup, isoprenoid sidechain, properties, and biological functions, including the shuttling of electrons between membrane-bound protein complexes within the electron transport chain. Lipoquinones are highly hydrophobic molecules that are soluble in organic solvents and insoluble in aqueous solution, causing obstacles in water-based assays that measure their chemical properties, enzyme activities and effects on cell growth. Little is known about the location and ultimately movement of lipoquinones in the membrane, and these properties are topics described in this review. Computational studies are particularly abundant in the recent years in this area, and there is far less experimental evidence to verify the often conflicting interpretations and conclusions that result from computational studies of very different membrane model systems. Some recent experimental studies have described using truncated lipoquinone derivatives, such as ubiquinone-2 (UQ-2) and menaquinone-2 (MK-2), to investigate their conformation, their location in the membrane, and their biological function. Truncated lipoquinone derivatives are soluble in water-based assays, and hence can serve as excellent analogs for study even though they are more mobile in the membrane than the longer chain counterparts. In this review, we will discuss the properties, location in the membrane, and syntheses of three main classes of lipoquinones including truncated derivatives. Our goal is to highlight the importance of bridging the gap between experimental and computational methods and to incorporate properties-focused considerations when proposing future studies relating to the function of lipoquinones in membranes.

A 4H+/4e- Electron-Coupled-Proton Buffer Based on a Mononuclear Cu Complex

In this research article, we describe a 4H⁺/4e^- electron-coupled-proton buffer (ECPB) based on Cu and a redox-active ligand. The protonated/reduced ECPB (complex 1: [Cu(8H⁺/14e^-)]¹⁺), consisting of Cu^I with 2 equiv of the ligand (^catLH₄: 1,1'-(4,5-dimethoxy-1,2-phenylene)bis(3-(tert-butyl)urea)), reacted with H⁺/e^- acceptors such as O₂ to generate the deprotonated/oxidized ECPB. The resulting compound, (complex 5: [Cu(4H⁺/10e^-)]¹⁺), was characterized by X-ray diffraction analysis, nuclear magnetic resonance (¹H-NMR), and density functional theory, and it is electronically described as a cuprous bis(benzoquinonediimine) species. The stoichiometric 4H⁺/4e^- reduction of 5 was carried out with H⁺/e^- donors to generate 1 (Cu^I and 2 equiv of ^catLH₄) and the corresponding oxidation products. The 1/5 ECPB system catalyzed the 4H⁺/4e^- reduction of O₂ to H₂O and the dehydrogenation of organic substrates in a decoupled (oxidations and reductions are separated in time and space) and a coupled fashion (oxidations and reductions coincide in time and space). Mechanistic analysis revealed that upon reductive protonation of 5 and oxidative deprotonation of 1, fast disproportionation reactions regenerate complexes 5 and 1 in a stoichiometric fashion to maintain the ECPB equilibrium.

Concerning the enigmatic cytochrome b-559 of oxygenic photosynthesis

Although there is an extensive literature on the properties and possible electron transfer pathways of cytochrome b-559, which is a prominent subunit of the multi-subunit photosystem II complex which functions in oxygenic photosynthesis, there is presently no consensus on the function of b-559 in the photosynthetic electron transport chain. The inability in earlier times to define a redox-linked function of this cytochrome was, to a large extent, a consequence of an absence of biochemical and structure information to complement an extensive array of spectrophotometric studies of the cytochrome in situ. Based on the location of hetero-dimeric b-559 in the photosystem II reaction center complex, derived from crystal crystallographic structure analysis, and the absence of a necessary redox function for the cytochrome in PSII, it is proposed that the main function of cytochrome b-559 is linked to its role as a structure component in the PSII reaction center complex. This function resides in the association of b-559 through its heme histidine residues in the trans-membrane domains of the PsbE and PsbF subunits of the PSII reaction center. These subunits, along with PsbJ, are inferred, from the analysis of structure, to define the intra-membrane portal in the PSII reaction center for plastoquinol (PQH₂) export which, through the PSII complex, provides the redox link to the cytochrome b₆f complex in the electron transfer chain.

Conformational Changes and H-Bond Rearrangements during Quinone Release in Photosystem II

In photosystem II (PSII) and photosynthetic reaction centers from purple bacteria (PbRC), the electron released from the electronically excited chlorophyll is transferred to the terminal electron acceptor quinone, Q_B. Q_B accepts two electrons and two protons before leaving the protein. We investigated the molecular mechanism of quinone exchange in PSII, conducting molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations. MD simulations suggest that the release of Q_B leads to the transformation of the short helix (D1-Phe260 to D1-Ser264), which is adjacent to the stromal helix de (D1-Asn247 to D1-Ile259), into a loop and to the formation of a water-intake channel. Water molecules enter the Q_B binding pocket via the channel and form an H-bond network. QM/MM calculations indicate that the H-bond network serves as a proton-transfer pathway for the reprotonation of D1-His215, the proton donor during Q_BH^-/Q_BH₂ conversion. Together with the absence of the corresponding short helix but the presence of Glu-L212 in PbRC, it seems likely that the two type-II reaction centers undergo quinone exchange via different mechanisms.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Tyr244 of the D2 Protein Is Required for Correct Assembly and Operation of the Quinone-Iron-Bicarbonate Acceptor Complex of Photosystem II

Two plastoquinone electron acceptors, Q_A and Q_B, are present in Photosystem II (PS II) with their binding sites formed by the D2 and D1 proteins, respectively. A hexacoordinate non-heme iron is bound between Q_A and Q_B by D2 and D1, each providing two histidine ligands, and a bicarbonate that is stabilized via hydrogen bonds with D2-Tyr244 and D1-Tyr246. Both tyrosines and bicarbonate are conserved in oxygenic photosynthetic organisms but absent from the corresponding quinone-iron electron acceptor complex of anoxygenic photosynthetic bacteria. We investigated the role of D2-Tyr244 by introducing mutations in the cyanobacterium Synechocystis sp. PCC 6803. Alanine, histidine, and phenylalanine substitutions were introduced creating the Y244A, Y244H, and Y244F mutants. Electron transfer between Q_A and Q_B was impaired, the back-reaction with the S2 state of the oxygen-evolving complex was modified, and PS II assembly was disrupted, with the Y244A strain being more affected than the Y244F and Y244H mutants. The strains were also highly susceptible to photodamage in the presence of PS II-specific electron acceptors. Thermoluminescence and chlorophyll a fluorescence decay measurements indicated that the redox potential of the Q_A/Q_A^- couple became more positive in the Y244F and Y244H mutants, consistent with bicarbonate binding being impacted. The replacement of Tyr244 by alanine also led to an insertion of two amino acid repeats from Gln239 to Ala249 within the DE loop of D2, resulting in an inactive PS II complex that lacked PS II-specific variable fluorescence. The 66 bp insertion giving rise to the inserted amino acids therefore created an obligate photoheterotrophic mutant.

New Structural and Mechanistic Insights Into Functional Roles of Cytochrome b 559 in Photosystem II

Cytochrome (Cyt) b ₅₅₉ is a key component of the photosystem II (PSII) complex for its assembly and proper function. Previous studies have suggested that Cytb ₅₅₉ has functional roles in early assembly of PSII and in secondary electron transfer pathways that protect PSII against photoinhibition. In addition, the Cytb ₅₅₉ in various PSII preparations exhibited multiple different redox potential forms. However, the precise functional roles of Cytb ₅₅₉ in PSII remain unclear. Recent site-directed mutagenesis studies combined with functional genomics and biochemical analysis, as well as high-resolution x-ray crystallography and cryo-electron microscopy studies on native, inactive, and assembly intermediates of PSII have provided important new structural and mechanistic insights into the functional roles of Cytb ₅₅₉. This mini-review gives an overview of new exciting results and their significance for understanding the structural and functional roles of Cytb ₅₅₉ in PSII.

Properties of Photosystem II lacking the PsbJ subunit

Photosystem II (PSII), the oxygen-evolving enzyme, consists of 17 trans-membrane and 3 extrinsic membrane proteins. Other subunits bind to PSII during assembly, like Psb27, Psb28, and Tsl0063. The presence of Psb27 has been proposed (Zabret et al. in Nat Plants 7:524-538, 2021; Huang et al. Proc Natl Acad Sci USA 118:e2018053118, 2021; Xiao et al. in Nat Plants 7:1132-1142, 2021) to prevent the binding of PsbJ, a single transmembrane α-helix close to the quinone Q_B binding site. Consequently, a PSII rid of Psb27, Psb28, and Tsl0034 prior to the binding of PsbJ would logically correspond to an assembly intermediate. The present work describes experiments aiming at further characterizing such a ∆PsbJ-PSII, purified from the thermophilic Thermosynechococcus elongatus, by means of MALDI-TOF spectroscopy, thermoluminescence, EPR spectroscopy, and UV-visible time-resolved spectroscopy. In the purified ∆PsbJ-PSII, an active Mn₄CaO₅ cluster is present in 60-70% of the centers. In these centers, although the forward electron transfer seems not affected, the Em of the Q_B/Q_B^- couple increases by ≥ 120 mV , thus disfavoring the electron coming back on Q_A. The increase of the energy gap between Q_A/Q_A^- and Q_B/Q_B^- could contribute in a protection against the charge recombination between the donor side and Q_B^-, identified at the origin of photoinhibition under low light (Keren et al. in Proc Natl Acad Sci USA 94:1579-1584, 1997), and possibly during the slow photoactivation process.

More than just lipid balls: quantitative analysis of plastoglobule attributes and their stress-related responses

Plastoglobules are ubiquitous under non-stress conditions and their morphology, closely related to their composition, changes differently depending on the specific stress that the plant undergoes. Plastoglobules are lipoprotein structures attached to thylakoid membranes, which participate in chloroplast metabolism and stress responses. Their structure contains a coating lipid monolayer and a hydrophobic core that differ in composition. Their function in chloroplasts has been studied focussing on their composition. However, we currently lack a comprehensive study that quantitatively evaluates the occurrence and morphology of plastoglobules. Following a literature search strategy, we quantified the main morphological attributes of plastoglobules from photosynthetic chloroplasts of more than 1000 TEM images published over the last 53 years, covering more than 100 taxa and 15 stress types. The analysis shows that plastoglobules under non-stress conditions are spherical, with an average diameter of 100-200 nm and cover less than 3% of the chloroplast cross-section area. This percentage rises under almost every type of stress, particularly in senescence. Interestingly, an apparent trade-off between increasing either the number or the diameter of plastoglobules governs this response. Our results show that plastoglobules are ubiquitous in chloroplasts of higher plants under non-stress conditions. Besides, provided the specific molecular composition of the core and coat of plastoglobules, we conclude that specific stress-related variation in plastoglobules attributes may allow inferring precise responses of the chloroplast metabolism.

PsbX maintains efficient electron transport in Photosystem II and reduces susceptibility to high light in Synechocystis sp. PCC 6803

PsbX is a 4.1 kDa intrinsic Photosystem II (PS II) protein, found together with the low-molecular-weight proteins, PsbY and PsbJ, in proximity to cytochrome b₅₅₉. The function of PsbX is not yet fully characterized but PsbX may play a role in the exchange of the secondary plastoquinone electron acceptor Q_B with the quinone pool in the thylakoid membrane. To study the role of PsbX, we have constructed a PsbX-lacking strain of Synechocystis sp. PCC 6803. Our studies indicate that the absence of PsbX causes sensitivity to high light and impairs electron transport within PS II. In addition to a change in the Q_B-binding pocket, PsbX-lacking cells exhibited sensitivity to sodium formate, suggesting altered binding of the bicarbonate ligand to the non-heme iron between the sequential plastoquinone electron acceptors Q_A and Q_B. Experiments using ³⁵S-methionine revealed high-light-treated PsbX-lacking cells restore PS II activity during recovery under low light by an increase in the turnover of PS II-associated core proteins. These labeling experiments indicate the recovery after exposure to high light requires both selective removal and replacement of the D1 protein and de novo PS II assembly.

Automatic multi-objective optimization of coarse-grained lipid force fields using SwarmCG

The development of coarse-grained (CG) molecular models typically requires a time-consuming iterative tuning of parameters in order to have the approximated CG models behave correctly and consistently with, e.g., available higher-resolution simulation data and/or experimental observables. Automatic data-driven approaches are increasingly used to develop accurate models for molecular dynamics simulations. However, the parameters obtained via such automatic methods often make use of specifically designed interaction potentials and are typically poorly transferable to molecular systems or conditions other than those used for training them. Using a multi-objective approach in combination with an automatic optimization engine (SwarmCG), here, we show that it is possible to optimize CG models that are also transferable, obtaining optimized CG force fields (FFs). As a proof of concept, here, we use lipids for which we can avail reference experimental data (area per lipid and bilayer thickness) and reliable atomistic simulations to guide the optimization. Once the resolution of the CG models (mapping) is set as an input, SwarmCG optimizes the parameters of the CG lipid models iteratively and simultaneously against higher-resolution simulations (bottom-up) and experimental data (top-down references). Including different types of lipid bilayers in the training set in a parallel optimization guarantees the transferability of the optimized lipid FF parameters. We demonstrate that SwarmCG can reach satisfactory agreement with experimental data for different resolution CG FFs. We also obtain stimulating insights into the precision-resolution balance of the FFs. The approach is general and can be effectively used to develop new FFs and to improve the existing ones.

Resolving Chemical Dynamics in Biological Energy Conversion: Long-Range Proton-Coupled Electron Transfer in Respiratory Complex I

Biological energy conversion is catalyzed by membrane-bound proteins that transduce chemical or light energy into energy forms that power endergonic processes in the cell. At a molecular level, these catalytic processes involve elementary electron-, proton-, charge-, and energy-transfer reactions that take place in the intricate molecular machineries of cell respiration and photosynthesis. Recent developments in structural biology, particularly cryo-electron microscopy (cryoEM), have resolved the molecular architecture of several energy transducing proteins, but detailed mechanistic principles of their charge transfer reactions still remain poorly understood and a major challenge for modern biochemical research. To this end, multiscale molecular simulations provide a powerful approach to probe mechanistic principles on a broad range of time scales (femtoseconds to milliseconds) and spatial resolutions (10¹–10⁶ atoms), although technical challenges also require balancing between the computational accuracy, cost, and approximations introduced within the model. Here we discuss how the combination of atomistic (aMD) and hybrid quantum/classical molecular dynamics (QM/MM MD) simulations with free energy (FE) sampling methods can be used to probe mechanistic principles of enzymes responsible for biological energy conversion. We present mechanistic explorations of long-range proton-coupled electron transfer (PCET) dynamics in the highly intricate respiratory chain enzyme Complex I, which functions as a redox-driven proton pump in bacterial and mitochondrial respiratory chains by catalyzing a 300 Å fully reversible PCET process. This process is initiated by a hydride (H^–) transfer between NADH and FMN, followed by long-range (>100 Å) electron transfer along a wire of 8 FeS centers leading to a quinone biding site. The reduction of the quinone to quinol initiates dissociation of the latter to a second membrane-bound binding site, and triggers proton pumping across the membrane domain of complex I, in subunits up to 200 Å away from the active site. Our simulations across different size and time scales suggest that transient charge transfer reactions lead to changes in the internal hydration state of key regions, local electric fields, and the conformation of conserved ion pairs, which in turn modulate the dynamics of functional steps along the reaction cycle. Similar functional principles, which operate on much shorter length scales, are also found in some unrelated proteins, suggesting that enzymes may employ conserved principles in the catalysis of biological energy transduction processes.

Does oxidatively damaged DNA drive amyloid-β generation in Alzheimer's disease? A hypothesis

In Alzheimer's disease (AD), amyloid-β (Aβ) generation and upstream β-secretase 1 (BACE1) expression appear to be driven by oxidative stress via c-Jun N-terminal kinase (JNK), p38, and Interferon-Induced, Double-Stranded RNA-Activated Protein Kinase (PKR). In addition, inflammatory molecules, including lipopolysaccharide (LPS), induce genes central to Aβ genesis, such as BACE1, via nuclear factor-κB (NFκB). However, additional triggers of Aβ generation remain poorly understood and might represent novel opportunities for therapeutic intervention. Based on mechanistic studies and elevated ectopic oxidatively damaged DNA (oxoDNA) levels in preclinical AD, mild cognitive impairment, and AD patients, we hypothesize oxoDNA contributes to β-amyloidosis starting from the earliest stages of AD through multiple pathways. OxoDNA induces mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4), thereby sensitizing the brain to oxidative stress-induced JNK activation and BACE1 transcription. It also induces myeloid differentiation primary response 88 (MyD88) and activates protein kinase CK2, thereby increasing NFκB activation and BACE1 induction. OxoDNA increases oxidative stress via nuclear factor erythroid 2-related factor 2 (Nrf2) ectopic localization, likely augmenting JNK-mediated BACE1 induction. OxoDNA likely also promotes β-amyloidosis via absent in melanoma 2 (AIM2) induction. Falsifiable predictions of this hypothesis include that deoxyribonuclease treatment should decrease Aβ and possibly slow cognitive decline in AD patients. While formal testing of this hypothesis remains to be performed, a case report has found deoxyribonuclease I treatment improved a severely demented AD patient's Mini-Mental Status Exam score from 3 to 18 at 2 months. There is preliminary preclinical and clinical evidence suggesting that ectopic oxidatively damaged DNA may act as an inflammatory damage-associated molecular pattern contributing to Aβ generation in AD, and deoxyribonuclease I should be formally evaluated to test whether it can decrease Aβ levels and slow cognitive decline in AD patients.

Martini 3: a general purpose force field for coarse-grained molecular dynamics

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Targeting mitochondrial dysfunction with small molecules in intervertebral disc aging and degeneration

The prevalence of rheumatic and musculoskeletal diseases (RMDs) including osteoarthritis (OA) and low back pain (LBP) in aging societies present significant cost burdens to health and social care systems. Intervertebral disc (IVD) degeneration, which is characterized by disc dehydration, anatomical alterations, and extensive changes in extracellular matrix (ECM) composition, is an important contributor to LBP. IVD cell homeostasis can be disrupted by mitochondrial dysfunction. Mitochondria are the main source of energy supply in IVD cells and a major contributor to the production of reactive oxygen species (ROS). Therefore, mitochondria represent a double-edged sword in IVD cells. Mitochondrial dysfunction results in oxidative stress, cell death, and premature cell senescence, which are all implicated in IVD degeneration. Considering the importance of optimal mitochondrial function for the preservation of IVD cell homeostasis, extensive studies have been done in recent years to evaluate the efficacy of small molecules targeting mitochondrial dysfunction. In this article, we review the pathogenesis of mitochondrial dysfunction, aiming to highlight the role of small molecules and a selected number of biological growth factors that regulate mitochondrial function and maintain IVD cell homeostasis. Furthermore, molecules that target mitochondria and their mechanisms of action and potential for IVD regeneration are identified. Finally, we discuss mitophagy as a key mediator of many cellular events and the small molecules regulating its function.

Coarse-Grained Parameterization of Nucleotide Cofactors and Metabolites: Protonation Constants, Partition Coefficients, and Model Topologies

Nucleotides are structural units relevant not only in nucleic acids but also as substrates or cofactors in key biochemical reactions. The size- and timescales of such nucleotide-protein interactions fall well within the scope of coarse-grained molecular dynamics, which holds promise of important mechanistic insight. However, the lack of specific parameters has prevented accurate coarse-grained simulations of protein interactions with most nucleotide compounds. In this work, we comprehensively develop coarse-grained parameters for key metabolites/cofactors (FAD, FMN, riboflavin, NAD, NADP, ATP, ADP, AMP, and thiamine pyrophosphate) in different oxidation and protonation states as well as for smaller molecules derived from them (among others, nicotinamide, adenosine, adenine, ribose, thiamine, and lumiflavin), summing up a total of 79 different molecules. In line with the Martini parameterization methodology, parameters were tuned to reproduce octanol-water partition coefficients. Given the lack of existing data, we set out to experimentally determine these partition coefficients, developing two methodological approaches, based on ³¹P-NMR and fluorescence spectroscopy, specifically tailored to the strong hydrophilicity of most of the parameterized compounds. To distinguish the partition of each relevant protonation species, we further potentiometrically characterized the protonation constants of key molecules. This work successfully builds a comprehensive and relevant set of computational models that will boost the biochemical application of coarse-grained simulations. It does so based on the measurement of partition and acid-base physicochemical data that, in turn, covers important gaps in nucleotide characterization.

Plastoquinone In and Beyond Photosynthesis

Plastoquinone-9 (PQ-9) is an essential component of photosynthesis that carries electrons in the linear and alternative electron transport chains, and is also a redox sensor that regulates state transitions and gene expression. However, a large fraction of the PQ pool is located outside the thylakoid membranes, in the plastoglobules and the chloroplast envelopes, reflecting a wider range of functions beyond electron transport. This review describes new functions of PQ in photoprotection, as a potent antioxidant, and in chloroplast metabolism as a cofactor in the biosynthesis of chloroplast metabolites. It also focuses on the essential need for tight environmental control of PQ biosynthesis and for active exchange of this compound between the thylakoid membranes and the plastoglobules. Through its multiple functions, PQ connects photosynthesis with metabolism, light acclimation, and stress tolerance.

Synergetic suppression effects upon the combination of UV-C irradiation and berberine on Microcystis aeruginosa and Scenedesmus obliquus in reclaimed water: Effectiveness and mechanisms

The risk of harmful algal blooms (HABs) in the water recharged with reclaimed water is a bottleneck for water reuse. The suppression effects and mechanisms of the combination of UV-C and berberine on Microcystis aeruginosa and Scenedesmus obliquus in reclaimed water were investigated. Mono UV-C irradiation at 75 mJ cm^-2 could suppress the growth of M. aeruginosa for 7 d and that at 90 mJ cm^-2 could suppress the growth of S. obliquus for 5 d. UV-C irradiation combined with 0.2-2 mg L^-1 berberine lengthened the inhibition period of M. aeruginosa to 10- > 22 d and that of S. obliquus to 7- > 22 d and induced more rapid lethal effects on the harmful microalgal cells, in significant synergetic patterns. The combination of UV-C and berberine suppressed total, intracellular and extracellular microcystin-LR (MC-LR) more effectively and decreased the MC-LR quota significantly, which further reduced the risks of microcystin production and release. Furthermore, synergetic mechanisms of the combined treatments were systematically investigated from the aspects of photosynthetic system (photosynthetic activity and pigments), metabolic activity (ATP and membrane potential), oxidation stress (reactive oxygen species (ROS) and glutathione (GSH)), and apoptosis-like cell death (phosphatidylserine (PS) ectropion, caspase-3 activity and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) positive rate). The combination treatment provided a joint attack of UV-C and berberine on photosynthetic transport chain of photosynthetic system II (PS II), and a synergetic pathway to achieve more severe disruptions in energy metabolism as well as aggravated oxidative stress. The accumulated ROS enhanced increases in programmed cell death (PCD) indicators of both microalgal species, which contributed to the enhancement effects on growth suppression. The results showed that the combination treatment achieved lower dose requirements of both UV-C irradiation and berberine for inducing the same inhibition effects on microalgal cells, which was promising to be applied in the HABs control of reclaimed water.

Light harvesting in oxygenic photosynthesis: Structural biology meets spectroscopy

Oxygenic photosynthesis is the main process that drives life on earth. It starts with the harvesting of solar photons that, after transformation into electronic excitations, lead to charge separation in the reaction centers of photosystems I and II (PSI and PSII). These photosystems are large, modular pigment-protein complexes that work in series to fuel the formation of carbohydrates, concomitantly producing molecular oxygen. Recent advances in cryo-electron microscopy have enabled the determination of PSI and PSII structures in complex with light-harvesting components called "supercomplexes" from different organisms at near-atomic resolution. Here, we review the structural and spectroscopic aspects of PSI and PSII from plants and algae that directly relate to their light-harvesting properties, with special attention paid to the pathways and efficiency of excitation energy transfer and the regulatory aspects.

Unique organization of photosystem II supercomplexes and megacomplexes in Norway spruce

Photosystem II (PSII) complexes are organized into large supercomplexes with variable amounts of light-harvesting proteins (Lhcb). A typical PSII supercomplex in plants is formed by four trimers of Lhcb proteins (LHCII trimers), which are bound to the PSII core dimer via monomeric antenna proteins. However, the architecture of PSII supercomplexes in Norway spruce[Picea abies (L.) Karst.] is different, most likely due to a lack of two Lhcb proteins, Lhcb6 and Lhcb3. Interestingly, the spruce PSII supercomplex shares similar structural features with its counterpart in the green alga Chlamydomonas reinhardtii [Kouřil et al. (2016) New Phytol. 210, 808-814]. Here we present a single-particle electron microscopy study of isolated PSII supercomplexes from Norway spruce that revealed binding of a variable amount of LHCII trimers to the PSII core dimer at positions that have never been observed in any other plant species so far. The largest spruce PSII supercomplex, which was found to bind eight LHCII trimers, is even larger than the current largest known PSII supercomplex from C. reinhardtii. We have also shown that the spruce PSII supercomplexes can form various types of PSII megacomplexes, which were also identified in intact grana membranes. Some of these large PSII supercomplexes and megacomplexes were identified also in Pinus sylvestris, another representative of the Pinaceae family. The structural variability and complexity of LHCII organization in Pinaceae seems to be related to the absence of Lhcb6 and Lhcb3 in this family, and may be beneficial for the optimization of light-harvesting under varying environmental conditions.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

Titratable Martini model for constant pH simulations

In this work, we deliver a proof of concept for a fast method that introduces pH effects into classical coarse-grained (CG) molecular dynamics simulations. Our approach is based upon the latest version of the popular Martini CG model to which explicit proton mimicking particles are added. We verify our approach against experimental data involving several different molecules and different environmental conditions. In particular, we compute titration curves, pH dependent free energies of transfer, and lipid bilayer membrane affinities as a function of pH. Using oleic acid as an example compound, we further illustrate that our method can be used to study passive translocation in lipid bilayers via protonation. Finally, our model reproduces qualitatively the expansion of the macromolecule dendrimer poly(propylene imine) as well as the associated pKa shift of its different generations. This example demonstrates that our model is able to pick up collective interactions between titratable sites in large molecules comprising many titratable functional groups.

Charge Transfer and Chemo-Mechanical Coupling in Respiratory Complex I

The mitochondrial respiratory chain, formed by five protein complexes, utilizes energy from catabolic processes to synthesize ATP. Complex I, the first and the largest protein complex of the chain, harvests electrons from NADH to reduce quinone, while pumping protons across the mitochondrial membrane. Detailed knowledge of the working principle of such coupled charge-transfer processes remains, however, fragmentary due to bottlenecks in understanding redox-driven conformational transitions and their interplay with the hydrated proton pathways. Complex I from Thermus thermophilus encases 16 subunits with nine iron-sulfur clusters, reduced by electrons from NADH. Here, employing the latest crystal structure of T. thermophilus complex I, we have used microsecond-scale molecular dynamics simulations to study the chemo-mechanical coupling between redox changes of the iron-sulfur clusters and conformational transitions across complex I. First, we identify the redox switches within complex I, which allosterically couple the dynamics of the quinone binding pocket to the site of NADH reduction. Second, our free-energy calculations reveal that the affinity of the quinone, specifically menaquinone, for the binding-site is higher than that of its reduced, menaquinol form-a design essential for menaquinol release. Remarkably, the barriers to diffusive menaquinone dynamics are lesser than that of the more ubiquitous ubiquinone, and the naphthoquinone headgroup of the former furnishes stronger binding interactions with the pocket, favoring menaquinone for charge transport in T. thermophilus. Our computations are consistent with experimentally validated mutations and hierarchize the key residues into three functional classes, identifying new mutation targets. Third, long-range hydrogen-bond networks connecting the quinone-binding site to the transmembrane subunits are found to be responsible for proton pumping. Put together, the simulations reveal the molecular design principles linking redox reactions to quinone turnover to proton translocation in complex I.

Molecular dynamics simulations in photosynthesis

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

Effect of different light intensity on physiology, antioxidant capacity and photosynthetic characteristics on wheat seedlings under high CO2 concentration in a closed artificial ecosystem

The growth of plants under high carbon dioxide (CO₂) concentrations (≥ 1000 ppm) is explored for the climate change and the bioregenerative life support system (BLSS) environment of long-duration space missions. Wheat (Triticum aestivum L.) is a grass cultivated for cereal grain-a global staple food including astronauts. Light and CO₂ are both indispensable conditions for wheat seedlings. This study provides insights on the physiology, antioxidant capacity and photosynthetic characteristics of wheat seedlings under a range of photosynthetic photon flux densities in a closed system simulating BLSS with high CO₂ concentration. We found that the F_v/F_m, F_v/F₀, chlorophyll content, intrinsic water use efficiencies (WUE_i), membrane stability index (MSI), and malondialdehyde (MDA) of wheat seedlings grown under an intermediate light intensity of 600 μmol m^-2 s^-1 environment were all largest. Interestingly, the high light intensity of 1200 mol m^-2 s^-1 treatment group exhibits the highest net photosynthetic rate but the lowest MDA content. The stomatal conductance and F₀ of high light intensity of 1000 μmol m^-2 s^-1 treatment group were both significantly higher than that of other groups. Our study provides basic knowledge on the wheat growth in different environments, especially in a closed ecosystem with artificial lights.

The study of conformational changes in photosystem II during a charge separation

Photosystem II (PSII) is a multi-subunit pigment-protein complex and is one of several protein assemblies that function cooperatively in photosynthesis in plants and cyanobacteria. As more structural data on PSII become available, new questions arise concerning the nature of the charge separation in PSII reaction center (RC). The crystal structure of PSII RC from cyanobacteria Thermosynechococcus vulcanus was selected for the computational study of conformational changes in photosystem II associated to the charge separation process. The parameterization of cofactors and lipids for classical MD simulation with Amber force field was performed. The parametrized complex of PSII was embedded in the lipid membrane for MD simulation with Amber in Gromacs. The conformational behavior of protein and the cofactors directly involved in the charge separation were studied by MD simulations and QM/MM calculations. This study identified the most likely mechanism of the proton-coupled reduction of plastoquinone Q_B. After the charge separation and the first electron transfer to Q_B, the system undergoes conformational change allowing the first proton transfer to Q_B^- mediated via Ser264. After the second electron transfer to Q_BH, the system again adopts conformation allowing the second proton transfer to Q_BH^-. The reduced Q_BH₂ would then leave the binding pocket.

Plastoquinone homoeostasis by Arabidopsis proton gradient regulation 6 is essential for photosynthetic efficiency

Photosynthesis produces organic carbon via a light-driven electron flow from H2O to CO2 that passes through a pool of plastoquinone molecules. These molecules are either present in the photosynthetic thylakoid membranes, participating in photochemistry (photoactive pool), or stored (non-photoactive pool) in thylakoid-attached lipid droplets, the plastoglobules. The photoactive pool acts also as a signal of photosynthetic activity allowing the adaptation to changes in light condition. Here we show that, in Arabidopsis thaliana, proton gradient regulation 6 (PGR6), a predicted atypical kinase located at plastoglobules, is required for plastoquinone homoeostasis, i.e. to maintain the photoactive plastoquinone pool. In a pgr6 mutant, the photoactive pool is depleted and becomes limiting under high light, affecting short-term acclimation and photosynthetic efficiency. In the long term, pgr6 seedlings fail to adapt to high light and develop a conditional variegated leaf phenotype. Therefore, PGR6 activity, by regulating plastoquinone homoeostasis, is required to cope with high light.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Ceramides bind VDAC2 to trigger mitochondrial apoptosis

Ceramides draw wide attention as tumor suppressor lipids that act directly on mitochondria to trigger apoptotic cell death. However, molecular details of the underlying mechanism are largely unknown. Using a photoactivatable ceramide probe, we here identify the voltage-dependent anion channels VDAC1 and VDAC2 as mitochondrial ceramide binding proteins. Coarse-grain molecular dynamics simulations reveal that both channels harbor a ceramide binding site on one side of the barrel wall. This site includes a membrane-buried glutamate that mediates direct contact with the ceramide head group. Substitution or chemical modification of this residue abolishes photolabeling of both channels with the ceramide probe. Unlike VDAC1 removal, loss of VDAC2 or replacing its membrane-facing glutamate with glutamine renders human colon cancer cells largely resistant to ceramide-induced apoptosis. Collectively, our data support a role of VDAC2 as direct effector of ceramide-mediated cell death, providing a molecular framework for how ceramides exert their anti-neoplastic activity.