Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of Vo complex

Theoretical insights into rotary mechanism of MotAB in the bacterial flagellar motor

Many bacteria enable locomotion by rotating their flagellum. It has been suggested that this rotation is realized by the rotary motion of the stator unit, MotAB, which is driven by proton transfer across the membrane. Recent cryo-electron microscopy studies have revealed a 5:2 MotAB configuration, in which a MotB dimer is encircled by a ring-shaped MotA pentamer. Although the structure implicates the rotary motion of the MotA wheel around the MotB axle, the molecular mechanisms of rotary motion and how they are coupled with proton transfer across the membrane remain elusive. In this study, we built a structure-based computational model for Campylobacter jejuni MotAB, conducted comprehensive protonation-state-dependent molecular dynamics simulations, and revealed a plausible proton-transfer-coupled rotation pathway. The model assumes rotation-dependent proton transfer, in which proton uptake from the periplasmic side to the conserved aspartic acid in MotB is followed by proton hopping to the MotA proton-carrying site, followed by proton export to the CP. We suggest that, by maintaining two of the proton-carrying sites of MotA in the deprotonated state, the MotA pentamer robustly rotates by ∼36° per proton transfer across the membrane. Our results provide a structure-based mechanistic model of the rotary motion of MotAB in bacterial flagellar motors and provide insights into various ion-driven rotary molecular motors.

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores.

Molecular mechanism of Oxr1p mediated disassembly of yeast V-ATPase

The eukaryotic vacuolar H+-ATPase (V-ATPase) is regulated by reversible disassembly into autoinhibited V1-ATPase and Vo proton channel subcomplexes. We recently reported that the TLDc protein Oxr1p induces V-ATPase disassembly in vitro. Whether and how Oxr1p is involved in enzyme disassembly in vivo, however, is not known. Here, using yeast genetics and fluorescence microscopy, we show that Oxr1p is essential for efficient V-ATPase disassembly in the cell. Supporting biochemical and biophysical in vitro experiments show that whereas Oxr1p-driven holoenzyme disassembly can occur in the absence of nucleotides, the presence of ATP greatly accelerates the process. ATP hydrolysis is needed, however, for subsequent release of Oxr1p so that the free V1 can adopt the autoinhibited conformation. Overall, our study unravels the molecular mechanism of Oxr1p-induced disassembly that occurs in vivo as part of the canonical V-ATPase regulation by reversible disassembly.

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.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

Mechanism of proton-powered c-ring rotation in a mitochondrial ATP synthase

Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.

The effect of lipid composition on the thermal stability of nanodiscs

Discoidal lipid nanoparticles (LNPs) called Nanodiscs (NDs) are derived from human high-density lipoprotein (HDL). Such biomimetics are ideally suited for the stabilization and delivery of pharmaceuticals, including chemicals, bio-active proteins and vaccines. The stability and circulation lifetimes of reconstituted HDL nanoparticles, including NDs, are variable. Lipids found in thermophilic archaea and bacteria are prime candidates for the stabilization of LNPs. We report the thermal stability of NDs prepared with lipids that differ in saturation, have either ether- or ester linkages between the fatty acid and glycerol backbone or contain isoprenoid fatty acid tails (phytanyl lipids). NDs with two saturated fatty acids show a much greater long-term thermostability than NDs with an unsaturated fatty acid. Ether fatty acid linkages, commonly found in thermophiles, did not improve stability of NDs compared to ester fatty acid linkages when using saturated lipids. NDs containing phytanyl and saturated alkyl fatty acids show similar stability at 37 °C. NDs assembled with phytanyl lipids contain three copies of the membrane scaffolding protein as opposed to the canonical dimer found in conventional NDs. The findings present a strong basis for the production of thermostable NDs through the selection of appropriate lipids and are likely broadly applicable to LNP development.

Fostering discoveries in the era of exascale computing: How the next generation of supercomputers empowers computational and experimental biophysics alike

Over a century ago, physicists started broadly relying on theoretical models to guide new experiments. Soon thereafter, chemists began doing the same. Now, biological research enters a new era when experiment and theory walk hand in hand. Novel software and specialized hardware became essential to understand experimental data and propose new models. In fact, current petascale computing resources already allow researchers to reach unprecedented levels of simulation throughput to connect in silico and in vitro experiments. The reduction in cost and improved access allowed a large number of research groups to adopt supercomputing resources and techniques. Here, we outline how large-scale computing has evolved to expand decades-old research, spark new research efforts, and continuously connect simulation and observation. For instance, multiple publicly and privately funded groups have dedicated extensive resources to develop artificial intelligence tools for computational biophysics, from accelerating quantum chemistry calculations to proposing protein structure models. Moreover, advances in computer hardware have accelerated data processing from single-molecule experimental observations and simulations of chemical reactions occurring throughout entire cells. The combination of software and hardware has opened the way for exascale computing and the production of the first public exascale supercomputer, Frontier, inaugurated by the Oak Ridge National Laboratory in 2022. Ultimately, the popularization and development of computational techniques and the training of researchers to use them will only accelerate the diversification of tools and learning resources for future generations.

