Search

The photosynthetic bacterial reaction centers from purple non-sulfur bacteria use light energy to drive the transfer of electrons from cytochrome c to ubiquinone. Ubiquinone bound in the QA site cycles between quinone, QA, and anionic semiquinone, QA·−, being reduced once and never binding protons. In the QB site, ubiquinone is reduced twice by QA·−, binds two protons and is released into the membrane as the quinol, QH2. The network of hydrogen bonds formed in a molecular dynamics trajectory was drawn to investigate proton transfer pathways from the cytoplasm to each quinone binding site. QA is isolated with no path for protons to enter from the surface. In contrast, there is a complex and tangled network requiring residues and waters that can bring protons to QB. There are three entries from clusters of surface residues centered around HisH126, GluH224, and HisH68. The network is in good agreement with earlier studies, Mutation of key nodes in the network, such as SerL223, were previously shown to slow proton delivery. Mutational studies had also shown that double mutations of residues such as AspM17 and AspL210 along multiple paths in the network presented here slow the reaction, while single mutations do not. Likewise, mutation of both HisH126 and HisH128, which are at the entry to two paths reduce the rate of proton uptake. [ABSTRACT FROM AUTHOR]

Water molecules play a key role in all biochemical processes. They help define the shape of proteins, and they are reactant or product in many reactions and are released as ligands are bound. They facilitate the transfer of protons through transmembrane proton channel, pump and transporter proteins. Continuum electrostatics (CE) force fields used by program Multiconformation CE (MCCE) capture electrostatic interactions in biomolecules with an implicit solvent, which captures the averaged solvent water equilibrium properties. Hybrid CE methods can use explicit water molecules within the protein surrounded by implicit solvent. These hybrid methods permit the study of explicit hydrogen bond networks within the protein and allow analysis of processes such as proton transfer reactions. Yet hybrid CE methods have not been rigorously tested. Here, we present an explicit treatment of water molecules in the Gramicidin A (gA) channel using MCCE and compare the resulting distributions of water molecules and key hydration features against those obtained with explicit solvent Molecular Dynamics (MD) simulations with the nonpolarizable CHARMM36 and polarizable Drude force fields. CHARMM36 leads to an aligned water wire in the channel characterized by a large absolute net water dipole moment; the MCCE and Drude analysis lead to a small net dipole moment as the water molecules change orientation within the channel. The correct orientation is not as yet known, so these calculations identify an open question. • Water is the primary cellular solvent, yet is challenging to simulate computationally. Here we simulate water molecules in the Gramicidin A channel comparing Monte Carlo (MC) sampling with a continuum electrostatics and Molecular Dynamics (MD) calculations with the non-polarizable CHARMM36 and polarizable Drude force fields. • These give different water properties, with classical MD yielding well oriented water wires, while the Drude or continuum electrostatics force fields lead to more disordered water molecules, often changing orientation in the middle of the channel. [ABSTRACT FROM AUTHOR]

The oxidation of water to O2 is catalyzed by the Oxygen Evolving Complex (OEC), a Mn4CaO5 complex in Photosystem II (PSII). The OEC is sequentially oxidized from state S0 to S4. The S2 state, (MnIII)(MnIV)3, coexists in two redox isomers: S2,g=2, where Mn4 is MnIV and S2,g=4.1, where Mn1 is MnIV. Mn4 has two terminal water ligands, whose proton affinity is affected by the Mn oxidation state. The relative energy of the two S2 redox isomers and the protonation state of the terminal water ligands are analyzed using classical multi-conformer continuum electrostatics (MCCE). The Monte Carlo simulations are done on QM/MM optimized S1 and S2 structures docked back into the complete PSII, keeping the protonation state of the protein at equilibrium with the OEC redox and protonation states. Wild-type PSII, chloride-depleted PSII, PSII in the presence of oxidized YZ/protonated D1-H190, and the PSII mutants D2-K317A, D1-D61A, and D1-S169A are studied at pH 6. The wild-type PSII at pH 8 is also described. In qualitative agreement with experiment, in wild-type PSII, the S2,g=2 redox isomer is the lower energy state; while chloride depletion or pH 8 stabilizes the S2,g=4.1 state and the mutants D2-K317A, D1-D61A, and D1-S169A favor the S2,g=2 state. The protonation states of D1-E329, D1-E65, D1-H337, D1-D61, and the terminal waters on Mn4 (W1 and W2) are affected by the OEC oxidation state. The terminal W2 on Mn4 is a mixture of water and hydroxyl in the S2,g=2 state, indicating the two water protonation states have similar energy, while it remains neutral in the S1 and S2,g=4.1 states. In wild-type PSII, advancement to S2 leads to negligible proton loss and so there is an accumulation of positive charge. In the analyzed mutations and Cl− depleted PSII, additional deprotonation is found upon formation of S2 state. [ABSTRACT FROM AUTHOR]

pH conditions are central to the functioning of all biomolecules. However, implications of pH changes are nontrivial on a molecular scale. Though a rigorous microscopic definition of pH exists, its implementation in classical molecular dynamics (MD) simulations is cumbersome, and more so in large integral membrane systems. In this chapter, an integrative pipeline is described that combines Multi-Conformation Continuum Electrostatics (MCCE) computations with MD simulations to capture the effect of transient protonation states on the coupled conformational changes in transmembrane proteins. The core methodologies are explained, and all the software required to set up this pipeline are outlined with their key parameters. All associated analyses of structure and function are provided using two case studies, namely those of bioenergetic complexes: NADH dehydrogenase (complex I) and V o domain of V-type ATPase. The hybrid MCCE-MD pipeline has allowed the discovery of hydrogen bond networks, ligand binding pathways, and disease-causing mutations., (© 2021. Springer Science+Business Media, LLC, part of Springer Nature.)