All-atom molecular dynamics simulation of a photosystem i/detergent complex

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

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.

Lipid shape determination of detergent solubilization in mixed-lipid liposomes

The effects of lipid charge and head group size on liposome partitioning by detergents is an important consideration for applications such as liposomal drug delivery or proteoliposome formation. Yet, the solubilization of mixed-lipid liposomes, those containing multiple types of lipids, by detergents has received insufficient attention. This study examines the incorporation into and subsequent dissolution of mixed-lipid liposomes comprised of both egg phosphatidylcholine (ePC) and egg phosphatidic acid (ePA) by the detergent Triton-X100 (TX). Liposomes were prepared with mixtures of the two lipids, ePC and ePA, at molar ratios from 0 to 1, then step-wise solubilized with TX. Changes in turbidity, size distribution, and molar heat power at constant temperature throughout the solubilization process were assessed. The data suggest that the difference in lipid shapes (shape factors = 0.74 and 1.4 [1,2]) affects packing in membranes, and hence influences how much TX can be incorporated before disruption. As such, liposomes containing the observed ratios of ePA incorporated higher concentrations of TX before initiating dissolution into detergent and lipid mixed-micelles. The cause was concluded to be increased mismatching in the bilayer from the conical shape of ePA compared to the cylindrical shape of ePC. Additionally, the degree to which ePA is approximated as conical versus cylindrical was modulated with pH. It was confirmed that less conical ePA behaved more similarly to ePC than more conical ePA. The understanding gained here on lipid shape in liposome incorporation of TX enables research to use in vitro liposomes that more closely mimic native membranes.

Electron Transfer through the Acceptor Side of Photosystem I: Interaction with Exogenous Acceptors and Molecular Oxygen

This review considers the state-of-the-art on mechanisms and alternative pathways of electron transfer in photosynthetic electron transport chains of chloroplasts and cyanobacteria. The mechanisms of electron transport control between photosystems (PS) I and II and the Calvin-Benson cycle are considered. The redistribution of electron fluxes between the noncyclic, cyclic, and pseudocyclic pathways plays an important role in the regulation of photosynthesis. Mathematical modeling of light-induced electron transport processes is considered. Particular attention is given to the electron transfer reactions on the acceptor side of PS I and to interactions of PS I with exogenous acceptors, including molecular oxygen. A kinetic model of PS I and its interaction with exogenous electron acceptors has been developed. This model is based on experimental kinetics of charge recombination in isolated PS I. Kinetic and thermodynamic parameters of the electron transfer reactions in PS I are scrutinized. The free energies of electron transfer between quinone acceptors A_1A/A_1B in the symmetric redox cofactor branches of PS I and iron-sulfur clusters F_X, F_A, and F_B have been estimated. The second-order rate constants of electron transfer from PS I to external acceptors have been determined. The data suggest that byproduct formation of superoxide radical in PS I due to the reduction of molecular oxygen in the A₁ site (Mehler reaction) can exceed 0.3% of the total electron flux in PS I.

From light-harvesting to photoprotection: structural basis of the dynamic switch of the major antenna complex of plants (LHCII)

Light-Harvesting Complex II (LHCII) is largely responsible for light absorption and excitation energy transfer in plants in light-limiting conditions, while in high-light it participates in photoprotection. It is generally believed that LHCII can change its function by switching between different conformations. However, the underlying molecular picture has not been elucidated yet. The available crystal structures represent the quenched form of the complex, while solubilized LHCII has the properties of the unquenched state. To determine the structural changes involved in the switch and to identify potential quenching sites, we have explored the structural dynamics of LHCII, by performing a series of microsecond Molecular Dynamics simulations. We show that LHCII in the membrane differs substantially from the crystal and has the signatures that were experimentally associated with the light-harvesting state. Local conformational changes at the N-terminus and at the xanthophyll neoxanthin are found to strongly correlate with changes in the interactions energies of two putative quenching sites. In particular conformational disorder is observed at the terminal emitter resulting in large variations of the excitonic coupling strength of this chlorophyll pair. Our results strongly support the hypothesis that light-harvesting regulation in LHCII is coupled with structural changes.

Monte Carlo Modeling of Spectral Diffusion Employing Multiwell Protein Energy Landscapes: Application to Pigment-Protein Complexes Involved in Photosynthesis

We are reporting development and initial applications of the light-induced and thermally induced spectral diffusion modeling software, covering nonphotochemical spectral hole burning (NPHB), hole recovery, and single-molecule spectroscopy and involving random generation of the multiwell protein energy landscapes. The model includes tunneling and activated barrier-hopping in both ground and excited states of a protein-chromophore system. Evolution of such a system is predicted by solving the system of rate equations. Using the barrier parameters from the range typical for the energy landscapes of the pigment-protein complexes involved in photosynthesis, we (a) show that realistic cooling of the sample must result in proteins quite far from thermodynamic equilibrium, (b) demonstrate hole evolution in the cases of burning, fixed-temperature recovery and thermocycling that mostly agrees with the experiment and modeling based on the NPHB master equation, and (c) explore the effects of different protein energy landscapes on the antihole shape. Introducing the multiwell energy landscapes and starting the hole burning experiments in realistic nonequilibrium conditions are not sufficient to explain all experimental observations even qualitatively. Therefore, for instance, one is required to invoke the modified NPHB mechanism where a complex interplay of several small conformational changes is poising the energy landscape of the pigment-protein system for downhill tunneling.

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

Sortase-mediated ligation of PsaE-modified photosystem I from Synechocystis sp. PCC 6803 to a conductive surface for enhanced photocurrent production on a gold electrode

Sortase-mediated ligation was used to attach the photosystem I (PSI) complex from Synechocystis sp. PCC 6803 in a preferential orientation to enhance photoinduced electron flow to a conductive gold surface. Ideally, this method can result in a uniform monolayer of protein, covalently bound unidirectionally to the electrode surface. The exposed C-termini of the psaE subunits of the PSI trimer were targeted to contain an LPETG-sortase recognition sequence to increase noncompeting electron transfer by uniformly orienting the PSI stromal side proximal to the surface. Surface characterization with atomic force microscopy suggested that monolayer formation and optimal surface coverage occurred when the gold surfaces were incubated with peptide at 100 to 500 μM concentrations. When photochronoamperometry with potassium ferrocyanide and ferricyanide as redox mediators was used, photocurrents in the range of 100 to 200 nA/cm(2) were produced, which is an improvement over other attachment techniques for photosystem monolayers that produce approximately 100 nA/cm(2) or less. This work demonstrated that sortase-mediated ligation aided in the control of PSI orientation on modified gold surfaces with a distribution of 94% stromal side proximal and 6% lumenal side proximal to the surface for current-producing PSI.