Tender love and disassembly: How a TLDc domain protein breaks the V-ATPase

Vacuolar ATPases (V-ATPases, V1 Vo -ATPases) are rotary motor proton pumps that acidify intracellular compartments, and, when localized to the plasma membrane, the extracellular space. V-ATPase is regulated by a unique process referred to as reversible disassembly, wherein V1 -ATPase disengages from Vo proton channel in response to diverse environmental signals. Whereas the disassembly step of this process is ATP dependent, the (re)assembly step is not, but requires the action of a heterotrimeric chaperone referred to as the RAVE complex. Recently, an alternative pathway of holoenzyme disassembly was discovered that involves binding of Oxidation Resistance 1 (Oxr1p), a poorly characterized protein implicated in oxidative stress response. Unlike conventional reversible disassembly, which depends on enzyme activity, Oxr1p induced dissociation can occur in absence of ATP. Yeast Oxr1p belongs to the family of TLDc domain containing proteins that are conserved from yeast to mammals, and have been implicated in V-ATPase function in a variety of tissues. This brief perspective summarizes what we know about the molecular mechanisms governing both reversible (ATP dependent) and Oxr1p driven (ATP independent) V-ATPase dissociation into autoinhibited V1 and Vo subcomplexes.

The cytosolic N-terminal domain of V-ATPase a-subunits is a regulatory hub targeted by multiple signals

Vacuolar H+-ATPases (V-ATPases) acidify several organelles in all eukaryotic cells and export protons across the plasma membrane in a subset of cell types. V-ATPases are multisubunit enzymes consisting of a peripheral subcomplex, V1, that is exposed to the cytosol and an integral membrane subcomplex, Vo, that contains the proton pore. The Vo a-subunit is the largest membrane subunit and consists of two domains. The N-terminal domain of the a-subunit (aNT) interacts with several V1 and Vo subunits and serves to bridge the V1 and Vo subcomplexes, while the C-terminal domain contains eight transmembrane helices, two of which are directly involved in proton transport. Although there can be multiple isoforms of several V-ATPase subunits, the a-subunit is encoded by the largest number of isoforms in most organisms. For example, the human genome encodes four a-subunit isoforms that exhibit a tissue- and organelle-specific distribution. In the yeast S. cerevisiae, the two a-subunit isoforms, Golgi-enriched Stv1 and vacuolar Vph1, are the only V-ATPase subunit isoforms. Current structural information indicates that a-subunit isoforms adopt a similar backbone structure but sequence variations allow for specific interactions during trafficking and in response to cellular signals. V-ATPases are subject to several types of environmental regulation that serve to tune their activity to their cellular location and environmental demands. The position of the aNT domain in the complex makes it an ideal target for modulating V1-Vo interactions and regulating enzyme activity. The yeast a-subunit isoforms have served as a paradigm for dissecting interactions of regulatory inputs with subunit isoforms. Importantly, structures of yeast V-ATPases containing each a-subunit isoform are available. Chimeric a-subunits combining elements of Stv1NT and Vph1NT have provided insights into how regulatory inputs can be integrated to allow V-ATPases to support cell growth under different stress conditions. Although the function and distribution of the four mammalian a-subunit isoforms present additional complexity, it is clear that the aNT domains of these isoforms are also subject to multiple regulatory interactions. Regulatory mechanisms that target mammalian a-subunit isoforms, and specifically the aNT domains, will be described. Altered V-ATPase function is associated with multiple diseases in humans. The possibility of regulating V-ATPase subpopulations via their isoform-specific regulatory interactions are discussed.

Structural and Functional Diversity of Two ATP-Driven Plant Proton Pumps

