Molecular dynamics simulations of ubiquinone; a survey over torsional potentials and hydrogen bonds

Chemistry of Lipoquinones: Properties, Synthesis, and Membrane Location of Ubiquinones, Plastoquinones, and Menaquinones

Lipoquinones are the topic of this review and are a class of hydrophobic lipid molecules with key biological functions that are linked to their structure, properties, and location within a biological membrane. Ubiquinones, plastoquinones, and menaquinones vary regarding their quinone headgroup, isoprenoid sidechain, properties, and biological functions, including the shuttling of electrons between membrane-bound protein complexes within the electron transport chain. Lipoquinones are highly hydrophobic molecules that are soluble in organic solvents and insoluble in aqueous solution, causing obstacles in water-based assays that measure their chemical properties, enzyme activities and effects on cell growth. Little is known about the location and ultimately movement of lipoquinones in the membrane, and these properties are topics described in this review. Computational studies are particularly abundant in the recent years in this area, and there is far less experimental evidence to verify the often conflicting interpretations and conclusions that result from computational studies of very different membrane model systems. Some recent experimental studies have described using truncated lipoquinone derivatives, such as ubiquinone-2 (UQ-2) and menaquinone-2 (MK-2), to investigate their conformation, their location in the membrane, and their biological function. Truncated lipoquinone derivatives are soluble in water-based assays, and hence can serve as excellent analogs for study even though they are more mobile in the membrane than the longer chain counterparts. In this review, we will discuss the properties, location in the membrane, and syntheses of three main classes of lipoquinones including truncated derivatives. Our goal is to highlight the importance of bridging the gap between experimental and computational methods and to incorporate properties-focused considerations when proposing future studies relating to the function of lipoquinones in membranes.

Electron Transport Lipids Fold Within Membrane-Like Interfaces

Lipoquinones, such as ubiquinones (UQ) and menaquinones (MK), function as essential lipid components of the electron transport system (ETS) by shuttling electrons and protons to facilitate the production of ATP in eukaryotes and prokaryotes. Lipoquinone function in membrane systems has been widely studied, but the exact location and conformation within membranes remains controversial. Lipoquinones, such as Coenzyme Q (UQ-10), are generally depicted simply as "Q" in life science diagrams or in extended conformations in primary literature even though specific conformations are important for function in the ETS. In this study, our goal was to determine the location, orientation, and conformation of UQ-2, a truncated analog of UQ-10, in model membrane systems and to compare our results to previously studied MK-2. Herein, we first carried out a six-step synthesis to yield UQ-2 and then demonstrated that UQ-2 adopts a folded conformation in organic solvents using 1H-1H 2D NOESY and ROESY NMR spectroscopic studies. Similarly, using 1H-1H 2D NOESY NMR spectroscopic studies, UQ-2 was found to adopt a folded, U-shaped conformation within the interface of an AOT reverse micelle model membrane system. UQ-2 was located slightly closer to the surfactant-water interface compared to the more hydrophobic MK-2. In addition, Langmuir monolayer studies determined UQ-2 resided within the monolayer water-phospholipid interface causing expansion, whereas MK-2 was more likely to be compressed out and reside within the phospholipid tails. All together these results support the model that lipoquinones fold regardless of the headgroup structure but that the polarity of the headgroup influences lipoquinone location within the membrane interface. These results have implications regarding the redox activity near the interface as quinone vs. quinol forms may facilitate locomotion of lipoquinones within the membrane. The location, orientation, and conformation of lipoquinones are critical for their function in generating cellular energy within membrane ETS, and the studies described herein shed light on the behavior of lipoquinones within membrane-like environments.

1,2H hyperfine spectroscopy and DFT modeling unveil the demethylmenasemiquinone binding mode to E. coli nitrate reductase A (NarGHI)

