Molecular Dynamics Simulation Study of Chlorophyll a in Different Organic Solvents

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

Thylakoid Composition Facilitates Chlorophyll a Dimerization through Stronger Interlipid Interactions

Plant thylakoid membrane serves as a crucial matrix for the aggregation of chlororophyll a (CLA) pigments, essential for light harvesting. To understand the role of lipid compositions in the stability of CLA aggregates, dimerization of chlorophyll a molecules (CLA) is studied in the presence of the thylakoid and the bilayers comprising either the least or the highest unsaturated lipids by using coarse-grained molecular dynamics simulations. The thylakoid membrane enhances the stability of the CLA dimer compared with other membranes due to very strong lipid-lipid interactions. The thylakoid exhibits a distinct distribution of lipids around the CLA dimer. Less unsaturated lipids reside in close proximity to the dimer, promoting increased order and efficient packing. Conversely, higher unsaturated lipids are depleted from the dimer, imparting flexibility to the membrane. The combination of tight packing near the dimer and membrane flexibility away from the dimer enhances the stability of the dimer in the thylakoid membrane. Our results suggest that lipid mixing, rather than lipid unsaturation, plays a critical role in facilitating CLA dimerization by modulating the membrane microenvironment through stronger lipid-lipid interactions. These insights will be useful in understanding how lipid compositions affect efficient light absorption and energy transfer during photosynthesis in the future.

The role of chirality and plastic crystallinity in the optical and mechanical properties of chlorosomes

The most efficient light-harvesting antennae found in nature, chlorosomes, are molecular tubular aggregates (TMAs) assembled by pigments without protein scaffolds. Here, we discuss a classification of chlorosomes as a unique tubular plastic crystal and we attribute the robust energy transfer in chlorosomes to this unique nature. To systematically study the role of supramolecular tube chirality by molecular simulation, a role that has remained unresolved, we share a protocol for generating realistic tubes at atomic resolution. We find that both the optical and the mechanical behavior are strongly dependent on chirality. The optical-chirality relation enables a direct interpretation of experimental spectra in terms of overall tube chirality. The mechanical response shows that the overall chirality regulates the hardness of the tube and provides a new characteristic for relating chlorosomes to distinct chirality. Our protocol also applies to other TMA systems and will inspire other systematic studies beyond lattice models.

•

Classifies chlorosomes in terms of a tubular plastic crystal phase

•

Clarifies the unique strategy of chlorosomes for harvesting and transporting energy

•

Presents a protocol for building atom-resolved helical tube structures

•

Maps tube chirality directly to measurable optical and mechanical responses

Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model

Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.

Interactions Determining the Structural Integrity of the Trimer of Plant Light Harvesting Complex in Lipid Membranes

The structural basis for the stability of the trimeric form of the light harvesting complex (LHCII), a pigmented protein from green plants pivotal for photosynthesis, remains elusive till date. The protein embedded in a dipalmitoylphosphatidylcholine (DPPC) lipid membrane is investigated using all-atom molecular dynamics simulations to find out the interactions responsible for the structural integrity of the trimer and its relation to antenna function. Central association of chlorophyll a (CLA) molecules near the LHCII chains is attributed to a conserved coordination between the Mg of CLA and the oxygen of a specific residue of the first helix of a chain. The residue forms a salt-bridge with the fourth helix of the same chain of the trimer, not of the monomer. In an earlier experiment, three residues (WYR) at each chain of the trimer have been found indispensable for the trimerization and referred to as trimerization motif. We find that the residues of the trimerization motif are connected to the lipids or pigments by a chain of interactions rather than a direct contact. Synergistic effects of sequentially located hydrogen bonds and salt-bridges within monomers of the trimer keep the trimer conformation stable in association with the pigments or the lipids. These interactions are exclusively present in the pigmented trimer and not present in the monomer or in the unpigmented trimer. Thus, our results provide a molecular basis for the inherent stability of the LHCII trimer in a lipid membrane and explain many pre-existing experimental data.

