Direct evidence for a deprotonated lysine serving as a H-bond "acceptor" in a photoreceptor protein.

Deprotonation or suppression of the pKa of the amino group of a lysine sidechain is a widely recognized phenomenon whereby the sidechain amino group transiently can act as a nucleophile at the active site of enzymatic reactions. However, a deprotonated lysine and its molecular interactions have not been directly experimentally detected. Here, we demonstrate a deprotonated lysine stably serving as an "acceptor" in a H-bond between the photosensor protein RcaE and its chromophore. Signal splitting and trans-H-bond J coupling observed by NMR spectroscopy provide direct evidence that Lys261 is deprotonated and serves as a H-bond acceptor for the chromophore NH group. Quantum mechanical/molecular mechanical calculations also indicate that this H-bond exists stably. Interestingly, the sidechain amino group of the lysine can act as both donor and acceptor. The remarkable shift in the H-bond characteristics arises from a decrease in solvation, triggered by photoisomerization. Our results provide insights into the dual role of this lysine. This mechanism has broad implications for other biological reactions in which lysine plays a role.

Substrate specificity controlled by the exit site of human P4-ATPases, revealed by de novo point mutations in neurological disorders.

The maintenance of lipid asymmetry on the plasma membrane is regulated by flippases, such as ATP8A2, ATP11A, and ATP11C, which translocate phosphatidylserine and phosphatidylethanolamine from the outer leaflet to the inner leaflet. We previously identified a patient-derived point mutation (Q84E) in ATP11A at the phospholipid entry site, which acquired the ability to flip phosphatidylcholine (PtdCho). This mutation led to elevated levels of sphingomyelin (SM) in the outer leaflet of the plasma membrane. We herein present two de novo ATP11A dominant mutations (E114G and S399L) in heterozygous patients exhibiting neurological and developmental disorders. These mutations, situated near the predicted phospholipid exit site, similarly confer the ability for ATP11A to recognize PtdCho as a substrate, resulting in its internalization into cells. Cells expressing these mutants had increased SM levels on their surface, attributed to the up-regulated expression of the sphingomyelin synthase-1 gene, rendering them more susceptible to SM phosphodiesterase-mediated cell lysis. Corresponding mutations in ATP11C and ATP8A2, paralogs of ATP11A, exerted similar effects on PtdCho-flipping activity and increased SM levels on the cell surface. Molecular dynamics simulations, based on the ATP11C structure, suggest that the E114G and S399L mutations enhance ATP11C's affinity toward PtdCho. These findings underscore the importance of the well-conserved exit and entry sites in determining phospholipid substrate specificity and indicate that aberrant flipping of PtdCho contributes to neurological disorders.

ZBTB2 links p53 deficiency to HIF-1-mediated hypoxia signaling to promote cancer aggressiveness.

Aberrant activation of the hypoxia-inducible transcription factor HIF-1 and dysfunction of the tumor suppressor p53 have been reported to induce malignant phenotypes and therapy resistance of cancers. However, their mechanistic and functional relationship remains largely unknown. Here, we reveal a mechanism by which p53 deficiency triggers the activation of HIF-1-dependent hypoxia signaling and identify zinc finger and BTB domain-containing protein 2 (ZBTB2) as an important mediator. ZBTB2 forms homodimers via its N-terminus region and increases the transactivation activity of HIF-1 only when functional p53 is absent. The ZBTB2 homodimer facilitates invasion, distant metastasis, and growth of p53-deficient, but not p53-proficient, cancers. The intratumoral expression levels of ZBTB2 are associated with poor prognosis in lung cancer patients. ZBTB2 N-terminus-mimetic polypeptides competitively inhibit ZBTB2 homodimerization and significantly suppress the ZBTB2-HIF-1 axis, leading to antitumor effects. Our data reveal an important link between aberrant activation of hypoxia signaling and loss of a tumor suppressor and provide a rationale for targeting a key mediator, ZBTB2, to suppress cancer aggressiveness.

