Occupancy Analysis of Water Molecules inside Channels within 25 Å Radius of the Oxygen-Evolving Center of Photosystem II in Molecular Dynamics Simulations.

At room temperature and neutral pH, the oxygen-evolving center (OEC) of photosystem II (PSII) catalyzes water oxidation. During this process, oxygen is released from the OEC, while substrate waters are delivered to the OEC and protons are passed from the OEC to the lumen through water channels known as the narrow or the O4 channel, broad or the Cl1 channel, and large or the O1 channel. Protein residues lining the surfaces of these channels play a critical role in stabilizing the hydrogen-bonding networks that assist in the process. We carried out an occupancy analysis to better understand the structural and possible substrate water dynamics in full PSII monomer molecular dynamics (MD) trajectories in both the S1 and S2 states. We find that the equilibrated positions of water molecules derived from MD-derived electron density maps largely match the experimentally observed positions in crystallography. Furthermore, the occupancy reduction in MD simulations of some water molecules inside the single-filed narrow channel also correlates well with the crystallographic data during a structural transition when the S1 state of the OEC advances to the S2 state. The overall reduced occupancies of water molecules are the source of their "vacancy-hopping" dynamic nature inside these channels, unlike water molecules inside an ice lattice where all water molecules have a fixed unit occupancy. We propose on the basis of findings in our structural and molecular dynamics analysis that the water molecule occupying a pocket formed by D1-D61, D1-S169, and O4 of the OEC could be the last steppingstone to enter into the OEC and that the broad channel may be favored for proton transfer.

Structural Insights into Binding of Remdesivir Triphosphate within the Replication-Transcription Complex of SARS-CoV-2.

Remdesivir is an adenosine analogue that has a cyano substitution in the C1' position of the ribosyl moiety and a modified base structure to stabilize the linkage of the base to the C1' atom with its strong electron-withdrawing cyano group. Within the replication-transcription complex (RTC) of SARS-CoV-2, the RNA-dependent RNA polymerase nsp12 selects remdesivir monophosphate (RMP) over adenosine monophosphate (AMP) for nucleotide incorporation but noticeably slows primer extension after the added RMP of the RNA duplex product is translocated by three base pairs. Cryo-EM structures have been determined for the RTC with RMP at the nucleotide-insertion (i) site or at the i + 1, i + 2, or i + 3 sites after product translocation to provide a structural basis for a delayed-inhibition mechanism by remdesivir. In this study, we applied molecular dynamics (MD) simulations to extend the resolution of structures to the measurable maximum that is intrinsically limited by MD properties of these complexes. Our MD simulations provide (i) a structural basis for nucleotide selectivity of the incoming substrates of remdesivir triphosphate over adenosine triphosphate and of ribonucleotide over deoxyribonucleotide, (ii) new detailed information on hydrogen atoms involved in H-bonding interactions between the enzyme and remdesivir, and (iii) direct information on the catalytically active complex that is not easily captured by experimental methods. Our improved resolution of interatomic interactions at the nucleotide-binding pocket between remedesivir and the polymerase could help to design a new class of anti-SARS-CoV-2 inhibitors.

Glycerol binding at the narrow channel of photosystem II stabilizes the low-spin S2 state of the oxygen-evolving complex.

The oxygen-evolving complex (OEC) of photosystem II (PSII) cycles through redox intermediate states Si (i = 0-4) during the photochemical oxidation of water. The S2 state involves an equilibrium of two isomers including the low-spin S2 (LS-S2) state with its characteristic electron paramagnetic resonance (EPR) multiline signal centered at g = 2.0, and a high-spin S2 (HS-S2) state with its g = 4.1 EPR signal. The relative intensities of the two EPR signals change under experimental conditions that shift the HS-S2/LS-S2 state equilibrium. Here, we analyze the effect of glycerol on the relative stability of the LS-S2 and HS-S2 states when bound at the narrow channel of PSII, as reported in an X-ray crystal structure of cyanobacterial PSII. Our quantum mechanics/molecular mechanics (QM/MM) hybrid models of cyanobacterial PSII show that the glycerol molecule perturbs the hydrogen-bond network in the narrow channel, increasing the pKa of D1-Asp61 and stabilizing the LS-S2 state relative to the HS-S2 state. The reported results are consistent with the absence of the HS-S2 state EPR signal in native cyanobacterial PSII EPR spectra and suggest that the narrow water channel hydrogen-bond network regulates the relative stability of OEC catalytic intermediates during water oxidation.

Insights into Binding of Single-Stranded Viral RNA Template to the Replication-Transcription Complex of SARS-CoV-2 for the Priming Reaction from Molecular Dynamics Simulations.