Unravelling the Role of Water in Ultrafast Excited-State Relaxation Dynamics within Nano-Architectures of Chlorophyll a

Water plays a pivotal role in structural stability of supramolecular pigment assemblies designed for natural light harvesting (for example, chlorosome antenna complex) as well as their artificial analogs. However, the dynamic role of water in the context of excite-state relaxation has not been explored till date, which we report here. Using femtosecond transient absorption spectroscopy, we investigate the excited-state dynamics of two types of nano-scale assemblies of chlorophyll a with different structural motifs, rod-shaped and micellar assemblies, that depend on the water content. We show how water participates in excess energy dissipation by vibrational cooling of the non-thermally populated Q_y band at different rates in different types of clusters but exhibits no polar solvation dynamics. For the micelles, we observe a bifurcation of stimulated emission line shape, whereas a positive-to-negative switching of differential absorption is observed for the rods; both these observations are correlated with their specific structural aspects. Density functional theory calculations reveal two possible stable ground state geometries of dimers, accounting for the bifurcation of line shape in micelles. Thus, our study elucidates water-mediated structure-function relationship within these pigment assemblies.

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.

Solvation of cholesterol in different solvents: a molecular dynamics simulation study

To the best of our knowledge, molecular dynamics simulations of an isolated cholesterol immersed in four different solvents of varying polarity, such as water, methanol, dimethyl sulfoxide and benzene, were reported for the first time to gain insights into the structural and dynamical properties. The study was mainly focused on the evaluation of solvation of cholesterol with respect to its hydrophilic and hydrophobic structural components in the form of respective functional groups interacting with the solvents. Structural evaluations suggested that both hydrophilic and hydrophobic groups of cholesterol were interacting with the solvents, in particular methanol and dimethyl sulfoxide, which presented both types of interactions that are polar and non-polar. On the other hand, the highly polar water and non-polar benzene demonstrated extreme solvation behavior, since water was involved only in hydrogen bonding to the solute hydroxyl group and non-polar benzene formed strong van der Waals interactions only. Furthermore, the hydrophobic effect of cholesterol was also analyzed mainly in polar solvents, as the effect was more pronounced in the polar environment thereby preventing the solvent mobility in the solvation layer(s). The dynamical properties in terms of lateral diffusion and hydrogen bond dynamics as well as free energies of solvation also corroborated the findings based on the structural data and the hydrophobic character of cholesterol was later quantified by the computation of the averaged solvent accessible surface area. The polarity effect of the solvents on the aggregation property of cholesterol was further investigated, which is of big concern from the clinical point of view due to its major role in cardiovascular ailments. It was another major finding of the present study that aggregation was shown to be facilitated by highly polar solvents like water.

Vibronic structure of photosynthetic pigments probed by polarized two-dimensional electronic spectroscopy and ab initio calculations

Bacteriochlorophyll a (Bchl a) and chlorophyll a (Chl a) play important roles as light absorbers in photosynthetic antennae and participate in the initial charge-separation steps in photosynthetic reaction centers. Despite decades of study, questions remain about the interplay of electronic and vibrational states within the Q-band and its effect on the photoexcited dynamics. Here we report results of polarized two-dimensional electronic spectroscopic measurements, performed on penta-coordinated Bchl a and Chl a and their interpretation based on state-of-the-art time-dependent density functional theory calculations and vibrational mode analysis for spectral shapes. We find that the Q-band of Bchl a is comprised of two independent bands, that are assigned following the Gouterman model to Q x and Q y states with orthogonal transition dipole moments. However, we measure the angle to be ∼75°, a finding that is confirmed by ab initio calculations. The internal conversion rate constant from Q x to Q y is found to be 11 ps-1. Unlike Bchl a, the Q-band of Chl a contains three distinct peaks with different polarizations. Ab initio calculations trace these features back to a spectral overlap between two electronic transitions and their vibrational replicas. The smaller energy gap and the mixing of vibronic states result in faster internal conversion rate constants of 38-50 ps-1. We analyze the spectra of penta-coordinated Bchl a and Chl a to highlight the interplay between low-lying vibronic states and their relationship to photoinduced relaxation. Our findings shed new light on the photoexcited dynamics in photosynthetic systems where these chromophores are primary pigments.

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.

