Folding-based molecular simulations reveal mechanisms of the rotary motor F1-ATPase

Switching Go̅-Martini for Investigating Protein Conformational Transitions and Associated Protein-Lipid Interactions

Proteins are dynamic biomolecules that can transform between different conformational states when exerting physiological functions, which is difficult to simulate using all-atom methods. Coarse-grained (CG) Go̅-like models are widely used to investigate large-scale conformational transitions, which usually adopt implicit solvent models and therefore cannot explicitly capture the interaction between proteins and surrounding molecules, such as water and lipid molecules. Here, we present a new method, named Switching Go̅-Martini, to simulate large-scale protein conformational transitions between different states, based on the switching Go̅ method and the CG Martini 3 force field. The method is straightforward and efficient, as demonstrated by the benchmarking applications for multiple protein systems, including glutamine binding protein (GlnBP), adenylate kinase (AdK), and β2-adrenergic receptor (β2AR). Moreover, by employing the Switching Go̅-Martini method, we can not only unveil the conformational transition from the E2Pi-PL state to E1 state of the type 4 P-type ATPase (P4-ATPase) flippase ATP8A1-CDC50 but also provide insights into the intricate details of lipid transport.

Simulating the conformational dynamics of the ATPase complex on proteasome using its free-energy landscape

The AAA+ ATPase complex on proteasome powers its functions through a series of intricate conformational transitions. Here, we describe a procedure to simulate the conformational dynamics of the proteasomal ATPase complex. We first empirically determined the free-energy landscape (FEL) of proteasome and then simulated proteasome's conformational changes as stochastic transitions on its FEL. We compared the FEL-predicted proteasomal behaviors with experimental measurements and analyzed the map of the ATPase's global dynamics to gain mechanistic insights into proteasomal degradation. For complete details on the use and execution of this protocol, please refer to Fang et al. (2022).1.

Physical pictures of rotation mechanisms of F1- and V1-ATPases: Leading roles of translational, configurational entropy of water

We aim to develop a theory based on a concept other than the chemo-mechanical coupling (transduction of chemical free energy of ATP to mechanical work) for an ATP-driven protein complex. Experimental results conflicting with the chemo-mechanical coupling have recently emerged. We claim that the system comprises not only the protein complex but also the aqueous solution in which the protein complex is immersed and the system performs essentially no mechanical work. We perform statistical-mechanical analyses on V₁-ATPase (the A₃B₃DF complex) for which crystal structures in more different states are experimentally known than for F₁-ATPase (the α₃β₃γ complex). Molecular and atomistic models are employed for water and the structure of V₁-ATPase, respectively. The entropy originating from the translational displacement of water molecules in the system is treated as a pivotal factor. We find that the packing structure of the catalytic dwell state of V₁-ATPase is constructed by the interplay of ATP bindings to two of the A subunits and incorporation of the DF subunit. The packing structure represents the nonuniformity with respect to the closeness of packing of the atoms in constituent proteins and protein interfaces. The physical picture of rotation mechanism of F₁-ATPase recently constructed by Kinoshita is examined, and common points and differences between F₁- and V₁-ATPases are revealed. An ATP hydrolysis cycle comprises binding of ATP to the protein complex, hydrolysis of ATP into ADP and Pi in it, and dissociation of ADP and Pi from it. During each cycle, the chemical compounds bound to the three A or β subunits and the packing structure of the A₃B₃ or α₃β₃ complex are sequentially changed, which induces the unidirectional rotation of the central shaft for retaining the packing structure of the A₃B₃DF or α₃β₃γ complex stabilized for almost maximizing the water entropy. The torque driving the rotation is generated by water with no input of chemical free energy. The presence of ATP is indispensable as a trigger of the torque generation. The ATP hydrolysis or synthesis reaction is tightly coupled to the rotation of the central shaft in the normal or inverse direction through the water-entropy effect.

Landscape-Based View on the Stepping Movement of Myosin VI

Myosin VI dimer walks toward the minus end of the actin filament with a large and variable step size of 25–36 nm. Two competing models have been put forward to explain this large step size. The Spudich model assumes that the myosin VI dimer associates at a distal tail near the cargo-binding domain, which makes two full-length single α-helix (SAH) domains serve as long legs. In contrast, the Houdusse–Sweeney model assumes that the association occurs in the middle (between residues 913 and 940) of the SAH domain and that the three-helix bundles unfold to ensure the large step size. Their consistency with the observation of stepping motion with a large and variable step size has not been examined in detail. To compare the two proposed models of myosin VI, we computationally characterized the free energy landscape experienced by the leading head during the stepping movement along the actin filament using the elastic network model of two heads and an implicit model of the SAH domains. Our results showed that the Spudich model is more consistent with the 25–36 nm step size than the Houdusse–Sweeney model. The unfolding of the three-helix bundles gives rise to the free energy bias toward a shorter distance between two heads. Besides, the stiffness of the SAH domain is a key factor for giving strong energetic bias toward the longer distance of stepping. Free energy analysis of the stepping motion complements the visual inspection of static structures and enables a deeper understanding of underlying mechanisms of molecular motors.

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.

Simulating a chemically fueled molecular motor with nonequilibrium molecular dynamics

Most computer simulations of molecular dynamics take place under equilibrium conditions-in a closed, isolated system, or perhaps one held at constant temperature or pressure. Sometimes, extra tensions, shears, or temperature gradients are introduced to those simulations to probe one type of nonequilibrium response to external forces. Catalysts and molecular motors, however, function based on the nonequilibrium dynamics induced by a chemical reaction's thermodynamic driving force. In this scenario, simulations require chemostats capable of preserving the chemical concentrations of the nonequilibrium steady state. We develop such a dynamic scheme and use it to observe cycles of a particle-based classical model of a catenane-like molecular motor. Molecular motors are frequently modeled with detailed-balance-breaking Markov models, and we explicitly construct such a picture by coarse graining the microscopic dynamics of our simulations in order to extract rates. This work identifies inter-particle interactions that tune those rates to create a functional motor, thereby yielding a computational playground to investigate the interplay between directional bias, current generation, and coupling strength in molecular information ratchets.

Multiscale Coarse-Grained Model for the Stepping of Molecular Motors with Application to Kinesin