A minimal replication-transcription complex (RTC) of SARS-CoV-2 for synthesis of viral RNAs includes the nsp12 RNA-dependent RNA polymerase and two nsp8 RNA primase subunits for de novo primer synthesis, one nsp8 in complex with its accessory nsp7 subunit and the other without it. The RTC is responsible for faithfully copying the entire (+) sense viral genome from its first 5'-end to the last 3'-end nucleotides through a replication-intermediate (RI) template. The single-stranded (ss) RNA template for the RI is its 33-nucleotide 3'-poly(A) tail adjacent to a well-characterized secondary structure. The ssRNA template for viral transcription is a 5'-UUUAU-3' next to stem-loop (SL) 1'. We analyze the electrostatic potential distribution of the nsp8 subunit within the RTC around the template strand of the primer/template (P/T) RNA duplex in recently published cryo-EM structures to address the priming reaction using the viral poly(A) template. We carried out molecular dynamics (MD) simulations with a P/T RNA duplex, the viral poly(A) template, or a generic ssRNA template. We find evidence that the viral poly(A) template binds similarly to the template strand of the P/T RNA duplex within the RTC, mainly through electrostatic interactions, providing new insights into the priming reaction by the nsp8 subunit within the RTC, which differs significantly from the existing proposal of the nsp7/nsp8 oligomer formed outside the RTC. High-order oligomerization of nsp8 and nsp7 for SARS-CoV observed outside the RTC of SARS-CoV-2 is not found in the RTC and not likely to be relevant to the priming reaction.

Binding of the substrate analog methanol in the oxygen-evolving complex of photosystem II in the D1-N87A genetic variant of cyanobacteria.

The solar water-splitting protein complex, photosystem II (PSII), catalyzes one of the most energetically demanding reactions in nature by using light energy to drive a catalyst capable of oxidizing water. The water oxidation reaction is catalyzed at the Mn4Ca-oxo cluster in the oxygen-evolving complex (OEC), which cycles through five light-driven S-state intermediates (S0-S4). A detailed mechanism of the reaction remains elusive as it requires knowledge of the delivery and binding of substrate water in the higher S-state intermediates. In this study, we use two-dimensional (2D) hyperfine sublevel correlation spectroscopy, in conjunction with quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT), to probe the binding of the substrate analog, methanol, in the S2 state of the D1-N87A variant of PSII from Synechocystis sp. PCC 6803. The results indicate that the size and specificity of the "narrow" channel is altered in D1-N87A PSII, allowing for the binding of deprotonated 13C-labeled methanol at the Mn4(IV) ion of the catalytic cluster in the S2 state. This has important implications on the mechanistic models for water oxidation in PSII.

High-resolution cryo-electron microscopy structure of photosystem II from the mesophilic cyanobacterium, Synechocystis sp. PCC 6803.

Photosystem II (PSII) enables global-scale, light-driven water oxidation. Genetic manipulation of PSII from the mesophilic cyanobacterium Synechocystis sp. PCC 6803 has provided insights into the mechanism of water oxidation; however, the lack of a high-resolution structure of oxygen-evolving PSII from this organism has limited the interpretation of biophysical data to models based on structures of thermophilic cyanobacterial PSII. Here, we report the cryo-electron microscopy structure of PSII from Synechocystis sp. PCC 6803 at 1.93-Å resolution. A number of differences are observed relative to thermophilic PSII structures, including the following: the extrinsic subunit PsbQ is maintained, the C terminus of the D1 subunit is flexible, some waters near the active site are partially occupied, and differences in the PsbV subunit block the Large (O1) water channel. These features strongly influence the structural picture of PSII, especially as it pertains to the mechanism of water oxidation.

Heterogeneous Composition of Oxygen-Evolving Complexes in Crystal Structures of Dark-Adapted Photosystem II.

Photosystem II (PSII) is a homodimeric protein complex that catalyzes water oxidation at the oxygen-evolving complex (OEC), a heterocubanoid calcium-tetramanganese cluster. Here, we analyze the omit electron density peaks of the OEC's metal ions in five X-ray free-electron laser PSII structures at resolutions between 2.15 and 1.95 Å. The omit peaks can be described by the total number of electrons and approximated by the variance of electron density distribution when the distributions are spherically symmetric. We show that the number of electrons of metal centers is different in the two OECs of PSII dimers, implying that the oxidation states and/or occupancies of individual metal ions are different in the two monomers. In either case, we find that the two OECs of dark-adapted PSII dimers in crystals are not fully synchronized in the S1 state. Differences in redox states of the OEC in PSII only partially account for the observation that the electron densities integrate to a smaller number of electrons than expected. Differences between the determined and expected relative electron numbers are much larger than the estimated errors, indicating heterogeneity in the OEC composition.