Dimerization Behavior of Methyl Chlorophyllide a as the Model of Chlorophyll a in the Presence of Water Molecules-Theoretical Study

A dimerization of methyl chlorophyllide a molecules and a role of water in stabilization and properties of methyl chlorophyllide a dimers were studied by means of symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), density-functional theory (DFT), and time-dependent DFT approaches. The quantification of various types of interactions, such as π-π stacking, coordinative, and hydrogen bonding by applying the F-SAPT energy decomposition scheme shows the major role of the magnesium atom and the pheophytin macrocycle in the stability of the complex. The examination of interaction energy components with respect to a mutual orientation of monomers and in the presence or absence of water molecules reveals that the dispersion energy is the main binding factor of the interaction, while water molecules tend to weaken the attraction between methyl chlorophyllide a species. The dimerization can be seen in computed UV-vis spectra, and results in a doubling of the lowest peaks, as compared to the monomer spectrum, and in an intensity rise of the lowest 1.8 and 2.4 eV peaks at a cost of the 3.5 eV peaks for the majority of dimer configurations. The complexation of water has little effect on the peaks' position; however, it affects the overall shape of simulated spectra through changes in peak intensities, which is strongly dependent on the structure of the complex. The VCD spectra for the dimers show several characteristic features attributed to the interaction of substituting groups and/or water ligand attached to macrocycle groups belonging to different monomers. VCD is sensitive to the type of the formed dimer, but not to the number of water molecules it contains. This and several other features, as well as the differential UV-vis spectra, may serve as the indicator of the presence of a given dimer structure in the experiment.

Contrasting Modes of Self-Assembly and Hydrogen-Bonding Heterogeneity in Chlorosomes of Chlorobaculum tepidum

Chlorosome antennae form an interesting class of materials for studying the role of structural motifs and dynamics in nonadiabatic energy transfer. They perform robust and highly quantum-efficient transfer of excitonic energy while allowing for compositional variation and completely lacking the usual regulatory proteins. Here, we first cast the geometrical analysis for ideal tubular scaffolding models into a formal framework, to relate effective helical properties of the assembly structures to established characterization data for various types of chlorosomes. This analysis shows that helicity is uniquely defined for chlorosomes composed of bacteriochlorophyll (BChl) d and that three chiral angles are consistent with the nuclear magnetic resonance (NMR) and electron microscope data for BChl c, including two novel ones that are at variance with current interpretations of optical data based on perfect cylindrical symmetry. We use this information as a starting point for investigating dynamic and static heterogeneity at the molecular level by unconstrained molecular dynamics. We first identify a rotational degree of freedom, along the Mg-OH coordination bond, that alternates along the syn-anti stacks and underlies the (flexible) curvature on a larger scale. Because rotation directly relates to the formation or breaking of interstack hydrogen bonds of the O-H···O=C structural motif along the syn-anti stacks, we analyzed the relative fractions of hydrogen-bonded and the nonbonded regions, forming stripe domains in otherwise spectroscopically homogeneous curved slabs. The ratios 7:3 for BChl c and 9:1 for BChl d for the two distinct structural components agree well with the signal intensities determined by NMR. In addition, rotation with curvature-independent formation of stripe domains offers a viable explanation for the localization and dispersion of exciton states over two fractions, as observed in single chlorosome fluorescence decay studies.

The Photosystem II Subunit S under Stress