Two ATP-dependent proton pumps function in plant cells. Plasma membrane H⁺-ATPase (PM H⁺-ATPase) transfers protons from the cytoplasm to the apoplast, while vacuolar H⁺-ATPase (V-ATPase), located in tonoplasts and other endomembranes, is responsible for proton pumping into the organelle lumen. Both enzymes belong to two different families of proteins and, therefore, differ significantly in their structure and mechanism of action. The plasma membrane H⁺-ATPase is a member of the P-ATPases that undergo conformational changes, associated with two distinct E1 and E2 states, and autophosphorylation during the catalytic cycle. The vacuolar H⁺-ATPase represents rotary enzymes functioning as a molecular motor. The plant V-ATPase consists of thirteen different subunits organized into two subcomplexes, the peripheral V₁ and the membrane-embedded V₀, in which the stator and rotor parts have been distinguished. In contrast, the plant plasma membrane proton pump is a functional single polypeptide chain. However, when the enzyme is active, it transforms into a large twelve-protein complex of six H⁺-ATPase molecules and six 14-3-3 proteins. Despite these differences, both proton pumps can be regulated by the same mechanisms (such as reversible phosphorylation) and, in some processes, such as cytosolic pH regulation, may act in a coordinated way.

Structural basis of V-ATPase VO region assembly by Vma12p, 21p, and 22p

Vacuolar-type adenosine triphosphatases (V-ATPases) are rotary proton pumps that acidify specific intracellular compartments in almost all eukaryotic cells. These multi-subunit enzymes consist of a soluble catalytic V₁ region and a membrane-embedded proton-translocating V_O region. V_O is assembled in the endoplasmic reticulum (ER) membrane, and V₁ is assembled in the cytosol. However, V₁ binds V_O only after V_O is transported to the Golgi membrane, thereby preventing acidification of the ER. We isolated V_O complexes and subcomplexes from Saccharomyces cerevisiae bound to V-ATPase assembly factors Vma12p, Vma21p, and Vma22p. Electron cryomicroscopy shows how the Vma12-22p complex recruits subunits a, e, and f to the rotor ring of V_O while blocking premature binding of V₁. Vma21p, which contains an ER-retrieval motif, binds the V_O:Vma12-22p complex, "mature" V_O, and a complex that appears to contain a ring of loosely packed rotor subunits and the proteins YAR027W and YAR028W. The structures suggest that Vma21p binds assembly intermediates that contain a rotor ring and that activation of proton pumping following assembly of V₁ with V_O removes Vma21p, allowing V-ATPase to remain in the Golgi. Together, these structures show how Vma12-22p and Vma21p function in V-ATPase assembly and quality control, ensuring the enzyme acidifies only its intended cellular targets.

Prediction of hydrophilic and hydrophobic hydration structure of protein by neural network optimized using experimental data

The hydration structures of proteins, which are necessary for their folding, stability, and functions, were visualized using X-ray and neutron crystallography and transmission electron microscopy. However, complete visualization of hydration structures over the entire protein surface remains difficult. To compensate for this incompleteness, we developed a three-dimensional convolutional neural network to predict the probability distribution of hydration water molecules on the hydrophilic and hydrophobic surfaces, and in the cavities of proteins. The neural network was optimized using the distribution patterns of protein atoms around the hydration water molecules identified in the high-resolution X-ray crystal structures. We examined the feasibility of the neural network using water sites in the protein crystal structures that were not included in the datasets. The predicted distribution covered most of the experimentally identified hydration sites, with local maxima appearing in their vicinity. This computational approach will help to highlight the relevance of hydration structures to the biological functions and dynamics of proteins.

Near-Atomic Resolution Cryo-EM Image Reconstruction of RNA

The rapid development of cryogenic electron microscopy (cryo-EM) enables the structure determination of macromolecules without the need for crystallization. Protein, protein-lipid, and protein-nucleic acid complexes can now be routinely resolved by cryo-EM single-particle analysis (SPA) to near-atomic or atomic resolution. Here we describe the structure determination of pure RNAs by SPA, from cryo-specimen preparation to data collection and 3D reconstruction. This protocol is useful to yield many cryo-EM structures of RNA, here exemplified by the Tetrahymena L-21 ScaI ribozyme at 3.1-Å resolution.

Bioelectric Fields at the Beginnings of Life

The consensus on the origins of life is that it involved organization of prebiotic chemicals according to the underlying principles of thermodynamics to dissipate energy derived from photochemical and/or geochemical sources. Leading theories tend to be chemistry-centric, revolving around either metabolism or information-containing polymers first. However, experimental data also suggest that bioelectricity and quantum effects play an important role in biology, which might suggest that a further factor is required to explain how life began. Intriguingly, in the early part of 20th century, the concept of the "morphogenetic field" was proposed by Gurwitsch to explain how the shape of an organism was determined, while a role for quantum mechanics in biology was suggested by Bohr and Schrödinger, among others. This raises the question as to the potential of these phenomena, especially bioelectric fields, to have been involved in the origin of life. It points to the possibility that as bioelectricity is universally prevalent in biological systems today, it represents a more complex echo of an electromagnetic skeleton which helped shape life into being. It could be argued that as a flow of ions creates an electric field, this could have been pivotal in the formation of an energy dissipating structure, for instance, in deep sea thermal vents. Moreover, a field theory might also hint at the potential involvement of nontrivial quantum effects in life. Not only might this perspective help indicate the origins of morphogenetic fields, but also perhaps suggest where life may have started, and whether metabolism or information came first. It might also help to provide an insight into aging, cancer, consciousness, and, perhaps, how we might identify life beyond our planet. In short, when thinking about life, not only do we have to consider the accepted chemistry, but also the fields that must also shape it. In effect, to fully understand life, as well as the yin of accepted particle-based chemistry, there is a yang of field-based interaction and an ethereal skeleton.

