The effect of protein conformational flexibility on the electronic properties of a chromophore

Development of artificial photosystems based on designed proteins for mechanistic insights into photosynthesis

This review aims to provide an overview of the progress in protein-based artificial photosystem design and their potential to uncover the underlying principles governing light-harvesting in photosynthesis. While significant advances have been made in this area, a gap persists in reviewing these advances. This review provides a perspective of the field, pinpointing knowledge gaps and unresolved challenges that warrant further inquiry. In particular, it delves into the key considerations when designing photosystems based on the chromophore and protein scaffold characteristics, presents the established strategies for artificial photosystems engineering with their advantages and disadvantages, and underscores the recent breakthroughs in understanding the molecular mechanisms governing light-harvesting, charge separation, and the role of the protein motions in the chromophore's excited state relaxation. By disseminating this knowledge, this article provides a foundational resource for defining the field of bio-hybrid photosystems and aims to inspire the continued exploration of artificial photosystems using protein design.

An exploration of the binding prediction of anatoxin-a and atropine to acetylcholinesterase enzyme using multi-level computer simulations

Acetylcholinesterase (AChE) is crucial for the breakdown of acetylcholine to acetate and choline, while the inhibition of AChE by anatoxin-a (ATX-a) results in severe health complications. This study explores the structural characteristics of ATX-a and its interactions with AChE, comparing to the reference molecule atropine for binding mechanisms. Molecular docking simulations reveal strong binding affinity of both ATX-a and atropine to AChE, interacting effectively with specific amino acids in the binding site as potential inhibitors. Quantitative assessment using the MM-PBSA method demonstrates a significantly negative binding free energy of -81.659 kJ mol-1for ATX-a, indicating robust binding, while atropine exhibits a stronger binding affinity with a free energy of -127.565 kJ mol-1. Umbrella sampling calculates the ΔGbindvalues to evaluate binding free energies, showing a favorable ΔGbindof -36.432 kJ mol-1for ATX-a and a slightly lower value of -30.12 kJ mol-1for atropine. This study reveals the dual functionality of ATX-a, acting as both a nicotinic acetylcholine receptor agonist and an AChE inhibitor. Remarkably, stable complexes form between ATX-a and atropine with AChE at its active site, exhibiting remarkable binding free energies. These findings provide valuable insights into the potential use of ATX-a and atropine as promising candidates for modulating AChE activity.

An exploration of binding of Hesperidin, Rutin, and Thymoquinone to acetylcholinesterase enzyme using multi-level computational approaches

Alzheimer's disease, an intricate neurological disorder, is impacting an ever-increasing number of individuals globally, particularly among the aging population. For several decades phytochemicals were used as Ayurveda to treat both communicable and non-communicable diseases. Acetylcholinesterase (AChE) is a widely chosen therapeutic target for the development of early prevention and effective management of neurodegenerative diseases. The primary objective of the present study was to investigate the binding potential between Rutin Thymoquinone, Hesperidin and the FDA-approved drug Donepezil with AChE. Additionally, a comparative analysis was conducted. These phytochemicals were docked with the binding site of the AChE experimental complex. The molecular dockings demonstrated that the Hesperidinh showed a better binding affinity of -22.0631 kcal/mol. The ADME/T investigations revealed that the selected phytochemicals are non-toxic and drug-like candidates. Molecular dynamics simulations were implemented to determine the conformational changes of Rutin, hesperidin, Thymoquinone, and Donepezil complexed with AChE. Hesperidin and Donepezil were more stable than Rutin, Thymoquinone complexed with AChE. Next, essential dynamics and defining the secondary structure of protein were to determine the conformational changes in AChE complexed with selected phytochemicals during simulations. Overall, the MD Simulations demonstrated that all complexes in this study achieved stability until 100 ns of the simulation period was performed thrice. The structural analysis of AChE was done using multiple search engines to explore the molecular functions, biological processes, and pathways in which AChE proteins are involved and to identify potential drug targets for various diseases. This present study concludes that Hesperidin was found to be a more potent AChE inhibitors than Rutin, and further experiments are required to determine the effectivity of Hesperidin against neurodegenerative diseases.Communicated by Ramaswamy H. Sarma.

Molecular modeling and simulations of some antiviral drugs, benzylisoquinoline alkaloid, and coumarin molecules to investigate the effects on Mpro main viral protease inhibition

Background

