Profile Hidden Markov Models for Analyzing Similarities and Dissimilarities in the Bacterial Reaction Center and Photosystem II

The bacterial photosynthetic reaction center is the evolutionary ancestor of the Photosystem II reaction center. These proteins share the same fold and perform the same biological function. Nevertheless, the details of their molecular reaction mechanism differ. It is of significant biological and biochemical interest to determine which functional characteristics are conserved at the level of the protein sequences. Since the level of sequence identity between the bacterial photosynthetic reaction center and Photosystem II is low, a progressive multiple-sequence alignment leads to errors in identifying the conserved residues. In such a situation, profile hidden Markov models (pHMM) can be used to obtain reliable multiple-sequence alignments. We therefore constructed the pHMM with the help of a sequence alignment based on a structural superposition of both proteins. To validate the multiple-sequence alignments obtained with the pHMM, the conservation of residues with known functional importance was examined. Having confirmed the correctness of the multiple-sequence alignments, we analyzed the conservation of residues involved in hydrogen bonding and redox potential tuning of the cofactors. Our analysis reveals similarities and dissimilarities between the bacterial photosynthetic reaction center and Photosystem II at the protein sequence level, hinting at different charge separation and charge transfer mechanisms. The conservation analysis that we perform in this paper can be considered as a model for analyzing the conservation in proteins with a low level of sequence identity

Profile Hidden Markov Models for Analyzing Similarities and Dissimilarities in the Bacterial Reaction Center and Photosystem II

The bacterial photosynthetic reaction center is the evolutionary ancestor of the Photosystem II reaction center. These proteins share the same fold and perform the same biological function. Nevertheless, the details of their molecular reaction mechanism differ. It is of significant biological and biochemical interest to determine which functional characteristics are conserved at the level of the protein sequences. Since the level of sequence identity between the bacterial photosynthetic reaction center and Photosystem II is low, a progressive multiple-sequence alignment leads to errors in identifying the conserved residues. In such a situation, profile hidden Markov models (pHMM) can be used to obtain reliable multiple-sequence alignments. We therefore constructed the pHMM with the help of a sequence alignment based on a structural superposition of both proteins. To validate the multiple-sequence alignments obtained with the pHMM, the conservation of residues with known functional importance was examined. Having confirmed the correctness of the multiple-sequence alignments, we analyzed the conservation of residues involved in hydrogen bonding and redox potential tuning of the cofactors. Our analysis reveals similarities and dissimilarities between the bacterial photosynthetic reaction center and Photosystem II at the protein sequence level, hinting at different charge separation and charge transfer mechanisms. The conservation analysis that we perform in this paper can be considered as a model for analyzing the conservation in proteins with a low level of sequence identity

Frontier Residues Lining Globin Internal Cavities Present Specific Mechanical Properties

The internal cavity matrix of globins plays a key role in their biological function. Previous studies have already highlighted the plasticity of this inner network, which can fluctuate with the proteins breathing motion, and the importance of a few key residues for the regulation of ligand diffusion within the protein. In this Article, we combine all-atom molecular dynamics and coarse-grain Brownian dynamics to establish a complete mechanical landscape for six different globins chain (myoglobin, neuroglobin, cytoglobin, truncated hemoglobin, and chains α and β of hemoglobin). We show that the rigidity profiles of these proteins can fluctuate along time, and how a limited set of residues present specific mechanical properties that are related to their position at the frontier between internal cavities. Eventually, we postulate the existence of conserved positions within the globin fold, which form a mechanical nucleus located at the center of the cavity network, and whose constituent residues are essential for controlling ligand migration in globins

Relating the Diffusion of Small Ligands in Human Neuroglobin to Its Structural and Mechanical Properties

Neuroglobin (Ngb), a recently discovered member of the globin family, is overexpressed in the brain tissues over oxygen deprivation. Unlike more classical globins, such as myoglobin and hemoglobin, it is characterized by a hexacoordinated heme, and its physiological role is still unknown, despite the numerous investigations made on the protein in recent years. Another important specific feature of human Ngb is the presence of two cysteine residues (Cys46 and Cys55), which are known to form an intramolecular disulfide bridge. Since previous work on human Ngb reported that its ligand binding properties could be controlled by the coordination state of the Fe2+ atom (in the heme moiety) and the redox state of the thiol groups, we choose to develop a simulation approach combining coarse-grain Brownian dynamics and all-atom molecular dynamics and metadynamics. We have studied the diffusion of small ligands (CO, NO, and O2) in the globin internal cavity network for various states of human Ngb. Our results show how the structural and mechanical properties of the protein can be related to the ligand migration pathway, which can be extensively modified when changing the thiol’s redox state and the iron’s coordination state. We suggest that ligand binding is favored in the pentacoordinated species bearing an internal disulfide bridge

Role of Ionic Strength and pH in Modulating Thermodynamic Profiles Associated with CO Escape from Rice Nonsymbiotic Hemoglobin 1

Type 1 nonsymbiotic hemoglobins are found in a wide variety of land plants and exhibit very high affinities for exogenous gaseous ligands. These proteins are presumed to have a role in protecting plant cells from oxidative stress under etiolated/hypoxic conditions through NO dioxygenase activity. In this study we have employed photoacoustic calorimetry, time-resolved absorption spectroscopy, and classical molecular dynamics simulations in order to elucidate thermodynamics, kinetics, and ligand migration pathways upon CO photodissociation from WT and a H73L mutant of type 1 nonsymbiotic hemoglobin from <i>Oryza sativa</i> (rice). We observe a temperature dependence of the resolved thermodynamic parameters for CO photodissociation from CO-rHb1 which we attribute to temperature dependent formation of a network of electrostatic interactions in the vicinity of the heme propionate groups. We also observe slower ligand escape from the protein matrix under mildly acidic conditions in both the WT and H73L mutant (τ = 134 ± 19 and 90 ± 15 ns). Visualization of transient hydrophobic channels within our classical molecular dynamics trajectories allows us to attribute this phenomenon to a change in the ligand migration pathway which occurs upon protonation of the distal His73, His117, and His152. Protonation of these residues may be relevant to the functioning of the protein in vivo given that etiolation/hypoxia can cause a decrease in intracellular pH in plant cells

A Time-Resolved Iron-Specific X-ray Absorption Experiment Yields No Evidence for an Fe²⁺ → Fe³⁺ Transition during Q_A^- → Q_B Electron Transfer in the Photosynthetic Reaction Center^†

Previous time-resolved FTIR measurements suggested the involvement of an intermediary component in the electron transfer step QA- → QB in the photosynthetic reaction center (RC) from Rhodobacter sphaeroides [Remy and Gerwert (2003) Nat. Struct. Biol. 10, 637]. By a kinetic X-ray absorption experiment at the Fe K-edge we investigated whether oxidation occurs at the ferric non-heme iron located between the two quinones. In isolated reaction centers with a high content of functional QB, at a time resolution of 30 μs and at room temperature, no evidence for transient oxidation of Fe was obtained. However, small X-ray transients occurred, in a similar micro- to millisecond time range as in the IR experiments, which may point to changes in the Fe ligand environment due to the charges on QA- and QB-. In addition, VIS measurements agree with the IR data and do not exclude an intermediate in the QA- → QB transition