Conventional kinesin, a motor protein that transports cargo within cells, walks by taking multiple steps toward the plus end of the microtubule (MT). While significant progress has been made in understanding the details of the walking mechanism of kinesin, there are many unresolved issues. From a computational perspective, a central challenge is the large size of the system, which limits the scope of time scales accessible in standard computer simulations. Here, we create a general multiscale coarse-grained model for motors that enables us to simulate the stepping process of motors on polar tracks (actin and MT) with a focus on kinesin. Our approach greatly shortens the computation times without a significant loss in detail, thus allowing us to better describe the molecular basis of the stepping kinetics. The small number of parameters, which are determined by fitting to experimental data, allows us to develop an accurate method that may be adopted to simulate stepping in other molecular motors. The model enables us to simulate a large number of steps, which was not possible previously. We show in agreement with experiments that due to the docking of the neck linker (NL) of kinesin, sometimes deemed as the power stroke, the space explored diffusively by the tethered head is severely restricted, allowing the step to be completed in tens of microseconds. We predict that increasing the interaction strength between the NL and the motor head, achievable by mutations in the NL, decreases the stepping time but reaches a saturation value. Furthermore, the full three-dimensional dynamics of the cargo are fully resolved in our model, contributing to the predictive power and allowing us to study the important aspects of cargo-motor interactions.

Simulating the directional translocation of a substrate by the AAA+ motor in the 26S proteasome

This work explored the molecular origin of substrate translocation by the AAA+ motor of the 26S proteasome. This exploration was performed by combining different simulation approaches including calculations of binding free energies, coarse-grained simulations, and considerations of the ATP hydrolysis energy. The simulations were used to construct the free energy landscape for the translocation process. This included the evaluation of the conformational barriers in different translocation steps. Our simulation reveals that the substrate translocation by the AAA+ motor is guided in part by electrostatic interactions. We also validated the experimental observation that bulkier residues in pore loop 1 are responsible for substrate translocation. However, our calculation also reveals that the lysine residues prior to the bulkier residues (conserved along pore loop 1) are also important for the translocation process. We believe that this computational study can help in guiding the ongoing research of the proteasome.

Conformational state switching and pathways of chromosome dynamics in cell cycle

The cell cycle is a process and function of a cell with different phases essential for cell growth, proliferation, and replication. It depends on the structure and dynamics of the underlying DNA molecule, which underpins the genome function. A microscopic structural-level understanding of how a genome or its functional module chromosome performs the cell cycle in terms of large-scale conformational transformation between different phases, such as the interphase and the mitotic phase, is still challenging. Here, we develop a non-equilibrium, excitation-relaxation energy landscape-switching model to quantify the underlying chromosome conformational transitions through (de-)condensation for a complete microscopic understanding of the cell cycle. We show that the chromosome conformational transition mechanism from the interphase to the mitotic phase follows a two-stage scenario, in good agreement with the experiments. In contrast, the mitotic exit pathways show the existence of an over-expanded chromosome that recapitulates the chromosome in the experimentally identified intermediate state at the telophase. We find the conformational pathways are heterogeneous and irreversible as a result of the non-equilibrium dynamics of the cell cycle from both structural and kinetic perspectives. We suggest that the irreversibility is mainly due to the distinct participation of the ATP-dependent structural maintenance of chromosomal protein complexes during the cell cycle. Our findings provide crucial insights into the microscopic molecular structural and dynamical physical mechanism for the cell cycle beyond the previous more macroscopic descriptions. Our non-equilibrium landscape framework is general and applicable to study diverse non-equilibrium physical and biological processes such as active matter, differentiation/development, and cancer.

Exploring Successful Parameter Region for Coarse-Grained Simulation of Biomolecules by Bayesian Optimization and Active Learning

Accompanied with an increase of revealed biomolecular structures owing to advancements in structural biology, the molecular dynamics (MD) approach, especially coarse-grained (CG) MD suitable for macromolecules, is becoming increasingly important for elucidating their dynamics and behavior. In fact, CG-MD simulation has succeeded in qualitatively reproducing numerous biological processes for various biomolecules such as conformational changes and protein folding with reasonable calculation costs. However, CG-MD simulations strongly depend on various parameters, and selecting an appropriate parameter set is necessary to reproduce a particular biological process. Because exhaustive examination of all candidate parameters is inefficient, it is important to identify successful parameters. Furthermore, the successful region, in which the desired process is reproducible, is essential for describing the detailed mechanics of functional processes and environmental sensitivity and robustness. We propose an efficient search method for identifying the successful region by using two machine learning techniques, Bayesian optimization and active learning. We evaluated its performance using F1-ATPase, a biological rotary motor, with CG-MD simulations. We successfully identified the successful region with lower computational costs (12.3% in the best case) without sacrificing accuracy compared to exhaustive search. This method can accelerate not only parameter search but also biological discussion of the detailed mechanics of functional processes and environmental sensitivity based on MD simulation studies.

Flexible Fitting of Biomolecular Structures to Atomic Force Microscopy Images via Biased Molecular Simulations

High-speed (HS) atomic force microscopy (AFM) is a prominent imaging technology that observes large-scale structural dynamics of biomolecules near the physiological condition, but the AFM data are limited to the surface shape of specimens. Rigid-body fitting methods were developed to obtain molecular structures that fit to an AFM image, without accounting for conformational changes. Here, we developed a method to fit flexibly a three-dimensional (3D) biomolecular structure into an AFM image. First, we describe a method to produce a pseudo-AFM image from a given 3D structure in a differentiable form. Then, using a correlation function between the experimental AFM image and the computational pseudo-AFM image, we developed a flexible fitting molecular dynamics (MD) simulation method by which we obtain protein structures that well fit to the given AFM image. We first test it with a twin experiment; using an AFM image produced from a protein structure different from its native conformation as a reference, we performed the flexible fitting MD simulations to sample conformations that fit well the reference AFM image, and the method was confirmed to work well. Then, parameter dependence in the protocol was discussed. Finally, we applied the method to a real experimental HS-AFM image for a flagellar protein FlhA, demonstrating its applicability. We also test the rigid-body fitting of a molecular structure to an AFM image. Our method will be a general tool for dynamic structure modeling based on HS-AFM images and is publicly available through the CafeMol software.