Molecular origins of absorption wavelength variation among phycocyanobilin-binding proteins.

Phycocyanobilin (PCB)-binding proteins, including cyanobacteriochromes and phytochromes, function as photoreceptors and exhibit a wide range of absorption maximum wavelengths. To elucidate the color-tuning mechanisms among these proteins, we investigated seven crystal structures of six PCB-binding proteins: Anacy_2551g3, AnPixJg2, phosphorylation-responsive photosensitive histidine kinase, RcaE, Sb.phyB(PG)-PCB, and Slr1393g3. Employing a quantum chemical/molecular mechanical approach combined with a polarizable continuum model, our analysis revealed that differences in absorption wavelengths among PCB-binding proteins primarily arise from variations in the shape of the PCB molecule itself, accounting for a ∼150 nm difference. Remarkably, calculated excitation energies sufficiently reproduced the absorption wavelengths of these proteins spanning ∼200 nm, including 728 nm for Anacy_2551g3. However, assuming the hypothesized lactim conformation resulted in a significant deviation from the experimentally measured absorption wavelength for Anacy_2551g3. The significantly red-shifted absorption wavelength of Anacy_2551g3 can unambiguously be explained by the significant overlap of molecular orbitals between the two pyrrole rings at both edges of the PCB chromophore without the need to hypothesize lactim formation.

Modulation of Electron Transfer Branches by Atrazine and Triazine Herbicides in Photosynthetic Reaction Centers.

Quinone analogue molecules, functioning as herbicides, bind to the secondary quinone site, QB, in type-II photosynthetic reaction centers, including those from purple bacteria (PbRC). Here, we investigated the impact of herbicide binding on electron transfer branches, using herbicide-bound PbRC crystal structures and employing the linear Poisson-Boltzmann equation. In contrast to urea and phenolic herbicides [Fufezan, C. Biochemistry 2005, 44, 12780-12789], binding of atrazine and triazine did not cause significant changes in the redox-potential (Em) values of the primary quinone (QA) in these crystal structures. However, a slight Em difference at the bacteriopheophytin in the electron transfer inactive branch (HM) was observed between the S(-)- and R(+)-triazine-bound PbRC structures. This discrepancy is linked to variations in the protonation pattern of the tightly coupled Glu-L212 and Glu-H177 pairs, crucial components of the proton uptake pathway in native PbRC. These findings suggest the existence of a QB-mediated link between the electron transfer inactive HM and the proton uptake pathway in PbRCs.

Energetics of the H-Bond Network in Exiguobacterium sibiricum Rhodopsin.

Exiguobacterium sibiricum rhodopsin (ESR) functions as a light-driven proton pump utilizing Lys96 for proton uptake and maintaining its activity over a wide pH range. Using a combination of methodologies including the linear Poisson-Boltzmann equation and a quantum mechanical/molecular mechanical approach with a polarizable continuum model, we explore the microscopic mechanisms underlying its pumping activity. Lys96, the primary proton uptake site, remains deprotonated owing to the loss of solvation in the ESR protein environment. Asp85, serving as a proton acceptor group for Lys96, does not form a low-barrier H-bond with His57. Instead, deprotonated Asp85 forms a salt-bridge with protonated His57, and the proton is predominantly located at the His57 moiety. Glu214, the only acidic residue at the end of the H-bond network exhibits a pKa value of ∼6, slightly elevated due to solvation loss. It seems likely that the H-bond network [Asp85···His57···H2O···Glu214] serves as a proton-conducting pathway toward the protein bulk surface.

ZBTB7A forms a heterodimer with ZBTB2 and inhibits ZBTB2 homodimerization required for full activation of HIF-1.

