Search

The functional units of natural photosynthetic systems control the process of converting sunlight into chemical energy. In this article, we explore a series of chemically and structurally modified bacteriochlorophyll and chlorophyll pigments through computational chemistry to evaluate their electronic spectroscopy properties. More specifically, we use multiconfigurational and time-dependent density functional theory methods, along with molecular dynamics simulations, to compute the models' energetics both in an implicit and explicit solvent environment. Structural modifications aimed at reducing the planarity of the macrocycle through alkyl-bridge anchoring reveal the significant role of the curvature in fine-tuning spectral properties, which mimics protein scaffold effects on naturally occurring pigments. Furthermore, chemical substitutions with a carbonyl group show potential for expanding absorption spectra toward the blue region, while incorporating an additional double bond decreases absorption efficiency. These insights lay the groundwork to design novel synthetic pigments, with potential applications in artificial light-harvesting systems and more efficient photovoltaic devices.

Molecular excitons play a major role within dye aggregates and hold significant potential for (opto)electronic and photovoltaic applications. Numerous studies have documented alterations in the spectral properties of dye homoaggregates, but only limited work has been reported for heteroaggregates. In this article, dimeric dye stacks were constructed from azobenzene-like dyes with identical or distinct structures, and their excitonic features were computationally investigated. Our results show that strong exciton coupling is not limited to identical chromophores, as often assumed, based on a recently made available Frenkel Exciton Hamiltonian and multiconfigurational plus second-order perturbation theory energetics methodology. Heteroaggregate stacks were found to exhibit different absorption features from the corresponding interacting monomers, indicating considerable coupling interactions between units. We analyzed how such coupling may vary according to various aspects, such as the relative positions of the interacting monomers or the differences in their energetics. Such qualitative and semiquantitative analyses allow the evaluation of the excitonic behavior of these dye aggregates to encourage further efforts toward a deeper understanding of the excitonic properties of tailored dye heteroaggregate systems.

Oxygenic photosynthesis begins in the reaction center (RC) of the protein complex photosystem II (PSII). PSII has an intriguing, nearly symmetrical arrangement of cofactors within its RC. Despite this symmetry, evolution has favored only one of the two branches of PSII for efficient electron transfer. Current spectroscopic experiments explore the electronic dynamics during the picoseconds after energy has entered the RC and until the electron transfers to the pheophytin of the first branch. We present state-of-the-art multiconfigurational multireference calculations of the excitation energies or site energies of the four chlorophyll pigments of the RC without protein environment considerations. We see a significant variation that breaks the apparent symmetry of the RC. The inner chlorophyll of the productive RC branch possessed the lowest excitation energy of the four central chlorophylls. Our computational method used here is expensive; thus, geometry optimization of the crystal structure is currently not possible. In future work, charge and energy dynamics within the RC will be included as well as a dynamic description of the protein environment and its coupling to the RC. Other state-of-the-art studies of the RC, at lower levels of electronic structure, include a static treatment of the protein environment. These almost unanimously report that the outer chlorophyll of the active branch had the lowest excitation energy. Future work is needed to reconcile this discrepancy., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter., (© 2022. The Author(s).)

Light-harvesting systems 2 and 3 (LH2 and LH3) act as antennas for the initial light capture by photosynthetic purple bacteria, thus initiating the conversion of solar energy into chemical energy. The main absorbers are carotenoids and bacteriochlorophylls (BChls), which harvest different parts of the solar spectrum. The first two optical transitions in BChl produce the Q y and Q x absorption bands. The large size of BChl molecules has prevented accurate computational determination of the electronic structures for the relevant states, until we recently succeeded in obtaining the excitation energies and transition dipole moments of the first (Q y ) transition for all BChls in LH2 and LH3 using multi-state multiconfigurational second-order perturbation theory calculations. In this work, we go one step further, compute the corresponding values for the Q x transition, in line with previous work [ J. Am. Chem. Soc . 2017 , 139 , 7558 - 7567 ], and compare and assess our data against excitation energies obtained through time-dependent density functional theory methods. Interestingly, we find that the two transitions respond differently to BChls' geometrical factors, such as the macrocycle ring curvature and the dihedral torsion of the acetyl moiety. These findings will aid the unraveling of structure-function relationships for absorption and energy transfer processes in purple bacteria, and once again this demonstrates the viability of multireference quantum chemical methods as computational tools for the photophysics of biomolecules.

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

