Identification of the betaTP site in the x-ray structure of F1-ATPase as the high-affinity catalytic site

Nonlethal deleterious mutation-induced stress accelerates bacterial aging

Random mutagenesis, including when it leads to loss of gene function, is a key mechanism enabling microorganisms' long-term adaptation to new environments. However, loss-of-function mutations are often deleterious, triggering, in turn, cellular stress and complex homeostatic stress responses, called "allostasis," to promote cell survival. Here, we characterize the differential impacts of 65 nonlethal, deleterious single-gene deletions on Escherichia coli growth in three different growth environments. Further assessments of select mutants, namely, those bearing single adenosine triphosphate (ATP) synthase subunit deletions, reveal that mutants display reorganized transcriptome profiles that reflect both the environment and the specific gene deletion. We also find that ATP synthase α-subunit deleted (ΔatpA) cells exhibit elevated metabolic rates while having slower growth compared to wild-type (wt) E. coli cells. At the single-cell level, compared to wt cells, individual ΔatpA cells display near normal proliferation profiles but enter a postreplicative state earlier and exhibit a distinct senescence phenotype. These results highlight the complex interplay between genomic diversity, adaptation, and stress response and uncover an "aging cost" to individual bacterial cells for maintaining population-level resilience to environmental and genetic stress; they also suggest potential bacteriostatic antibiotic targets and -as select human genetic diseases display highly similar phenotypes, - a bacterial origin of some human diseases.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

A review of TNP-ATP in protein binding studies: benefits and pitfalls

We review 50 years of use of 2',3'-O-trinitrophenyl (TNP)-ATP, a fluorescently tagged ATP analog. It has been extensively used to detect binding interactions of ATP to proteins and to measure parameters of those interactions such as the dissociation constant, K_d, or inhibitor dissociation constant, K_i. TNP-ATP has also found use in other applications, for example, as a fluorescence marker in microscopy, as a FRET pair, or as an antagonist (e.g., of P2X receptors). However, its use in protein binding studies has limitations because the TNP moiety often enhances binding affinity, and the fluorescence changes that occur with binding can be masked or mimicked in unexpected ways. The goal of this review is to provide a clear perspective of the pros and cons of using TNP-ATP to allow for better experimental design and less ambiguous data in future experiments using TNP-ATP and other TNP nucleotides.

The nucleotide binding affinities of two critical conformations of Escherichia coli ATP synthase

ATP synthase is essential in aerobic energy metabolism, and the rotary catalytic mechanism is one of the core concepts to understand the energetic functions of ATP synthase. Disulfide bonds formed by oxidizing a pair of cysteine mutations halted the rotation of the γ subunit in two critical conformations, the ATP-waiting dwell (αE284C/γQ274C) and the catalytic dwell (αE284C/γL276C). Tryptophan fluorescence was used to measure the nucleotide binding affinities for MgATP, MgADP and MgADP-AlF₄ (a transition state analog) to wild-type and mutant F₁ under reducing and oxidizing conditions. In the reduced state, αE284C/γL276C F₁ showed a wild-type-like nucleotide binding pattern; after oxidation to lock the enzyme in the catalytic dwell state, the nucleotide binding parameters remained unchanged. In contrast, αE284C/γQ274C F₁ showed significant differences in the affinities of the oxidized versus the reduced state. Locking the enzyme in the ATP-waiting dwell reduced nucleotide binding affinities of all three catalytic sites. Most importantly, the affinity of the low affinity site was reduced to such an extent that it could no longer be detected in the binding assay (K_d > 5 mM). The results of the present study allow to present a model for the catalytic mechanism of ATP synthase under consideration of the nucleotide affinity changes during a 360° cycle of the rotor.

Structural Asymmetry and Kinetic Limping of Single Rotary F-ATP Synthases