Nonphotochemical quenching is the protective mechanism against overexcitation of photosystem II, triggered by excess ΔpH in photosynthetic membranes. The light-harvesting complexes (LHCs), the de-epoxidation of violaxanthin to zeaxanthin, and the photosystem II subunit S (PsbS) work in synergy for an optimized multilevel response. Understanding the fine details of this synergy has proven challenging to scientific research. Here, we employ large-scale, all-atom molecular simulations and beyond experimental insight, we proceed a step further in identifying the PsbS dynamics that could possibly be associated with this synergy. For the first time, to our knowledge, we probe the distinct behavior of PsbS under ΔpH that probes the details of the potential dimer-to-monomer transition, and in a violaxanthin/zeaxanthin-rich membrane, at an all-atom resolution. We propose that the lumen-exposed residues, threonine 162 and glutamic acid 173, form stabilizing hydrogen bonds between the PsbS monomers only at high lumen pH, whereas at low pH (excess ΔpH) this interaction is lost, and leads to higher flexibility of the protein and potentially to the dimer-to-monomer transition. Lastly, we discuss how conformational changes under the presence of ΔpH/zeaxanthin are related to the PsbS role in the current nonphotochemical quenching model in the literature. For the latter, we probe a PsbS-monomeric LHCII association. The association is proposed to potentially alter the monomeric LHCII sensitivity to ΔpH by changing the pKa values of interacting LHCII residues. This serves as an example where protonation-ligation events enhance protein-protein interactions fundamental to many life processes.

Dust evolution, a global view I. Nanoparticles, nascence, nitrogen and natural selection … joining the dots

The role and importance of nanoparticles for interstellar chemistry and beyond is explored within the framework of The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS), focusing on their active surface chemistry, the effects of nitrogen doping and the natural selection of interesting nanoparticle sub-structures. Nanoparticle-driven chemistry, and in particular the role of intrinsic epoxide-type structures, could provide a viable route to the observed gas phase OH in tenuous interstellar clouds en route to becoming molecular clouds. The aromatic-rich moieties present in asphaltenes probably provide a viable model for the structures present within aromatic-rich interstellar carbonaceous grains. The observed doping of such nanoparticle structures with nitrogen, if also prevalent in interstellar dust, could perhaps have important and observable consequences for surface chemistry and the formation of precursor pre-biotic species.

Energizing the light harvesting antenna: Insight from CP29

How do plants cope with excess light energy? Crop health and stress tolerance are governed by molecular photoprotective mechanisms. Protective exciton quenching in plants is activated by membrane energization, via unclear conformational changes in proteins called antennas. Here we show that pH and salt gradients stimulate the response of such an antenna under low and high energization by all-atom Molecular Dynamics Simulations. Novel insight establishes that helix-5 (H5) conformation in CP29 from spinach is regulated by chemiosmotic factors. This is selectively correlated with the chl-614 macrocycle deformation and interactions with nearby pigments, that could suggest a role in plant photoprotection. Adding to the significance of our findings, H5 domain is conserved among five antennas (LHCB1-5). These results suggest that light harvesting complexes of Photosystem II, one of the most abundant proteins on earth, can sense chemiosmotic gradients via their H5 domains in an upgraded role from a solar detector to also a chemiosmotic sensor.

Optical absorption and electronic spectra of chlorophylls a and b

We report optical and electronic properties of the two main chlorophylls in green plants, namely, chlorophylls a and b. We estimate the electric moments of these molecules and study absorption spectra of the chlorophylls.

Derivation of coarse-grained simulation models of chlorophyll molecules in lipid bilayers for applications in light harvesting systems

The correct interplay of interactions between protein, pigment and lipid molecules is highly relevant for our understanding of the association behavior of the light harvesting complex (LHCII) of green plants. To cover the relevant time and length scales in this multicomponent system, a multi-scale simulation ansatz is employed that subsequently uses a classical all atomistic (AA) model to derive a suitable coarse grained (CG) model which can be backmapped into the AA resolution, aiming for a seamless conversion between two scales. Such an approach requires a faithful description of not only the protein and lipid components, but also the interaction functions for the indispensable pigment molecules, chlorophyll b and chlorophyll a (referred to as chl b/chl a). In this paper we develop a CG model for chl b and chl a in a dipalmitoylphosphatidyl choline (DPPC) bilayer system. The structural properties and the distribution behavior of chl within the lipid bilayer in the CG simulations are consistent with those of AA reference simulations. The non-bonded potentials are parameterized such that they fit to the thermodynamics based MARTINI force-field for the lipid bilayer and the protein. The CG simulation shows chl aggregation in the lipid bilayer which is supported by fluorescence quenching experiments. It is shown that the derived chl model is well suited for CG simulations of stable, structurally consistent, trimeric LHCII and can in the future be used to study their large scale aggregation behavior.

