Docking of photosystem I subunit C using a constrained geometric simulation

Competition between intra-protein charge recombination and electron transfer outside photosystem I complexes used for photovoltaic applications

Photosystem I (PSI) complexes isolated from three different species were electrodeposited on FTO conducting glass, forming a photoactive multilayer of the photo-electrode, for investigation of intricate electron transfer (ET) properties in such green hybrid nanosystems. The internal quantum efficiency of photo-electrochemical cells (PEC) containing the PSI-based photo-electrodes did not exceed ~ 0.5%. To reveal the reason for such a low efficiency of photocurrent generation, the temporal evolution of the transient concentration of the photo-oxidized primary electron donor, P⁺, was studied in aqueous suspensions of the PSI complexes by time-resolved absorption spectroscopy. The results of these measurements provided the information on: (1) completeness of charge separation in PSI reaction centers (RCs), (2) dynamics of internal charge recombination, and (3) efficiency of electron transfer from PSI to the electrolyte, which is the reaction competing with the internal charge recombination in the PSI RC. The efficiency of the full charge separation in the PSI complexes used for functionalization of the electrodes was ~ 90%, indicating that incomplete charge separation was not the main reason for the small yield of photocurrents. For the PSI particles isolated from a green alga Chlamydomonas reinhardtii, the probability of ET outside PSI was ~ 30-40%, whereas for their counterparts isolated from a cyanobacterium Synechocystis sp. PCC 6803 and a red alga Cyanidioschyzon merolae, it represented a mere ~ 4%. We conclude from the transient absorption data for the PSI biocatalysts in solution that the observed small photocurrent efficiency of ~ 0.5% for all the PECs analyzed in this study is likely due to: (1) limited efficiency of ET outside PSI, particularly in the case of PECs based on PSI from Synechocystis and C. merolae, and (2) the electrolyte-mediated electric short-circuiting in PSI particles forming the photoactive layer, particularly in the case of the C. reinhardtii PEC.

Flexibility and mobility of SARS-CoV-2-related protein structures

The worldwide CoVid-19 pandemic has led to an unprecedented push across the whole of the scientific community to develop a potent antiviral drug and vaccine as soon as possible. Existing academic, governmental and industrial institutions and companies have engaged in large-scale screening of existing drugs, in vitro, in vivo and in silico. Here, we are using in silico modelling of possible SARS-CoV-2 drug targets, as deposited on the Protein Databank (PDB), and ascertain their dynamics, flexibility and rigidity. For example, for the SARS-CoV-2 spike protein-using its complete homo-trimer configuration with 2905 residues-our method identifies a large-scale opening and closing of the S1 subunit through movement of the S[Formula: see text] domain. We compute the full structural information of this process, allowing for docking studies with possible drug structures. In a dedicated database, we present similarly detailed results for the further, nearly 300, thus far resolved SARS-CoV-2-related protein structures in the PDB.

The flexibility and dynamics of protein disulfide isomerase

We have studied the mobility of the multidomain folding catalyst, protein disulfide isomerase (PDI), by a coarse‐graining approach based on flexibility. We analyze our simulations of yeast PDI (yPDI) using measures of backbone movement, relative positions and orientations of domains, and distances between functional sites. We find that there is interdomain flexibility at every interdomain junction but these show very different characteristics. The extent of interdomain flexibility is such that yPDI's two active sites can approach much more closely than is found in crystal structures—and indeed hinge motion to bring these sites into proximity is the lowest energy normal mode of motion of the protein. The flexibility predicted for yPDI (based on one structure) includes the other known conformation of yPDI and is consistent with (i) the mobility observed experimentally for mammalian PDI and (ii) molecular dynamics. We also observe intradomain flexibility and clear differences between the domains in their propensity for internal motion. Our results suggest that PDI flexibility enables it to interact with many different partner molecules of widely different sizes and shapes, and highlights considerable similarities of yPDI and mammalian PDI. Proteins 2016; 84:1776–1785. © 2016 Wiley Periodicals, Inc.

Residue Geometry Networks: A Rigidity-Based Approach to the Amino Acid Network and Evolutionary Rate Analysis