Direct observation of stepping rotation of V-ATPase reveals rigid component in coupling between Vo and V1 motors

V-ATPases are rotary motor proteins that convert the chemical energy of ATP into the electrochemical potential of ions across cell membranes. V-ATPases consist of two rotary motors, Vo and V1, and Enterococcus hirae V-ATPase (EhVoV1) actively transports Na+ in Vo (EhVo) by using torque generated by ATP hydrolysis in V1 (EhV1). Here, we observed ATP-driven stepping rotation of detergent-solubilized EhVoV1 wild-type, aE634A, and BR350K mutants under various Na+ and ATP concentrations ([Na+] and [ATP], respectively) by using a 40-nm gold nanoparticle as a low-load probe. When [Na+] was low and [ATP] was high, under the condition that only Na+ binding to EhVo is rate limiting, wild-type and aE634A exhibited 10 pausing positions reflecting 10-fold symmetry of the EhVo rotor and almost no backward steps. Duration time before the forward steps was inversely proportional to [Na+], confirming that Na+ binding triggers the steps. When both [ATP] and [Na+] were low, under the condition that both Na+ and ATP bindings are rate limiting, aE634A exhibited 13 pausing positions reflecting 10- and 3-fold symmetries of EhVo and EhV1, respectively. The distribution of duration time before the forward step was fitted well by the sum of two exponential decay functions with distinct time constants. Furthermore, occasional backward steps smaller than 36° were observed. Small backward steps were also observed during three long ATP cleavage pauses of BR350K. These results indicate that EhVo and EhV1 do not share pausing positions, Na+ and ATP bindings occur at different angles, and the coupling between EhVo and EhV1 has a rigid component.

F1FO ATP synthase molecular motor mechanisms

The F-ATP synthase, consisting of F1 and FO motors connected by a central rotor and the stators, is the enzyme responsible for synthesizing the majority of ATP in all organisms. The F1 (αβ)3 ring stator contains three catalytic sites. Single-molecule F1 rotation studies revealed that ATP hydrolysis at each catalytic site (0°) precedes a power-stroke that rotates subunit-γ 120° with angular velocities that vary with rotational position. Catalytic site conformations vary relative to subunit-γ position (βE, empty; βD, ADP bound; βT, ATP-bound). During a power stroke, βE binds ATP (0°–60°) and βD releases ADP (60°–120°). Årrhenius analysis of the power stroke revealed that elastic energy powers rotation via unwinding the γ-subunit coiled-coil. Energy from ATP binding at 34° closes βE upon subunit-γ to drive rotation to 120° and forcing the subunit-γ to exchange its tether from βE to βD, which changes catalytic site conformations. In F1FO, the membrane-bound FO complex contains a ring of c-subunits that is attached to subunit-γ. This c-ring rotates relative to the subunit-a stator in response to transmembrane proton flow driven by a pH gradient, which drives subunit-γ rotation in the opposite direction to force ATP synthesis in F1. Single-molecule studies of F1FO embedded in lipid bilayer nanodisks showed that the c-ring transiently stopped F1-ATPase-driven rotation every 36° (at each c-subunit in the c10-ring of E. coli F1FO) and was able to rotate 11° in the direction of ATP synthesis. Protonation and deprotonation of the conserved carboxyl group on each c-subunit is facilitated by separate groups of subunit-a residues, which were determined to have different pKa’s. Mutations of any of any residue from either group changed both pKa values, which changed the occurrence of the 11° rotation proportionately. This supports a Grotthuss mechanism for proton translocation and indicates that proton translocation occurs during the 11° steps. This is consistent with a mechanism in which each 36° of rotation the c-ring during ATP synthesis involves a proton translocation-dependent 11° rotation of the c-ring, followed by a 25° rotation driven by electrostatic interaction of the negatively charged unprotonated carboxyl group to the positively charged essential arginine in subunit-a.

The Plant V-ATPase