Building a macro-mixing dual-basin Gō model using the Multistate Bennett Acceptance Ratio

The dual-basin Gō-model is a structural-based coarsegrained model for simulating a conformational transition between two known structures of a protein. Two parameters are required to produce a dual-basin potential mixed using two single-basin potentials, although the determination of mixing parameters is usually not straightforward. Here, we have developed an efficient scheme to determine the mixing parameters using the Multistate Bennett Acceptance Ratio (MBAR) method after short simulations with a set of parameters. In the scheme, MBAR allows us to predict observables at various unsimulated conditions, which are useful to improve the mixing parameters in the next round of iterative simulations. The number of iterations that are necessary for obtaining the converged mixing parameters are significantly reduced in the scheme. We applied the scheme to two proteins, the glutamine binding protein and the ribose binding protein, for showing the effectiveness in the parameter determination. After obtaining the converged parameters, both proteins show frequent conformational transitions between open and closed states, providing the theoretical basis to investigate structure-dynamics-function relationships of the proteins.

Gō model revisited

This review discusses Gō models broadly used in biomolecular simulations. I start with a brief description of the original lattice model study by Nobuhiro Gō. Then, the theory of protein folding behind Gō model, free energy approaches, and off-lattice Gō models are reviewed. I also mention a stringent test for the assumption in Gō models given from all-atom molecular dynamics simulations. Subsequently, I move to application of Gō models to protein dynamical functions. Various extension of Gō models is also reviewed. Finally, some publicly available tools to use Gō models are listed.

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.

Is F1-ATPase a Rotary Motor with Nearly 100% Efficiency? Quantitative Analysis of Chemomechanical Coupling and Mechanical Slip

We present a chemomechanical network model of the rotary molecular motor F₁-ATPase which quantitatively describes not only the rotary motor dynamics driven by ATP hydrolysis but also the ATP synthesis caused by forced reverse rotations. We observe a high reversibility of F₁-ATPase, that is, the main cycle of ATP synthesis corresponds to the reversal of the main cycle in the hydrolysis-driven motor rotation. However, our quantitative analysis indicates that torque-induced mechanical slip without chemomechanical coupling occurs under high external torque and reduces the maximal efficiency of the free energy transduction to 40-80% below the optimal efficiency. Heat irreversibly dissipates not only through the viscous friction of the probe but also directly from the motor due to torque-induced mechanical slip. Such irreversible heat dissipation is a crucial limitation for achieving a 100% free-energy transduction efficiency with biological nanomachines because biomolecules are easily deformed by external torque.

Chromatin remodelers couple inchworm motion with twist-defect formation to slide nucleosomal DNA

ATP-dependent chromatin remodelers are molecular machines that control genome organization by repositioning, ejecting, or editing nucleosomes, activities that confer them essential regulatory roles on gene expression and DNA replication. Here, we investigate the molecular mechanism of active nucleosome sliding by means of molecular dynamics simulations of the Snf2 remodeler translocase in complex with a nucleosome. During its inchworm motion driven by ATP consumption, the translocase overwrites the original nucleosome energy landscape via steric and electrostatic interactions to induce sliding of nucleosomal DNA unidirectionally. The sliding is initiated at the remodeler binding location via the generation of a pair of twist defects, which then spontaneously propagate to complete sliding throughout the entire nucleosome. We also reveal how remodeler mutations and DNA sequence control active nucleosome repositioning, explaining several past experimental observations. These results offer a detailed mechanistic picture of remodeling important for the complete understanding of these key biological processes.

Nucleosomes are the protein-DNA complexes underlying Eukaryotic genome organization, and serve as regulators of gene expression by occluding DNA to other proteins. This regulation requires the precise positioning of nucleosomes along DNA. Chromatin remodelers are the molecular machines that consume ATP to slide nucleosome at their correct locations, but the mechanisms of remodeling are still unclear. Based on the static structural information of a remodeler bound on nucleosome, we performed molecular dynamics computer simulations revealing the details of how remodelers slide nucleosomal DNA: the inchworm-like motion of remodelers create small DNA deformations called twist defects, which then spontaneously propagate throughout the nucleosome to induce sliding. These simulations explain several past experimental findings and are important for our understanding of genome organization.

Evidence for a Partially Stalled γ Rotor in F1-ATPase from Hydrogen-Deuterium Exchange Experiments and Molecular Dynamics Simulations

F₁-ATPase uses ATP hydrolysis to drive rotation of the γ subunit. The γ C-terminal helix constitutes the rotor tip that is seated in an apical bearing formed by α₃β₃. It remains uncertain to what extent the γ conformation during rotation differs from that seen in rigid crystal structures. Existing models assume that the entire γ subunit participates in every rotation. Here we interrogated E. coli F₁-ATPase by hydrogen-deuterium exchange (HDX) mass spectrometry. Rotation of γ caused greatly enhanced deuteration in the γ C-terminal helix. The HDX kinetics implied that most F₁ complexes operate with an intact rotor at any given time, but that the rotor tip is prone to occasional unfolding. A molecular dynamics (MD) strategy was developed to model the off-axis forces acting on γ. MD runs showed stalling of the rotor tip and unfolding of the γ C-terminal helix. MD-predicted H-bond opening events coincided with experimental HDX patterns. Our data suggest that in vitro operation of F₁-ATPase is associated with significant rotational resistance in the apical bearing. These conditions cause the γ C-terminal helix to get "stuck" (and unfold) sporadically while the remainder of γ continues to rotate. This scenario contrasts the traditional "greasy bearing" model that envisions smooth rotation of the γ C-terminal helix. The fragility of the apical rotor tip in F₁-ATPase is attributed to the absence of a c₁₀ ring that stabilizes the rotation axis in intact F_oF₁. Overall, the MD/HDX strategy introduced here appears well suited for interrogating the inner workings of molecular motors.

Deciphering Intrinsic Inter-subunit Couplings that Lead to Sequential Hydrolysis of F1-ATPase Ring