Hypoxia-inducible factor 1 (HIF-1), recognized as a master transcription factor for adaptation to hypoxia, is associated with malignant characteristics and therapy resistance in cancers. It has become clear that cofactors such as ZBTB2 are critical for the full activation of HIF-1; however, the mechanisms downregulating the ZBTB2-HIF-1 axis remain poorly understood. In this study, we identified ZBTB7A as a negative regulator of ZBTB2 by analyzing protein sequences and structures. We found that ZBTB7A forms a heterodimer with ZBTB2, inhibits ZBTB2 homodimerization necessary for the full expression of ZBTB2-HIF-1 downstream genes, and ultimately delays the proliferation of cancer cells under hypoxic conditions. The Cancer Genome Atlas (TCGA) analyses revealed that overall survival is better in patients with high ZBTB7A expression in their tumor tissues. These findings highlight the potential of targeting the ZBTB7A-ZBTB2 interaction as a novel therapeutic strategy to inhibit HIF-1 activity and improve treatment outcomes in hypoxia-related cancers.

Mechanism of the formation of proton transfer pathways in photosynthetic reaction centers.

In photosynthetic reaction centers from purple bacteria (PbRCs) from Rhodobacter sphaeroides, the secondary quinone QB accepts two electrons and two protons via electron-coupled proton transfer (PT). Here, we identify PT pathways that proceed toward the QB binding site, using a quantum mechanical/molecular mechanical approach. As the first electron is transferred to QB, the formation of the Grotthuss-like pre-PT H-bond network is observed along Asp-L213, Ser-L223, and the distal QB carbonyl O site. As the second electron is transferred, the formation of a low-barrier H-bond is observed between His-L190 at Fe and the proximal QB carbonyl O site, which facilitates the second PT. As QBH2 leaves PbRC, a chain of water molecules connects protonated Glu-L212 and deprotonated His-L190 forms, which serves as a pathway for the His-L190 reprotonation. The findings of the second pathway, which does not involve Glu-L212, and the third pathway, which proceeds from Glu-L212 to His-L190, provide a mechanism for PT commonly used among PbRCs.

Stretching vibrational frequencies and pKa differences in H-bond networks of protein environments.

The experimentally measured stretching vibrational frequencies of O-D [νO-D(donor)] and C=O [νC=O(donor)] H-bond donor groups can provide valuable information about the H-bonds in proteins. Here, using a quantum mechanical/molecular mechanical approach, the relationship between these vibrational frequencies and the difference in pKa values between H-bond donor and acceptor groups [ΔpKa(donor … acceptor)] in bacteriorhodopsin and photoactive yellow protein environments was investigated. The results show that νO-D(donor) is correlated with ΔpKa(donor … acceptor), regardless of the specific protein environment. νC=O(donor) is also correlated with ΔpKa(donor … acceptor), although the correlation is weak because the C=O bond does not have a proton. Importantly, the shifts in νO-D(donor) and νC=O(donor) are not caused by changes in pKa(donor) alone, but rather by changes in ΔpKa(donor … acceptor). Specifically, a decrease in ΔpKa(donor … acceptor) can lead to proton release from the H-bond donor group toward the acceptor group, resulting in shifts in the vibrational frequencies of the protein environment. These findings suggest that changes in the stretching vibrational frequencies, in particular νO-D(donor), can be used to monitor proton transfer in protein environments.

Quantum mechanical analysis of excitation energy transfer couplings in photosystem II.

We evaluated excitation energy transfer (EET) coupling (J) between all pairs of chlorophylls (Chls) and pheophytins (Pheos) in the protein environment of photosystem II based on the time-dependent density functional theory with a quantum mechanical/molecular mechanics approach. In the reaction center, the EET coupling between Chls PD1 and PD2 is weaker (|J(PD1/PD2)| = 79 cm-1), irrespective of a short edge-to-edge distance of 3.6 Å (Mg-to-Mg distance of 8.1 Å), than the couplings between PD1 and the accessory ChlD1 (|J(PD1/ChlD2)| = 104 cm-1) and between PD2 and ChlD2 (|J(PD2/ChlD1)| = 101 cm-1), suggesting that PD1 and PD2 are two monomeric Chls rather than a "special pair". There exist strongly coupled Chl pairs (|J| > âˆ¼100 cm-1) in the CP47 and CP43 core antennas, which may be candidates for the red-shifted Chls observed in spectroscopic studies. In CP47 and CP43, Chls ligated to CP47-His26 and CP43-His56, which are located in the middle layer of the thylakoid membrane, play a role in the "hub" that mediates the EET from the lumenal to stromal layers. In the stromal layer, Chls ligated to CP47-His466, CP43-His441, and CP43-His444 mediate the EET from CP47 to ChlD2/PheoD2 and from CP43 to ChlD1/PheoD1 in the reaction center. Thus, the excitation energy from both CP47 and CP43 can always be utilized for the charge-separation reaction in the reaction center.

