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.

Mutation-induced shift of the photosystem II active site reveals insight into conserved water channels.

Photosystem II (PSII) is the water-plastoquinone photo-oxidoreductase central to oxygenic photosynthesis. PSII has been extensively studied for its ability to catalyze light-driven water oxidation at a Mn4CaO5 cluster called the oxygen-evolving complex (OEC). Despite these efforts, the complete reaction mechanism for water oxidation by PSII is still heavily debated. Previous mutagenesis studies have investigated the roles of conserved amino acids, but these studies have lacked a direct structural basis that would allow for a more meaningful interpretation. Here, we report a 2.14-Å resolution cryo-EM structure of a PSII complex containing the substitution Asp170Glu on the D1 subunit. This mutation directly perturbs a bridging carboxylate ligand of the OEC, which alters the spectroscopic properties of the OEC without fully abolishing water oxidation. The structure reveals that the mutation shifts the position of the OEC within the active site without markedly distorting the Mn4CaO5 cluster metal-metal geometry, instead shifting the OEC as a rigid body. This shift disturbs the hydrogen-bonding network of structured waters near the OEC, causing disorder in the conserved water channels. This mutation-induced disorder appears consistent with previous FTIR spectroscopic data. We further show using quantum mechanics/molecular mechanics methods that the mutation-induced structural changes can affect the magnetic properties of the OEC by altering the axes of the Jahn-Teller distortion of the Mn(III) ion coordinated to D1-170. These results offer new perspectives on the conserved water channels, the rigid body property of the OEC, and the role of D1-Asp170 in the enzymatic water oxidation mechanism.

Mechanism of Action of KL-50, a Candidate Imidazotetrazine for the Treatment of Drug-Resistant Brain Cancers.

Aberrant DNA repair is a hallmark of cancer, and many tumors display reduced DNA repair capacities that sensitize them to genotoxins. Here, we demonstrate that the differential DNA repair capacities of healthy and transformed tissue may be exploited to obtain highly selective chemotherapies. We show that the novel N3-(2-fluoroethyl)imidazotetrazine "KL-50" is a selective toxin toward tumors that lack the DNA repair protein O6-methylguanine-DNA-methyltransferase (MGMT), which reverses the formation of O6-alkylguanine lesions. We establish that KL-50 generates DNA interstrand cross-links (ICLs) by a multistep process comprising DNA alkylation to generate an O6-(2-fluoroethyl)guanine (O6FEtG) lesion, slow unimolecular displacement of fluoride to form an N1,O6-ethanoguanine (N1,O6EtG) intermediate, and ring-opening by the adjacent cytidine. The slow rate of N1,O6EtG formation allows healthy cells expressing MGMT to reverse the initial O6FEtG lesion before it evolves to N1,O6EtG, thereby suppressing the formation of toxic DNA-MGMT cross-links and reducing the amount of DNA ICLs generated in healthy cells. In contrast, O6-(2-chloroethyl)guanine lesions produced by agents such as lomustine and the N3-(2-chloroethyl)imidazotetrazine mitozolomide rapidly evolve to N1,O6EtG, resulting in the formation of DNA-MGMT cross-links and DNA ICLs in healthy tissue. These studies suggest that careful consideration of the rates of chemical DNA modification and biochemical DNA repair may lead to the identification of other tumor-specific genotoxic agents.

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.

Drug Binding to BamA Targets Its Lateral Gate.

BamA, the core component of the ß-barrel assembly machinery (BAM) complex, is an outer-membrane protein (OMP) in Gram-negative bacteria. Its function is to insert and fold substrate OMPs into the outer membrane (OM). Evidence suggests that BamA follows the asymmetric hybrid-barrel model where the first and last strands of BamA separate, a process known as lateral gate opening, to allow nascent substrate OMP ß-strands to sequentially insert and fold through ß-augmentation. Recently, multiple lead compounds that interfere with BamA's function have been identified. We modeled and then docked one of these compounds into either the extracellular loops of BamA or the open lateral gate. With the compound docked in the loops, we found that the lateral gate remains closed during 5 µs molecular dynamics simulations. The same compound when docked in the open lateral gate stays bound to the ß16 strand of BamA during the simulation, which would prevent substrate OMP folding. In addition, we simulated mutants of BamA that are resistant to one or more of the identified lead compounds. In these simulations, we observed a differing degree and/or frequency of opening of BamA's lateral gate compared to BamA-apo, suggesting that the mutations grant resistance by altering the dynamics at the gate. We conclude that the compounds act by inhibiting BamA lateral gate opening and/or binding of substrate, thus preventing subsequent OMP folding and insertion.