Rotary sequential hydrolysis of the metabolic machine F₁-ATPase is a prominent manifestation of high coordination among multiple chemical sites in ring-shaped molecular machines, and it is also functionally essential for F₁ to tightly couple chemical reactions and central γ-shaft rotation. High-speed AFM experiments have identified that sequential hydrolysis is maintained in the F₁ stator ring even in the absence of the γ-rotor. To explore the origins of intrinsic sequential performance, we computationally investigated essential inter-subunit couplings on the hexameric ring of mitochondrial and bacterial F₁. We first reproduced in stochastic Monte Carlo simulations the experimentally determined sequential hydrolysis schemes by kinetically imposing inter-subunit couplings and following subsequent tri-site ATP hydrolysis cycles on the F₁ ring. We found that the key couplings to support the sequential hydrolysis are those that accelerate neighbor-site ADP and Pi release upon a certain ATP binding or hydrolysis reaction. The kinetically identified couplings were then examined in atomistic molecular dynamics simulations at a coarse-grained level to reveal the underlying structural mechanisms. To do that, we enforced targeted conformational changes of ATP binding or hydrolysis to one chemical site on the F₁ ring and monitored the ensuing conformational responses of the neighboring sites using structure-based simulations. Notably, we found asymmetrical neighbor-site opening that facilitates ADP release upon enforced ATP binding. We also captured a complete charge-hopping process of the Pi release subsequent to enforced ATP hydrolysis in the neighbor site, confirming recent single-molecule analyses with regard to the role of ATP hydrolysis in F₁. Our studies therefore elucidate both the coordinated chemical kinetics and structural dynamics mechanisms underpinning the sequential operation of the F₁ ring.

Allosteric conformational change cascade in cytoplasmic dynein revealed by structure-based molecular simulations

Cytoplasmic dynein is a giant ATP-driven molecular motor that proceeds to the minus end of the microtubule (MT). Dynein hydrolyzes ATP in a ring-like structure, containing 6 AAA+ (ATPases associated with diverse cellular activities) modules, which is ~15 nm away from the MT binding domain (MTBD). This architecture implies that long-distance allosteric couplings exist between the AAA+ ring and the MTBD in order for dynein to move on the MT, although little is known about the mechanisms involved. Here, we have performed comprehensive molecular simulations of the dynein motor domain based on pre- and post- power-stroke structural information and in doing so we address the allosteric conformational changes that occur during the power-stroke and recovery-stroke processes. In the power-stroke process, the N-terminal linker movement was the prerequisite to the nucleotide-dependent AAA1 transition, from which a transition cascade propagated, on average, in a circular manner on the AAA+ ring until it reached the AAA6/C-terminal module. The recovery-stroke process was initiated by the transition of the AAA6/C-terminal, from which the transition cascade split into the two directions of the AAA+ ring, occurring both clockwise and anti-clockwise. In both processes, the MTBD conformational change was regulated by the AAA4 module and the AAA5/Strut module.

The linear molecular motor dynein is an intriguing allosteric model protein. ATP hydrolysis, catalyzed by modules in the AAA+ ring, regulates the binding to the rail molecule, microtubule, which is ~15 nm away from the AAA+ ring. The molecular mechanisms underpinning this long-distance communication are unclear. Based on recently solved pre- and post- power-stroke crystal structure information, we performed, for the first time to our knowledge, molecular simulations of complete conformational changes between the two structures. The simulation revealed that module-by-module allosteric conformational changes occur. Interestingly, the transition cascade from the pre- to the post-power-stroke states propagated in a circular manner around the AAA+ ring, while that of the recovery transitions propagated in a bi-directional manner around the ring.

Role of non-native electrostatic interactions in the coupled folding and binding of PUMA with Mcl-1

PUMA, which belongs to the BH3-only protein family, is an intrinsically disordered protein (IDP). It binds to its cellular partner Mcl-1 through its BH3 motif, which folds upon binding into an α helix. We have applied a structure-based coarse-grained model, with an explicit Debye-Hückel charge model, to probe the importance of electrostatic interactions both in the early and the later stages of this model coupled folding and binding process. This model was carefully calibrated with the experimental data on helical content and affinity, and shown to be consistent with previously published experimental data on binding rate changes with respect to ionic strength. We find that intramolecular electrostatic interactions influence the unbound states of PUMA only marginally. Our results further suggest that intermolecular electrostatic interactions, and in particular non-native electrostatic interactions, are involved in formation of the initial encounter complex. We are able to reveal the binding mechanism in more detail than is possible using experimental data alone however, and in particular we uncover the role of non-native electrostatic interactions. We highlight the potential importance of such electrostatic interactions for describing the binding reactions of IDPs. Such approaches could be used to provide predictions for the results of mutational studies.

Geometric somersaults of a polymer chain through cyclic twisting motions

This study explores the significance of geometric angle shifts, which we call geometric somersaults, arising from cyclic twisting motions of a polymer chain. A five-bead polymer chain serves as a concise and minimal model of a molecular shaft throughout this study. We first show that this polymer chain can change its orientation about its longitudinal axis largely, e.g., 120^{∘}, under conditions of zero total angular momentum by changing the two dihedral angles in a cyclic manner. This phenomenon is an example of the so-called "falling cat" phenomenon, where a falling cat undergoes a geometric somersault by changing its body shape under conditions of zero total angular momentum. We then extend the geometric somersault of the polymer chain to a noisy and viscous environment, where the polymer chain is steered by external driving forces. This extension shows that the polymer chain can achieve an orientation change keeping its total angular momentum and total external torque fluctuating around zero in a noisy and viscous environment. As an application, we argue that the geometric somersault of the polymer chain by 120^{∘} may serve as a prototypical and coarse-grained model for the rotary motion of the central shaft of ATP synthase (F_{O}F_{1}-ATPase). This geometric somersault is in clear contrast to the standard picture for the rotary motion of the central shaft as a rigid body, which generally incurs nonzero total angular momentum and nonzero total external torque. The power profile of the geometric somersault implies a preliminary mechanism for elastic power transmission. The results of this study may be of fundamental interest in twisting and rotary motions of biomolecules.

Rotation of artificial rotor axles in rotary molecular motors

