Biophysical comparison of ATP synthesis mechanisms shows a kinetic advantage for the rotary process

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

WESTPA 2.0: High-Performance Upgrades for Weighted Ensemble Simulations and Analysis of Longer-Timescale Applications

The weighted ensemble (WE) family of methods is one of several statistical mechanics-based path sampling strategies that can provide estimates of key observables (rate constants and pathways) using a fraction of the time required by direct simulation methods such as molecular dynamics or discrete-state stochastic algorithms. WE methods oversee numerous parallel trajectories using intermittent overhead operations at fixed time intervals, enabling facile interoperability with any dynamics engine. Here, we report on the major upgrades to the WESTPA software package, an open-source, high-performance framework that implements both basic and recently developed WE methods. These upgrades offer substantial improvements over traditional WE methods. The key features of the new WESTPA 2.0 software enhance the efficiency and ease of use: an adaptive binning scheme for more efficient surmounting of large free energy barriers, streamlined handling of large simulation data sets, exponentially improved analysis of kinetics, and developer-friendly tools for creating new WE methods, including a Python API and resampler module for implementing both binned and "binless" WE strategies.

From OCR and ECAR to energy: Perspectives on the design and interpretation of bioenergetics studies

Biological energy transduction underlies all physiological phenomena in cells. The metabolic systems that support energy transduction have been of great interest due to their association with numerous pathologies including diabetes, cancer, rare genetic diseases, and aberrant cell death. Commercially available bioenergetics technologies (e.g., extracellular flux analysis, high-resolution respirometry, fluorescent dye kits, etc.) have made practical assessment of metabolic parameters widely accessible. This has facilitated an explosion in the number of studies exploring, in particular, the biological implications of oxygen consumption rate (OCR) and substrate level phosphorylation via glycolysis (i.e., via extracellular acidification rate (ECAR)). Though these technologies have demonstrated substantial utility and broad applicability to cell biology research, they are also susceptible to historical assumptions, experimental limitations, and other caveats that have led to premature and/or erroneous interpretations. This review enumerates various important considerations for designing and interpreting cellular and mitochondrial bioenergetics experiments, some common challenges and pitfalls in data interpretation, and some potential "next steps" to be taken that can address these highlighted challenges.

Poor Person's pH Simulation of Membrane Proteins

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.

Enhanced hydrogen production of Rhodobacter sphaeroides promoted by extracellular H+ of Halobacterium salinarum

Purpose

This study utilized the principle that the bacteriorhodopsin (BR) produced by Halobacterium salinarum could increase the hydrogen production of Rhodobacter sphaeroides. H. salinarum are co-cultured with R. sphaeroides to determine the impact of purple membrane fragments (PM) on R. sphaeroides and improve its hydrogen production capacity.

Methods

In this study, low-salinity in 14 % NaCl domesticates H salinarum. Then, 0–160 nmol of different concentration gradient groups of bacteriorhodopsin (BR) and R. sphaeroides was co-cultivated, and the hydrogen production and pH are measured; then, R. sphaeroides and immobilized BR of different concentrations are used to produce hydrogen to detect the amount of hydrogen. Two-chamber microbial hydrogen production system with proton exchange membrane-assisted proton flow was established, and the system was operated. As additional electricity added under 0.3 V, the hydrogen production rate increased with voltages in the coupled system.

Results

H salinarum can still grow well after low salt in 14% NaCl domestication. When the BR concentration is 80 nmol, the highest hydrogen production reached 217 mL per hour. Both immobilized PC (packed cells) and immobilized PM (purple membrane) of H. salinarum could promote hydrogen production of R. sphaeroides to some extent. The highest production of hydrogen was obtained by the coupled system with 40 nmol BR of immobilized PC, which increased from 127 to 232 mL, and the maximum H2 production rate was 18.2 mL−1 h−1 L culture. In the 192 h experiment time, when the potential is 0.3 V, the hydrogen production amount can reach 920 mL, which is 50.3% higher than the control group.

Conclusions

The stability of the system greatly improved after PC was immobilized, and the time for hydrogen production of R. sphaeroides significantly extended on same condition. As additional electricity added under 0.3 V, the hydrogen production rate increased with voltages in the coupled system. These results are helpful to build a hydrogen production-coupled system by nitrogenase of R. sphaeroides and proton pump of H. salinarum.