V-ATPase is the dominant proton pump in plant cells. It contributes to cytosolic pH homeostasis and energizes transport processes across endomembranes of the secretory pathway. Its localization in the trans Golgi network/early endosomes is essential for vesicle transport, for instance for the delivery of cell wall components. Furthermore, it is crucial for response to abiotic and biotic stresses. The V-ATPase's rather complex structure and multiple subunit isoforms enable high structural flexibility with respect to requirements for different organs, developmental stages, and organelles. This complexity further demands a sophisticated assembly machinery and transport routes in cells, a process that is still not fully understood. Regulation of V-ATPase is a target of phosphorylation and redox-modifications but also involves interactions with regulatory proteins like 14-3-3 proteins and the lipid environment. Regulation by reversible assembly, as reported for yeast and the mammalian enzyme, has not be proven in plants but seems to be absent in autotrophic cells. Addressing the regulation of V-ATPase is a promising approach to adjust its activity for improved stress resistance or higher crop yield.

Rotational Mechanism of FO Motor in the F-Type ATP Synthase Driven by the Proton Motive Force

In FOF1 ATP synthase, driven by the proton motive force across the membrane, the FO motor rotates the central rotor and induces conformational changes in the F1 motor, resulting in ATP synthesis. Recently, many near-atomic resolution structural models have been obtained using cryo-electron microscopy. Despite high resolution, however, static information alone cannot elucidate how and where the protons pass through the FO and how proton passage is coupled to FO rotation. Here, we review theoretical and computational studies based on FO structure models. All-atom molecular dynamics (MD) simulations elucidated changes in the protonation/deprotonation of glutamate-the protein-carrier residue-during rotation and revealed the protonation states that form the "water wire" required for long-range proton hopping. Coarse-grained MD simulations unveiled a free energy surface based on the protonation state and rotational angle of the rotor. Hybrid Monte Carlo and MD simulations showed how proton transfer is coupled to rotation.

Molecular basis of mEAK7-mediated human V-ATPase regulation

The activity of V-ATPase is well-known to be regulated by reversible dissociation of its V1 and Vo domains in response to growth factor stimulation, nutrient sensing, and cellular differentiation. The molecular basis of its regulation by an endogenous modulator without affecting V-ATPase assembly remains unclear. Here, we discover that a lysosome-anchored protein termed (mammalian Enhancer-of-Akt-1-7 (mEAK7)) binds to intact V-ATPase. We determine cryo-EM structure of human mEAK7 in complex with human V-ATPase in native lipid-containing nanodiscs. The structure reveals that the TLDc domain of mEAK7 engages with subunits A, B, and E, while its C-terminal domain binds to subunit D, presumably blocking V1-Vo torque transmission. Our functional studies suggest that mEAK7, which may act as a V-ATPase inhibitor, does not affect the activity of V-ATPase in vitro. However, overexpression of mEAK7 in HCT116 cells that stably express subunit a4 of V-ATPase represses the phosphorylation of ribosomal protein S6. Thus, this finding suggests that mEAK7 potentially links mTOR signaling with V-ATPase activity.

Structural insights into ClpP protease side exit pore-opening by a pH drop coupled with substrate hydrolysis

The ClpP serine peptidase is a tetradecameric degradation molecular machine involved in many physiological processes. It becomes a competent ATP-dependent protease when coupled with Clp-ATPases. Small chemical compounds, acyldepsipeptides (ADEPs), are known to cause the dysregulation and activation of ClpP without ATPases and have potential as novel antibiotics. Previously, structural studies of ClpP from various species revealed its structural details, conformational changes, and activation mechanism. Although product release through side exit pores has been proposed, the detailed driving force for product release remains elusive. Herein, we report crystal structures of ClpP from Bacillus subtilis (BsClpP) in unforeseen ADEP-bound states. Cryo-electron microscopy structures of BsClpP revealed various conformational states under different pH conditions. To understand the conformational change required for product release, we investigated the relationship between substrate hydrolysis and the pH-lowering process. The production of hydrolyzed peptides from acidic and basic substrates by proteinase K and BsClpP lowered the pH values. Our data, together with those of previous findings, provide insight into the molecular mechanism of product release by the ClpP self-compartmentalizing protease.

A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls

The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.

Oxidative stress protein Oxr1 promotes V-ATPase holoenzyme disassembly in catalytic activity-independent manner