Insights into substrate transport and water permeation in the mycobacterial transporter MmpL3.

Mycobacteria, such as Mycobacterium tuberculosis, are characterized by a uniquely thick and waxy cell envelope that consists of two membranes, with a variety of mycolates comprising their outer membrane (OM). The protein Mycobacterial membrane protein Large 3 (MmpL3) is responsible for the transport of a primary OM component, trehalose monomycolate (TMM), from the inner (cytoplasmic) membrane (IM) to the periplasmic space, a process driven by the proton gradient. Although multiple structures of MmpL3 with bound substrates have been solved, the exact pathway(s) for TMM or proton transport remains elusive. Here, employing molecular dynamics simulations we investigate putative pathways for either transport species. We hypothesized that MmpL3 will cycle through similar conformational states as the related transporter AcrB, which we used as targets for modeling the conformation of MmpL3. A continuous water pathway through the transmembrane region was found in one of these states, illustrating a putative pathway for protons. Additional equilibrium simulations revealed that TMM can diffuse from the membrane into a binding pocket in MmpL3 spontaneously. We also found that acetylation of TMM, which is required for transport, makes it more stable within MmpL3's periplasmic cavity compared with the unacetylated form.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers.

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the QB site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to QB. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with QB•- are compared with earlier results obtained with QB [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

How to correct relative voxel scale factors for calculations of vector-difference Fourier maps in cryo-EM.

The atomic coordinates derived from cryo-electron microscopy (cryo-EM) maps can be inaccurate when the voxel scaling factors are not properly calibrated. Here, we describe a method for correcting relative voxel scaling factors between pairs of cryo-EM maps for the same or similar structures that are expanded or contracted relative to each other. We find that the correction of scaling factors reduces the amplitude differences of Fourier-inverted structure factors from voxel-rescaled maps by up to 20-30%, as shown by two cryo-EM maps of the SARS-CoV-2 spike protein measured at pH 4.0 and pH 8.0. This allows for the calculation of the difference map after properly scaling, revealing differences between the two structures for individual amino acid residues. Unexpectedly, the analysis uncovers two previously overlooked differences of amino acid residues in structures and their local structural changes. Furthermore, we demonstrate the method as applied to two cryo-EM maps of monomeric apo-photosystem II from the cyanobacteria Synechocystis sp. PCC 6803 and Thermosynechococcus elongatus. The resulting difference maps reveal many changes in the peripheral transmembrane PsbX subunit between the two species.

Proton transfer activity of the reconstituted Mycobacterium tuberculosis MmpL3 is modulated by substrate mimics and inhibitors.

Transporters belonging to the Resistance-Nodulation-cell Division (RND) superfamily of proteins such as Mycobacterium tuberculosis MmpL3 and its analogs are the focus of intense investigations due to their importance in the physiology of Corynebacterium-Mycobacterium-Nocardia species and antimycobacterial drug discovery. These transporters deliver trehalose monomycolates, the precursors of major lipids of the outer membrane, to the periplasm by a proton motive force-dependent mechanism. In this study, we successfully purified, from native membranes, the full-length and the C-terminal truncated M. tuberculosis MmpL3 and Corynebacterium glutamicum CmpL1 proteins and reconstituted them into proteoliposomes. We also generated a series of substrate mimics and inhibitors specific to these transporters, analyzed their activities in the reconstituted proteoliposomes, and carried out molecular dynamics simulations of the model MmpL3 transporter at different pH. We found that all reconstituted proteins facilitate proton translocation across a phospholipid bilayer, but MmpL3 and CmpL1 differ dramatically in their responses to pH and interactions with substrate mimics and indole-2-carboxamide inhibitors. Our results further suggest that some inhibitors abolish the transport activity of MmpL3 and CmpL1 by inhibition of proton translocation.

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.

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.

Identification and mechanistic basis of non-ACE2 blocking neutralizing antibodies from COVID-19 patients with deep RNA sequencing and molecular dynamics simulations.

Variants of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continue to cause disease and impair the effectiveness of treatments. The therapeutic potential of convergent neutralizing antibodies (NAbs) from fully recovered patients has been explored in several early stages of novel drugs. Here, we identified initially elicited NAbs (Ig Heavy, Ig lambda, Ig kappa) in response to COVID-19 infection in patients admitted to the intensive care unit at a single center with deep RNA sequencing (>100 million reads) of peripheral blood as a diagnostic tool for predicting the severity of the disease and as a means to pinpoint specific compensatory NAb treatments. Clinical data were prospectively collected at multiple time points during ICU admission, and amino acid sequences for the NAb CDR3 segments were identified. Patients who survived severe COVID-19 had significantly more of a Class 3 antibody (C135) to SARS-CoV-2 compared to non-survivors (15059.4 vs. 1412.7, p = 0.016). In addition to highlighting the utility of RNA sequencing in revealing unique NAb profiles in COVID-19 patients with different outcomes, we provided a physical basis for our findings via atomistic modeling combined with molecular dynamics simulations. We established the interactions of the Class 3 NAb C135 with the SARS-CoV-2 spike protein, proposing a mechanistic basis for inhibition via multiple conformations that can effectively prevent ACE2 from binding to the spike protein, despite C135 not directly blocking the ACE2 binding motif. Overall, we demonstrate that deep RNA sequencing combined with structural modeling offers the new potential to identify and understand novel therapeutic(s) NAbs in individuals lacking certain immune responses due to their poor endogenous production. Our results suggest a possible window of opportunity for administration of such NAbs when their full sequence becomes available. A method involving rapid deep RNA sequencing of patients infected with SARS-CoV-2 or its variants at the earliest infection time could help to develop personalized treatments using the identified specific NAbs.

Preparing Membrane Proteins for Simulation Using CHARMM-GUI.

Molecular dynamics simulations of membrane proteins have grown dramatically in the last 20 years. Running these simulations first requires embedding the protein's three-dimensional structure in a lipid bilayer of a suitable composition, one that resembles its native environment. This step is far from trivial, especially for modeling heterogeneous mixtures of lipids. CHARMM-GUI, a webserver for simulation system preparation greatly simplifies this step, allowing for the construction of complex heterogeneous and/or asymmetric membranes. Here, we demonstrate how to use CHARMM-GUI to build the membrane for the outer-membrane protein BamA.

Membrane thinning and lateral gating are consistent features of BamA across multiple species.

In Gram-negative bacteria, the folding and insertion of ß-barrel outer membrane proteins (OMPs) to the outer membrane are mediated by the ß-barrel assembly machinery (BAM) complex. Two leading models of this process have been put forth: the hybrid barrel model, which claims that a lateral gate in BamA's ß-barrel can serve as a template for incoming OMPs, and the passive model, which claims that a thinned membrane near the lateral gate of BamA accelerates spontaneous OMP insertion. To examine the key elements of these two models, we have carried out 45.5 µs of equilibrium molecular dynamics simulations of BamA with and without POTRA domains from Escherichia coli, Salmonella enterica, Haemophilus ducreyi and Neisseria gonorrhoeae, together with BamA's homolog, TamA from E. coli, in their native, species-specific outer membranes. In these equilibrium simulations, we consistently observe membrane thinning near the lateral gate for all proteins. We also see occasional spontaneous lateral gate opening and sliding of the ß-strands at the gate interface for N. gonorrhoeae, indicating that the gate is dynamic. An additional 14 µs of free-energy calculations shows that the energy necessary to open the lateral gate in BamA/TamA varies by species, but is always lower than the Omp85 homolog, FhaC. Our combined results suggest OMP insertion utilizes aspects of both the hybrid barrel and passive models.