The quinol oxidation site Q_D in E. coli respiratory nitrate reductase A (EcNarGHI) reacts with the three isoprenoid quinones naturally synthesized by the bacterium, i.e. ubiquinones (UQ), menaquinones (MK) and demethylmenaquinones (DMK). The binding mode of the demethylmenasemiquinone (DMSK) intermediate to the EcNarGHI Q_D quinol oxidation site is analyzed in detail using ^1,2H hyperfine (hf) spectroscopy in combination with H₂O/D₂O exchange experiments and DFT modeling, and compared to the menasemiquinone one bound to the Q_D site (MSK_D) previously studied by us. DMSK_D and MSK_D are shown to bind in a similar and strongly asymmetric manner through a short (~1.7 Å) H-bond. The origin of the specific hf pattern resolved on the DMSK_D field-swept EPR spectrum is unambiguously ascribed to slightly inequivalent contributions from two β-methylene protons of the isoprenoid side chain. DFT calculations show that their large isotropic hf coupling constants (A_iso ~12 and 15 MHz) are consistent with both (i) a specific highly asymmetric binding mode of DMSK_D and (ii) a near in-plane orientation of its isoprenyl chain at Cβ relative to the aromatic ring, which differs by ~90° to that predicted for free or NarGHI-bound MSK. Our results provide new insights into how the conformation and the redox properties of different natural quinones are selectively fine-tuned by the protein environment at a single Q site. Such a fine-tuning most likely contributes to render NarGHI as an efficient and flexible respiratory enzyme to be used upon rapid variations of the Q-pool content.

Distribution and dynamics of quinones in the lipid bilayer mimicking the inner membrane of mitochondria

Quinone and its analogues (Q) constitute an important class of compounds that perform key electron transfer reactions in oxidative- and photo-phosphorylation. In the inner membrane of mitochondria, ubiquinone molecules undergo continuous redox transitions enabling electron transfer between the respiratory complexes. In such a dynamic system undergoing continuous turnover for ATP synthesis, an uninterrupted supply of substrate molecules is absolutely necessary. In the current work, we have performed atomistic molecular dynamics simulations and free energy calculations to assess the structure, dynamics, and localization of quinone and its analogues in a lipid bilayer, whose composition mimics the one in the inner mitochondrial membrane. The results show that there is a strong tendency of both quinone and quinol molecules to localize in the vicinity of the lipids' acyl groups, right under the lipid head group region. Additionally, we observe a second location in the middle of the bilayer where quinone molecules tend to stabilize. Translocation of quinone through a lipid bilayer is very fast and occurs in 10-100ns time scale, whereas the translocation of quinol is at least an order of magnitude slower. We suggest that this has important mechanistic implications given that the localization of Q ensures maximal occupancy of the Q-binding sites or Q-entry points in electron transport chain complexes, thereby maintaining an optimal turnover rate for ATP synthesis.

Impact of Antioxidants on Cardiolipin Oxidation in Liposomes: Why Mitochondrial Cardiolipin Serves as an Apoptotic Signal?

Molecules of mitochondrial cardiolipin (CL) get selectively oxidized upon oxidative stress, which triggers the intrinsic apoptotic pathway. In a chemical model most closely resembling the mitochondrial membrane-liposomes of pure bovine heart CL-we compared ubiquinol-10, ubiquinol-6, and alpha-tocopherol, the most widespread naturally occurring antioxidants, with man-made, quinol-based amphiphilic antioxidants. Lipid peroxidation was induced by addition of an azo initiator in the absence and presence of diverse antioxidants, respectively. The kinetics of CL oxidation was monitored via formation of conjugated dienes at 234 nm. We found that natural ubiquinols and ubiquinol-based amphiphilic antioxidants were equally efficient in protecting CL liposomes from peroxidation; the chromanol-based antioxidants, including alpha-tocopherol, were 2-3 times less efficient. Amphiphilic antioxidants, but not natural ubiquinols and alpha-tocopherol, were able, additionally, to protect the CL bilayer from oxidation by acting from the water phase. We suggest that the previously reported therapeutic efficiency of mitochondrially targeted amphiphilic antioxidants is owing to their ability to protect those CL molecules that are inaccessible to natural hydrophobic antioxidants, being trapped within respiratory supercomplexes. The high susceptibility of such occluded CL molecules to oxidation may have prompted their recruitment as apoptotic signaling molecules by nature.

Partition, orientation and mobility of ubiquinones in a lipid bilayer