Computational insights into the membrane fusion mechanism of SARS-CoV-2 at the cellular level.

The membrane fusion mechanism of SARS-CoV-2 offers an attractive target for the development of small molecule antiviral inhibitors. Fusion involves an initial binding of the crown-like trimeric spike glycoproteins of SARS-CoV-2 to the receptor angiotensin II-converting enzyme 2 (ACE2) on the permissive host cellular membrane and a prefusion to post-fusion conversion of the spike trimer. During this conversion, the fusion peptides of the spike trimer are inserted into the host membrane to bring together the host and viral membranes for membrane fusion in highly choreographic events. However, geometric constraints due to interactions with the membranes remain poorly understood. In this study, we build structural models of super-complexes of spike trimer/ACE2 dimers based on the molecular structures of the ACE2/neutral amino acid transporter B(0)AT heterodimer. We determine the conformational constraints due to the membrane geometry on the enzymatic activity of ACE2 and on the viral fusion process. Furthermore, we find that binding three ACE2 dimers per spike trimer is essential to open the central pore as necessary for triggering productive membrane fusion through an elongation of the central stalk. The reported findings thus provide valuable insights for targeting the membrane fusion mechanism for drug design at the molecular level.

A structurally preserved allosteric site in the MIF superfamily affects enzymatic activity and CD74 activation in D-dopachrome tautomerase.

The macrophage migration inhibitory factor (MIF) family of cytokines contains multiple ligand-binding sites and mediates immunomodulatory processes through an undefined mechanism(s). Previously, we reported a dynamic relay connecting the MIF catalytic site to an allosteric site at its solvent channel. Despite structural and functional similarity, the MIF homolog D-dopachrome tautomerase (also called MIF-2) has low sequence identity (35%), prompting the question of whether this dynamic regulatory network is conserved. Here, we establish the structural basis of an allosteric site in MIF-2, showing with solution NMR that dynamic communication is preserved in MIF-2 despite differences in the primary sequence. X-ray crystallography and NMR detail the structural consequences of perturbing residues in this pathway, which include conformational changes surrounding the allosteric site, despite global preservation of the MIF-2 fold. Molecular simulations reveal MIF-2 to contain a comparable hydrogen bond network to that of MIF, which was previously hypothesized to influence catalytic activity by modulating the strength of allosteric coupling. Disruption of the allosteric relay by mutagenesis also attenuates MIF-2 enzymatic activity in vitro and the activation of the cluster of differentiation 74 receptor in vivo, highlighting a conserved point of control for nonoverlapping functions in the MIF superfamily.

Mechanism of Inhibition of the Reproduction of SARS-CoV-2 and Ebola Viruses by Remdesivir.

Remdesivir is an antiviral drug initially designed against the Ebola virus. The results obtained with it both in biochemical studies in vitro and in cell line assays in vivo were very promising, but it proved to be ineffective in clinical trials. Remdesivir exhibited far better efficacy when repurposed against SARS-CoV-2. The chemistry that accounts for this difference is the subject of this study. Here, we examine the hypothesis that remdesivir monophosphate (RMP)-containing RNA functions as a template at the polymerase site for the second run of RNA synthesis, and as mRNA at the decoding center for protein synthesis. Our hypothesis is supported by the observation that RMP can be incorporated into RNA by the RNA-dependent RNA polymerases (RdRps) of both viruses, although some of the incorporated RMPs are subsequently removed by exoribonucleases. Furthermore, our hypothesis is consistent with the fact that RdRp of SARS-CoV-2 selects RMP for incorporation over AMP by 3-fold in vitro, and that RMP-added RNA can be rapidly extended, even though primer extension is often paused when the added RMP is translocated at the i + 3 position (with i the nascent base pair at an initial insertion site of RMP) or when the concentrations of the subsequent nucleoside triphosphates (NTPs) are below their physiological concentrations. These observations have led to the hypothesis that remdesivir might be a delayed chain terminator. However, that hypothesis is challenged under physiological concentrations of NTPs by the observation that approximately three-quarters of RNA products efficiently overrun the pause.

Proton exit pathways surrounding the oxygen evolving complex of photosystem II.