The vacuolar ATPase (V-ATPase) is a rotary motor proton pump that is regulated by an assembly equilibrium between active holoenzyme and autoinhibited V1 -ATPase and Vo proton channel subcomplexes. Here, we report cryo-EM structures of yeast V-ATPase assembled in vitro from lipid nanodisc reconstituted Vo and mutant V1 . Our analysis identified holoenzymes in three active rotary states, indicating that binding of V1 to Vo provides sufficient free energy to overcome Vo autoinhibition. Moreover, the structures suggest that the unequal spacing of Vo 's proton-carrying glutamic acid residues serves to alleviate the symmetry mismatch between V1 and Vo motors, a notion that is supported by mutagenesis experiments. We also uncover a structure of free V1 bound to Oxr1, a conserved but poorly characterized factor involved in the oxidative stress response. Biochemical experiments show that Oxr1 inhibits V1 -ATPase and causes disassembly of the holoenzyme, suggesting that Oxr1 plays a direct role in V-ATPase regulation.

Determinants of Directionality and Efficiency of the ATP Synthase Fo Motor at Atomic Resolution

F_o subcomplex of ATP synthase is a membrane-embedded rotary motor that converts proton motive force into mechanical energy. Despite a rapid increase in the number of high-resolution structures, the mechanism of tight coupling between proton transport and motion of the rotary c-ring remains elusive. Here, using extensive all-atom free energy simulations, we show how the motor's directionality naturally arises from the interplay between intraprotein interactions and energetics of protonation of the c-ring. Notably, our calculations reveal that the strictly conserved arginine in the a-subunit (R176) serves as a jack-of-all-trades: it dictates the direction of rotation, controls the protonation state of the proton-release site, and separates the two proton-access half-channels. Therefore, arginine is necessary to avoid slippage between the proton flux and the mechanical output and guarantees highly efficient energy conversion. We also provide mechanistic explanations for the reported defective mutations of R176, reconciling the structural information on the F_o motor with previous functional and single-molecule data.

pH-dependent 11° F1FO ATP synthase sub-steps reveal insight into the FO torque generating mechanism

Most cellular ATP is made by rotary F₁F_O ATP synthases using proton translocation-generated clockwise torque on the F_O c-ring rotor, while F₁-ATP hydrolysis can force counterclockwise rotation and proton pumping. The F_O torque-generating mechanism remains elusive even though the F_O interface of stator subunit-a, which contains the transmembrane proton half-channels, and the c-ring is known from recent F₁F_O structures. Here, single-molecule F₁F_O rotation studies determined that the pKa values of the half-channels differ, show that mutations of residues in these channels change the pKa values of both half-channels, and reveal the ability of F_O to undergo single c-subunit rotational stepping. These experiments provide evidence to support the hypothesis that proton translocation through F_O operates via a Grotthuss mechanism involving a column of single water molecules in each half-channel linked by proton translocation-dependent c-ring rotation. We also observed pH-dependent 11° ATP synthase-direction sub-steps of the Escherichia coli c₁₀-ring of F₁F_O against the torque of F₁-ATPase-dependent rotation that result from H⁺ transfer events from F_O subunit-a groups with a low pKa to one c-subunit in the c-ring, and from an adjacent c-subunit to stator groups with a high pKa. These results support a mechanism in which alternating proton translocation-dependent 11° and 25° synthase-direction rotational sub-steps of the c₁₀-ring occur to sustain F₁F_O ATP synthesis.

Structure of the ATP synthase from Mycobacterium smegmatis provides targets for treating tuberculosis

The structure has been determined by electron cryomicroscopy of the adenosine triphosphate (ATP) synthase from Mycobacterium smegmatis This analysis confirms features in a prior description of the structure of the enzyme, but it also describes other highly significant attributes not recognized before that are crucial for understanding the mechanism and regulation of the mycobacterial enzyme. First, we resolved not only the three main states in the catalytic cycle described before but also eight substates that portray structural and mechanistic changes occurring during a 360° catalytic cycle. Second, a mechanism of auto-inhibition of ATP hydrolysis involves not only the engagement of the C-terminal region of an α-subunit in a loop in the γ-subunit, as proposed before, but also a "fail-safe" mechanism involving the b'-subunit in the peripheral stalk that enhances engagement. A third unreported characteristic is that the fused bδ-subunit contains a duplicated domain in its N-terminal region where the two copies of the domain participate in similar modes of attachment of the two of three N-terminal regions of the α-subunits. The auto-inhibitory plus the associated "fail-safe" mechanisms and the modes of attachment of the α-subunits provide targets for development of innovative antitubercular drugs. The structure also provides support for an observation made in the bovine ATP synthase that the transmembrane proton-motive force that provides the energy to drive the rotary mechanism is delivered directly and tangentially to the rotor via a Grotthuss water chain in a polar L-shaped tunnel.

Multiscale simulations of large complexes in conjunction with cryo-EM analysis

