Search

ATP synthases play a crucial role in energy production by utilizing the proton motive force (pmf) across the membrane to rotate their membrane-embedded rotor c-ring, and thus driving ATP synthesis in the hydrophilic catalytic hexamer. However, the mechanism of how pmf converts into c-ring rotation remains unclear. This study presents a 2.8 Å cryo-EM structure of the Vo domain of V/A-ATPase from Thermus thermophilus, revealing precise orientations of glutamate (Glu) residues in the c12-ring. Three Glu residues face a water channel, with one forming a salt bridge with the Arginine in the stator (a/Arg). Molecular dynamics (MD) simulations show that protonation of specific Glu residues triggers unidirectional Brownian motion of the c12-ring towards ATP synthesis. When the key Glu remains unprotonated, the salt bridge persists, blocking rotation. These findings suggest that asymmetry in the protonation of c/Glu residues biases c12-ring movement, facilitating rotation and ATP synthesis. Here the authors determine a 2.8 Å cryo-EM structure of the Vo domain of V/A-ATPase, showing how protonation of Glu residues in the c12-ring drives ATP synthesis. Molecular dynamics simulations highlight that asymmetry in protonation of Glu residues biases c12-ring rotation. [ABSTRACT FROM AUTHOR]

A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing the product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases. [ABSTRACT FROM AUTHOR]

We propose a cross-entropy minimization method for finding the reaction coordinate from a large number of collective variables in complex molecular systems. This method is an extension of the likelihood maximization approach describing the committor function with a sigmoid. By design, the reaction coordinate as a function of various collective variables is optimized such that the distribution of the committor p B * values generated from molecular dynamics simulations can be described in a sigmoidal manner. We also introduce the L2-norm regularization used in the machine learning field to prevent overfitting when the number of considered collective variables is large. The current method is applied to study the isomerization of alanine dipeptide in vacuum, where 45 dihedral angles are used as candidate variables. The regularization parameter is determined by cross-validation using training and test datasets. It is demonstrated that the optimal reaction coordinate involves important dihedral angles, which are consistent with the previously reported results. Furthermore, the points with p B * ∼ 0.5 clearly indicate a separatrix distinguishing reactant and product states on the potential of mean force using the extracted dihedral angles. [ABSTRACT FROM AUTHOR]

Multidomain proteins can exhibit sophisticated functions based on cooperative interactions and allosteric regulation through spatial rearrangements of the multiple domains. This study explored the potential of using multidomain proteins as a basis for Förster resonance energy transfer (FRET) biosensors, focusing on protein disulfide isomerase (PDI) as a representative example. PDI, a well-studied multidomain protein, undergoes redox-dependent conformational changes, enabling the exposure of a hydrophobic surface extending across the b' and a' domains that serves as the primary binding site for substrates. Taking advantage of the dynamic domain rearrangements of PDI, we developed FRET-based biosensors by fusing the b' and a' domains of thermophilic fungal PDI with fluorescent proteins as the FRET acceptor and donor, respectively. Both experimental and computational approaches were used to characterize FRET efficiency in different redox states. In vitro and in vivo evaluations demonstrated higher FRET efficiency of this biosensor in the oxidized form, reflecting the domain rearrangement and its responsiveness to intracellular redox environments. This novel approach of exploiting redox-dependent domain dynamics in multidomain proteins offers promising opportunities for designing innovative FRET-based biosensors with potential applications in studying cellular redox regulation and beyond. [ABSTRACT FROM AUTHOR]

We combine molecular simulations and mechanical modeling to explore the mechanism of energy conversion in the coupled rotary motors of FoF1-ATP synthase. A torsional viscoelastic model with frictional dissipation quantitatively reproduces the dynamics and energetics seen in atomistic molecular dynamics simulations of torque-driven γ-subunit rotation in the F1-ATPase rotary motor. The torsional elastic coefficients determined from the simulations agree with results from independent single-molecule experiments probing different segments of the γ-subunit, which resolves a long-lasting controversy. At steady rotational speeds of ~1 kHz corresponding to experimental turnover, the calculated frictional dissipation of less than kBT per rotation is consistent with the high thermodynamic efficiency of the fully reversible motor. Without load, the maximum rotational speed during transitions between dwells is reached at ~1 MHz. Energetic constraints dictate a unique pathway for the coupled rotations of the Fo and F1 rotary motors in ATP synthase, and explain the need for the finer stepping of the F1 motor in the mammalian system, as seen in recent experiments. Compensating for incommensurate eightfold and threefold rotational symmetries in Fo and F1, respectively, a significant fraction of the external mechanical work is transiently stored as elastic energy in the γ-subunit. The general framework developed here should be applicable to other molecular machines. [ABSTRACT FROM AUTHOR]

