Water Ligands Regulate the Redox Leveling Mechanism of the Oxygen-Evolving Complex of the Photosystem II.

Understanding how water ligands regulate the conformational changes and functionality of the oxygen-evolving complex (OEC) in photosystem II (PSII) throughout the catalytic cycle of oxygen evolution remains a highly intriguing and unresolved challenge. In this study, we investigate the effect of water insertion (WI) on the redox state of the OEC by using the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) hybrid methods. We find that water binding significantly reduces the free energy change for proton-coupled electron transfer (PCET) from Mn to YZ•, underscoring the important regulatory role of water binding, which is essential for enabling the OEC redox-leveling mechanism along the catalytic cycle. We propose a water binding mechanism in which WI is thermodynamically favored by the closed-cubane form of the OEC, with water delivery mediated by Ca2+ ligand exchange. Isomerization from the closed- to open-cubane conformation at three post-WI states highlights the importance of the location of the MnIII center in the OEC and the orientation of its Jahn-Teller axis to conformational changes of the OEC, which might be critical for the formation of the O-O bond. These findings reveal a complex interplay between conformational changes in the OEC and the ligand environment during the activation of the OEC by YZ•. Analogous regulatory effects due to water ligand binding are expected to be important for a wide range of catalysts activated by redox state transitions in aqueous environments.

A Resilient Platform for the Discrete Functionalization of Gold Surfaces Based on N-Heterocyclic Carbene Self-Assembled Monolayers.

We describe the synthesis and characterization of a versatile platform for gold functionalization, based on self-assembled monolayers (SAMs) of distal-pyridine-functionalized N-heterocyclic carbenes (NHC) derived from bis(NHC) Au(I) complexes. The SAMs are characterized using polarization-modulation infrared reflectance-absorption spectroscopy, surface-enhanced Raman spectroscopy, and X-ray photoelectron spectroscopy. The binding mode is examined computationally using density functional theory, including calculations of vibrational spectra and direct comparisons to the experimental spectroscopic signatures of the monolayers. Our joint computational and experimental analyses provide structural information about the SAM binding geometries under ambient conditions. Additionally, we examine the reactivity of the pyridine-functionalized SAMs toward H2SO4 and W(CO)5(THF) and verify the preservation of the introduced functionality at the interface. Our results demonstrate the versatility of N-heterocyclic carbenes as robust platforms for on-surface acid-base and ligand exchange reactions.

A quantitative assessment of (bacterio)chlorophyll assignments in the cryo-EM structure of the Chloracidobacterium thermophilum reaction center.

Chlorophylls and bacteriochlorophylls are the primary pigments used by photosynthetic organisms for light harvesting, energy transfer, and electron transfer. Many molecular structures of (bacterio)chlorophyll-containing protein complexes are available, some of which contain mixtures of different (bacterio)chlorophyll types. Differentiating these, which sometimes are structurally similar, is challenging but is required for leveraging structural data to gain functional insight. The reaction center complex from Chloroacidobacterium thermophilum has a hybrid (bacterio)chlorophyll antenna system containing both chlorophyll a and bacteriochlorophyll a molecules. The recent availability of its cryogenic electron microscopy (cryo-EM) structure provides an opportunity for a quantitative analysis of their identities and chemical environments. Here, we describe a theoretical basis for differentiating chlorophyll a and bacteriochlorophyll a in a cryo-EM map, and apply the approach to the experimental cryo-EM maps of the (bacterio)chlorophyll sites of the chloroacidobacterial reaction center. The comparison reveals that at ~ 2.2-Å resolution, chlorophyll a and bacteriochlorophyll a are easily distinguishable, but the orientation of the bacteriochlorophyll a acetyl moiety is not; however, the latter can confidently be assigned by identifying a hydrogen bond donor from the protein environment. This study reveals the opportunities and challenges in assigning (bacterio)chlorophyll types in structural biology, the accuracy of which is vital for downstream investigations.

Hydroxide Diffusion in Functionalized Cylindrical Nanopores as Idealized Models of Anion Exchange Membrane Environments: An Ab Initio Molecular Dynamics Study.

Anion exchange membranes (AEMs) have attracted significant interest for their applications in fuel cells and other electrochemical devices in recent years. Understanding water distributions and hydroxide transport mechanisms within AEMs is critical to improving their performance as concerns hydroxide conductivity. Recently, nanoconfined environments have been used to mimic AEM environments. Following this approach, we construct nanoconfined cylindrical pore structures using graphane nanotubes (GNs) functionalized with trimethylammonium cations as models of local AEM morphology. These structures were then used to investigate hydroxide transport using ab initio molecular dynamics (AIMD). The simulations showed that hydroxide transport is suppressed in these confined environments relative to the bulk solution although the mechanism is dominated by structural diffusion. One factor causing the suppressed hydroxide transport is the reduced proton transfer (PT) rates due to changes in hydroxide and water solvation patterns under confinement compared to bulk solution as well as strong interactions between hydroxide ions and the tethered cation groups.