Energetic Diversity in the Electron-Transfer Pathways of Type I Photosynthetic Reaction Centers.

Photosynthetic reaction centers from heliobacteria (HbRC) and green sulfur bacteria (GsbRC) are homodimeric proteins and share a common ancestor with photosystem I (PSI), classified as type I reaction centers. Using the HbRC crystal structure, we calculated the redox potential (Em) values in the electron-transfer branches, solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable sites in the entire protein-pigment complex. Em(A-1) for bacteriochlorophyll g at the secondary site in HbRC (-1157 mV) is as low as Em(A-1) for chlorophyll a in PSI (-1173 mV). Em(A0/HbRC) is at the same level as Em(A0/GsbRC) and is 200 mV higher than Em(A0/PSI) due to the replacement of PsaA-Trp697/PsaB-Trp677 in PSI with PshA-Arg554 in HbRC. In contrast, Em(FX) for the Fe4S4 cluster in HbRC (-420 mV) is significantly higher than Em(FX) in GsbRC (-719 mV) and PSI (-705 mV) due to the absence of acidic residues that correspond to PscA-Asp634 in GsbRC and PsaB-Asp575 in PSI. It seems likely that type I reaction centers have evolved, adopting (bacterio)chlorophylls suitable for their light environments while maintaining electron-transfer cascades.

Mechanism of Asparagine-Mediated Proton Transfer in Photosynthetic Reaction Centers.

In photosynthetic reaction centers from purple bacteria (PbRCs), light-induced charge separation leads to the reduction of the terminal electron acceptor quinone, QB. The reduction of QB to QB•- is followed by protonation via Asp-L213 and Ser-L223 in PbRC from Rhodobacter sphaeroides. However, Asp-L213 is replaced with nontitratable Asn-L222 and Asn-L213 in PbRCs from Thermochromatium tepidum and Blastochloris viridis, respectively. Here, we investigated the energetics of proton transfer along the asparagine-involved H-bond network using a quantum mechanical/molecular mechanical approach. The potential energy profile for the H-bond between H3O+ and the carbonyl O site of Asn-L222 shows that the proton is predominantly localized at the Asn-L222 moiety in the T. tepidum PbRC protein environment, easily forming the enol species. The release of the proton from the amide -NH2 site toward Ser-L232 via tautomerization suffers from the energy barrier. Upon reorientation of Asn-L222, the enol -OH site forms a short low-barrier H-bond with Ser-L232, facilitating protonation of QB•- in a Grotthuss-like mechanism. This is a basis of how asparagine or glutamine side chains function as acceptors/donors in proton transfer pathways.

pH-Dependent Binding and Releasing Mechanism of Acetate in the Inner Water Cavity of Heliorhodopsin.

The high-resolution structure of heliorhodopsin crystallized at low pH reveals the presence of a planar triangle molecule, acetate, in the inner water cavity. Here, we investigate how the acetate molecule is stabilized at the counterion Glu107 moiety, using molecular dynamics (MD) simulations and a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations indicate that the density is best described as acetate among triangle acids, including nitric acid and bicarbonate. The calculated protonation state indicates that protonated acetate donates an H-bond to deprotonated Glu107 in the low-pH crystal structure. The observed red-shift of ∼30 nm in the absorption wavelength with pKa ≈ 4 is likely due to the His23/His80 protonation, rather than the Glu107 protonation. MD simulations also show that acetate can exist at the Glu107 moiety only when it is protonated. When ionized, acetate is released from the Glu107 moiety via Asn101 at the channel bottleneck and Arg91 on the intracellular protein surface. These observations could explain how acetate binds at low pH and releases at high pH.