F1-ATPase, the catalytic domain of ATP synthase, synthesizes most of the ATP in living organisms. Running in reverse powered by ATP hydrolysis, this hexameric ring-shaped molecular motor formed by three αβ-dimers creates torque on its central γ-subunit. This reverse operation enables detailed explorations of the mechanochemical coupling mechanisms in experiment and simulation. Here, we use molecular dynamics simulations to construct a first atomistic conformation of the intermediate state following the 40° substep of rotary motion, and to study the timing and molecular mechanism of inorganic phosphate (Pi) release coupled to the rotation. In response to torque-driven rotation of the γ-subunit in the hydrolysis direction, the nucleotide-free αβE interface forming the "empty" E site loosens and singly charged Pi readily escapes to the P loop. By contrast, the interface stays closed with doubly charged Pi. The γ-rotation tightens the ATP-bound αβTP interface, as required for hydrolysis. The calculated rate for the outward release of doubly charged Pi from the αβE interface 120° after ATP hydrolysis closely matches the ~1-ms functional timescale. Conversely, Pi release from the ADP-bound αβDP interface postulated in earlier models would occur through a kinetically infeasible inward-directed pathway. Our simulations help reconcile conflicting interpretations of single-molecule experiments and crystallographic studies by clarifying the timing of Pi exit, its pathway and kinetics, associated changes in Pi protonation, and changes of the F1-ATPase structure in the 40° substep. Important elements of the molecular mechanism of Pi release emerging from our simulations appear to be conserved in myosin despite the different functional motions. [ABSTRACT FROM AUTHOR]

Association of protein molecules constitutes the basis for the interaction network in a cell. Despite its fundamental importance, the thermodynamic aspect of protein-protein binding, particularly the issues relating to the entropy change upon binding, remains elusive, The binding of actin and myosin, which are vital proteins in motility, is a typical example, in which two different binding mechanisms have been argued: the binding affinity increases with increasing temperature and with decreasing salt-concentration, indicating the entropy-driven binding and the enthalpy-driven binding, respectively. How can these thermodynamically different binding mechanisms coexist? To address this question, which is of general importance in understanding protein-protein bindings, we conducted an in silico titration of the actin-myosin system by molecular dynamics simulation using a residue-level coarse-grained model, with particular focus on the role of the electrostatic interaction. We found a good agreement between in silico and in vitro experiments on the salt-concentration dependence and the temperature dependence of the binding aftinity. We then figured out how the two binding mechanisms can coexist: the enthalpy (due to electrostatic interaction between actin and myosin) provides the basal binding affinity, and the entropy (due to the orientational disorder of water molecules) enhances it at higher temperatures. In addition, we analyzed the actin-myosin complex structures observed during the simulation and obtained a variety of weak-binding complex structures, among which were found an unusual binding mode suggested by an earlier experiment and precursor structures of the strong-binding complex proposed by electron microscopy. These results collectively indicate the potential capability of a residue-level coarse-grained model to simulate the association-dissociation dynamics (particularly for transient weak-bindings) exhibited by larger and more complicated systems, as in a cell. [ABSTRACT FROM AUTHOR]

Summary: F1-ATPase, a rotary motor powered by adenosine triphosphate hydrolysis, has been extensively studied by various methods. Here, we performed a systematic comparison of 29 X-ray crystal structures of F1-complexes, finding fine interplay among enzyme structures, catalysis, and rotations. First, analyzing the 87 structures of enzymatic αβ-subunits, we confirmed that the two modes, the hinge motion of β-subunit and the loose/tight motion of the αβ-interface, dominate the variations. The structural ensemble was nearly contiguous bridging three clusters, αβTP, αβDP, and αβE. Second, the catalytic site analysis suggested the correlation between the phosphate binding and the tightening of the αβ-interface. Third, addressing correlations of enzymatic structures with the orientations of the central stalk γ, we found that the γ rotation highly correlates with loosening of αβE-interface and βDP hinge motions. Finally, calculating the helix 6 angle of β, we identified the recently observed partially closed conformation being consistent with βHC. [ABSTRACT FROM AUTHOR]