Natural light-harvesting is performed by pigment-protein complexes, which collect and funnel the solar energy at the start of photosynthesis. The identity and arrangement of pigments largely define the absorption spectrum of the antenna complex, which is further regulated by a palette of structural factors. Small alterations are induced by pigment-protein interactions. In light-harvesting systems 2 and 3 from Rhodoblastus acidophilus , the pigments are arranged identically, yet the former has an absorption peak at 850 nm that is blue-shifted to 820 nm in the latter. While the shift has previously been attributed to the removal of hydrogen bonds, which brings changes in the acetyl moiety of the bacteriochlorophyll, recent work has shown that other mechanisms are also present. Using computational and modeling tools on the corresponding crystal structures, we reach a different conclusion: The most critical factor for the shift is the curvature of the macrocycle ring. The bending of the planar part of the pigment is identified as the second-most important design principle for the function of pigment-protein complexes-a finding that can inspire the design of novel artificial systems., Competing Interests: The authors declare no conflict of interest.

Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.

Over the last decade, we have investigated and exploited the photophysical properties of triangulenium dyes. Azadioxatriangulenium (ADOTA) and diazaoxatriangulenium (DAOTA), in particular, have features that make them useful in various fluorescence-based technologies (e.g., bioimaging). Through our work with ADOTA and DAOTA, we became aware that the reported fluorescence quantum yields (ϕ fl ) for these dyes are lower than their actual values. We thus set out to further investigate the fundamental structure-property relationships in these unique conjugated cationic systems. The nonradiative processes in the systems were explored using transient absorption spectroscopy and time-resolved emission spectroscopy in combination with computational chemistry. The influence of molecular oxygen on the fluorescence properties was explored, and the singlet oxygen sensitization efficiencies of ADOTA and DAOTA were determined. We conclude that, for these dyes, the amount of nonradiative deactivation of the first excited singlet state (S 1 ) of the azaoxa-triangulenium fluorophores is low, that the rate of such deactivation is slower than what is observed in common cationic dyes, that there are no observable radiative transitions occurring from the first excited triplet state (T 1 ) of these dyes, and that the efficiency of sensitized singlet oxygen production is low (ϕ Δ ≤ 10%). These photophysical results provide a solid base upon which technological applications of these fluorescent dyes can be built., Competing Interests: The authors declare no competing financial interest.

Benign prostatic hyperplasia is a common disease in men aged over 50 years old, with an incidence increasing to more than 80% over the age of 70, that is increasingly going to attract pharmaceutical interest. Within conventional therapies, such as α -adrenoreceptor antagonists and 5 α -reductase inhibitor, there is a large requirement for treatments with less adverse events on, e.g., blood pressure and sexual function: phytotherapy may be the right way to fill this need. Serenoa repens standardized extract has been widely studied and its ability to reduce lower urinary tract symptoms related to benign prostatic hyperplasia is comprehensively described in literature. An innovative investigation on the mechanism of inhibition of 5 α -reductase by Serenoa repens extract active principles is proposed in this work through computational methods, performing molecular docking simulations on the crystal structure of human liver 5 β -reductase. The results confirm that both sterols and fatty acids can play a role in the inhibition of the enzyme, thus, suggesting a competitive mechanism of inhibition. This work proposes a further confirmation for the rational use of herbal products in the management of benign prostatic hyperplasia, and suggests computational methods as an innovative, low cost, and non-invasive process for the study of phytocomplex activity toward proteic targets., Competing Interests: The authors declare there are no competing interests.