Substitution of Ca2+ and changes in the H-bond network near the oxygen-evolving complex of photosystem II.

Ca2+, which provides binding sites for ligand water molecules W3 and W4 in the Mn4CaO5 cluster, is a prerequisite for O2 evolution in photosystem II (PSII). We report structural changes in the H-bond network and the catalytic cluster itself upon the replacement of Ca2+ with other alkaline earth metals, using a quantum mechanical/molecular mechanical approach. The small radius of Mg2+ makes W3 donate an H-bond to D1-Glu189 in Mg2+-PSII. If an additional water molecule binds at the large surface of Ba2+, it donates H-bonds to D1-Glu189 and the ligand water molecule at the dangling Mn, altering the H-bond network. The potential energy profiles of the H-bond between D1-Tyr161 (TyrZ) and D1-His190 and the interconversion between the open- and closed-cubane S2 conformations remain substantially unaltered upon the replacement of Ca2+. Remarkably, the O5⋯Ca2+ distance is shortest among all O5⋯metal distances irrespective of the radius being larger than that of Mg2+. Furthermore, Ca2+ is the only alkaline earth metal that equalizes the O5⋯metal and O2⋯metal distances and facilitates the formation of the symmetric cubane structure.

Acquirement of water-splitting ability and alteration of the charge-separation mechanism in photosynthetic reaction centers.

In photosynthetic reaction centers from purple bacteria (PbRC) and the water-oxidizing enzyme, photosystem II (PSII), charge separation occurs along one of the two symmetrical electron-transfer branches. Here we report the microscopic origin of the unidirectional charge separation, fully considering electron-hole interaction, electronic coupling of the pigments, and electrostatic interaction with the polarizable entire protein environments. The electronic coupling between the pair of bacteriochlorophylls is large in PbRC, forming a delocalized excited state with the lowest excitation energy (i.e., the special pair). The charge-separated state in the active branch is stabilized by uncharged polar residues in the transmembrane region and charged residues on the cytochrome c2 binding surface. In contrast, the accessory chlorophyll in the D1 protein (ChlD1) has the lowest excitation energy in PSII. The charge-separated state involves ChlD1•+ and is stabilized predominantly by charged residues near the Mn4CaO5 cluster and the proceeding proton-transfer pathway. It seems likely that the acquirement of water-splitting ability makes ChlD1 the initial electron donor in PSII.

Energetics of the Electron Transfer Pathways in the Homodimeric Photosynthetic Reaction Center.

Photosynthetic reaction centers from a green sulfur bacterium (GsbRC), the PscA/PscA proteins, and photosystem I (PSI), PsaA/PsaB proteins, share structural similarities. Here, we report the redox potential (Em) values of GsbRC by solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable sites in the entire GsbRC protein and identify the factors that shift the Em values with respect to PSI. The Em values for one-electron reduction of the accessory (A-1) and adjacent (A0) chlorophylls in GsbRC are 100-250 mV higher than those in PSI, whereas the Em values for the Fe4S4 cluster (FX) are at the same level. The PsaA-Trp697/PsaB-Trp677 pair in PSI, which forms the A1-quinone binding site, is replaced with PscA-Arg638 in GsbRC. PsaB-Asp575 in PSI, which is responsible for the Em difference between A1A and A1B quinones in PSI, is absent in GsbRC. These discrepancies also contribute to the upshift in Em(A-1) and Em(A0) in GsbRC with respect to PSI. It seems likely that the upshifted Em for chlorophylls in GsbRC ultimately originates from the characteristics of the electrostatic environment that corresponds to the A1 site of PSI.

Conformational Changes and H-Bond Rearrangements during Quinone Release in Photosystem II.