Atomistic and Coarse Grain Topologies for the Cofactors Associated with the Photosystem II Core Complex

Electron transfers within and between protein complexes are core processes of the electron transport chains occurring in thylakoid (chloroplast), mitochondrial, and bacterial membranes. These electron transfers involve a number of cofactors. Here we describe the derivation of molecular mechanics parameters for the cofactors associated with the function of the photosystem II core complex: plastoquinone, plastoquinol, heme b, chlorophyll A, pheophytin, and β-carotene. Parameters were also obtained for ubiquinol and ubiquinone, related cofactors involved in the respiratory chain. Parameters were derived at both atomistic and coarse grain (CG) resolutions, compatible with the building blocks of the GROMOS united-atom and Martini CG force fields, respectively. Structural and thermodynamic properties of the cofactors were compared to experimental values when available. The topologies were further tested in molecular dynamics simulations of the cofactors in their physiological environment, e.g., either in a lipid membrane environment or in complex with the heme binding protein bacterioferritin.

Aggregation/disaggregation of chlorophyll a in model phospholipid-detergent vesicles and micelles

The photosynthetic pigments of higher plants exist in complex oligomeric states, which are difficult to study in vivo. To investigate aggregation processes of chlorophyll a (Chl a), we used an in vitro reconstitution procedure, with this pigment incorporated into liposomes of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), micelles and pre-micelle media of the detergent n-dodecyltrimethylammonium chloride (DTAC), and mixed, spontaneous, DMPC-DTAC vesicles and micelles. Chl a oligomers were characterized by UV-visible absorption, steady-state and time-resolved fluorescence, and fluorescence lifetime imaging microscopy. Equivalent diameters of the colloidal structures were obtained by fluorescence correlation spectroscopy. In DMPC liposomes and DMPC-DTAC vesicles and micelles, three fluorescence lifetimes indicated the coexistence of Chl a monomers (≈5 ns) and oligomers (≈1-2 to ≈0.1 ns). The increase in DTAC amount, in the mixed system, induces a progressive solubilization of DMPC liposomes (from vesicles to micelles) and simultaneous disruption of Chl a aggregates; in pure DTAC micelles, mostly monomers were found. The present work aims for a better understanding of chlorophyll-chlorophyll (Chl-Chl), Chl-lipid, and Chl-detergent interactions in spontaneous colloidal micro- and nanostructures.

Symmetry-adapted perturbation theory study for some magnesium complexes

A systematic theoretical study on Mg–ligand interactions has been carried out employing both ab initio correlated wave function and density functional methods. The interactions of the Mg(CH3N2)2moiety with BF, CO, N2, NH3, and H2O ligands have been investigated by performing calculations at the B3LYP, MP2, MP4, and CCSD(T)/6–311++G(3df,3pd) levels of theory. Results indicate that the interaction energies of the Mg(CH3N2)2–L complexes increase in the order NH3 > H2O > BF > CO > N2. Symmetry-adapted perturbation theory (SAPT) analysis has been carried out to understand the nature of the forces involved in the bonding. The SAPT results indicate that the stabilities of the Mg–L interactions are attributed mainly to electrostatic effects, while induction and dispersion forces also play a significant role. The evaluated SAPT interaction energies for the Mg(CH3N2)2–L complexes are generally in good agreement with those obtained using the supermolecule CCSD(T) methods, suggesting that SAPT is a proper method to study the intermolecular interactions in these complexes. The results also suggest an explanation for the unique role of Mg2+as a carrier of water molecules that mediate enzymatic hydrolysis reactions.