Amino acid networks (AANs) abstract the protein structure by recording the amino acid contacts and can provide insight into protein function. Herein, we describe a novel AAN construction technique that employs the rigidity analysis tool, FIRST, to build the AAN, which we refer to as the residue geometry network (RGN). We show that this new construction can be combined with network theory methods to include the effects of allowed conformal motions and local chemical environments. Importantly, this is done without costly molecular dynamics simulations required by other AAN-related methods, which allows us to analyse large proteins and/or data sets. We have calculated the centrality of the residues belonging to 795 proteins. The results display a strong, negative correlation between residue centrality and the evolutionary rate. Furthermore, among residues with high closeness, those with low degree were particularly strongly conserved. Random walk simulations using the RGN were also successful in identifying allosteric residues in proteins involved in GPCR signalling. The dynamic function of these residues largely remain hidden in the traditional distance-cutoff construction technique. Despite being constructed from only the crystal structure, the results in this paper suggests that the RGN can identify residues that fulfil a dynamical function.

Constrained geometric dynamics of the Fenna-Matthews-Olson complex: the role of correlated motion in reducing uncertainty in excitation energy transfer

The trimeric Fenna-Mathews-Olson (FMO) complex of green sulphur bacteria is a well-studied example of a photosynthetic pigment-protein complex, in which the electronic properties of the pigments are modified by the protein environment to promote efficient excitonic energy transfer from antenna complexes to the reaction centres. By a range of simulation methods, many of the electronic properties of the FMO complex can be extracted from knowledge of the static crystal structure. However, the recent observation and analysis of long-lasting quantum dynamics in the FMO complex point to protein dynamics as a key factor in protecting and generating quantum coherence under laboratory conditions. While fast inter- and intra-molecular vibrations have been investigated extensively, the slow, conformational dynamics which effectively determine the optical inhomogeneous broadening of experimental ensembles has received less attention. The following study employs constrained geometric dynamics to study the flexibility in the protein network by efficiently generating the accessible conformational states from the published crystal structure. Statistical and principle component analyses reveal highly correlated low frequency motions between functionally relevant elements, including strong correlations between pigments that are excitonically coupled. Our analysis reveals a hierarchy of structural interactions which enforce these correlated motions, from the level of monomer-monomer interfaces right down to the α-helices, β-sheets and pigments. In addition to inducing strong spatial correlations across the conformational ensemble, we find that the overall rigidity of the FMO complex is exceptionally high. We suggest that these observations support the idea of highly correlated inhomogeneous disorder of the electronic excited states, which is further supported by the remarkably low variance (typically <5%) of the excitonic couplings of the conformational ensemble.

Constrained geometric simulation of the nicotinic acetylcholine receptor

Constrained geometric simulations have been performed for the recently published closed-channel state of the nicotinic acetylcholine receptor. These simulations support the theory that correlated motion in the flexible β-sheet structure of the extracellular domain helps to communicate a "conformational wave", spreading from the acetylcholine binding pocket. Furthermore, we have identified key residues that act at the interface between subunits and between domains that could potentially facilitate rapid communication between the binding site and the transmembrane gate.

Protein flexibility is key to cisplatin crosslinking in calmodulin

Chemical crosslinking in combination with Fourier transform ion cyclotron resonance mass spectrometry (FTICR MS) has significant potential for studying protein structures and protein-protein interactions. Previously, cisplatin has been shown to be a crosslinker and crosslinks multiple methionine (Met) residues in apo-calmodulin (apo-CaM). However, the inter-residue distances obtained from nuclear magnetic resonance structures are inconsistent with the measured distance constraints by crosslinking. Met residues lie too far apart to be crosslinked by cisplatin. Here, by combining FTICR MS with a novel computational flexibility analysis, the flexible nature of the CaM structure is found to be key to cisplatin crosslinking in CaM. It is found that the side chains of Met residues can be brought together by flexible motions in both apo-CaM and calcium-bound CaM (Ca₄-CaM). The possibility of cisplatin crosslinking Ca₄-CaM is then confirmed by MS data. Therefore, flexibility analysis as a fast and low-cost computational method can be a useful tool for predicting crosslinking pairs in protein crosslinking analysis and facilitating MS data analysis. Finally, flexibility analysis also indicates that the crosslinking of platinum to pairs of Met residues will effectively close the nonpolar groove and thus will likely interfere with the binding of CaM to its protein targets, as was proved by comparing assays for cisplatin-modified/unmodified CaM binding to melittin. Collectively, these results suggest that cisplatin crosslinking of apo-CaM or Ca₄-CaM can inhibit the ability of CaM to recognize its target proteins, which may have important implications for understanding the mechanism of tumor resistance to platinum anticancer drugs.