In photosystem II (PSII) and photosynthetic reaction centers from purple bacteria (PbRC), the electron released from the electronically excited chlorophyll is transferred to the terminal electron acceptor quinone, QB. QB accepts two electrons and two protons before leaving the protein. We investigated the molecular mechanism of quinone exchange in PSII, conducting molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations. MD simulations suggest that the release of QB leads to the transformation of the short helix (D1-Phe260 to D1-Ser264), which is adjacent to the stromal helix de (D1-Asn247 to D1-Ile259), into a loop and to the formation of a water-intake channel. Water molecules enter the QB binding pocket via the channel and form an H-bond network. QM/MM calculations indicate that the H-bond network serves as a proton-transfer pathway for the reprotonation of D1-His215, the proton donor during QBH-/QBH2 conversion. Together with the absence of the corresponding short helix but the presence of Glu-L212 in PbRC, it seems likely that the two type-II reaction centers undergo quinone exchange via different mechanisms.

Identification of intermediate conformations in the photocycle of the light-driven sodium-pumping rhodopsin KR2.

The light-driven rhodopsin KR2 transports Na+via the M- and O-states. However, the mechanisms by which the retinal regulates Na+ pumping is unknown, in part because KR2 adopts both pentamer and monomer forms in crystal structures and in part because these structures show differences in the protein conformation near the Schiff base, even when they are of the same intermediate state within the photocycle. A particular open question is the nature of the H-bond networks and protonation state in the active site, including Asp116. Here, we analyze the protonation state and the absorption wavelength for each crystal structure, using a quantum mechanical/molecular mechanical approach. In the pentamer ground state, the calculated absorption wavelength reproduces the experimentally measured absorption wavelength (530 nm). The analysis also shows that ionized Asp116 is stabilized by the H-bond donations of both Ser70 and a cluster of water molecules. The absorption wavelength of 400 nm in the M-state can be best reproduced when the two O atoms of Asp116 interact strongly with the Schiff base, as reported in one of the previous monomer ground state structures. The absorption wavelengths calculated for the two Na+-incorporated O-state structures are consistent with the measured absorption wavelength (∼600 nm), which suggests that two conformations represent the O-state. These results may provide a key to designing enhanced tools in optogenetics.

A self-assembled coordination cage enhances the reactivity of confined amides via mechanical bond-twisting.

Self-assembled coordination cages composed of metal cations and ligands can enhance the hydrolysis of non-covalently trapped amides in mild conditions as demonstrated in recent experiments. Here, we reveal the mechanism that accelerates base-catalyzed amide hydrolysis inside the octahedral coordination cage, by means of a quantum mechanics/molecular mechanics/polarizable continuum model. The calculated activation barrier of the nucleophilic OH- addition to a planar diaryl amide drastically decreases in the cage because of mechanical bond-twisting due to host-guest π-stacking. By contrast, the OH- addition to an N-acylindole, which possesses a twisted amide bond in bulk water, is not enhanced in the cage. Even though the cage hinders OH- collisions with the confined amide, the cage can twist the dihedral angle of the planar amide so as to mimic the transition state of OH- addition.

Role of redox-inactive metals in controlling the redox potential of heterometallic manganese-oxido clusters.

Photosystem II (PSII) contains Ca2+, which is essential to the oxygen-evolving activity of the catalytic Mn4CaO5 complex. Replacement of Ca2+ with other redox-inactive metals results in a loss/decrease of oxygen-evolving activity. To investigate the role of Ca2+ in this catalytic reaction, we investigate artificial Mn3[M]O2 clusters redox-inactive metals  [M] ([M] = Mg2+, Ca2+, Zn2+, Sr2+, and Y3+), which were synthesized by Tsui et al. (Nat Chem 5:293, 2013). The experimentally measured redox potentials (Em) of these clusters are best described by the energy of their highest occupied molecular orbitals. Quantum chemical calculations showed that the valence of metals predominantly affects Em(MnIII/IV), whereas the ionic radius of metals affects Em(MnIII/IV) only slightly.