F₁- and V₁-ATPase are rotary molecular motors that convert chemical energy released upon ATP hydrolysis into torque to rotate a central rotor axle against the surrounding catalytic stator cylinder with high efficiency. How conformational change occurring in the stator is coupled to the rotary motion of the axle is the key unknown in the mechanism of rotary motors. Here, we generated chimeric motor proteins by inserting an exogenous rod protein, FliJ, into the stator ring of F₁ or of V₁ and tested the rotation properties of these chimeric motors. Both motors showed unidirectional and continuous rotation, despite no obvious homology in amino acid sequence between FliJ and the intrinsic rotor subunit of F₁ or V₁ These results showed that any residue-specific interactions between the stator and rotor are not a prerequisite for unidirectional rotation of both F₁ and V₁ The torque of chimeric motors estimated from viscous friction of the rotation probe against medium revealed that whereas the F₁-FliJ chimera generates only 10% of WT F₁, the V₁-FliJ chimera generates torque comparable to that of V₁ with the native axle protein that is structurally more similar to FliJ than the native rotor of F₁ This suggests that the gross structural mismatch hinders smooth rotation of FliJ accompanied with the stator ring of F₁.

F1-ATPase conformational cycle from simultaneous single-molecule FRET and rotation measurements

Despite extensive studies, the structural basis for the mechanochemical coupling in the rotary molecular motor F1-ATPase (F1) is still incomplete. We performed single-molecule FRET measurements to monitor conformational changes in the stator ring-α3β3, while simultaneously monitoring rotations of the central shaft-γ. In the ATP waiting dwell, two of three β-subunits simultaneously adopt low FRET nonclosed forms. By contrast, in the catalytic intermediate dwell, two β-subunits are simultaneously in a high FRET closed form. These differences allow us to assign crystal structures directly to both major dwell states, thus resolving a long-standing issue and establishing a firm connection between F1 structure and the rotation angle of the motor. Remarkably, a structure of F1 in an ε-inhibited state is consistent with the unique FRET signature of the ATP waiting dwell, while most crystal structures capture the structure in the catalytic dwell. Principal component analysis of the available crystal structures further clarifies the five-step conformational transitions of the αβ-dimer in the ATPase cycle, highlighting the two dominant modes: the opening/closing motions of β and the loosening/tightening motions at the αβ-interface. These results provide a new view of tripartite coupling among chemical reactions, stator conformations, and rotary angles in F1-ATPase.

Inherent conformational flexibility of F1-ATPase α-subunit

The core of F1-ATPase consists of three catalytic (β) and three noncatalytic (α) subunits, forming a hexameric ring in alternating positions. A wealth of experimental and theoretical data has provided a detailed picture of the complex role played by catalytic subunits. Although major conformational changes have only been seen in β-subunits, it is clear that α-subunits have to respond to these changes in order to be able to transmit information during the rotary mechanism. However, the conformational behavior of α-subunits has not been explored in detail. Here, we have combined unbiased molecular dynamics (MD) simulations and calorimetrically measured thermodynamic signatures to investigate the conformational flexibility of isolated α-subunits, as a step toward deepening our understanding of its function inside the α3β3 ring. The simulations indicate that the open-to-closed conformational transition of the α-subunit is essentially barrierless, which is ideal to accompany and transmit the movement of the catalytic subunits. Calorimetric measurements of the recombinant α-subunit from Geobacillus kaustophilus indicate that the isolated subunit undergoes no significant conformational changes upon nucleotide binding. Simulations confirm that the nucleotide-free and nucleotide-bound subunits show average conformations similar to that observed in the F1 crystal structure, but they reveal an increased conformational flexibility of the isolated α-subunit upon MgATP binding, which might explain the evolutionary conserved capacity of α-subunits to recognize nucleotides with considerable strength. Furthermore, we elucidate the different dependencies that α- and β-subunits show on Mg(II) for recognizing ATP.

Modeling Structural Dynamics of Biomolecular Complexes by Coarse-Grained Molecular Simulations

Due to hierarchic nature of biomolecular systems, their computational modeling calls for multiscale approaches, in which coarse-grained (CG) simulations are used to address long-time dynamics of large systems. Here, we review recent developments and applications of CG modeling methods, focusing on our methods primarily for proteins, DNA, and their complexes. These methods have been implemented in the CG biomolecular simulator, CafeMol. Our CG model has resolution such that ∼10 non-hydrogen atoms are grouped into one CG particle on average. For proteins, each amino acid is represented by one CG particle. For DNA, one nucleotide is simplified by three CG particles, representing sugar, phosphate, and base. The protein modeling is based on the idea that proteins have a globally funnel-like energy landscape, which is encoded in the structure-based potential energy function. We first describe two representative minimal models of proteins, called the elastic network model and the classic Go̅ model. We then present a more elaborate protein model, which extends the minimal model to incorporate sequence and context dependent local flexibility and nonlocal contacts. For DNA, we describe a model developed by de Pablo's group that was tuned to well reproduce sequence-dependent structural and thermodynamic experimental data for single- and double-stranded DNAs. Protein-DNA interactions are modeled either by the structure-based term for specific cases or by electrostatic and excluded volume terms for nonspecific cases. We also discuss the time scale mapping in CG molecular dynamics simulations. While the apparent single time step of our CGMD is about 10 times larger than that in the fully atomistic molecular dynamics for small-scale dynamics, large-scale motions can be further accelerated by two-orders of magnitude with the use of CG model and a low friction constant in Langevin dynamics. Next, we present four examples of applications. First, the classic Go̅ model was used to emulate one ATP cycle of a molecular motor, kinesin. Second, nonspecific protein-DNA binding was studied by a combination of elaborate protein and DNA models. Third, a transcription factor, p53, that contains highly fluctuating regions was simulated on two perpendicularly arranged DNA segments, addressing intersegmental transfer of p53. Fourth, we simulated structural dynamics of dinucleosomes connected by a linker DNA finding distinct types of internucleosome docking and salt-concentration-dependent compaction. Finally, we discuss many of limitations in the current approaches and future directions. Especially, more accurate electrostatic treatment and a phospholipid model that matches our CG resolutions are of immediate importance.

Domain Motion Enhanced (DoME) Model for Efficient Conformational Sampling of Multidomain Proteins