Elucidating the Proton Transport Pathways in Liquid Imidazole with First-Principles Molecular Dynamics.

Imidazole is a promising anhydrous proton conductor with a high conductivity comparable to that of water at a similar temperature relative to its melting point. Previous theoretical studies of the mechanism of proton transport in imidazole have relied either on empirical models or on ab initio trajectories that have been too short to draw significant conclusions. Here, we present the results of multiple time-step ab initio molecular dynamics simulations of an excess proton in liquid imidazole reaching 1 ns in total simulation time. We find that the proton transport is dominated by structural diffusion, with the diffusion constant of the proton defect being ∼8 times higher than that of self-diffusion of the imidazole molecules. By using correlation function analysis, we decompose the mechanism for proton transport into a series of first-order processes and show that the proton transport mechanism occurs over three distinct time and length scales. Although the mechanism at intermediate times is dominated by hopping along pseudo-one-dimensional chains, at longer times the overall rate of diffusion is limited by the re-formation of these chains. These results provide a more complete picture of the traditional idealized Grotthuss structural diffusion mechanism.

Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline.

A range of techniques including physical property measurements, neutron scattering experiments, ab initio molecular dynamics, and classical molecular dynamics simulations are used to probe the structural, thermodynamic, and transport properties of a deep eutectic solvent comprised of a 1:2 molar ratio of choline chloride and ethylene glycol. This mixture, known as Ethaline, has many desirable properties for use in a range of applications, and therefore, understanding its liquid structure and transport properties is of interest. Simulation results are able to capture experimental densities, diffusivities, viscosities, and structure factors extremely well. The solvation environment is dynamic and dominated by different hydrogen bonding interactions. Dynamic heterogeneities resulting from hydrogen bonding interactions are quantified. Rotational dynamics of molecular dipole moments of choline and ethylene glycol are computed and found to exhibit a fast and slow mode.

Structure of water confined between two parallel graphene plates.

We study, in this paper, the physical properties of water confined between two parallel graphene plates with different slit widths to understand the effects of confinement on the water structure and how bulk properties are reached as the water layer thickens. It was found that the microscopic structures of the interfacial liquid layer close to graphene vary with the slit width. Water tends to locate at the center of the six-membered ring of graphene planes to form triangular patterns, as found by others. The narrower the slit width is, the more pronounced this pattern is, except for the slit width of 9.5 Å, for which a well-defined two-layer structure of water forms. On the other hand, squared structures can be clearly seen in single snapshots at small (6.5 Å and 7.5 Å) but not large slit widths. Even at small slit widths, the square-like geometry is observed only when an average is taken for a short trajectory, and averaging over a long time yields a triangular pattern dictated by the graphene geometry. We estimate the length of time needed to observe two patterns, respectively. We also used the two-phase thermodynamic model to study the variation of entropy of confined water and found that at 8.5 Å, the entropy of confined water is larger than that of bulk water. The rotational entropy of confined water is higher than that of bulk water for all slit widths due to the reduction of the hydrogen bond in the confined space.

Single Mutations Reshape the Structural Correlation Network of the DMXAA-Human STING Complex.

Subtle changes in protein sequences are able to alter ligand-protein interactions. Unraveling the mechanism of such phenomena is important for understanding ligand-protein interactions, including the DMXAA-STING interaction. DMXAA specifically binds to mouse STING instead of human STING. However, the S162A mutation and a newly discovered E260I mutation endow human STINGAQ with DMXAA sensitivity. Through molecular dynamics simulations, we revealed how these single mutations alter the DMXAA-STING interaction. Compared to mutated systems, structural correlations in the interaction of STINGAQ with DMXAA are stronger, and the correlations are cross-protomers in the dimeric protein. Analyses on correlation coefficients lead to the identification of two key interactions that mediate the strong cross-protomer correlation in the DMXAA-STINGAQ interaction network: DMXAA-267T-162S* and 238R-260E*. These two interactions are partially and totally interrupted by the S162A and E260I mutations, respectively. Moreover, a smaller number of water molecules are displaced upon DMXAA binding to STINGAQ than that on binding to its mutants, leading to a larger entropic penalty for the former. Considering the sensitivity of STINGAQ and two of its mutants to DMXAA, a strong structural correlation appears to discourage DMXAA-STING binding. Such an observation suggests that DMXAA derivatives, which are deprived of hydrogen-bond interaction with both 162S* and 267T, are potential agonists of human STING.