F-ATP synthases use proton flow through the FO domain to synthesize ATP in the F₁ domain. In Escherichia coli, the enzyme consists of rotor subunits γεc10 and stator subunits (αβ)₃δab₂. Subunits c10 or (αβ)₃ alone are rotationally symmetric. However, symmetry is broken by the b₂ homodimer, which together with subunit δa, forms a single eccentric stalk connecting the membrane embedded FO domain with the soluble F₁ domain, and the central rotating and curved stalk composed of subunit γε. Although each of the three catalytic binding sites in (αβ)₃ catalyzes the same set of partial reactions in the time average, they might not be fully equivalent at any moment, because the structural symmetry is broken by contact with b₂δ in F₁ and with b₂a in FO. We monitored the enzyme's rotary progression during ATP hydrolysis by three single-molecule techniques: fluorescence video-microscopy with attached actin filaments, Förster resonance energy transfer between pairs of fluorescence probes, and a polarization assay using gold nanorods. We found that one dwell in the three-stepped rotary progression lasting longer than the other two by a factor of up to 1.6. This effect of the structural asymmetry is small due to the internal elastic coupling.

Identification of two segments of the γ subunit of ATP synthase responsible for the different affinities of the catalytic nucleotide-binding sites

ATP synthase uses a rotary mechanism to couple transmembrane proton translocation to ATP synthesis and hydrolysis, which occur at the catalytic sites in the β subunits. In the presence of Mg²⁺, the three catalytic sites of ATP synthase have vastly different affinities for nucleotides, and the position of the central γ subunit determines which site has high, medium, or low affinity. Affinity differences and their changes as rotation progresses underpin the ATP synthase catalytic mechanism. Here, we used a series of variants with up to 45- and 60-residue-long truncations of the N- and C-terminal helices of the γ subunit, respectively, to identify the segment(s) responsible for the affinity differences of the catalytic sites. We found that each helix carries an affinity-determining segment of ∼10 residues. Our findings suggest that the affinity regulation by these segments is transmitted to the catalytic sites by the DELSEED loop in the C-terminal domain of the β subunits. For the N-terminal truncation variants, presence of the affinity-determining segment and therefore emergence of a high-affinity binding site resulted in WT-like catalytic activity. At the C terminus, additional residues outside of the affinity-determining segment were required for optimal enzymatic activity. Alanine substitutions revealed that the affinity changes of the catalytic sites required no specific interactions between amino acid side chains in the γ and α₃β₃ subunits but were caused by the presence of the helices themselves. Our findings help unravel the molecular basis for the affinity changes of the catalytic sites during ATP synthase rotation.

The regulatory subunit ε in Escherichia coli FOF1-ATP synthase

F-type ATP synthases are extraordinary multisubunit proteins that operate as nanomotors. The Escherichia coli (E. coli) enzyme uses the proton motive force (pmf) across the bacterial plasma membrane to drive rotation of the central rotor subunits within a stator subunit complex. Through this mechanical rotation, the rotor coordinates three nucleotide binding sites that sequentially catalyze the synthesis of ATP. Moreover, the enzyme can hydrolyze ATP to turn the rotor in the opposite direction and generate pmf. The direction of net catalysis, i.e. synthesis or hydrolysis of ATP, depends on the cell's bioenergetic conditions. Different control mechanisms have been found for ATP synthases in mitochondria, chloroplasts and bacteria. This review discusses the auto-inhibitory behavior of subunit ε found in F_OF₁-ATP synthases of many bacteria. We focus on E. coli F_OF₁-ATP synthase, with insights into the regulatory mechanism of subunit ε arising from structural and biochemical studies complemented by single-molecule microscopy experiments.

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.

Mechanism of the αβ conformational change in F1-ATPase after ATP hydrolysis: free-energy simulations