Differential degradation of photosystem I subunits under iron deficiency in rice

Rice (Oryza sativa) is one of the staple foods of the world. Iron (Fe) deficiency is a major abiotic stress factor that contributes world-wide to losses in crop yield and decline in nutritional quality. As cofactor for many enzymes and proteins, iron is an essential element. It plays a pivotal role in chlorophyll (Chl) biosynthesis, and iron deficiency may result in decreased Chl production and, thus, reduced photosynthetic capacity. Photosystem I (PSI) is a prime target of iron deficiency because of its high iron content (12 Fe per PS). To understand the protein level changes in the light-harvesting complex (LHC) of PSI (LHCI) under iron deficiency, rice seedlings were grown in Hoagland's nutrient medium with and without Fe. Chlorophyll content and photosynthetic efficiency decreased under iron deficiency. Protein gel blots probed with antibodies against the PSI core and Lhca 1-4 proteins revealed that the core subunits PsaA and PsaB remained stable under iron deficiency, whereas PsaC and PsaD decreased by about 50%, and PsaE was completely degraded. Among the LHCI subunits, Lhca1 and Lhca2 decreased by 40 and 50%, respectively, whereas Lhca3 and Lhca4 were completely degraded. We propose that the dissociation of LHCI subunits may be due to increased levels of reactive oxygen species, which is suggested by the increased activity of superoxide dismutase.

Rapid simulation of protein motion: merging flexibility, rigidity and normal mode analyses

Protein function frequently involves conformational changes with large amplitude on timescales which are difficult and computationally expensive to access using molecular dynamics. In this paper, we report on the combination of three computationally inexpensive simulation methods--normal mode analysis using the elastic network model, rigidity analysis using the pebble game algorithm, and geometric simulation of protein motion--to explore conformational change along normal mode eigenvectors. Using a combination of ElNemo and First/Froda software, large-amplitude motions in proteins with hundreds or thousands of residues can be rapidly explored within minutes using desktop computing resources. We apply the method to a representative set of six proteins covering a range of sizes and structural characteristics and show that the method identifies specific types of motion in each case and determines their amplitude limits.

A normal mode-based geometric simulation approach for exploring biologically relevant conformational transitions in proteins

A three-step approach for multiscale modeling of protein conformational changes is presented that incorporates information about preferred directions of protein motions into a geometric simulation algorithm. The first two steps are based on a rigid cluster normal-mode analysis (RCNMA). Low-frequency normal modes are used in the third step (NMSim) to extend the recently introduced idea of constrained geometric simulations of diffusive motions in proteins by biasing backbone motions of the protein, whereas side-chain motions are biased toward favorable rotamer states. The generated structures are iteratively corrected regarding steric clashes and stereochemical constraint violations. The approach allows performing three simulation types: unbiased exploration of conformational space; pathway generation by a targeted simulation; and radius of gyration-guided simulation. When applied to a data set of proteins with experimentally observed conformational changes, conformational variabilities are reproduced very well for 4 out of 5 proteins that show domain motions, with correlation coefficients r > 0.70 and as high as r = 0.92 in the case of adenylate kinase. In 7 out of 8 cases, NMSim simulations starting from unbound structures are able to sample conformations that are similar (root-mean-square deviation = 1.0-3.1 Å) to ligand bound conformations. An NMSim generated pathway of conformational change of adenylate kinase correctly describes the sequence of domain closing. The NMSim approach is a computationally efficient alternative to molecular dynamics simulations for conformational sampling of proteins. The generated conformations and pathways of conformational transitions can serve as input to docking approaches or as starting points for more sophisticated sampling techniques.

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 A. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.

Structural and functional changes of PSI-LHCI supercomplexes of Chlamydomonas reinhardtii cells grown under high salt conditions