Large conformational changes of multidomain proteins are difficult to simulate using all-atom molecular dynamics (MD) due to the slow time scale. We show that a simple modification of the structure-based coarse-grained (CG) model enables a stable and efficient MD simulation of those proteins. "Motion Tree", a tree diagram that describes conformational changes between two structures in a protein, provides information on rigid structural units (domains) and the magnitudes of domain motions. In our new CG model, which we call the DoME (domain motion enhanced) model, interdomain interactions are defined as being inversely proportional to the magnitude of the domain motions in the diagram, whereas intradomain interactions are kept constant. We applied the DoME model in combination with the Go model to simulations of adenylate kinase (AdK). The results of the DoME-Go simulation are consistent with an all-atom MD simulation for 10 μs as well as known experimental data. Unlike the conventional Go model, the DoME-Go model yields stable simulation trajectories against temperature changes and conformational transitions are easily sampled despite domain rigidity. Evidently, identification of domains and their interfaces is useful approach for CG modeling of multidomain proteins.

Elasticity, friction, and pathway of γ-subunit rotation in FoF1-ATP synthase

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 k(B)T 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.

Rotational mechanism of Enterococcus hirae V1-ATPase by crystal-structure and single-molecule analyses

In ion-transporting rotary ATPases, the mechanical rotation of inner rotor subunits against other stator subunits in the complex mediates conversion of chemical free energy from ATP hydrolysis into electrochemical potential by pumping ions across the cell membrane. To fully understand the rotational mechanism of energy conversion, it is essential to analyze a target sample by multiple advanced methods that differ in spatiotemporal resolutions and sample environments. Here, we describe such a strategy applied to the water-soluble V1 moiety of Enterococcus hirae V-ATPase; this strategy involves integration of crystal structure studies and single-molecule analysis of rotary dynamics and torque generation. In addition, we describe our current model of the chemo-mechanical coupling scheme obtained by this approach, as well as future prospects.

Nucleotide-dependent structural fluctuations and regulation of microtubule-binding affinity of KIF1A

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single-headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō-like model. We found that the helix α4 at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix α4. These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix α4 during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state.

Dissecting the role of the γ-subunit in the rotary-chemical coupling and torque generation of F1-ATPase

Unraveling the molecular nature of the conversion of chemical energy (ATP hydrolysis in the α/β-subunits) to mechanical energy and torque (rotation of the γ-subunit) in F1-ATPase is very challenging. A major part of the challenge involves understanding the rotary-chemical coupling by a nonphenomenological structure-energy description, while accounting for the observed torque generated on the γ-subunit and its change due to mutation of this unit. Here we extend our previous study that used a coarse-grained model of the F1-ATPase to generate a structure-based free energy landscape of the rotary-chemical process. Our quantitative analysis of the landscape reproduced the observed torque for the wild-type enzyme. In doing so, we found that there are several possibilities of torque generation from landscapes with various shapes and demonstrated that a downhill slope along the chemical coordinate could still result in negligible torque, due to ineffective coupling of the chemistry to the γ-subunit rotation. We then explored the relationship between the functionality and the underlying sequence through systematic examination of the effect of various parts of the γ-subunit on free energy surfaces of F1-ATPase. Furthermore, by constructing several types of γ-deletion systems and calculating the corresponding torque generation, we gained previously unknown insights into the molecular nature of the F1-ATPase rotary motor. Significantly, our results are in excellent agreement with recent experimental findings and indicate that the rotary-chemical coupling is primarily established through electrostatic effects, although specific contacts through γ-ionizable residue side chains are not essential for establishing the basic features of the coupling.

None of the rotor residues of F1-ATPase are essential for torque generation

F₁-ATPase is a powerful rotary molecular motor that can rotate an object several hundred times as large as the motor itself against the viscous friction of water. Forced reverse rotation has been shown to lead to ATP synthesis, implying that the mechanical work against the motor’s high torque can be converted into the chemical energy of ATP. The minimal composition of the motor protein is α₃β₃γ subunits, where the central rotor subunit γ turns inside a stator cylinder made of alternately arranged α₃β₃ subunits using the energy derived from ATP hydrolysis. The rotor consists of an axle, a coiled coil of the amino- and carboxyl-terminal α-helices of γ, which deeply penetrates the stator cylinder, and a globular protrusion that juts out from the stator. Previous work has shown that, for a thermophilic F₁, significant portions of the axle can be truncated and the motor still rotates a submicron sized bead duplex, indicating generation of up to half the wild-type (WT) torque. Here, we inquire if any specific interactions between the stator and the rest of the rotor are needed for the generation of a sizable torque. We truncated the protruding portion of the rotor and replaced part of the remaining axle residues such that every residue of the rotor has been deleted or replaced in this or previous truncation mutants. This protrusionless construct showed an unloaded rotary speed about a quarter of the WT, and generated one-third to one-half of the WT torque. No residue-specific interactions are needed for this much performance. F₁ is so designed that the basic rotor-stator interactions for torque generation and control of catalysis rely solely upon the shape and size of the rotor at very low resolution. Additional tailored interactions augment the torque to allow ATP synthesis under physiological conditions.

Trapping the ATP binding state leads to a detailed understanding of the F1-ATPase mechanism

The rotary motor enzyme FoF1-ATP synthase uses the proton-motive force across a membrane to synthesize ATP from ADP and Pi (H2PO4(-)) under cellular conditions that favor the hydrolysis reaction by a factor of 2 × 10(5). This remarkable ability to drive a reaction away from equilibrium by harnessing an external force differentiates it from an ordinary enzyme, which increases the rate of reaction without shifting the equilibrium. Hydrolysis takes place in the neighborhood of one conformation of the catalytic moiety F1-ATPase, whose structure is known from crystallography. By use of molecular dynamics simulations we trap a second structure, which is rotated by 40° from the catalytic dwell conformation and represents the state associated with ATP binding, in accord with single-molecule experiments. Using the two structures, we show why Pi is not released immediately after ATP hydrolysis, but only after a subsequent 120° rotation, in agreement with experiment. A concerted conformational change of the α3β3 crown is shown to induce the 40° rotation of the γ-subunit only when the βE subunit is empty, whereas with Pi bound, βE serves as a latch to prevent the rotation of γ. The present results provide a rationalization of how F1-ATPase achieves the coupling between the small changes in the active site of βDP and the 40° rotation of γ.