Ubiquinone is the universal mobile charge carrier involved in biological electron transfer processes. Its redox properties and biological function depend on the molecular partition and lateral diffusion over biological membranes. However, ubiquinone localization and dynamics within lipid bilayers are long debated and still uncertain. Here we present molecular dynamics simulations of several ubiquinone homologs with variable isoprenoid tail lengths complexed to phosphatidylcholine bilayers. Initially, a new force-field parametrization for ubiquinone is derived from and compared to high level quantum chemical data. Free energy profiles for ubiquinone insertion in the lipid bilayer are obtained with the new force-field. The profiles allow for the determination of the equilibrium location of ubiquinone in the membrane as well as for the validation of the simulation model by direct comparison with experimental partition coefficients. A detailed analysis of structural properties and interactions shows that the ubiquinone polar head group is localized at the water-bilayer interface at the same depth of the lipid glycerol groups and oriented normal to the membrane plane. Both the localization and orientation of ubiquinone head groups do not change significantly when increasing the number of isoprenoid units. The isoprenoid tail is extended and packed with the lipid acyl chains. For ubiquinones with long tails, the terminal isoprenoid units have high flexibility. Calculated ubiquinone diffusion coefficients are similar to that found for the phosphatidylcholine lipid. These results may have further implications for the mechanisms of ubiquinone transport and binding to respiratory and photosynthetic protein complexes.

Redox potential tuning through differential quinone binding in the photosynthetic reaction center of Rhodobacter sphaeroides

Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, we have investigated and characterized the interactions of the protein with the quinones in the Q(A) and Q(B) sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the Q(B) site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the Q(A) and Q(B) sites. Disruption of this interaction leads to weaker binding in a ubiquinone analogue that lacks a 2-methoxy group, a finding supported by reverse electron transfer electron paramagnetic resonance experiments of the Q(A)⁻Q(B)⁻ biradical and competitive binding assays.

p-Benzoquinone in aqueous solution: Stark shifts in spectra, asymmetry in solvent structure

Results from a simulation of p-benzoquinone (PBQ) in water is presented. An explicit solvent representation is used together with a multiconfigurational ab initio quantum chemical method. The electronic n --> pi* transitions are studied in aqueous solution and the two such transitions are both blue-shifted but to different degree. Both non-equilibrium and many-body effects are found to have decisive influence on the solvation: despite stronger hydrogen bonding between solute and solvent in an excited state than in the ground state, there is a blue-shift, and the solvent structure around the non-polar PBQ is asymmetric, which is argued to come from special many-body effects. The unusual result of strengthened hydrogen bonds in the excited state that follows from an excitation of a non-bonding electron on a proton acceptor, is explained by the near-linear Stark shift that is present in the transition.

Extended Car-Parrinello molecular dynamics and electronic g-tensors study of benzosemiquinone radical anion

Car-Parrinello molecular dynamics simulations of benzoquinone and benzosemiquinone radical anion in both aqueous solution and the gas phase have been carried out at ambient conditions. Hydrogen bonding is considerably more extensive to the anionic than to the neutral aqueous system. In addition to the conventional hydrogen bonding to the carbonyl oxygen atoms, T-stacked hydrogen bonding to the pi-system is statistically and energetically significant for the semiquinone anion but not for the neutral quinone. EPR g-tensors have been calculated at DFT level for snapshots taken at regular intervals from the gas-phase and solution semiquinone anion trajectories. Different criteria for extraction of semiquinone/water clusters from the solution trajectory give insight into the effects of different interactions on the g-tensor, as does correlation of the g-tensor with statistically significant hydrogen-bond configurations identified along the trajectories. Comparison of gas-phase and solution results indicates opposite directions of direct electronic and indirect structural influences of hydrogen bonding on g-tensors. Short-time oscillations in g(x) along the trajectory are due mainly to C-O bond vibrations.

Ab Initio Molecular Dynamics Simulations and g-Tensor Calculations of Aqueous Benzosemiquinone Radical Anion: Effects of Regular and “T-Stacked” Hydrogen Bonds

Car-Parrinello molecular dynamics (CP-MD) simulations of the benzosemiquinone radical anion in aqueous solution have been performed at ambient conditions. Analysis of the trajectory shows not only extensive hydrogen bonding to the carbonyl oxygen atoms (ca. 4-5.6 water molecules depending on distance criteria), but also relatively long-lived "T-stacked" hydrogen bonds to the semiquinone pi-system. These results are discussed in the context of recent findings on semiquinone-protein interactions in photosynthetic reaction centers, and of EPR and vibration spectroscopical data for the aqueous system. Snapshots from the CP-MD trajectory are used for the first quantum chemical analyses of dynamical effects on electronic g-tensors, using cluster models and a recently developed density functional method. In particular, the effects of intermolecular hydrogen-bond dynamics on the g-tensor components are examined, in comparison with recent EPR and ENDOR studies.