The eVect of high salt concentration (100 mM NaCl) on the organization of photosystem I-light harvesting complex I supercomplexes (PSI-LHCI) of Chlamydomonas reinhardtii was studied. The electron transfer activity was reduced by 39% in isolated PSI-LHCI supercomplexes. The visible circular dichroism (CD) spectra associated with strongly coupled chlorophyll (Chl) dimers were reduced in intensity, indicating that pigment-pigment interactions were disrupted. This data is consistent with results from Xuorescence streak camera spectroscopy, which suggest that red-shifted pigments in the PSI-LHCI antenna had been lost. Denaturing gel electrophoresis and immunoblot analysis reveals that levels of the PSI reaction center proteins PsaD, PsaE and PsaF were reduced due to salt stress. PsaE is almost completely absent under high salt conditions. It is known that the membrane-extrinsic subunits PsaD and E form the ferredoxin-docking site. Our results indicate that the PSI-LHCI supercomplex is damaged by reactive oxygen species at high salt concentration, with particular impact on the ferredoxin-docking site and the PSILHCI interface.

Unconstrained motions, dynamic heterogeneities, and relaxation in disordered solids

A disordered network of bonds with a fixed configuration can relax via a variety of unconstrained motions. These motions can be directly inferred from the topological arrangement of constraints without any geometrical information. We use the pebble game algorithm of Jacobs and Thorpe [D. J. Jacobs and M. F. Thorpe, Phys. Rev. Lett. 75, 4051 (1995)] to decompose the system into separate rigid clusters and identify the remaining degrees of freedom. Unconstrained motions can then be resolved in the form of translations and rotations of isolated groups of bonds and the internal motion within bond groups. We show that each motion can be assigned a characteristic thermal velocity and hence a relaxation time scale. We use this information to construct a relaxation function and also examine the spatial distribution of relaxation time scales. We investigate the sensitivity of the relaxation time scales and their spatial distribution when making individual bond changes in the system, and we consider the dependence of these time scales on the underlying structure.

Rigidity percolation and the spatial heterogeneity of soft modes in disordered materials

We show how the spatial character of unconstrained motion in a network of bonds can be directly inferred from the topological arrangement of constraints. Relaxation time scales of these soft modes are determined, and from this information we generate spatial maps of the heterogeneous distribution of relaxation times in the disordered network. We show that the nature of the dynamic heterogeneity and its sensitivity to changes in bond configuration depends dramatically on the proximity of the system to the rigidity percolation point.

Understanding of the binding interface between PsaC and the PsaA/PsaB heterodimer in photosystem I

The PsaC subunit of Photosystem I (PS I) is tightly bound to the PsaA/PsaB heterodimer via an extensive network of ionic and hydrogen bonds. To improve our understanding of the design of the PsaC-PsaA/PsaB binding interface, variants of PsaC were generated, each lacking a key binding contact with the PsaA/PsaB heterodimer. The characteristics of the reconstituted, variant PS I complexes were monitored by time-resolved optical spectroscopy, low-temperature EPR spectroscopy, and electron transfer throughput measurements. In the absence of the ionic bond forming contacts R52(C) or R65(C), a markedly slower charge recombination occurs between P(700)(+) and [F(A)/F(B)](-). The addition of PsaD leads to the restoration of native recombination kinetics in a fraction of the PS I complexes reconstituted with R52A(C), but not with R65A(C). Contrary to expectation, the absence of Y80(C), which forms two symmetry-breaking H-bonds with PsaB, does not significantly affect the binding of PsaC as judged by the rate of charge recombination between P(700)(+) and [F(A)/F(B)](-). However, the removal of the entire C-terminus results in a dramatic decrease in the rate of charge recombination. Low-temperature EPR spectra of the variant PS I complexes indicate that the magnetic environments of F(A) and F(B) are altered when compared to that of native PS I. The slowing of the rate of charge recombination in the variant PS I complexes could be due to an increase in the distance between F(X) and F(A)/F(B) as the result of non-native binding or to an altered reduction potential of the iron-sulfur clusters, which would result in a different rate of thermalization up the electron acceptor chain. The most significant finding is that the variant PS I complexes support lower rates of light-induced flavodoxin reduction and that the rates deteriorate rapidly on exposure to dioxygen due to the degradation of F(A) and F(B). We suggest that the extensive set of ionic bonds and H-bonds between PsaC and the PsaA/PsaB heterodimer has evolved to ensure an exceedingly tight binding interface, thereby rendering the [4Fe-4S] clusters in PsaC inaccessible to dioxygen at the onset of oxygenic photosynthesis.