Photosystem II allows water to be the primary electron source for the photosynthetic electron transfer chain. Water is oxidized to dioxygen at the Oxygen Evolving Complex (OEC), a Mn4CaO5 inorganic core embedded on the lumenal side of PSII. Water-filled channels surrounding the OEC must bring in substrate water molecules, remove the product protons to the lumen, and may transport the product oxygen. Three water-filled channels, denoted large, narrow, and broad, extend from the OEC towards the aqueous surface more than 15 Å away. However, the role of each pathway in the transport in and out of the OEC is yet to be established. Here, we combine Molecular Dynamics (MD), Multi Conformation Continuum Electrostatics (MCCE) and Network Analysis to compare and contrast the three potential proton transfer paths. Hydrogen bond network analysis shows that near the OEC the waters are highly interconnected with similar free energy for hydronium at all locations. The paths diverge as they move towards the lumen. The water chain in the broad channel is better connected than in the narrow and large channels, where disruptions in the network are observed approximately 10 Å from the OEC. In addition, the barrier for hydronium translocation is lower in the broad channel. Thus, a proton released from any location on the OEC can access all paths, but the likely exit to the lumen passes through PsbO via the broad channel.

Quantitative assessment of chlorophyll types in cryo-EM maps of photosystem I acclimated to far-red light.

Chlorophyll cofactors are vital for the metabolism of photosynthetic organisms. Cryo-electron microscopy (cryo-EM) has been used to elucidate molecular structures of pigment-protein complexes, but the minor structural differences between multiple types of chlorophylls make them difficult to distinguish in cryo-EM maps. This is exemplified by inconsistencies in the assignments of chlorophyll f molecules in structures of photosystem I acclimated to far-red light (FRL-PSI). A quantitative assessment of chlorophyll substituents in cryo-EM maps was used to identify chlorophyll f-binding sites in structures of FRL-PSI from two cyanobacteria. The two cryo-EM maps provide direct evidence for chlorophyll f-binding at two and three binding sites, respectively, and three more sites in each structure exhibit strong indirect evidence for chlorophyll f-binding. Common themes in chlorophyll f-binding are described that clarify the current understanding of the molecular basis for FRL photoacclimation in photosystems.

Corrigendum to Quantitative assessment of chlorophyll types in cryo-EM maps of photosystem I acclimated to far-red light BBA Advances 1 (2021) 100019.

[This corrects the article DOI: 10.1016/j.bbadva.2021.100019.].

D1-S169A substitution of photosystem II reveals a novel S2-state structure.

In photosystem II (PSII), photosynthetic water oxidation occurs at the O2-evolving complex (OEC), a tetramanganese-calcium cluster that cycles through light-induced redox intermediates (S0-S4) to produce oxygen from two substrate water molecules. The OEC is surrounded by a hydrogen-bonded network of amino-acid residues that plays a crucial role in proton transfer and substrate water delivery. Previously, we found that D1-S169 was crucial for water oxidation and its mutation to alanine perturbed the hydrogen-bonding network. In this study, we demonstrate that the activation energy for the S2 to S1 transition of D1-S169A PSII is higher than wild-type PSII with a ~1.7-2.7× slower rate of charge recombination with QA- relative to wild-type PSII. Arrhenius analysis of the decay kinetics shows an Ea of 5.87 ± 1.15 kcal mol-1 for decay back to the S1 state, compared to 0.80 ± 0.13 kcal mol-1 for the wild-type S2 state. In addition, we find that ammonia does not affect the S2-state EPR signal, indicating that ammonia does not bind to the Mn cluster in D1-S169A PSII. Finally, a QM/MM analysis indicates that an additional water molecule binds to the Mn4 ion in place of an oxo ligand O5 in the S2 state of D1-S169A PSII. The altered S2 state of D1-S169A PSII provides insight into the S2➔S3 state transition.

Electric field stimulates production of highly conductive microbial OmcZ nanowires.

Multifunctional living materials are attractive due to their powerful ability to self-repair and replicate. However, most natural materials lack electronic functionality. Here we show that an electric field, applied to electricity-producing Geobacter sulfurreducens biofilms, stimulates production of cytochrome OmcZ nanowires with 1,000-fold higher conductivity (30 S cm-1) and threefold higher stiffness (1.5 GPa) than the cytochrome OmcS nanowires that are important in natural environments. Using chemical imaging-based multimodal nanospectroscopy, we correlate protein structure with function and observe pH-induced conformational switching to ß-sheets in individual nanowires, which increases their stiffness and conductivity by 100-fold due to enhanced π-stacking of heme groups; this was further confirmed by computational modeling and bulk spectroscopic studies. These nanowires can transduce mechanical and chemical stimuli into electrical signals to perform sensing, synthesis and energy production. These findings of biologically produced, highly conductive protein nanowires may help to guide the development of seamless, bidirectional interfaces between biological and electronic systems.

