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.

The site energies and spectral densities of the major light-harvesting complex LHCII have been determined using QM/MM molecular dynamics simulations. In particular, the spectral densities show a remarkable agreement with experimental findings.

The structure of photosystem II in the presence of subunit S (PsbS) under photoprotection: implications for the energy flow.

The ability to harvest light effectively in a changing environment is necessary to ensure efficient photosynthesis and crop growth. One mechanism, known as qE, protects photosystem II (PSII) and regulates electron transfer through the harmless dissipation of excess absorbed photons as heat. This process involves reversible clustering of the major light-harvesting complexes of PSII (LHCII) in the thylakoid membrane and relies upon the ΔpH gradient and the allosteric modulator protein PsbS. To date, the exact role of PsbS in the qE mechanism has remained elusive. Here, we show that PsbS induces hydrophobic mismatch in the thylakoid membrane through dynamic rearrangement of lipids around LHCII leading to observed membrane thinning. We found that upon illumination, the thylakoid membrane reversibly shrinks from around 4.3 to 3.2 nm, without PsbS, this response is eliminated. Furthermore, we show that the lipid digalactosyldiacylglycerol (DGDG) is repelled from the LHCII-PsbS complex due to an increase in both the pK_a of lumenal residues and in the dipole moment of LHCII, which allows for further conformational change and clustering in the membrane. Our results suggest a mechanistic role for PsbS as a facilitator of a hydrophobic mismatch-mediated phase transition between LHCII-PsbS and its environment. This could act as the driving force to sort LHCII into photoprotective nanodomains in the thylakoid membrane. This work shows an example of the key role of the hydrophobic mismatch process in regulating membrane protein function in plants.

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.