Graphical abstract

A Cerium Vanadate Nanozyme with Specific Superoxide Dismutase Activity Regulates Mitochondrial Function and ATP Synthesis in Neuronal Cells

Nanoparticles that functionally mimic the activity of metal-containing enzymes (metallo-nanozymes) are of therapeutic importance for treating various diseases. However, it is still not clear whether such nanozymes can completely substitute the function of natural enzymes in living cells. In this work, we show for the first time that a cerium vanadate (CeVO₄ ) nanozyme can substitute the function of superoxide dismutase 1 and 2 (SOD1 and SOD2) in the neuronal cells even when the natural enzyme is down-regulated by specific gene silencing. The nanozyme prevents the mitochondrial damage in SOD1- and SOD2-depleted cells by regulating the superoxide levels and restores the physiological levels of the anti-apoptotic Bcl-2 family proteins. Furthermore, the nanozyme effectively prevents the mitochondrial depolarization, leading to a significant improvement in the cellular levels of ATP under oxidative stress.

Nonequilibrium Energy Transduction in Stochastic Strongly Coupled Rotary Motors

Living systems at the molecular scale are composed of many constituents with strong and heterogeneous interactions, operating far from equilibrium, and subject to strong fluctuations. These conditions pose significant challenges to efficient, precise, and rapid free energy transduction, yet nature has evolved numerous molecular machines that do just this. Using a simple model of the ingenious rotary machine F_oF₁-ATP synthase, we investigate the interplay between nonequilibrium driving forces, thermal fluctuations, and interactions between strongly coupled subsystems. This model reveals design principles for effective free energy transduction. Most notably, while tight coupling is intuitively appealing, we find that output power is maximized at intermediate-strength coupling, which permits lubrication by stochastic fluctuations with only minimal slippage.

Opposing Pressures of Speed and Efficiency Guide the Evolution of Molecular Machines

Many biomolecular machines need to be both fast and efficient. How has evolution optimized these machines along the tradeoff between speed and efficiency? We explore this question using optimizable dynamical models along coordinates that are plausible evolutionary degrees of freedom. Data on 11 motors and ion pumps are consistent with the hypothesis that evolution seeks an optimal balance of speed and efficiency, where any further small increase in one of these quantities would come at great expense to the other. For F_oF₁-ATPases in different species, we also find apparent optimization of the number of subunits in the c-ring, which determines the number of protons pumped per ATP synthesized. Interestingly, these ATPases appear to more optimized for efficiency than for speed, which can be rationalized through their key role as energy transducers in biology. The present modeling shows how the dynamical performance properties of biomolecular motors and pumps may have evolved to suit their corresponding biological actions.

Voltage-energized Calcium-sensitive ATP Production by Mitochondria

Regulation of ATP production by mitochondria, critical to multicellular life, is poorly understood. Here we investigate the molecular controls of this process in heart and provide a framework for its Ca²⁺-dependent regulation. We find that the entry of Ca²⁺ into the matrix through the mitochondrial calcium uniporter (MCU) in heart has neither an apparent cytosolic Ca²⁺ threshold nor gating function and guides ATP production by its influence on the inner mitochondrial membrane (IMM) potential, ΔΨ_m. This regulation occurs by matrix Ca²⁺-dependent modulation of pyruvate and glutamate dehydrogenase activity and not through any effect of Ca²⁺ on ATP Synthase or on Electron Transport Chain Complexes II, III or IV. Examining the ΔΨ_m dependence of ATP production over the range of -60 mV to -170 mV in detail reveals that cardiac ATP synthase has a voltage dependence that distinguishes it fundamentally from the previous standard, the bacterial ATP synthase. Cardiac ATP synthase operates with a different ΔΨ_m threshold for ATP production than bacterial ATP synthase and reveals a concave-upwards shape without saturation. Skeletal muscle MCU Ca²⁺ flux, while also having no apparent cytosolic Ca²⁺ threshold, is substantially different from the cardiac MCU, yet the ATP synthase voltage dependence in skeletal muscle is identical to that in the heart. These results suggest that while the conduction of cytosolic Ca²⁺ signals through the MCU appears to be tissue-dependent, as shown by earlier work¹, the control of ATP synthase by ΔΨ_m appears to be broadly consistent among tissues but is clearly different from bacteria.