One of the motive forces for F1-ATPase rotation is the conformational change of the catalytically active β subunit due to closing and opening motions caused by ATP binding and hydrolysis, respectively. The closing motion is accomplished in two steps: the hydrogen-bond network around ATP changes and then the entire structure changes via B-helix sliding, as shown in our previous study. Here, we investigated the opening motion induced by ATP hydrolysis using all-atom free-energy simulations, combining the nudged elastic band method and umbrella sampling molecular-dynamics simulations. Because hydrolysis requires residues in the α subunit, the simulations were performed with the αβ dimer. The results indicate that the large-scale opening motion is also achieved by the B-helix sliding (in the reverse direction). However, the sliding mechanism is different from that of ATP binding because sliding is triggered by separation of the hydrolysis products ADP and Pi. We also addressed several important issues: 1), the timing of the product Pi release; 2), the unresolved half-closed β structure; and 3), the ADP release mechanism. These issues are fundamental for motor function; thus, the rotational mechanism of the entire F1-ATPase is also elucidated through this αβ study. During the conformational change, conserved residues among the ATPase proteins play important roles, suggesting that the obtained mechanism may be shared with other ATPase proteins. When combined with our previous studies, these results provide a comprehensive view of the β-subunit conformational change that drives the ATPase.

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).

Principal role of the arginine finger in rotary catalysis of F1-ATPase

F(1)-ATPase (F(1)) is an ATP-driven rotary motor wherein the γ subunit rotates against the surrounding α(3)β(3) stator ring. The 3 catalytic sites of F(1) reside on the interface of the α and β subunits of the α(3)β(3) ring. While the catalytic residues predominantly reside on the β subunit, the α subunit has 1 catalytically critical arginine, termed the arginine finger, with stereogeometric similarities with the arginine finger of G-protein-activating proteins. However, the principal role of the arginine finger of F(1) remains controversial. We studied the role of the arginine finger by analyzing the rotation of a mutant F(1) with a lysine substitution of the arginine finger. The mutant showed a 350-fold longer catalytic pause than the wild-type; this pause was further lengthened by the slowly hydrolyzed ATP analog ATPγS. On the other hand, the mutant F(1) showed highly unidirectional rotation with a coupling ratio of 3 ATPs/turn, the same as wild-type, suggesting that cooperative torque generation by the 3 β subunits was not impaired. The hybrid F(1) carrying a single copy of the α mutant revealed that the reaction step slowed by the mutation occurs at +200° from the binding angle of the mutant subunit. Thus, the principal role of the arginine finger is not to mediate cooperativity among the catalytic sites, but to enhance the rate of the ATP cleavage by stabilizing the transition state of ATP hydrolysis. Lysine substitution also caused frequent pauses because of severe ADP inhibition, and a slight decrease in ATP-binding rate.

The beta subunit loop that couples catalysis and rotation in ATP synthase has a critical length

ATP synthase uses a unique rotational mechanism to convert chemical energy into mechanical energy and back into chemical energy. The helix-turn-helix structure in the C-terminal domain of the β subunit containing the conserved DELSEED motif, termed "DELSEED-loop," was suggested to be involved in coupling between catalysis and rotation. If this is indeed the role of the loop, it must have a critical length, the minimum length required to sustain its function. Here, the critical length of the DELSEED-loop was determined by functional analysis of mutants of Bacillus PS3 ATP synthase that had 7-14 amino acids within the loop deleted. A 10 residue deletion lost the ability to catalyze ATP synthesis, but was still an active ATPase. Deletion of 14 residues abolished any enzymatic activity. Modeling indicated that in both deletion mutants the DELSEED-loop was shortened by ∼10 Å; fluorescence resonance energy transfer experiments confirmed the modeling results. This appears to define the minimum length for DELSEED-loop required for coupling of catalysis and rotation. In addition, we could demonstrate that the loss of high-affinity binding to the catalytic site(s) that had been observed previously in two deletion mutants with 3-4 residues removed was not due to the loss of negative charged residues of the DELSEED motif in these mutants. An AALSAAA mutant in which all negative charges of the DELSEED motif were removed showed a normal pattern for MgATP binding to the catalytic sites, with a clearly present high-affinity site.

Domain motion of individual F1-ATPase β-subunits during unbiased molecular dynamics simulations