The cellular environment is highly crowded with most proteins and RNA/DNA forming homomeric and heteromeric complexes. Essential questions regarding how these complexes switch between functional, rest, and abnormal states with regulators or modifications remain challenging and complicated. Here, we review the recent progress integrating cryoelectron microscopy and multiscale molecular modeling to understand the dynamics and function-related mechanism in protein-RNA/DNA complexes, protein-protein complexes/assemblies, and membrane protein complexes. One future direction of multiscale simulations will be to interpret the large complex multibody regulation in assembly-induced function enhancement in conjunction with advanced atomic resolution structural-biology techniques and specialized computing architectures.

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.

Purification of active human vacuolar H+-ATPase in native lipid-containing nanodiscs

Vacuolar H+-ATPases (V-ATPases) are large, multisubunit proton pumps that acidify the lumen of organelles in virtually every eukaryotic cell and in specialized acid-secreting animal cells, the enzyme pumps protons into the extracellular space. In higher organisms, most of the subunits are expressed as multiple isoforms, with some enriched in specific compartments or tissues and others expressed ubiquitously. In mammals, subunit a is expressed as four isoforms (a1-4) that target the enzyme to distinct biological membranes. Mutations in a isoforms are known to give rise to tissue-specific disease, and some a isoforms are upregulated and mislocalized to the plasma membrane in invasive cancers. However, isoform complexity and low abundance greatly complicate purification of active human V-ATPase, a prerequisite for developing isoform-specific therapeutics. Here, we report the purification of an active human V-ATPase in native lipid nanodiscs from a cell line stably expressing affinity-tagged a isoform 4 (a4). We find that exogenous expression of this single subunit in HEK293F cells permits assembly of a functional V-ATPase by incorporation of endogenous subunits. The ATPase activity of the preparation is >95% sensitive to concanamycin A, indicating that the lipid nanodisc-reconstituted enzyme is functionally coupled. Moreover, this strategy permits purification of the enzyme's isolated membrane subcomplex together with biosynthetic assembly factors coiled-coil domain-containing protein 115, transmembrane protein 199, and vacuolar H+-ATPase assembly integral membrane protein 21. Our work thus lays the groundwork for biochemical characterization of active human V-ATPase in an a subunit isoform-specific manner and establishes a platform for the study of the assembly and regulation of the human holoenzyme.

The milk-derived lactoferrin inhibits V-ATPase activity by targeting its V1 domain

Lactoferrin (Lf), a bioactive milk protein, exhibits strong anticancer and antifungal activities. The search for Lf targets and mechanisms of action is of utmost importance to enhance its effective applications. A common feature among Lf-treated cancer and fungal cells is the inhibition of a proton pump called V-ATPase. Lf-driven V-ATPase inhibition leads to cytosolic acidification, ultimately causing cell death of cancer and fungal cells. Given that a detailed elucidation of how Lf and V-ATPase interact is still missing, herein we aimed to fill this gap by employing a five-stage computational approach. Molecular dynamics simulations of both proteins were performed to obtain a robust sampling of their conformational landscape, followed by clustering, which allowed retrieving representative structures, to then perform protein-protein docking. Subsequently, molecular dynamics simulations of the docked complexes and free binding energy calculations were carried out to evaluate the dynamic binding process and build a final ranking based on the binding affinities. Detailed atomist analysis of the top ranked complexes clearly indicates that Lf binds to the V₁ cytosolic domain of V-ATPase. Particularly, our data suggest that Lf binds to the interfaces between A/B subunits, where the ATP hydrolysis occurs, thus inhibiting this process. The free energy decomposition analysis further identified key binding residues that will certainly aid in the rational design of follow-up experimental studies, hence bridging computational and experimental biochemistry.

RAVE and Rabconnectin-3 Complexes as Signal Dependent Regulators of Organelle Acidification

The yeast RAVE (Regulator of H⁺-ATPase of Vacuolar and Endosomal membranes) complex and Rabconnectin-3 complexes of higher eukaryotes regulate acidification of organelles such as lysosomes and endosomes by catalyzing V-ATPase assembly. V-ATPases are highly conserved proton pumps consisting of a peripheral V₁ subcomplex that contains the sites of ATP hydrolysis, attached to an integral membrane V_o subcomplex that forms the transmembrane proton pore. Reversible disassembly of the V-ATPase is a conserved regulatory mechanism that occurs in response to multiple signals, serving to tune ATPase activity and compartment acidification to changing extracellular conditions. Signals such as glucose deprivation can induce release of V₁ from V_o, which inhibits both ATPase activity and proton transport. Reassembly of V₁ with V_o restores ATP-driven proton transport, but requires assistance of the RAVE or Rabconnectin-3 complexes. Glucose deprivation triggers V-ATPase disassembly in yeast and is accompanied by binding of RAVE to V₁ subcomplexes. Upon glucose readdition, RAVE catalyzes both recruitment of V₁ to the vacuolar membrane and its reassembly with V_o. The RAVE complex can be recruited to the vacuolar membrane by glucose in the absence of V₁ subunits, indicating that the interaction between RAVE and the V_o membrane domain is glucose-sensitive. Yeast RAVE complexes also distinguish between organelle-specific isoforms of the V_o a-subunit and thus regulate distinct V-ATPase subpopulations. Rabconnectin-3 complexes in higher eukaryotes appear to be functionally equivalent to yeast RAVE. Originally isolated as a two-subunit complex from rat brain, the Rabconnectin-3 complex has regions of homology with yeast RAVE and was shown to interact with V-ATPase subunits and promote endosomal acidification. Current understanding of the structure and function of RAVE and Rabconnectin-3 complexes, their interactions with the V-ATPase, their role in signal-dependent modulation of organelle acidification, and their impact on downstream pathways will be discussed.