Theory of Nonequilibrium Free Energy Transduction by Molecular Machines

Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared with more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

Mechanisms for achieving high speed and efficiency in biomolecular machines

How does a biomolecular machine achieve high speed at high efficiency? We explore optimization principles using a simple two-state dynamical model. With this model, we establish physical principles-such as the optimal way to distribute free-energy changes and barriers across the machine cycle-and connect them to biological mechanisms. We find that a machine can achieve high speed without sacrificing efficiency by varying its conformational free energy to directly link the downhill, chemical energy to the uphill, mechanical work and by splitting a large work step into more numerous, smaller substeps. Experimental evidence suggests that these mechanisms are commonly used by biomolecular machines. This model is useful for exploring questions of evolution and optimization in molecular machines.

Optimal control of rotary motors

Single-molecule experiments have found near-perfect thermodynamic efficiency in the rotary motor F_{1}-ATP synthase. To help elucidate the principles underlying nonequilibrium energetic efficiency in such stochastic machines, we investigate driving protocols that minimize dissipation near equilibrium in a simple model rotary mechanochemical motor, as determined by a generalized friction coefficient. Our simple model has a periodic friction coefficient that peaks near system energy barriers. This implies a minimum-dissipation protocol that proceeds rapidly when the system is overwhelmingly in a single macrostate but slows significantly near energy barriers, thereby harnessing thermal fluctuations to kick the system over energy barriers with minimal work input. This model also manifests a phenomenon not seen in otherwise similar nonperiodic systems: Sufficiently fast protocols can effectively lap the system. While this leads to a trade-off between accuracy of driving and energetic cost, we find that our designed protocols outperform naive protocols.

Allocating dissipation across a molecular machine cycle to maximize flux

Biomolecular machines consume free energy to break symmetry and make directed progress. Nonequilibrium ATP concentrations are the typical free energy source, with one cycle of a molecular machine consuming a certain number of ATP, providing a fixed free energy budget. Since evolution is expected to favor rapid-turnover machines that operate efficiently, we investigate how this free energy budget can be allocated to maximize flux. Unconstrained optimization eliminates intermediate metastable states, indicating that flux is enhanced in molecular machines with fewer states. When maintaining a set number of states, we show that-in contrast to previous findings-the flux-maximizing allocation of dissipation is not even. This result is consistent with the coexistence of both "irreversible" and reversible transitions in molecular machine models that successfully describe experimental data, which suggests that, in evolved machines, different transitions differ significantly in their dissipation.

Biophysical comparison of ATP-driven proton pumping mechanisms suggests a kinetic advantage for the rotary process depending on coupling ratio

ATP-driven proton pumps, which are critical to the operation of a cell, maintain cytosolic and organellar pH levels within a narrow functional range. These pumps employ two very different mechanisms: an elaborate rotary mechanism used by V-ATPase H⁺ pumps, and a simpler alternating access mechanism used by P-ATPase H⁺ pumps. Why are two different mechanisms used to perform the same function? Systematic analysis, without parameter fitting, of kinetic models of the rotary, alternating access and other possible mechanisms suggest that, when the ratio of protons transported per ATP hydrolyzed exceeds one, the one-at-a-time proton transport by the rotary mechanism is faster than other possible mechanisms across a wide range of driving conditions. When the ratio is one, there is no intrinsic difference in the free energy landscape between mechanisms, and therefore all mechanisms can exhibit the same kinetic performance. To our knowledge all known rotary pumps have an H⁺:ATP ratio greater than one, and all known alternating access ATP-driven proton pumps have a ratio of one. Our analysis suggests a possible explanation for this apparent relationship between coupling ratio and mechanism. When the conditions under which the pump must operate permit a coupling ratio greater than one, the rotary mechanism may have been selected for its kinetic advantage. On the other hand, when conditions require a coupling ratio of one or less, the alternating access mechanism may have been selected for other possible advantages resulting from its structural and functional simplicity.