F(1)-ATPase is the catalytic domain of F(1)F(o)-ATP synthase and consists of a hexameric arrangement of three noncatalytic α and three catalytic β subunits. We have used unbiased molecular dynamics simulations with a total simulation time of 900 ns to investigate the dynamic relaxation properties of isolated β-subunits as a step toward explaining the function of the integral F(1) unit. To this end, we simulated the open (β(E)) and the closed (β(TP)) conformations under unbiased conditions for up to 120 ns each using several samples. The simulations confirm that nucleotide-free β(E) retains its open configuration over the course of the simulations. The same is true when the neighboring α subunits are included. The nucleotide-depleted as well as the nucleotide-bound isolated β(TP) subunits show a significant trend toward the open conformation during our simulations, with one trajectory per case opening completely. Hence, our simulations suggest that the equilibrium conformation of a nucleotide-free β-subunit is the open conformation and that the transition from the closed to the open conformation can occur on a time scale of a few tens of nanoseconds.

Structural divergence of the rotary ATPases

The rotary ATPase family of membrane protein complexes may have only three members, but each one plays a fundamental role in biological energy conversion. The F1Fo-ATPase (F-ATPase) couples ATP synthesis to the electrochemical membrane potential in bacteria, mitochondria and chloroplasts, while the vacuolar H+-ATPase (V-ATPase) operates as an ATP-driven proton pump in eukaryotic membranes. In different species of archaea and bacteria, the A1Ao-ATPase (A-ATPase) can function as either an ATP synthase or an ion pump. All three of these multi-subunit complexes are rotary molecular motors, sharing a fundamentally similar mechanism in which rotational movement drives the energy conversion process. By analogy to macroscopic systems, individual subunits can be assigned to rotor, axle or stator functions. Recently, three-dimensional reconstructions from electron microscopy and single particle image processing have led to a significant step forward in understanding of the overall architecture of all three forms of these complexes and have allowed the organisation of subunits within the rotor and stator parts of the motors to be more clearly mapped out. This review describes the emerging consensus regarding the organisation of the rotor and stator components of V-, A- and F-ATPases, examining core similarities that point to a common evolutionary origin, and highlighting key differences. In particular, it discusses how newly revealed variation in the complexity of the inter-domain connections may impact on the mechanics and regulation of these molecular machines.

Structural fluctuation and concerted motions in F(1)-ATPase: A molecular dynamics study

F(1)-ATPase is an adenosine tri-phosphate (ATP)-driven rotary motor enzyme. We investigated the structural fluctuations and concerted motions of subunits in F(1)-ATPase using molecular dynamics (MD) simulations. An MD simulation for the alpha(3)beta(3)gamma complex was carried out for 30 ns. Although large fluctuations of the N-terminal domain observed in simulations of the isolated beta(E) subunit were suppressed in the complex simulation, the magnitude of fluctuations in the C-terminal domain was clearly different among the three beta subunits (beta(E), beta(TP), and beta(DP)). Despite fairly similar conformations of the beta(TP) and beta(DP) subunits, the beta(DP) subunit exhibits smaller fluctuations in the C-terminal domain than the beta(TP) subunit due to their dissimilar interface configurations. Compared with the beta(TP) subunit, the beta(DP) subunit stably interacts with both the adjacent alpha(DP) and alpha(E) subunits. This sandwiched configuration in the beta(DP) subunit leads to strongly correlated motions between the beta(DP) and adjacent alpha subunits. The beta(DP) subunit exhibits an extensive network of highly correlated motions with bound ATP and the gamma subunit, as well as with the adjacent alpha subunits, suggesting that the structural changes occurring in the catalytically active beta(DP) subunit can effectively induce movements of the gamma subunit.

Single molecule behavior of inhibited and active states of Escherichia coli ATP synthase F1 rotation