Molecular structure and rotary dynamics of Enterococcus hirae V₁-ATPase

V1-ATPase is a rotary molecular motor in which the mechanical rotation of the rotor DF subunits against the stator A3B3 ring is driven by the chemical free energy of ATP hydrolysis. Recently, using X-ray crystallography, we solved the high-resolution molecular structure of Enterococcus hirae V1-ATPase (EhV1) and revealed how the three catalytic sites in the stator A3B3 ring change their structure on nucleotide binding and interaction with the rotor DF subunits. Furthermore, recently, we also demonstrated directly the rotary catalysis of EhV1 by using single-molecule high-speed imaging and analyzed the properties of the rotary motion in detail. In this critical review, we introduce the molecular structure and rotary dynamics of EhV1 and discuss a possible model of its chemomechanical coupling scheme.

Structure-based simulations reveal concerted dynamics of GPCR activation

G protein-coupled receptors (GPCRs) are a vital class of proteins that transduce biological signals across the cell membrane. However, their allosteric activation mechanism is not fully understood; crystal structures of active and inactive receptors have been reported, but the functional pathway between these two states remains elusive. Here, we use structure-based (Gō-like) models to simulate activation of two GPCRs, rhodopsin and the β₂ adrenergic receptor (β₂AR). We used data-derived reaction coordinates that capture the activation mechanism for both proteins, showing that activation proceeds through quantitatively different paths in the two systems. Both reaction coordinates are determined from the dominant concerted motions in the simulations so the technique is broadly applicable. There were two surprising results. First, the main structural changes in the simulations were distributed throughout the transmembrane bundle, and not localized to the obvious areas of interest, such as the intracellular portion of Helix 6. Second, the activation (and deactivation) paths were distinctly nonmonotonic, populating states that were not simply interpolations between the inactive and active structures. These transitions also suggest a functional explanation for β₂AR's basal activity: it can proceed through a more broadly defined path during the observed transitions.

Single-molecule analysis of the rotation of F₁-ATPase under high hydrostatic pressure

F1-ATPase is the water-soluble part of ATP synthase and is an ATP-driven rotary molecular motor that rotates the rotary shaft against the surrounding stator ring, hydrolyzing ATP. Although the mechanochemical coupling mechanism of F1-ATPase has been well studied, the molecular details of individual reaction steps remain unclear. In this study, we conducted a single-molecule rotation assay of F1 from thermophilic bacteria under various pressures from 0.1 to 140 MPa. Even at 140 MPa, F1 actively rotated with regular 120° steps in a counterclockwise direction, showing high conformational stability and retention of native properties. Rotational torque was also not affected. However, high hydrostatic pressure induced a distinct intervening pause at the ATP-binding angles during continuous rotation. The pause was observed under both ATP-limiting and ATP-saturating conditions, suggesting that F1 has two pressure-sensitive reactions, one of which is evidently ATP binding. The rotation assay using a mutant F1(βE190D) suggested that the other pressure-sensitive reaction occurs at the same angle at which ATP binding occurs. The activation volumes were determined from the pressure dependence of the rate constants to be +100 Å(3) and +88 Å(3) for ATP binding and the other pressure-sensitive reaction, respectively. These results are discussed in relation to recent single-molecule studies of F1 and pressure-induced protein unfolding.

Anatomy of F1-ATPase powered rotation

F1-ATPase, the catalytic complex of the ATP synthase, is a molecular motor that can consume ATP to drive rotation of the γ-subunit inside the ring of three αβ-subunit heterodimers in 120° power strokes. To elucidate the mechanism of ATPase-powered rotation, we determined the angular velocity as a function of rotational position from single-molecule data collected at 200,000 frames per second with unprecedented signal-to-noise. Power stroke rotation is more complex than previously understood. This paper reports the unexpected discovery that a series of angular accelerations and decelerations occur during the power stroke. The decreases in angular velocity that occurred with the lower-affinity substrate ITP, which could not be explained by an increase in substrate-binding dwells, provides direct evidence that rotation depends on substrate binding affinity. The presence of elevated ADP concentrations not only increased dwells at 35° from the catalytic dwell consistent with competitive product inhibition but also decreased the angular velocity from 85° to 120°, indicating that ADP can remain bound to the catalytic site where product release occurs for the duration of the power stroke. The angular velocity profile also supports a model in which rotation is powered by Van der Waals repulsive forces during the final 85° of rotation, consistent with a transition from F1 structures 2HLD1 and 1H8E (Protein Data Bank).

Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics

Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.

Force-induced unzipping transitions in an athermal crowded environment

Using theoretical arguments and extensive Monte Carlo (MC) simulations of a coarse-grained three-dimensional off-lattice model of a β-hairpin, we demonstrate that the equilibrium critical force, Fc, needed to unfold the biopolymer increases nonlinearly with increasing volume fraction occupied by the spherical macromolecular crowding agent. Both scaling arguments and MC simulations show that the critical force increases as Fc ≈ φc(α). The exponent α is linked to the Flory exponent relating the size of the unfolded state of the biopolymer and the number of amino acids. The predicted power law dependence is confirmed in simulations of the dependence of the isothermal extensibility and the fraction of native contacts on φc. We also show using MC simulations that Fc is linearly dependent on the average osmotic pressure (P) exerted by the crowding agents on the β-hairpin. The highly significant linear correlation coefficient of 0.99657 between Fc and P makes it straightforward to predict the dependence of the critical force on the density of crowders. Our predictions are amenable to experimental verification using laser optical tweezers.

Comparing a simple theoretical model for protein folding with all-atom molecular dynamics simulations

Advances in computing have enabled microsecond all-atom molecular dynamics trajectories of protein folding that can be used to compare with and test critical assumptions of theoretical models. We show that recent simulations by the Shaw group (10, 11, 14, 15) are consistent with a key assumption of an Ising-like theoretical model that native structure grows in only a few regions of the amino acid sequence as folding progresses. The distribution of mechanisms predicted by simulating the master equation of this native-centric model for the benchmark villin subdomain, with only two adjustable thermodynamic parameters and one temperature-dependent kinetic parameter, is remarkably similar to the distribution in the molecular dynamics trajectories.

Phosphate release coupled to rotary motion of F1-ATPase

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.