SARS-CoV-2 is a deadly viral disease and uncounted deaths occurs since its first appearance in the year 2019. The antiviral drugs, benzylisoquinoline alkaloids, and coumarin molecules were searched using different online engines for drug repurposing with SARS-CoV-2 and to investigate the effects on main viral protease (Mpro) upon their bindings.

Methods

A database composed of antiviral drugs, benzylisoquinoline alkaloids, and Coumarin molecules was screened through a molecular docking strategy to uncover the interactions of collected molecules with SARS-CoV-2 Mpro. Further, molecular dynamics simulations (MDS) were implemented for 100 ns to calculate the stability of the best complexed molecular scaffold with Mpro. The conformations of the simulated complexes were investigated by using principal component analysis (PCA) and Gibbs energy landscape (FEL) and DSSP together. Next, free binding energy (ΔG_bind) was calculated using the mmpbsa method.

Results

Molecular docking simulations demonstrate 17 molecules exhibited better binding affinity out of 99 molecules present in the database with the viral protease Mpro, followed ADMET properties and were documented. The Coumarin-EM04 molecular scaffold exhibited interactions with catalytical dyad HIS41, CYS145, and neighboring amino acids SER165 and GLN189 in the catalytical site. The crucial factor RMSD was calculated to determine the orientations of Coumarin-EM04. The Coumarin-EM04 complexed with Mpro was found stable in the binding site during MDS. Furthermore, the free energy binding ΔG_bind of Coumarin-EM04 was found to be −187.471 ± 2.230 kJ/mol, and for Remdesivir ΔG_bind was −171.926 ± 2.237 kJ/mol with SARS-CoV-2 Mpro^.

Conclusion

In this study, we identify potent molecules that exhibit interactions with catalytical dyad HIS41 and CYS145 amino acids and unravel Coumarin-EM04 exhibited ΔG_bind higher than Remdesivir against Mpro and thus may serve better antiviral agent against SARS-CoV-2.

Protein photodegradation in the visible range? Insights into protein photooxidation with respect to protein concentration

Visible light (400-800 nm) can lead to photooxidation of protein formulations, which might impair protein integrity. However, the relevant mechanism of photooxidation upon visible light exposure is still unclear for therapeutic proteins, since proteinogenic structures do not absorb light in the visible range. Here, we show that exposure of monoclonal antibody formulations to visible light, lead to the formation of reactive oxygen species (ROS), which subsequently induce specific protein degradations. The formation of ROS and singlet oxygen upon visible light exposure is investigated using electron paramagnetic resonance (EPR) spectroscopy. We describe the initial formation of ROS, most likely after direct reaction of molecular oxygen with a triplet state photosensitizer, generated from intersystem crossing of the excited singlet state. Since these radicals affect the oxygen content in the headspace of the vial, we monitored photooxidation of these mAb formulations. With increasing protein concentrations, we found (i) a decreasing headspace oxygen content in the sample, (ii) a higher relative number of radicals in solution and (iii) a higher protein degradation. Thus, the protein concentration dependence indicates the presence of higher concentration of a currently unknown photosensitizer.

Structural insights into conformational stability and binding of thiazolo-[2,3-b] quinazolinone derivatives with EGFR-TKD and in-vitro study

Heterocyclic molecules are well-known drugs against various diseases including cancer. Many tyrosine kinase inhibitors including erlotinib, osimertinib, and sunitinib were developed and approved but caused adverse effects among treated patients. Which prevents them from being used as cancer therapeutics. In this study, we strategically developed heterocyclic thiazolo-[2,3-b]quinazolinone derivatives by an organic synthesis approach. These synthesized molecules were assessed against the epidermal growth factor receptor tyrosine kinase domain (EGFR-TKD) by in silico methods. Molecular docking simulations unravel derivative 17 showed better binding energy scores and followed Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties. The binding affinity displayed by synthetic congener and reference molecule erlotinib was found to be -8.26 ± 0.0033 kcal/mol and -7.54 ± 0.1411 kcal/mol with the kinase domain. Further, molecular dynamic simulations were conducted thrice to validate the molecular docking study and achieved significant results. Both synthetic derivative and reference molecule attained stability in the active site of the TKD. The synthetic congener and erlotinib showed free energy binding (ΔGbind) -102.975 ± 3.714 kJ/mol and -130.378 ± 0.355 kJ/mol computed by Molecular Mechanics Poison Boltzmann Surface Area (MM-PBSA) method. In addition, the motions of each sampled system including the Apo complex were determined by the principal component analysis and Gibbs energy landscape analysis. The in-vitro apoptosis study was performed using MCF-7 and H-1299 cancer cell lines. However, thiazolo-[2,3-]-quinazoline derivative 17 showed fair anti-proliferative activity against MCF-7 and H-1299. Further, the in-vivo study is necessary to determine the effectivity of the potent anti-proliferative, non-toxic molecule against TKD.