ATP hydrolysis-dependent rotation of the F(1) sector of the ATP synthase is a successive cycle of catalytic dwells (∼0.2 ms at 24 °C) and 120° rotation steps (∼0.6 ms) when observed under V(max) conditions using a low viscous drag 60-nm bead attached to the γ subunit (Sekiya, M., Nakamoto, R. K., Al-Shawi, M. K., Nakanishi-Matsui, M., and Futai, M. (2009) J. Biol. Chem. 284, 22401-22410). During the normal course of observation, the γ subunit pauses in a stochastic manner to a catalytically inhibited state that averages ∼1 s in duration. The rotation behavior with adenosine 5'-O-(3-thiotriphosphate) as the substrate or at a low ATP concentration (4 μM) indicates that the rotation is inhibited at the catalytic dwell when the bound ATP undergoes reversible hydrolysis/synthesis. The temperature dependence of rotation shows that F(1) requires ∼2-fold higher activation energy for the transition from the active to the inhibited state compared with that for normal steady-state rotation during the active state. Addition of superstoichiometric ε subunit, the inhibitor of F(1)-ATPase, decreases the rotation rate and at the same time increases the duration time of the inhibited state. Arrhenius analysis shows that the ε subunit has little effect on the transition between active and inhibited states. Rather, the ε subunit confers lower activation energy of steady-state rotation. These results suggest that the ε subunit plays a role in guiding the enzyme through the proper and efficient catalytic and transport rotational pathway but does not influence the transition to the inhibited state.

Phosphate release in F1-ATPase catalytic cycle follows ADP release

F(1)-ATPase is an ATP-driven rotary motor protein in which the γ-subunit rotates against the catalytic stator ring. Although the reaction scheme of F(1) has mostly been revealed, the timing of inorganic phosphate (P(i)) release remains controversial. Here we addressed this issue by verifying the reversibility of ATP hydrolysis on arrested F(1) with magnetic tweezers. ATP hydrolysis was found to be essentially reversible, implying that P(i) is released after the γ rotation and ADP release, although extremely slow P(i) release was found at the ATP hydrolysis angle as an uncoupling side reaction. On the basis of this finding, we deduced the chemomechanical coupling scheme of F(1). We found that the affinity for P(i) was strongly angle dependent, implying a large contribution by P(i) release to torque generation. These findings imply that under ATP synthesis conditions, P(i) binds to an empty catalytic site, preventing solution ATP (though not ADP) from binding. Thus, this supports the concept of selective ADP binding for efficient ATP synthesis.

Crystal structure of the Mg·ADP-inhibited state of the yeast F1c10-ATP synthase

The F(1)c(10) subcomplex of the yeast F(1)F(0)-ATP synthase includes the membrane rotor part c(10)-ring linked to a catalytic head, (αβ)(3), by a central stalk, γδε. The Saccharomyces cerevisiae yF(1)c(10)·ADP subcomplex was crystallized in the presence of Mg·ADP, dicyclohexylcarbodiimide (DCCD), and azide. The structure was solved by molecular replacement using a high resolution model of the yeast F(1) and a bacterial c-ring model with 10 copies of the c-subunit. The structure refined to 3.43-Å resolution displays new features compared with the original yF(1)c(10) and with the yF(1) inhibited by adenylyl imidodiphosphate (AMP-PNP) (yF(1)(I-III)). An ADP molecule was bound in both β(DP) and β(TP) catalytic sites. The α(DP)-β(DP) pair is slightly open and resembles the novel conformation identified in yF(1), whereas the α(TP)-β(TP) pair is very closed and resembles more a DP pair. yF(1)c(10)·ADP provides a model of a new Mg·ADP-inhibited state of the yeast F(1). As for the original yF(1) and yF(1)c(10) structures, the foot of the central stalk is rotated by ∼40 ° with respect to bovine structures. The assembly of the F(1) central stalk with the F(0) c-ring rotor is mainly provided by electrostatic interactions. On the rotor ring, the essential cGlu(59) carboxylate group is surrounded by hydrophobic residues and is not involved in hydrogen bonding.

Normal-mode-based modeling of allosteric couplings that underlie cyclic conformational transition in F(1) ATPase