Drug uptake pathways of multidrug transporter AcrB studied by molecular simulations and site-directed mutagenesis experiments

Multidrug resistance has been a critical issue in current chemotherapy. In Escherichia coli , a major efflux pump responsible for the multidrug resistance contains a transporter AcrB. Crystallographic studies and mutational assays of AcrB provided much of structural and overall functional insights, which led to the functionally rotating mechanism. However, the drug uptake pathways are somewhat controversial because at least two possible pathways, the vestibule and the cleft paths, were suggested. Here, combining molecular simulations and site-directed mutagenesis experiments, we addressed the uptake mechanism finding that the drug uptake pathways can be significantly different depending on the properties of drugs. First, in the computational free energy analysis of drug movements along AcrB tunnels, we found a ligand-dependent drug uptake mechanism. With the same molecular sizes, drugs that are both strongly hydrophobic and lipophilic were preferentially taken in via the vestibule path, while other drugs favored the cleft path. Second, direct simulations realized totally about 3500 events of drug uptake by AcrB for a broad range of drug property. These simulations confirmed the ligand-dependent drug uptake and further suggested that a smaller drug favors the vestibule path, while a larger one is taken in via the cleft path. Moreover, the direct simulations identified an alternative uptake path which is not visible in the crystal structure. Third, site-directed mutagenesis of AcrB in E. coli verified that mutations of residues located along the newly identified path significantly reduced the efflux efficiency, supporting its relevance in in vivo function.

Molecular dynamics simulations of yeast F1-ATPase before and after 16° rotation of the γ subunit

We have recently proposed the "packing exchange mechanism" for F1-ATPase, wherein the perturbation by a substrate binding/release or an ATP hydrolysis is followed by the reorganization of the asymmetric packing structure of the α3β3 complex, accompanying the γ subunit rotation. As part of a further investigation of this rotational mechanism, we performed all-atom molecular dynamics simulations for yeast F1-ATPase both before and after a 16° rotation of the γ subunit triggered by a Pi release. We analyzed the structural fluctuations, the subunit interface interactions, and the dynamics of the relative subunit arrangements before and after the rotation. We found that with the Pi release the αEβE subunit interface becomes looser, which also allosterically makes the αDPβDP subunit interface looser. This structural communication between these interfaces takes place through a tightening of the αTPβTP subunit interface. The γ subunit interacts less strongly with αDP and βDP and more strongly with αTP and βTP. After the Pi release, the tightly packed interfaces are reorganized from the interfaces around βDP to those around βTP, inducing the 16° rotation. These results, which are consistent with the packing exchange mechanism, allow us to deduce a view of the structural change during the 40° rotation.

Structure-based molecular simulations reveal the enhancement of biased Brownian motions in single-headed kinesin

Kinesin is a family of molecular motors that move unidirectionally along microtubules (MT) using ATP hydrolysis free energy. In the family, the conventional two-headed kinesin was experimentally characterized to move unidirectionally through “walking” in a hand-over-hand fashion by coordinated motions of the two heads. Interestingly a single-headed kinesin, a truncated KIF1A, still can generate a biased Brownian movement along MT, as observed by in vitro single molecule experiments. Thus, KIF1A must use a different mechanism from the conventional kinesin to achieve the unidirectional motions. Based on the energy landscape view of proteins, for the first time, we conducted a set of molecular simulations of the truncated KIF1A movements over an ATP hydrolysis cycle and found a mechanism exhibiting and enhancing stochastic forward-biased movements in a similar way to those in experiments. First, simulating stand-alone KIF1A, we did not find any biased movements, while we found that KIF1A with a large friction cargo-analog attached to the C-terminus can generate clearly biased Brownian movements upon an ATP hydrolysis cycle. The linked cargo-analog enhanced the detachment of the KIF1A from MT. Once detached, diffusion of the KIF1A head was restricted around the large cargo which was located in front of the head at the time of detachment, thus generating a forward bias of the diffusion. The cargo plays the role of a diffusional anchor, or cane, in KIF1A “walking.”

It is one of the major issues in biophysics how molecular motors such as conventional two-headed kinesin convert the chemical energy released at ATP hydrolysis into mechanical work. While most molecular motors move with more than one catalytic domain working in coordinated fashions, there are some motors that can move with only a single catalytic domain, which provides us a possibly simpler case to understand. A single-headed kinesin, KIF1A, with only one catalytic domain, has been characterized by in vitro single-molecule assay to generate a biased Brownian movement along the microtubule. Here, we conducted a set of structure-based coarse-grained molecular simulations for KIF1A system over an ATP hydrolysis cycle for the first time to our knowledge. Without cargo the simulated stand-alone KIF1A could not generate any directional movement, while a large-friction cargo-analog linked to the C-terminus of KIF1A clearly enhanced the forward-biased Brownian movement of KIF1A significantly. Interestingly, the cargo-analog here is not merely load but an important promoter to enable efficient movements of KIF1A.

Role of the DELSEED loop in torque transmission of F1-ATPase

F(1)-ATPase is an ATP-driven rotary motor that generates torque at the interface between the catalytic β-subunits and the rotor γ-subunit. The β-subunit inwardly rotates the C-terminal domain upon nucleotide binding/dissociation; hence, the region of the C-terminal domain that is in direct contact with γ-termed the DELSEED loop-is thought to play a critical role in torque transmission. We substituted all the DELSEED loop residues with alanine to diminish specific DELSEED loop-γ interactions and with glycine to disrupt the loop structure. All the mutants rotated unidirectionally with kinetic parameters comparable to those of the wild-type F(1), suggesting that the specific interactions between DELSEED loop and γ is not involved in cooperative interplays between the catalytic β-subunits. Glycine substitution mutants generated half the torque of the wild-type F(1), whereas the alanine mutant generated comparable torque. Fluctuation analyses of the glycine/alanine mutants revealed that the γ-subunit was less tightly held in the α(3)β(3)-stator ring of the glycine mutant than in the wild-type F(1) and the alanine mutant. Molecular dynamics simulation showed that the DELSEED loop was disordered by the glycine substitution, whereas it formed an α-helix in the alanine mutant. Our results emphasize the importance of loop rigidity for efficient torque transmissions.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.