Identification of a Na+-Binding Site near the Oxygen-Evolving Complex of Spinach Photosystem II.

The oxygen-evolving complex (OEC) of photosystem II (PSII) is an oxomanganese cluster composed of four redox-active Mn ions and one redox-inactive Ca2+ ion, with two nearby bound Cl- ions. Sodium is a common counterion of both chloride and hydroxide anions, and a sodium-specific binding site has not been identified near the OEC. Here, we find that the oxygen-evolution activity of spinach PSII increases with Na+ concentration, particularly at high pH. A Na+-specific binding site next to the OEC, becomes available after deprotonation of the D1-H337 amino acid residue, is suggested by the analysis of two recently published PSII cryo-electron microscopy maps in combination with quantum mechanical calculations and multiconformation continuum electrostatics simulations.

Regulation of MIF Enzymatic Activity by an Allosteric Site at the Central Solvent Channel.

In proteins with multiple functions, such as macrophage migration inhibitory factor (MIF), the study of its intramolecular dynamic network can offer a unique opportunity to understand how a single protein is able to carry out several nonoverlapping functions. A dynamic mechanism that controls the MIF-induced activation of CD74 was recently discovered. In this study, the regulation of tautomerase activity was explored. The catalytic base Pro1 is found to form dynamic communications with the same allosteric node that regulates CD74 activation. Signal transmission between the allosteric and catalytic sites take place through intramolecular aromatic interactions and a hydrogen bond network that involves residues and water molecules of the MIF solvent channel. Once thought to be a consequence of trimerization, a regulatory function for the solvent channel is now defined. These results provide mechanistic insights into the regulation of catalytic activity and the role of solvent channel water molecules in MIF catalysis.

Allosteric Impact of the Variable Insert Loop in Vaccinia H1-Related (VHR) Phosphatase.

Dynamics and conformational motions are important to the activity of enzymes, including protein tyrosine phosphatases. These motions often extend to regions outside the active site, called allosteric regions. In the tyrosine phosphatase Vaccinia H1-related (VHR) enzyme, we demonstrate the importance of the allosteric interaction between the variable insert region and the active-site loops in VHR. These studies include solution nuclear magnetic resonance, computation, steady-state, and rapid kinetic measurements. Overall, the data indicate concerted millisecond motions exist between the variable insert and the catalytic acid loop in wild-type (WT) VHR. The 150 ns computation studies show a flexible acid loop in WT VHR that opens during the simulation from its initial closed structure. Mutation of the variable insert residue, asparagine 74, to alanine results in a rigidification of the acid loop as observed by molecular dynamics simulations and a disruption of crucial active-site hydrogen bonds. Moreover, enzyme kinetic analysis shows a weakening of substrate affinity in the N74A mutant and a >2-fold decrease in substrate cleavage and hydrolysis rates. These data show that despite being nearly 20 Å from the active site, the variable insert region is linked to the acid loop by coupled millisecond motions, and that disruption of the communication between the variable insert and active site alters the normal catalytic function of VHR and perturbs the active-site environment.

Copper-mediated thiol potentiation and mutagenesis-guided modeling suggest a highly conserved copper-binding motif in human OR2M3.

Sulfur-containing compounds within a physiological relevant, natural odor space, such as the key food odorants, typically constitute the group of volatiles with the lowest odor thresholds. The observation that certain metals, such as copper, potentiate the smell of sulfur-containing, metal-coordinating odorants led to the hypothesis that their cognate receptors are metalloproteins. However, experimental evidence is sparse-so far, only one human odorant receptor, OR2T11, and a few mouse receptors, have been reported to be activated by sulfur-containing odorants in a copper-dependent way, while the activation of other receptors by sulfur-containing odorants did not depend on the presence of metals. Here we identified an evolutionary conserved putative copper interaction motif CC/CSSH, comprising two copper-binding sites in TMH5 and TMH6, together with the binding pocket for 3-mercapto-2-methylpentan-1-ol in the narrowly tuned human receptor OR2M3. To characterize the copper-binding motif, we combined homology modeling, docking studies, site-directed mutagenesis, and functional expression of recombinant ORs in a cell-based, real-time luminescence assay. Ligand activation of OR2M3 was potentiated in the presence of copper. This effect of copper was mimicked by ionic and colloidal silver. In two broadly tuned receptors, OR1A1 and OR2W1, which did not reveal a putative copper interaction motif, activation by their most potent, sulfur-containing key food odorants did not depend on the presence of copper. Our results suggest a highly conserved putative copper-binding motif to be necessary for a copper-modulated and thiol-specific function of members from three subfamilies of family 2 ORs.