F(1) ATPase, a rotary motor comprised of a central stalk (gamma subunit) enclosed by three alpha and beta subunits alternately arranged in a hexamer, features highly cooperative binding and hydrolysis of ATP. Despite steady progress in biophysical, biochemical, and computational studies of this fascinating motor, the structural basis for cooperative ATPase involving its three catalytic sites remains not fully understood. To illuminate this key mechanistic puzzle, we have employed a coarse-grained elastic network model to probe the allosteric couplings underlying the cyclic conformational transition in F(1) ATPase at a residue level of detail. We will elucidate how ATP binding and product (ADP and phosphate) release at two catalytic sites are coupled with the rotation of gamma subunit via various domain motions in alpha(3)beta(3) hexamer (including intrasubunit hinge-bending motions in beta subunits and intersubunit rigid-body rotations between adjacent alpha and beta subunits). To this end, we have used a normal-mode-based correlation analysis to quantify the allosteric couplings of these domain motions to local motions at catalytic sites and the rotation of gamma subunit. We have then identified key amino acid residues involved in the above couplings, some of which have been validated against past studies of mutated and gamma-truncated F(1) ATPase. Our finding strongly supports a binding change mechanism where ATP binding to the empty catalytic site triggers a series of intra- and intersubunit domain motions leading to ATP hydrolysis and product release at the other two closed catalytic sites.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

ATP synthase with its gamma subunit reduced to the N-terminal helix can still catalyze ATP synthesis

ATP synthase uses a unique rotary mechanism to couple ATP synthesis and hydrolysis to transmembrane proton translocation. As part of the synthesis mechanism, the torque of the rotor has to be converted into conformational rearrangements of the catalytic binding sites on the stator to allow synthesis and release of ATP. The gamma subunit of the rotor, which plays a central role in the energy conversion, consists of two long helices inside the central cavity of the stator cylinder plus a globular portion outside the cylinder. Here, we show that the N-terminal helix alone is able to fulfill the function of full-length gamma in ATP synthesis as long as it connects to the rest of the rotor. This connection can occur via the epsilon subunit. No direct contact between gamma and the c ring seems to be required. In addition, the results indicate that the epsilon subunit of the rotor exists in two different conformations during ATP synthesis and ATP hydrolysis.

Torque generation and elastic power transmission in the rotary F(O)F(1)-ATPase

Adenosine triphosphate (ATP), the universal fuel of the cell, is synthesized from adenosine diphosphate (ADP) and inorganic phosphate (P(i)) by 'ATP synthase' (F(O)F(1)-ATPase). During respiration or photosynthesis, an electrochemical potential difference of protons is set up across the respective membranes. This powers the enzyme's electrical rotary nanomotor (F(O)), which drives the chemical nanomotor (F(1)) by elastic mechanical-power transmission, producing ATP with high kinetic efficiency. Attempts to understand in detail the mechanisms of torque generation in this simple and robust system have been both aided and complicated by a wealth of sometimes conflicting data.

The role of the betaDELSEED-loop of ATP synthase

ATP synthase uses a unique rotational mechanism to convert chemical energy into mechanical energy and back into chemical energy. The helix-turn-helix motif, termed "DELSEED-loop," in the C-terminal domain of the beta subunit was suggested to be involved in coupling between catalysis and rotation. Here, the role of the DELSEED-loop was investigated by functional analysis of mutants of Bacillus PS3 ATP synthase that had 3-7 amino acids within the loop deleted. All mutants were able to catalyze ATP hydrolysis, some at rates several times higher than the wild-type enzyme. In most cases ATP hydrolysis in membrane vesicles generated a transmembrane proton gradient, indicating that hydrolysis occurred via the normal rotational mechanism. Except for two mutants that showed low activity and low abundance in the membrane preparations, the deletion mutants were able to catalyze ATP synthesis. In general, the mutants seemed less well coupled than the wild-type enzyme, to a varying degree. Arrhenius analysis demonstrated that in the mutants fewer bonds had to be rearranged during the rate-limiting catalytic step; the extent of this effect was dependent on the size of the deletion. The results support the idea of a significant involvement of the DELSEED-loop in mechanochemical coupling in ATP synthase. In addition, for two deletion mutants it was possible to prepare an alpha(3)beta(3)gamma subcomplex and measure nucleotide binding to the catalytic sites. Interestingly, both mutants showed a severely reduced affinity for MgATP at the high affinity site.

