Acclimation mechanism of microalgal photosynthetic apparatus under low atmospheric pressures - new astrobiological perspectives in a Mars-like atmosphere

This study reveals a new acclimation mechanism of the eukaryotic unicellular green alga Chlorella vulgaris in terms of the effect of varying atmospheric pressures on the structure and function of its photosynthetic apparatus using fluorescence induction measurements (JIP-test). The results indicate that low (400mbar) and extreme low (2 atmosphere (simulating the Mars atmosphere), reveals that the impact of extremely low atmospheric pressure on PQ mobility within the photosynthetic membrane, coupled with the low density of an almost 100% CO2 Mars-like atmosphere, results to a similar photosynthetic efficiency to that on Earth. These findings pave the way for the identification of novel functional acclimation mechanisms of microalgae to extreme environments that are vastly distinct from those found on Earth.

Unravelling the fluorescence kinetics of light-harvesting proteins with simulated measurements

The plant light-harvesting pigment-protein complex LHCII is the major antenna sub-unit of PSII and is generally (though not universally) accepted to play a role in photoprotective energy dissipation under high light conditions, a process known Non-Photochemical Quenching (NPQ). The underlying mechanisms of energy trapping and dissipation within LHCII are still debated. Various models have been proposed for the underlying molecular detail of NPQ, but they are often based on different interpretations of very similar transient absorption measurements of isolated complexes. Here we present a simulated measurement of the fluorescence decay kinetics of quenched LHCII aggregates to determine whether this relatively simple measurement can discriminate between different potential NPQ mechanisms. We simulate not just the underlying physics (excitation, energy migration, quenching and singlet-singlet annihilation) but also the signal detection and typical experimental data analysis. Comparing this to a selection of published fluorescence decay kinetics we find that: (1) Different proposed quenching mechanisms produce noticeably different fluorescence kinetics even at low (annihilation free) excitation density, though the degree of difference is dependent on pulse width. (2) Measured decay kinetics are consistent with most LHCII trimers becoming relatively slow excitation quenchers. A small sub-population of very fast quenchers produces kinetics which do not resemble any observed measurement. (3) It is necessary to consider at least two distinct quenching mechanisms in order to accurately reproduce experimental kinetics, supporting the idea that NPQ is not a simple binary switch.

Cryo-EM structures of LHCII in photo-active and photo-protecting states reveal allosteric regulation of light harvesting and excess energy dissipation

The major light-harvesting complex of photosystem II (LHCII) has a dual regulatory function in a process called non-photochemical quenching to avoid the formation of reactive oxygen. LHCII undergoes reversible conformation transitions to switch between a light-harvesting state for excited-state energy transfer and an energy-quenching state for dissipating excess energy under full sunshine. Here we report cryo-electron microscopy structures of LHCII in membrane nanodiscs, which mimic in vivo LHCII, and in detergent solution at pH 7.8 and 5.4, respectively. We found that, under low pH conditions, the salt bridges at the lumenal side of LHCII are broken, accompanied by the formation of two local α-helices on the lumen side. The formation of α-helices in turn triggers allosterically global protein conformational change, resulting in a smaller crossing angle between transmembrane helices. The fluorescence decay rates corresponding to different conformational states follow the Dexter energy transfer mechanism with a characteristic transition distance of 5.6 Å between Lut1 and Chl612. The experimental observations are consistent with the computed electronic coupling strengths using multistate density function theory.

A perspective on the major light-harvesting complex dynamics under the effect of pH, salts, and the photoprotective PsbS protein

The photosynthetic apparatus is a highly modular assembly of large pigment-binding proteins. Complexes called antennae can capture the sunlight and direct it from the periphery of two Photosystems (I, II) to the core reaction centers, where it is converted into chemical energy. The apparatus must cope with the natural light fluctuations that can become detrimental to the viability of the photosynthetic organism. Here we present an atomic scale view of the photoprotective mechanism that is activated on this line of defense by several photosynthetic organisms to avoid overexcitation upon excess illumination. We provide a complete macroscopic to microscopic picture with specific details on the conformations of the major antenna of Photosystem II that could be associated with the switch from the light-harvesting to the photoprotective state. This is achieved by combining insight from both experiments and all-atom simulations from our group and the literature in a perspective article.

The origin of pigment-binding differences in CP29 and LHCII: the role of protein structure and dynamics

The first step of photosynthesis in plants is performed by the light-harvesting complexes (LHC), a large family of pigment-binding proteins embedded in the photosynthetic membranes. These complexes are conserved across species, suggesting that each has a distinct role. However, they display a high degree of sequence homology and their static structures are almost identical. What are then the structural features that determine their different properties? In this work, we compared the two best-characterized LHCs of plants: LHCII and CP29. Using molecular dynamics simulations, we could rationalize the difference between them in terms of pigment-binding properties. The data also show that while the loops between the helices are very flexible, the structure of the transmembrane regions remains very similar in the crystal and the membranes. However, the small structural differences significantly affect the excitonic coupling between some pigment pairs. Finally, we analyzed in detail the structure of the long N-terminus of CP29, showing that it is structurally stable and it remains on top of the membrane even in the absence of other proteins. Although the structural changes upon phosphorylation are minor, they can explain the differences in the absorption properties of the pigments observed experimentally.

Recognition motifs for importin 4 [(L)PPRS(G/P)P] and importin 5 [KP(K/Y)LV] binding, identified by bio-informatic simulation and experimental in vitro validation

Nuclear translocation of large proteins is mediated through karyopherins, carrier proteins recognizing specific motifs of cargo proteins, known as nuclear localization signals (NLS). However, only few NLS signals have been reported until now. In the present work, NLS signals for Importins 4 and 5 were identified through an unsupervised in silico approach, followed by experimental in vitro validation. The sequences LPPRS(G/P)P and KP(K/Y)LV were identified and are proposed as recognition motifs for Importins 4 and 5 binding, respectively. They are involved in the trafficking of important proteins into the nucleus. These sequences were validated in the breast cancer cell line T47D, which expresses both Importins 4 and 5. Elucidating the complex relationships of the nuclear transporters and their cargo proteins is very important in better understanding the mechanism of nuclear transport of proteins and laying the foundation for the development of novel therapeutics, targeting specific importins.

Deciphering the QR Code of the CRISPR-Cas9 System: Synergy between Gln768 (Q) and Arg976 (R)

Markov state models (MSMs) and machine learning (ML) algorithms can extrapolate the long-time-scale behavior of large biomolecules from molecular dynamics (MD) trajectories. In this study, an MD-MSM-ML scheme has been applied to probe the large endonuclease (Cas9) in the bacterial adaptive immunity CRISPR-Cas9 system. CRISPR has become a programmable and state-of-the-art powerful genome editing tool that has already revolutionized life sciences. CRISPR-Cas9 is programmed to process specific DNA sequences in the genome. However, human/biomedical applications are compromised by off-target DNA damage. Characterization of Cas9 at the structural and biophysical levels is a prerequisite for the development of efficient and high-fidelity Cas9 variants. The Cas9 wild type and two variants (R63A-R66A-R70A, R69A-R71A-R74A-R78A) are studied herein. The configurational space of Cas9 is provided with a focus on the conformations of the side chains of two residues (Gln768 and Arg976). A model for the synergy between those two residues is proposed. The results are discussed within the context of experimental literature. The results and methodology can be exploited for the study of large biomolecules in general and for the engineering of more efficient and safer Cas9 variants for applications.

Conformation of Light-Harvesting Complex II Trimer Depends upon Its Binding Site

The light-harvesting complex II (LHCII) trimer in plants functions as a major antenna complex and a quencher to protect it from photooxidative damage. Theoretical studies on the structure of an LHCII trimer have demonstrated that excitation energy transfer between chlorophylls (Chls) in LHCII can be modulated by its exquisite conformational fluctuation. However, conformational changes depending on its binding location have not yet been investigated, even though reorganization of protein complexes occurs by physiological regulations. In this study, we investigated conformational differences in LHCII by comparing published structures of an identical LHCII trimer in the three different photosystem supercomplexes from the green alga Chlamydomonas reinhardtii. Our results revealed distinct differences in Chl configurations as well as polypeptide conformations of the LHCII trimers depending on its binding location. We propose that these configurational differences readily modulate the function of LHCII and possibly lead to a change in excitation-energy flow over the photosynthetic supercomplex.

From Traditional Ethnopharmacology to Modern Natural Drug Discovery: A Methodology Discussion and Specific Examples

Ethnopharmacology, through the description of the beneficial effects of plants, has provided an early framework for the therapeutic use of natural compounds. Natural products, either in their native form or after crude extraction of their active ingredients, have long been used by different populations and explored as invaluable sources for drug design. The transition from traditional ethnopharmacology to drug discovery has followed a straightforward path, assisted by the evolution of isolation and characterization methods, the increase in computational power, and the development of specific chemoinformatic methods. The deriving extensive exploitation of the natural product chemical space has led to the discovery of novel compounds with pharmaceutical properties, although this was not followed by an analogous increase in novel drugs. In this work, we discuss the evolution of ideas and methods, from traditional ethnopharmacology to in silico drug discovery, applied to natural products. We point out that, in the past, the starting point was the plant itself, identified by sustained ethnopharmacological research, with the active compound deriving after extensive analysis and testing. In contrast, in recent years, the active substance has been pinpointed by computational methods (in silico docking and molecular dynamics, network pharmacology), followed by the identification of the plant(s) containing the active ingredient, identified by existing or putative ethnopharmacological information. We further stress the potential pitfalls of recent in silico methods and discuss the absolute need for in vitro and in vivo validation as an absolute requirement. Finally, we present our contribution to natural products’ drug discovery by discussing specific examples, applying the whole continuum of this rapidly evolving field. In detail, we report the isolation of novel antiviral compounds, based on natural products active against influenza and SARS-CoV-2 and novel substances active on a specific GPCR, OXER1.

A fast method for electronic couplings in embedded multichromophoric systems

Electronic couplings are key to understanding exciton delocalization and transport in natural and artificial light harvesting processes. We develop a method to compute couplings in multichromophoric aggregates embedded in complex environments without running expensive quantum chemical calculations. We use a transition charge approximation to represent the quantum mechanical transition densities of the chromophores and an atomistic and polarizable classical model to describe the environment atoms. We extend our framework to estimate transition charges directly from the chromophore geometry, i.e., bypassing completely the quantum mechanical calculations using a regression approach. The method allows to rapidly compute accurate couplings for a large number of geometries along molecular dynamics trajectories.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

Protein dynamics and lipid affinity of monomeric, zeaxanthin-binding LHCII in thylakoid membranes

The xanthophyll cycle in the antenna of photosynthetic organisms under light stress is one of the most well-known processes in photosynthesis, but its role is not well understood. In the xanthophyll cycle, violaxanthin (Vio) is reversibly transformed to zeaxanthin (Zea) that occupies Vio binding sites of light-harvesting antenna proteins. Higher monomer/trimer ratios of the most abundant light-harvesting protein, the light-harvesting complex II (LHCII), usually occur in Zea accumulating membranes and have been observed in plants after prolonged illumination and during high-light acclimation. We present a combined NMR and coarse-grained simulation study on monomeric LHCII from the npq2 mutant that constitutively binds Zea in the Vio binding pocket. LHCII was isolated from ¹³C-enriched npq2 Chlamydomonas reinhardtii (Cr) cells and reconstituted in thylakoid lipid membranes. NMR results reveal selective changes in the fold and dynamics of npq2 LHCII compared with the trimeric, wild-type and show that npq2 LHCII contains multiple mono- or digalactosyl diacylglycerol lipids (MGDG and DGDG) that are strongly protein bound. Coarse-grained simulations on npq2 LHCII embedded in a thylakoid lipid membrane agree with these observations. The simulations show that LHCII monomers have more extensive lipid contacts than LHCII trimers and that protein-lipid contacts are influenced by Zea. We propose that both monomerization and Zea binding could have a functional role in modulating membrane fluidity and influence the aggregation and conformational dynamics of LHCII with a likely impact on photoprotection ability.

A different perspective for nonphotochemical quenching in plant antenna complexes

Light-harvesting complexes of plants exert a dual function of light-harvesting (LH) and photoprotection through processes collectively called nonphotochemical quenching (NPQ). While LH processes are relatively well characterized, those involved in NPQ are less understood. Here, we characterize the quenching mechanisms of CP29, a minor LHC of plants, through the integration of two complementary enhanced-sampling techniques, dimensionality reduction schemes, electronic calculations and the analysis of cryo-EM data in the light of the predicted conformational ensemble. Our study reveals that the switch between LH and quenching state is more complex than previously thought. Several conformations of the lumenal side of the protein occur and differently affect the pigments' relative geometries and interactions. Moreover, we show that a quenching mechanism localized on a single chlorophyll-carotenoid pair is not sufficient but many chlorophylls are simultaneously involved. In such a diffuse mechanism, short-range interactions between each carotenoid and different chlorophylls combined with a protein-mediated tuning of the carotenoid excitation energies have to be considered in addition to the commonly suggested Coulomb interactions.

Natural Polyphenols Inhibit the Dimerization of the SARS-CoV-2 Main Protease: The Case of Fortunellin and Its Structural Analogs

3CL-Pro is the SARS-CoV-2 main protease (MPro). It acts as a homodimer to cleave the large polyprotein 1ab transcript into proteins that are necessary for viral growth and replication. 3CL-Pro has been one of the most studied SARS-CoV-2 proteins and a main target of therapeutics. A number of drug candidates have been reported, including natural products. Here, we employ elaborate computational methods to explore the dimerization of the 3CL-Pro protein, and we formulate a computational context to identify potential inhibitors of this process. We report that fortunellin (acacetin 7-O-neohesperidoside), a natural flavonoid O-glycoside, and its structural analogs are potent inhibitors of 3CL-Pro dimerization, inhibiting viral plaque formation in vitro. We thus propose a novel basis for the search of pharmaceuticals as well as dietary supplements in the fight against SARS-CoV-2 and COVID-19.

Robust Strategy for Photoprotection in the Light-Harvesting Antenna of Diatoms: A Molecular Dynamics Study

Diatoms generate a large portion of the oxygen produced on earth due to their exceptional light-harvesting properties involving fucoxanthin and chlorophyll-binding proteins (FCP). At the same time, an efficient adaptation of these complexes to fluctuating light conditions is necessary to protect the diatoms against photodamage. So far, structural and dynamic data for the interaction between FCP and the photoprotective LHCX family of proteins in diatoms are lacking. In this computational study, we provide a structural basis for a remarkable pH-dependent adaptation at the molecular level. Upon binding of the LHCX1 protein to the FCP complex together with a change in pH, conformational changes within the FCP protein result in a variation of the electronic coupling in a specific chlorophyll-fucoxanthin pair, leading to a change in the exciton transfer rate by almost an order of magnitude. A common strategy for photoprotection between diatoms and higher plants is identified and discussed.

Time-dependent atomistic simulations of the CP29 light-harvesting complex

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

Hydrogen gas as a central on-off functional switch of reversible metabolic arrest - New perspectives for biotechnological applications

Metabolism is the sum of all chemical reactions that sustain life. There is an ongoing effort to control metabolic rate, which correlates with the maximum lifespan potential and constitutes one of the oldest scientific questions. Herein, we report on the complete reversible arrest of cellular metabolism and cell growth in a series of organisms, from microalgae to yeast upon exposure to a 100 % hydrogen atmosphere. We also report a tolerance of the microalgae under these conditions against extreme stress conditions, like high salt concentrations. The addition of oxygen or air almost completely restores the metabolic rate and cell growth. Molecular dynamics simulations are employed to decipher this phenomenon at atomic scale. Various proteins, including photosynthetic and respiratory complexes (LHCII, cytochrome c5) are probed in the interaction with hydrogen. Exposure to hydrogen, as opposed to oxygen, decreases the fluctuations of protein residues indicating thermostability. According to the above mechanism, an absolute hydrogen atmosphere can preserve biological products (e.g. fruits) for a long time without consuming any energy. By combining biological, chemical and computational methods, in this research we provide the basis for future innovative studies and advances in the field of biotechnology.

Uncovering the interactions driving carotenoid binding in light-harvesting complexes

Carotenoids are essential constituents of plant light-harvesting complexes (LHCs), being involved in protein stability, light harvesting, and photoprotection. Unlike chlorophylls, whose binding to LHCs is known to require coordination of the central magnesium, carotenoid binding relies on weaker intermolecular interactions (such as hydrogen bonds and van der Waals forces), whose character is far more elusive. Here we addressed the key interactions responsible for carotenoid binding to LHCs by combining molecular dynamics simulations and polarizable quantum mechanics/molecular mechanics calculations on the major LHC, LHCII. We found that carotenoid binding is mainly stabilized by van der Waals interactions with the surrounding chlorophyll macrocycles rather than by hydrogen bonds to the protein, the latter being more labile than predicted from structural data. Furthermore, the interaction network in the binding pockets is relatively insensitive to the chemical structure of the embedded carotenoid. Our results are consistent with a number of experimental data and challenge the role played by specific interactions in the assembly of pigment-protein complexes.

Trivial Excitation Energy Transfer to Carotenoids Is an Unlikely Mechanism for Non-photochemical Quenching in LHCII

Higher plants defend themselves from bursts of intense light via the mechanism of Non-Photochemical Quenching (NPQ). It involves the Photosystem II (PSII) antenna protein (LHCII) adopting a conformation that favors excitation quenching. In recent years several structural models have suggested that quenching proceeds via energy transfer to the optically forbidden and short-lived S ₁ states of a carotenoid. It was proposed that this pathway was controlled by subtle changes in the relative orientation of a small number of pigments. However, quantum chemical calculations of S ₁ properties are not trivial and therefore its energy, oscillator strength and lifetime are treated as rather loose parameters. Moreover, the models were based either on a single LHCII crystal structure or Molecular Dynamics (MD) trajectories about a single minimum. Here we try and address these limitations by parameterizing the vibronic structure and relaxation dynamics of lutein in terms of observable quantities, namely its linear absorption (LA), transient absorption (TA) and two-photon excitation (TPE) spectra. We also analyze a number of minima taken from an exhaustive meta-dynamical search of the LHCII free energy surface. We show that trivial, Coulomb-mediated energy transfer to S ₁ is an unlikely quenching mechanism, with pigment movements insufficiently pronounced to switch the system between quenched and unquenched states. Modulation of S ₁ energy level as a quenching switch is similarly unlikely. Moreover, the quenching predicted by previous models is possibly an artifact of quantum chemical over-estimation of S ₁ oscillator strength and the real mechanism likely involves short-range interaction and/or non-trivial inter-molecular states.