Violaxanthin and Zeaxanthin May Replace Lutein at the L1 Site of LHCII, Conserving the Interactions with Surrounding Chlorophylls and the Capability of Triplet-Triplet Energy Transfer

Carotenoids represent the first line of defence of photosystems against singlet oxygen (¹O₂) toxicity, because of their capacity to quench the chlorophyll triplet state (³Chl) through a physical mechanism based on the transfer of triplet excitation (triplet-triplet energy transfer, TTET). In previous works, we showed that the antenna LHCII is characterised by a robust photoprotective mechanism, able to adapt to the removal of individual chlorophylls while maintaining a remarkable capacity for ³Chl quenching. In this work, we investigated the effects on this quenching induced in LHCII by the replacement of the lutein bound at the L1 site with violaxanthin and zeaxanthin. We studied LHCII isolated from the Arabidopsis thaliana mutants lut2-in which lutein is replaced by violaxanthin-and lut2 npq2, in which all xanthophylls are replaced constitutively by zeaxanthin. We characterised the photophysics of these systems via optically detected magnetic resonance (ODMR) and time-resolved electron paramagnetic resonance (TR-EPR). We concluded that, in LHCII, lutein-binding sites have conserved characteristics, and ensure efficient TTET regardless of the identity of the carotenoid accommodated.

The unique photophysical properties of the Peridinin-Chlorophyll-α-Protein

Peridinin-Chlorophyll-a-Proteins (PCPs) are water-soluble light harvesting complexes from dinoflagellates. They have unique light-harvesting and energy transfer properties which have been studied in details in the last 15 years. This review aims to give an overview on all the main aspects of PCPs photophysics, with an emphasis on some aspects which have not been reviewed in details so far, such as vibrational spectroscopy studies, theoretical calculations, and magnetic resonance studies. A paragraph on the present development of PCPs towards technological applications is also included.

Environmental effects on vibrational properties of carotenoids: experiments and calculations on peridinin

Carotenoids are employed in light-harvesting complexes of dinoflagellates with the two-fold aim to extend the spectral range of the antenna and to protect it from radiation damage. We have studied the effect of the environment on the vibrational properties of the carotenoid peridinin in different solvents by means of vibrational spectroscopies and QM/MM molecular dynamics simulations. Three prototypical solvents were considered: cyclohexane (an apolar/aprotic solvent), deuterated acetonitrile (a polar/aprotic solvent) and methanol (a polar/protic solvent). Thanks to effective normal mode analysis, we were able to assign the experimental Raman and IR bands and to clarify the effect of the solvent on band shifts. In the 1500-1650 cm(-1) region, seven vibrational modes of the polyene chain were identified and assigned to specific molecular vibrations. In the 1700-1800 cm(-1) region a strong progressive down-shift of the lactonic carbonyl frequency is observed passing from cyclohexane to methanol solutions. This has been rationalized here in terms of solvent polarity and solute-solvent hydrogen bond interactions. On the basis of our data we propose a classification of non-equivalent peridinins in the Peridinin-Chlorophyll-Proteins, light-harvesting complexes of dinoflagellates.

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

On the effect of a point mutation on the reactivity of CuZn superoxide dismutase: a theoretical study

In this paper, we investigate the effects of a point mutation on the enzymatic activity of copper-zinc superoxide dismutase, which we recently studied in detail by means of a theoretical-computational procedure. Comparison of the reactivity of the initial catalytic steps in this mutant (G93A mutation far from the active site) with our previous data, reveals the beautiful mechanical-dynamical architecture of the enzyme, altered by such an apparently irrelevant mutation. Finally, our results suggest a possible atomic-molecular-based explanation for the mutant-pathology correlation, in line with the most recent experimental data.

Cyan fluorescent protein: molecular dynamics, simulations, and electronic absorption spectrum