Functional halt positions of rotary FOF1-ATPase correlated with crystal structures

The F(O)F(1)-ATPase is a rotary molecular motor. Driven by ATP-hydrolysis, its central shaft rotates in 80 degrees and 40 degrees steps, interrupted by catalytic and ATP-waiting dwells. We recorded rotations and halts by means of microvideography in laboratory coordinates. A correlation with molecular coordinates was established by using an engineered pair of cysteines that, under oxidizing conditions, formed zero-length cross-links between the rotor and the stator in an orientation as found in crystals. The fixed orientation coincided with that of the catalytic dwell, whereas the ATP waiting dwell was displaced from it by +40 degrees . In crystals, the convex side of the cranked central shaft faces an empty nucleotide binding site, as if holding it open for arriving ATP. Functional studies suggest that three sites are occupied during a catalytic dwell. Our data imply that the convex side faces a nucleotide-occupied rather than an empty site. The enzyme conformation in crystals seems to differ from the conformation during either dwell of the active enzyme. A revision of current schemes of the mechanism is proposed.

Structural frameworks for considering microbial protein- and nucleic acid-dependent motor ATPases

Many fundamental cellular processes depend on enzymes that utilize chemical energy to catalyse unfavourable reactions. Certain classes of ATPases provide a particularly vivid example of the process of energy conversion, employing cycles of nucleotide turnover to move and/or rearrange biological polymers such as proteins and nucleic acids. Four well-characterized classes of ATP-dependent protein/nucleic acid translocases and remodelling factors are found in all three domains of life (bacteria, archaea and eukarya): additional strand catalytic 'E' (ASCE) P-loop NTPases, GHL proteins, actin-fold enzymes and chaperonins. These unrelated protein superfamilies have each evolved the ability to couple ATP binding and hydrolysis to the generation of motion and force along or within their substrates. The past several years have witnessed the emergence of a wealth of structural data that help explain how such molecular engines link nucleotide turnover to conformational change. In this review, we highlight several recent advances to illustrate some of the mechanisms by which each family of ATP-dependent motors facilitates the rearrangement and movement of proteins, protein complexes and nucleic acids.

A rotor-stator cross-link in the F1-ATPase blocks the rate-limiting step of rotational catalysis

The F(0)F(1)-ATP synthase couples the functions of H(+) transport and ATP synthesis/hydrolysis through the efficient transmission of energy mediated by rotation of the centrally located gamma, epsilon, and c subunits. To understand the gamma subunit role in the catalytic mechanism, we previously determined the partial rate constants and devised a minimal kinetic model for the rotational hydrolytic mode of the F(1)-ATPase enzyme that uniquely fits the pre-steady state and steady state data ( Baylis Scanlon, J. A., Al-Shawi, M. K., Le, N. P., and Nakamoto, R. K. (2007) Biochemistry 46, 8785-8797 ). Here we directly test the model using two single cysteine mutants, betaD380C and betaE381C, which can be used to reversibly inhibit rotation upon formation of a cross-link with the conserved gammaCys-87. In the pre-steady state, the gamma-beta cross-linked enzyme at high Mg.ATP conditions retained the burst of hydrolysis but was not able to release P(i). These data show that the rate-limiting rotation step, k(gamma), occurs after hydrolysis and before P(i) release. This analysis provides additional insights into how the enzyme achieves efficient coupling and implicates the betaGlu-381 residue for proper formation of the rate-limiting transition state involving gamma subunit rotation.

A functionally important hydrogen-bonding network at the betaDP/alphaDP interface of ATP synthase