Amphipathic environments for determining the structure of membrane proteins by single-particle electron cryo-microscopy

Over the past decade, the structural biology of membrane proteins (MPs) has taken a new turn thanks to epoch-making technical progress in single-particle electron cryo-microscopy (cryo-EM) as well as to improvements in sample preparation. The present analysis provides an overview of the extent and modes of usage of the various types of surfactants for cryo-EM studies. Digitonin, dodecylmaltoside, protein-based nanodiscs, lauryl maltoside-neopentyl glycol, glyco-diosgenin, and amphipols (APols) are the most popular surfactants at the vitrification step. Surfactant exchange is frequently used between MP purification and grid preparation, requiring extensive optimization each time the study of a new MP is undertaken. The variety of both the surfactants and experimental approaches used over the past few years bears witness to the need to continue developing innovative surfactants and optimizing conditions for sample preparation. The possibilities offered by novel APols for EM applications are discussed.

Molecular basis of V-ATPase inhibition by bafilomycin A1

Pharmacological inhibition of vacuolar-type H+-ATPase (V-ATPase) by its specific inhibitor can abrogate tumor metastasis, prevent autophagy, and reduce cellular signaling responses. Bafilomycin A1, a member of macrolide antibiotics and an autophagy inhibitor, serves as a specific and potent V-ATPases inhibitor. Although there are many V-ATPase structures reported, the molecular basis of specific inhibitors on V-ATPase remains unknown. Here, we report the cryo-EM structure of bafilomycin A1 bound intact bovine V-ATPase at an overall resolution of 3.6-Å. The structure reveals six bafilomycin A1 molecules bound to the c-ring. One bafilomycin A1 molecule engages with two c subunits and disrupts the interactions between the c-ring and subunit a, thereby preventing proton translocation. Structural and sequence analyses demonstrate that the bafilomycin A1-binding residues are conserved in yeast and mammalian species and the 7'-hydroxyl group of bafilomycin A1 acts as a unique feature recognized by subunit c.

Emerging insights on the role of V-ATPase in human diseases: Therapeutic challenges and opportunities

The control of the intracellular pH is vital for the survival of all organisms. Membrane transporters, both at the plasma and intracellular membranes, are key players in maintaining a finely tuned pH balance between intra- and extracellular spaces, and therefore in cellular homeostasis. V-ATPase is a housekeeping ATP-driven proton pump highly conserved among prokaryotes and eukaryotes. This proton pump, which exhibits a complex multisubunit structure based on cell type-specific isoforms, is essential for pH regulation and for a multitude of ubiquitous and specialized functions. Thus, it is not surprising that V-ATPase aberrant overexpression, mislocalization, and mutations in V-ATPase subunit-encoding genes have been associated with several human diseases. However, the ubiquitous expression of this transporter and the high toxicity driven by its off-target inhibition, renders V-ATPase-directed therapies very challenging and increases the need for selective strategies. Here we review emerging evidence linking V-ATPase and both inherited and acquired human diseases, explore the therapeutic challenges and opportunities envisaged from recent data, and advance future research avenues. We highlight the importance of V-ATPases with unique subunit isoform molecular signatures and disease-associated isoforms to design selective V-ATPase-directed therapies. We also discuss the rational design of drug development pipelines and cutting-edge methodological approaches toward V-ATPase-centered drug discovery. Diseases like cancer, osteoporosis, and even fungal infections can benefit from V-ATPase-directed therapies.

A "Sugar-Coated" Proton Pump Comes into Focus: High-Resolution Structure of a Human V-ATPase

In this issue of Molecular Cell, Wang et al. (2020a) present near-atomic resolution cryoEM structures of a proton-pumping vacuolar ATPase from human cells, illuminating the glycosylation of integral subunits and identifying a novel and conserved glycolipid ligand.