The dynamics and electronic absorption spectrum of enhanced cyan fluorescent protein (ECFP), a mutant of green fluorescent protein (GFP), have been studied by means of a 1 ns molecular dynamics (MD) simulation. The two X-ray conformations A' and B' of ECFP were considered. The chromophore was assumed to be neutral, and all titratable residues were taken in their standard protonation state at neutral pH. The protein was embedded in a box of water molecules (and counterions). The first result is that the two conformations A' and B' are found to be stable all along the simulation. Then, an analysis of the hydrogen-bond networks shows strong differences between the two conformations in the surroundings of the nitrogen atom of the indolic part of the chromophore. This is partly due to the imperfection in the beta barrel near the His148 residue, which allows the access of one solvent molecule inside the protein in conformation A'. Finally, quantum mechanical calculations of the electronic transition energies of the chromophore in the charge cloud of the protein and solvent water molecules were performed using the TDDFT method on 160 snapshots extracted every 5 ps of the MD trajectories. It is found that conformations A' and B' exhibit very similar spectra despite different H-bond networks involving the chromophore. This similarity is related to the weak charge transfer involved in the electronic transition and the weak electrostatic field created by ECFP near the chromophore, within the hypotheses made in the present simulation.

Conformational and Electronic Properties of a Microperoxidase in Aqueous Solution: A Computational Study

A theoretical study of the conformational properties of a small heme peptide in aqueous solution is carried out by classical, long-timescale molecular dynamics simulations. The electronic properties of this species, that is, the relative energies of its excited electronic states and the redox potential, are reproduced and related to the conformational behavior using the perturbed matrix method and basic statistical mechanics. Our results show an interesting coupling between the conformational transitions and the electronic properties. These investigations, beyond the biophysically relevant results addressing the long-standing question of the actual role of the enzyme structure on the enzyme activity, are also of some methodological interest since they offer a further computational perspective for including the electronic degrees of freedom into the modeling of rather complex molecular systems.

Force field impact and spin-probe modeling in molecular dynamics simulations of spin-labeled T4 lysozyme

Several attempts have been made to compute electron paramagnetic resonance (EPR) spectra of biomolecules, using motional models or simulated trajectories to describe dynamics. Ideally, the simulated trajectories should capture "fast" (picosecond) snapshots of spin-probe rotations accurately, while being lengthy enough to ensure a proper Fourier integration of the time-domain signal. It is the interplay of the two criteria that poses computational challenges to the method. In this context, an analysis of the spin-probe and protein conformational sampling and equilibration, with different force fields and with explicit solvent, may be a useful attempt. The present work reports a comparative study of the effect of the molecular dynamics (MD) force field on conformational sampling and equilibration in two spin-labeled T4 lysozyme (T4L) variants, N40C and K48C. Ensembles of 10x 3 ns-trajectories per variant and per force field (OPLS/AMBER and AMBER99) are analyzed for a reliable assessment of convergence and sampling. It is found that subtle site-dependent differences in spin-probe rotations and torsions are more readily captured in the AMBER99 trajectories than in the OPLS/AMBER simulations. On the other hand, sampling and equilibration are found to be better with the OPLS/AMBER force field at equal trajectory lengths.

Theoretical modeling of vibroelectronic quantum states in complex molecular systems: Solvated carbon monoxide, a test case

In this paper we extend the perturbed matrix method by explicitly including the nuclear degrees of freedom, in order to treat quantum vibrational states in a perturbed molecule. In a previous paper we showed how to include, in a simple way, nuclear degrees of freedom for the calculation of molecular polarizability. In the present work we extend and generalize this approach to model vibroelectronic transitions, requiring a more sophisticated treatment.

Ground and excited electronic state thermodynamics of aqueous carbon monoxide: a theoretical study

By using the quasi Gaussian entropy theory in combination with molecular dynamics simulations and the perturbed matrix method, we investigate the ground and excited state thermodynamics of aqueous carbon monoxide. Results show that the model used is rather accurate and provides a great detail in the description of the excitation thermodynamics.

Conformational fluctuations and electronic properties in myoglobin

In this article we use the recently developed perturbed matrix method (PMM) to investigate the effect of conformational fluctuations on the electronic properties of heme in Myoglobin. This widely studied biomolecule has been chosen as a benchmark for evaluating the accuracy of PMM in a large and complex system. Using a long, 80-ns, molecular dynamics simulation and unperturbed Configuration Interaction (CISD) calculations in PMM, we reproduced the main spectroscopic features of deoxy-Myoglobin. Moreover, in line with our previous results on a photosensitive protein, this study reveals a clear dynamical coupling between electronic properties and conformational fluctuations, suggesting that this correlation could be a general feature of proteins.