ATP synthase uses a unique rotary mechanism to couple ATP synthesis and hydrolysis to transmembrane proton translocation. The F1 subcomplex has three catalytic nucleotide binding sites, one on each beta subunit, at the interface to the adjacent alpha subunit. In the x-ray structure of F1 (Abrahams, J. P., Leslie, A. G. W., Lutter, R., and Walker, J. E. (1994) Nature 370, 621-628), the three catalytic beta/alpha interfaces differ in the extent of inter-subunit interactions between the C termini of the beta and alpha subunits. At the closed betaDP/alphaDP interface, a hydrogen-bonding network is formed between both subunits, which is absent at the more open betaTP/alphaTP interface and at the wide open betaE/alphaE interface. The hydrogen-bonding network reaches from betaL328 (Escherichia coli numbering) and betaQ441 via alphaQ399, betaR398, and alphaE402 to betaR394, and ends in a cation/pi interaction between betaR394 and alphaF406. Using mutational analysis in E. coli ATP synthase, the functional importance of the betaDP/alphaDP hydrogen-bonding network is demonstrated. Its elimination results in a severely impaired enzyme but has no pronounced effect on the binding affinities of the catalytic sites. A possible role for the hydrogen-bonding network in coupling of ATP synthesis/hydrolysis and rotation will be discussed.

Molecular processes of inhibition and stimulation of ATP synthase caused by the phytotoxin tentoxin

F1-ATPase is the smallest mechanical motor known. Tentoxin, a cyclic peptide produced by phytopathogenic fungi, inactivates the F1 motor in sensitive plants at nanomolar to micromolar concentrations, whereas higher concentrations surpass the natural activity of the enzyme. Single molecule studies now have clarified the molecular steps involved in both processes. Inactivation delays the dwell time of a single step in the complete 360 degrees turn and results in an asymmetric rotation of the central rotor subunit. In contrast, rotation in the stimulated F1 particle is smooth and accompanied by strongly reduced ADP inhibition. Our study provides for the first time the direct observation of a noncompetitively inhibited state of the enzyme and directly visualizes the regulation of the molecular motor by an external natural compound. In addition, the ADP release step during catalysis was revealed by analysis of the single molecule rotation behavior. Hence, tentoxin is a sophisticated molecular tool to mark and control certain catalytic steps within the reaction pathway of the molecular F1 motor.

The rotary mechanism of the ATP synthase

The F0F1 ATP synthase is a large complex of at least 22 subunits, more than half of which are in the membranous F0 sector. This nearly ubiquitous transporter is responsible for the majority of ATP synthesis in oxidative and photo-phosphorylation, and its overall structure and mechanism have remained conserved throughout evolution. Most examples utilize the proton motive force to drive ATP synthesis except for a few bacteria, which use a sodium motive force. A remarkable feature of the complex is the rotary movement of an assembly of subunits that plays essential roles in both transport and catalytic mechanisms. This review addresses the role of rotation in catalysis of ATP synthesis/hydrolysis and the transport of protons or sodium.

How subunit coupling produces the gamma-subunit rotary motion in F1-ATPase

F(o)F(1)-ATP synthase manufactures the energy "currency," ATP, of living cells. The soluble F(1) portion, called F(1)-ATPase, can act as a rotary motor, with ATP binding, hydrolysis, and product release, inducing a torque on the gamma-subunit. A coarse-grained plastic network model is used to show at a residue level of detail how the conformational changes of the catalytic beta-subunits act on the gamma-subunit through repulsive van der Waals interactions to generate a torque that drives unidirectional rotation, as observed experimentally. The simulations suggest that the calculated 85 degrees substep rotation is driven primarily by ATP binding and that the subsequent 35 degrees substep rotation is produced by product release from one beta-subunit and a concomitant binding pocket expansion of another beta-subunit. The results of the simulation agree with single-molecule experiments [see, for example, Adachi K, et al. (2007) Cell 130:309-321] and support a tri-site rotary mechanism for F(1)-ATPase under physiological condition.