Intrinsic curvature properties of photosynthetic proteins in chromatophores

In purple bacteria, photosynthesis is carried out on large indentations of the bacterial plasma membrane termed chromatophores. Acting as primitive organelles, chromatophores are densely packed with the membrane proteins necessary for photosynthesis, including light harvesting complexes LH1 and LH2, reaction center (RC), and cytochrome bc(1). The shape of chromatophores is primarily dependent on species, and is typically spherical or flat. How these shapes arise from the protein-protein and protein-membrane interactions is still unknown. Now, using molecular dynamics simulations, we have observed the dynamic curvature of membranes caused by proteins in the chromatophore. A membrane-embedded array of LH2s was found to relax to a curved state, both for LH2 from Rps. acidophila and a homology-modeled LH2 from Rb. sphaeroides. A modeled LH1-RC-PufX dimer was found to develop a bend at the dimerizing interface resulting in a curved shape as well. In contrast, the bc(1) complex, which has not been imaged yet in native chromatophores, did not induce a preferred membrane curvature in simulation. Based on these results, a model for how the different photosynthetic proteins influence chromatophore shape is presented.

Structure of Photosystems I and II

Photosynthesis is the major process that converts solar energy into chemical energy on Earth. Two and a half billion years ago, the ancestors of cyanobacteria were able to use water as electron source for the photosynthetic process, thereby evolving oxygen and changing the atmosphere of our planet Earth. Two large membrane protein complexes, Photosystems I and II, catalyze the primary step in this energy conversion, the light-induced charge separation across the photosynthetic membrane. This chapter describes and compares the structure of two Photosystems and discusses their function in respect to the mechanism of light harvesting, electron transfer and water splitting.

Fitting low-resolution cryo-EM maps of proteins using constrained geometric simulations

Recent experimental advances in producing density maps from cryo-electron microscopy (cryo-EM) have challenged theorists to develop improved techniques to provide structural models that are consistent with the data and that preserve all the local stereochemistry associated with the biomolecule. We develop a new technique that maintains the local geometry and chemistry at each stage of the fitting procedure. A geometric simulation is used to drive the structure from some appropriate starting point (a nearby experimental structure or a modeled structure) toward the experimental density, via a set of small incremental motions. Structural motifs such as alpha-helices can be held rigid during the fitting procedure as the starting structure is brought into alignment with the experimental density. After validating this procedure on simulated data for adenylate kinase and lactoferrin, we show how cryo-EM data for two different GroEL structures can be fit using a starting x-ray crystal structure. We show that by incorporating the correct local stereochemistry in the modeling, structures can be obtained with effective resolution that is significantly higher than might be expected from the nominal cryo-EM resolution.

Molecular docking and spatial coarse graining simulations as tools to investigate substrate recognition, enhancer binding and conformational transitions in indoleamine-2,3-dioxygenase (IDO)

Indoleamine 2,3-dioxygenase (IDO) is an heme-containing enzyme involved in the regulation of important immunological responses and neurological processes. The enzyme catalyzes the oxidative cleavage of the pyrrole ring of the indole nucleus of tryptophan (Trp) to yield N-formylkynurenine, that is the initial and rate limiting step of the kynurenine pathway. Some indole derivatives have been reported to act as effectors of the enzyme by enhancing its catalytic activity. On the basis of the recent availability of the crystal structure of IDO, in this work we investigate substrate recognition and enhancer binding to IDO using molecular docking experiments. In addition, conformational transitions of IDO in response to substrate and enhancer binding are studied using coarse graining simulations with the program FIRST. The results enable us to identify (i) the binding site of enhancer modulators; (ii) the motion of an electrostatic gate that regulates the access of the substrate to the catalytic site of the enzyme; (iii) the movement of the anchoring region of a hairpin loop that may assist the shuttle of substrates/products to/from the catalytic site of IDO. These data, combined with available site-directed mutagenesis experiments, reveal that conformational transitions of IDO in response to substrate and enhancer binding are controlled by distinct combination of two conformational states (open and close) of the above structural motifs. On this basis, a molecular mechanism regarding substrate recognition and activity enhancement by indole derivatives is proposed.