In this study a number of chromophores based on boron subphthalocyanines are investigated for use in the future design of organic photovoltaic devices based on molecular triad systems. The computational study is performed at the TD-DFT CAM-B3LYP/6-311G(d) level of theory. The absorption spectra of these chromophores are simulated using TD-DFT and compared to experimental results. All investigated chromophores absorb light in the visible range and thus are suitable for absorption of sunlight in solar cell applications. On the basis of energy-level alignments, suitable combinations of moieties for a molecular triad system are proposed. The molecular triads will be used in future work as the functional part of organic photovoltaic devices, where the chromophore will be used both to absorb the incoming solar radiation and to increase the distance between the separated charges on donor and acceptor units to increase the lifetime of the charge-separated state.

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636-1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.

Our previously presented method for high throughput computational screening of mutant activity (Hediger et al., 2012) is benchmarked against experimentally measured amidase activity for 22 mutants of Candida antarctica lipase B (CalB). Using an appropriate cutoff criterion for the computed barriers, the qualitative activity of 15 out of 22 mutants is correctly predicted. The method identifies four of the six most active mutants with ≥3-fold wild type activity and seven out of the eight least active mutants with ≤0.5-fold wild type activity. The method is further used to screen all sterically possible (386) double-, triple- and quadruple-mutants constructed from the most active single mutants. Based on the benchmark test at least 20 new promising mutants are identified.

We present a semi-empirical (PM6-based) computational method for systematically estimating the effect of all possible single mutants, within a certain radius of the active site, on the barrier height of an enzymatic reaction. The intent of this method is not a quantitative prediction of the barrier heights, but rather to identify promising mutants for further computational or experimental study. The method is applied to identify promising single and double mutants of Bacillus circulans xylanase (BCX) with increased hydrolytic activity for the artificial substrate ortho-nitrophenyl β-xylobioside (ONPX2). The estimated reaction barrier for wild-type (WT) BCX is 18.5 kcal/mol, which is in good agreement with the experimental activation free energy value of 17.0 kcal/mol extracted from the observed k cat using transition state theory (Joshi et al., 2001). The PM6 reaction profiles for eight single point mutations are recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase in barrier height for all eight mutants while FMO predicts an increase for six of the eight mutants. Both methods predict that the largest change in barrier occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase, respectively. We thus conclude that PM6 is sufficiently accurate to identify promising mutants for further study. We prepared a set of all theoretically possible (342) single mutants in which every amino acid of the active site (except for the catalytically active residues E78 and E172) was mutated to every other amino acid. Based on results from the single mutants we construct a set of 111 double mutants consisting of all possible pairs of single mutants with the lowest barrier for a particular position and compute their reaction profile. None of the mutants have, to our knowledge, been prepared experimentally and therefore present experimentally testable predictions.

We present a fast computational method to efficiently screen enzyme activity. In the presented method, the effect of mutations on the barrier height of an enzyme-catalysed reaction can be computed within 24 hours on roughly 10 processors. The methodology is based on the PM6 and MOZYME methods as implemented in MOPAC2009, and is tested on the first step of the amide hydrolysis reaction catalyzed by the Candida Antarctica lipase B (CalB) enzyme. The barrier heights are estimated using adiabatic mapping and shown to give barrier heights to within 3 kcal/mol of B3LYP/6-31G(d)//RHF/3-21G results for a small model system. Relatively strict convergence criteria (0.5 kcal/(molÅ)), long NDDO cutoff distances within the MOZYME method (15 Å) and single point evaluations using conventional PM6 are needed for reliable results. The generation of mutant structures and subsequent setup of the semiempirical calculations are automated so that the effect on barrier heights can be estimated for hundreds of mutants in a matter of weeks using high performance computing.

We present a web interface which allows us to conveniently set up calculations based on the BioFET-SIM model. With the interface, the signal of a BioFET sensor can be calculated depending on its parameters, as well as the signal dependence on pH. As an illustration, two case studies are presented. In the first case, a generic peptide with opposite charges on both ends is inverted in orientation on a semiconducting nanowire surface leading to a corresponding change in sign of the computed sensitivity of the device. In the second case, the binding of an antibody/antigen complex on the nanowire surface is studied in terms of orientation and analyte/nanowire surface distance. We demonstrate how the BioFET-SIM web interface can aid in the understanding of experimental data and postulate alternative ways of antibody/antigen orientation on the nanowire surface.