High efficacy of the F-ATP synthase inhibitor TBAJ-5307 against nontuberculous mycobacteria in vitro and in vivo.

The F1FO-ATP synthase engine is essential for viability and growth of nontuberculous mycobacteria (NTM) by providing the biological energy ATP and keeping ATP homeostasis under hypoxic stress conditions. Here, we report the discovery of the diarylquinoline TBAJ-5307 as a broad spectrum anti-NTM inhibitor, targeting the FO domain of the engine and preventing rotation and proton translocation. TBAJ-5307 is active at low nanomolar concentrations against fast- and slow-growing NTM as well as clinical isolates by depleting intrabacterial ATP. As demonstrated for the fast grower Mycobacterium abscessus, the compound is potent in vitro and in vivo, without inducing toxicity. Combining TBAJ-5307 with anti-NTM antibiotics or the oral tebipenem-avibactam pair showed attractive potentiation. Furthermore, the TBAJ-5307-tebipenem-avibactam cocktail kills the pathogen, suggesting a novel oral combination for the treatment of NTM lung infections.

Variations of the Mycobacterium abscessus F-ATP synthase subunit a-c interface alter binding and potency of the anti-TB drug bedaquiline.

The anti-tuberculosis therapeutic bedaquiline (BDQ) is used against Mycobacterium abscessus. In M. abscessus BDQ is only bacteriostatic and less potent compared to M. tuberculosis or M. smegmatis. Here we demonstrate its reduced ATP synthesis inhibition against M. abscessus inside-out vesicles, including the F1FO-ATP synthase. Molecular dynamics simulations and binding free energy calculations highlight the differences in drug-binding of the M. abscessus and M. smegmatis FO-domain at the lagging site, where the drug deploys its mechanistic action, inhibiting ATP synthesis. These data pave the way for improved anti-M. abscessus BDQ analogs.

The Molecular Basis for Purine Binding Selectivity in the Bacterial ATP Synthase ϵ Subunit.

The ϵ subunit of ATP synthases has been proposed to regulate ATP hydrolysis in bacteria. Prevailing evidence supports the notion that when the ATP concentration falls below a certain threshold, the ϵ subunit changes its conformation from a non-inhibitory down-state to an extended up-state that then inhibits enzymatic ATP hydrolysis by binding to the catalytic domain. It has been demonstrated that the ϵ subunit from Bacillus PS3 is selective for ATP over other nucleotides, including GTP. In this study, the purine triphosphate selectivity is rationalized by using results from MD simulations and free energy calculations for the R103A/R115A mutant of the ϵ subunit from Bacillus PS3, which binds ATP more strongly than the wild-type protein. Our results are in good agreement with experimental data, and the elucidated molecular basis for selectivity could help to guide the design of novel GTP sensors.

How Ligand Binding Affects the Dynamical Transition Temperature in Proteins.

The biochemical functions of proteins are activated at the protein glass transition temperature, which has been proposed to be dependent upon protein-water interactions. However, at the molecular level it is unclear how ligand binding to well-defined binding sites can influence this transition temperature. We thus report molecular dynamics (MD) simulations of the ϵ subunit from thermophilic Bacillus PS3 in the ATP-free and ligand-bound states over a range of temperatures from 20 to 300 K, to study the influence of ligand association upon the transition temperature. We also measure the protein mean square displacement (MSD) in each state, which is well established as a means to quantify this dynamical temperature dependence. We find that the transition temperature is largely unaffected by ligand association, but the MSD beyond the transition temperature increases more rapidly in the ATP-free state. Our data suggests that ligands can effectively "shield" a binding site from solvent, and hence stabilize protein domains with increasing temperature.

Extending the Martini Coarse-Grained Force Field to N-Glycans.

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.

An Alternative HIV-1 Non-Nucleoside Reverse Transcriptase Inhibition Mechanism: Targeting the p51 Subunit.

The ongoing development of drug resistance in HIV continues to push for the need of alternative drug targets in inhibiting HIV. One such target is the Reverse transcriptase (RT) enzyme which is unique and critical in the viral life cycle-a rational target that is likely to have less off-target effects in humans. Serendipitously, we found two chemical scaffolds from the National Cancer Institute (NCI) Diversity Set V that inhibited HIV-1 RT catalytic activity. Computational structural analyses and subsequent experimental testing demonstrated that one of the two chemical scaffolds binds to a novel location in the HIV-1 RT p51 subunit, interacting with residue Y183, which has no known association with previously reported drug resistance. This finding supports the possibility of a novel druggable site on p51 for a new class of non-nucleoside RT inhibitors that may inhibit HIV-1 RT allosterically. Although inhibitory activity was shown experimentally to only be in the micromolar range, the scaffolds serve as a proof-of-concept of targeting the HIV RT p51 subunit, with the possibility of medical chemistry methods being applied to improve inhibitory activity towards more effective drugs.

On the ATP binding site of the ε subunit from bacterial F-type ATP synthases.

F-type ATP synthases are reversible machinery that not only synthesize adenosine triphosphate (ATP) using an electrochemical gradient across the membrane, but also can hydrolyze ATP to pump ions under certain conditions. To prevent wasteful ATP hydrolysis, subunit ε in bacterial ATP synthases changes its conformation from the non-inhibitory down- to the inhibitory up-state at a low cellular ATP concentration. Recently, a crystal structure of the ε subunit in complex with ATP was solved in a non-biologically relevant dimeric form. Here, to derive the functional ATP binding site motif, we carried out molecular dynamics simulations and free energy calculations. Our results suggest that the ATP binding site markedly differs from the experimental resolved one; we observe a reorientation of several residues, which bind to ATP in the crystal structure. In addition we find that an Mg(2+) ion is coordinated by ATP, replacing interactions of the second chain in the crystal structure. Thus we demonstrate more generally the influence of crystallization effects on ligand binding sites and their respective binding modes. Furthermore, we propose a role for two highly conserved residues to control the ATP binding/unbinding event, which have not been considered before. Additionally our results provide the basis for the rational development of new biosensors based on subunit ε, as shown previously for novel sensors measuring the ATP concentration in cells.

Insights into the ion-coupling mechanism in the MATE transporter NorM-VC.

Bacteria have developed a variety of different mechanisms to defend themselves from compounds that are toxic to them, such as antibiotics. One of these defence mechanisms is the expulsion of drugs or other noxious compounds by multidrug efflux pumps. Multidrug and toxic compound extrusion (MATE) transporters are efflux pumps that extrude metabolic waste and a variety of antibiotics out of the cell, using an ion gradient as energy source. They function via an alternating-access mechanism. When ions bind in the outward facing conformation, a large conformational change to the inward facing conformation is induced, from which the ion is released and the extruded chemical compound is bound. NorM proteins, which are usually coupled to a Na+ gradient, are members of the MATE family. However, for NorM-VC from Vibrio cholerae, it has been shown that this MATE transporter is additionally coupled to protons. How H+ and Na+ binding are coupled mechanistically to enable drug antiport is not well understood. In this study, we use molecular dynamics simulations to illuminate the sequence of ion binding events that enable efflux. Understanding this antiport mechanism is important to support the development of novel compounds that specifically inhibit the functional cycle of NorM transporters.

A role for loop G in the ß1 strand in GABAA receptor activation.

KEY POINTS: The role of the ß1 strand in GABAA receptor function is unclear. It lies anti-parallel to the ß2 strand, which is known to participate in receptor activation. Molecular dynamics simulation revealed solvent accessible residues within the ß1 strand of the GABAA ß3 homopentamer that might be amenable to analysis using the substituted Cys accessibility method. Cys substitutions from Asp43 to Thr47 in the GABAA α1 subunit showed that D43C and T47C reduced the apparent potency of GABA. F45C caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Cys43 and Cys47 were accessible to 2-aminoethyl methanethiosulphonate (MTSEA) modification, whereas Cys45 was not. Both GABA and the allosteric agonist propofol reduced MTSEA modification of Cys43 and Cys47. By contrast, modification of Cys64 in the ß2 strand loop D was impeded by GABA but unaffected by propofol. These data reveal movement of ß1 strand loop G residues during agonist activation of the GABAA receptor. ABSTRACT: The GABAA receptor α subunit ß1 strand runs anti-parallel to the ß2 strand, which contains loop D, known to participate in receptor activation and agonist binding. However, a role for the ß1 strand has yet to be established. We used molecular dynamics simulation to quantify the solvent accessible surface area (SASA) of ß1 strand residues in the GABAA ß3 homopentamer structure. Residues in the complementary interface equivalent to those between Asp43 and Thr47 in the α1 subunit have an alternating pattern of high and low SASA consistent with a ß strand structure. We investigated the functional role of these ß1 strand residues in the α1 subunit by individually replacing them with Cys residues. D43C and T47C substitutions reduced the apparent potency of GABA at α1ß2γ2 receptors by 50-fold and eight-fold, respectively, whereas the F45C substitution caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Receptors with D43C or T47C substitutions were sensitive to 2-aminoethyl methanethiosulphonate (MTSEA) modification. However, GABA-evoked currents mediated by α1(F45C)ß2γ2 receptors were unaffected by MTSEA, suggesting that this residue is inaccessible. Both GABA and the allosteric agonist propofol reduced MTSEA modification of α1(D43C)ß2γ2 and α1(T47C)ß2γ2 receptors, indicating movement of the ß1 strand even during allosteric activation. This is in contrast to α1(F64C)ß2γ2 receptors, where only GABA, but not propofol, reduced MTSEA modification. These findings provide the first functional evidence for movement of the ß1 strand during gating of the receptor and identify residues that are critical for maintaining GABAA receptor function.

On the Mg(2+) binding site of the ε subunit from bacterial F-type ATP synthases.

F-type ATP synthases, central energy conversion machines of the cell synthesize adenosine triphosphate (ATP) using an electrochemical gradient across the membrane and, reversely, can also hydrolyze ATP to pump ions across the membrane, depending on cellular conditions such as ATP concentration. To prevent wasteful ATP hydrolysis, mammalian and bacterial ATP synthases possess different regulatory mechanisms. In bacteria, a low ATP concentration induces a conformational change in the ε subunit from the down- to up-states, which inhibits ATP hydrolysis. Moreover, the conformational change of the ε subunit depends on Mg(2+) concentration in some bacteria such as Bacillus subtilis, but not in others. This diversity makes the ε subunit a potential target for antibiotics. Here, performing molecular dynamics simulations, we identify the Mg(2+) binding site in the ε subunit from B. subtilis as E59 and E86. The free energy analysis shows that the first-sphere bi-dentate coordination of the Mg(2+) ion by the two glutamates is the most stable state. In comparison, we also clarify the reason for the absence of Mg(2+) dependency in the ε subunit from thermophilic Bacillus PS3, despite the high homology to that from B. subtilis. Sequence alignment suggests that this Mg(2+) binding motif is present in the ε subunits of some pathogenic bacteria. In addition we discuss strategies to stabilize an isolated ε subunit carrying the Mg(2+) binding motif by site directed mutagenesis, which also can be used to crystallize Mg(2+) dependent ε subunits in future.

A new type of Na(+)-driven ATP synthase membrane rotor with a two-carboxylate ion-coupling motif.

The anaerobic bacterium Fusobacterium nucleatum uses glutamate decarboxylation to generate a transmembrane gradient of Na⁺. Here, we demonstrate that this ion-motive force is directly coupled to ATP synthesis, via an F1F0-ATP synthase with a novel Na⁺ recognition motif, shared by other human pathogens. Molecular modeling and free-energy simulations of the rotary element of the enzyme, the c-ring, indicate Na⁺ specificity in physiological settings. Consistently, activity measurements showed Na⁺ stimulation of the enzyme, either membrane-embedded or isolated, and ATP synthesis was sensitive to the Na⁺ ionophore monensin. Furthermore, Na⁺ has a protective effect against inhibitors targeting the ion-binding sites, both in the complete ATP synthase and the isolated c-ring. Definitive evidence of Na⁺ coupling is provided by two identical crystal structures of the c11 ring, solved by X-ray crystallography at 2.2 and 2.6 Šresolution, at pH 5.3 and 8.7, respectively. Na⁺ ions occupy all binding sites, each coordinated by four amino acids and a water molecule. Intriguingly, two carboxylates instead of one mediate ion binding. Simulations and experiments demonstrate that this motif implies that a proton is concurrently bound to all sites, although Na⁺ alone drives the rotary mechanism. The structure thus reveals a new mode of ion coupling in ATP synthases and provides a basis for drug-design efforts against this opportunistic pathogen.

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations.

Proteoglycans contain glycosaminoglycans (GAGs) which are negatively charged linear polymers made of repeating disaccharide units of uronic acid and hexosamine units. They play vital roles in numerous physiological and pathological processes, particularly in governing cellular communication and attachment. Depending on their sulfonation state, acetylation, and glycosidic linkages, GAGs belong to different families. The high molecular weight, heterogeneity, and flexibility of GAGs hamper their characterization at atomic resolution, but this may be circumvented via coarse-grained (CG) approaches. In this work, we report a CG model for a library of common GAG types in their isolated or proteoglycan-linked states compatible with version 2.2 (v2.2) of the widely popular CG Martini force field. The model reproduces conformational and thermodynamic properties for a wide variety of GAGs, as well as matching structural and binding data for selected proteoglycan test systems. The parameters developed here may thus be employed to study a range of GAG-containing biomolecular systems, thereby benefiting from the efficiency and broad applicability of the Martini framework.

Exploring bat-inspired cyclic tryptophan diketopiperazines as ABCB1 Inhibitors.

Chemotherapy-induced drug resistance remains a major cause of cancer recurrence and patient mortality. ATP binding cassette subfamily B member 1 (ABCB1) transporter overexpression in tumors contributes to resistance, yet current ABCB1 inhibitors have been unsuccessful in clinical trials. To address this challenge, we propose a new strategy using tryptophan as a lead molecule for developing ABCB1 inhibitors. Our idea stems from our studies on bat cells, as bats have low cancer incidences and high ABCB1 expression. We hypothesized that potential ABCB1 substrates in bats could act as competitive inhibitors in humans. By molecular simulations of ABCB1-substrate interactions, we generated a benzylated Cyclo-tryptophan (C3N-Dbn-Trp2) that inhibits ABCB1 activity with efficacy comparable to or better than the classical inhibitor, verapamil. C3N-Dbn-Trp2 restored chemotherapy sensitivity in drug-resistant human cancer cells with no adverse effect on cell proliferation. Our unique approach presents a promising lead toward developing effective ABCB1 inhibitors to treat drug-resistant cancers.

ATP binding by an F1Fo ATP synthase ε subunit is pH dependent, suggesting a diversity of ε subunit functional regulation in bacteria.

It is a conjecture that the ε subunit regulates ATP hydrolytic function of the F1Fo ATP synthase in bacteria. This has been proposed by the ε subunit taking an extended conformation, with a terminal helix probing into the central architecture of the hexameric catalytic domain, preventing ATP hydrolysis. The ε subunit takes a contracted conformation when bound to ATP, thus would not interfere with catalysis. A recent crystallographic study has disputed this; the Caldalkalibacillus thermarum TA2.A1 F1Fo ATP synthase cannot natively hydrolyse ATP, yet studies have demonstrated that the loss of the ε subunit terminal helix results in an ATP synthase capable of ATP hydrolysis, supporting ε subunit function. Analysis of sequence and crystallographic data of the C. thermarum F1Fo ATP synthase revealed two unique histidine residues. Molecular dynamics simulations suggested that the protonation state of these residues may influence ATP binding site stability. Yet these residues lie outside the ATP/Mg2+ binding site of the ε subunit. We then probed the effect of pH on the ATP binding affinity of the ε subunit from the C. thermarum F1Fo ATP synthase at various physiologically relevant pH values. We show that binding affinity changes 5.9 fold between pH 7.0, where binding is weakest, to pH 8.5 where it is strongest. Since the C. thermarum cytoplasm is pH 8.0 when it grows optimally, this correlates to the ε subunit being down due to ATP/Mg2+ affinity, and not being involved in blocking ATP hydrolysis. Here, we have experimentally correlated that the pH of the bacterial cytoplasm is of critical importance for ε subunit ATP affinity regulated by second-shell residues thus the function of the ε subunit changes with growth conditions.

In-depth biophysical analysis of interactions between therapeutic antibodies and the extracellular domain of the epidermal growth factor receptor.

Targeting of the epidermal growth factor receptor (EGFR) with monoclonal antibodies has become an established antitumor strategy in clinical use or in late stages of drug development. The mAbs effector mechanisms have been widely analyzed based on in vivo or cell studies. Hereby we intend to complement these functional studies by investigating the mAb-EGFR interactions on a molecular level. Surface plasmon resonance, isothermal titration calorimetry, and static light scattering were employed to characterize the interactions of matuzumab, cetuximab, and panitumumab with the extracellular soluble form ecEGFR. The kinetic and thermodynamic determinants dissected the differences in mAbs binding mechanism toward ecEGFR. The quantitative stoichiometric data clearly demonstrated the bivalent binding of the mAbs to two ecEGFR molecules. Our results complement earlier studies on simultaneous binding of cetuximab and matuzumab. The antibodies retain their bivalent binding mode achieving a 1:2:1 complex formation. Interestingly the binding parameters remain nearly constant for the individual antibodies in this ternary assembly. In contrast the binding of panitumumab is almost exclusive either by directly blocking the accessibility for the second antibody or by negative allosteric modulation. Overall we provide a comprehensive biophysical dataset on binding parameters, the complex assembly, and relative epitope accessibility for therapeutic anti-EGFR antibodies.

Microscopic rotary mechanism of ion translocation in the F(o) complex of ATP synthases.

The microscopic mechanism of coupled c-ring rotation and ion translocation in F(1)F(o)-ATP synthases is unknown. Here we present conclusive evidence supporting the notion that the ability of c-rings to rotate within the F(o) complex derives from the interplay between the ion-binding sites and their nonhomogenous microenvironment. This evidence rests on three atomic structures of the c(15) rotor from crystals grown at low pH, soaked at high pH and, after N,N'-dicyclohexylcarbodiimide (DCCD) modification, resolved at 1.8, 3.0 and 2.2 Å, respectively. Alongside a quantitative DCCD-labeling assay and free-energy molecular dynamics calculations, these data demonstrate how the thermodynamic stability of the so-called proton-locked state is maximized by the lipid membrane. By contrast, a hydrophilic environment at the a-subunit-c-ring interface appears to unlock the binding-site conformation and promotes proton exchange with the surrounding solution. Rotation thus occurs as c-subunits stochastically alternate between these environments, directionally biased by the electrochemical transmembrane gradient.

Binding properties of the anti-TB drugs bedaquiline and TBAJ-876 to a mycobacterial F-ATP synthase.

Tuberculosis (TB), the deadly disease caused by Mycobacterium tuberculosis (Mtb), kills more people worldwide than any other bacterial infectious disease. There has been a recent resurgence of TB drug discovery activities, resulting in the identification of a number of novel enzyme inhibitors. Many of these inhibitors target the electron transport chain complexes and the F1FO-ATP synthase; these enzymes represent new target spaces for drug discovery, since the generation of ATP is essential for the bacterial pathogen's physiology, persistence, and pathogenicity. The anti-TB drug bedaquiline (BDQ) targets the Mtb F-ATP synthase and is used as salvage therapy against this disease. Medicinal chemistry efforts to improve the physio-chemical properties of BDQ resulted in the discovery of 3,5-dialkoxypyridine (DARQ) analogs to which TBAJ-876 belongs. TBAJ-876, a clinical development candidate, shows attractive in vitro and in vivo antitubercular activity. Both BDQ and TBAJ-876 inhibit the mycobacterial F1FO-ATP synthase by stopping rotation of the c-ring turbine within the FO domain, thereby preventing proton translocation and ATP synthesis to occur. While structural data for the BDQ bound state are available, no structural information about TBAJ-876 binding have been described. In this study, we show how TBAJ-876 binds to the FO domain of the M. smegmatis F1FO-ATP synthase. We further calculate the binding free energy of both drugs bound to their target and predict an increased affinity of TBAJ-876 for the FO domain. This approach will be useful in future efforts to design new and highly potent DARQ analogs targeting F-ATP synthases of Mtb, nontuberculosis mycobacteria (NTM) as well as the M. leprosis complex.

Bridging the N-terminal and middle domains in FliG of the flagellar rotor.

Flagella are necessary for bacterial movement and contribute to various aspects of virulence. They are complex cylindrical structures built of multiple molecular rings with self-assembly properties. The flagellar rotor is composed of the MS-ring and the C-ring. The FliG protein of the C-ring is central to flagellar assembly and function due to its roles in linking the C-ring with the MS-ring and in torque transmission from stator to rotor. No high-resolution structure of an assembled C-ring has been resolved to date, and the conformation adopted by FliG within the ring is unclear due to variations in available crystallographic data. Here, we use molecular dynamics (MD) simulations to study the conformation and dynamics of FliG in different states of assembly, including both in physiologically relevant and crystallographic lattice environments. We conclude that the linker between the FliG N-terminal and middle domain likely adopts an extended helical conformation in vivo, in contrast with the contracted conformation observed in some previous X-ray studies. We further support our findings with integrative model building of full-length FliG and a FliG ring model that is compatible with cryo-electron tomography (cryo-ET) and electron microscopy (EM) densities of the C-ring. Collectively, our study contributes to a better mechanistic understanding of the flagellar rotor assembly and its function.

Structural and energetic basis for H+ versus Na+ binding selectivity in ATP synthase Fo rotors.

The functional mechanism of the F1Fo ATP synthase, like many membrane transporters and pumps, entails a conformational cycle that is coupled to the movement of H+ or Na+ ions across its transmembrane domain, down an electrochemical gradient. This coupling is an efficient means of energy transduction and regulation, provided that ion binding to the membrane domain, known as Fo, is appropriately selective. In this study we set out to establish the structural and energetic basis for the ion-binding selectivity of the membrane-embedded Fo rotors of two representative ATP synthases. First, we use a biochemical approach to demonstrate the inherent binding selectivity of these rotors, that is, independently from the rest of the enzyme. We then use atomically detailed computer simulations of wild-type and mutagenized rotors to calculate and rationalize their selectivity, on the basis of the structure, dynamics and coordination chemistry of the binding sites. We conclude that H+ selectivity is most likely a robust property of all Fo rotors, arising from the prominent presence of a conserved carboxylic acid and its intrinsic chemical propensity for protonation, as well as from the structural plasticity of the binding sites. In H+-coupled rotors, the incorporation of hydrophobic side chains to the binding sites enhances this inherent H+ selectivity. Size restriction may also favor H+ over Na+, but increasing size alone does not confer Na+ selectivity. Rather, the degree to which Fo rotors may exhibit Na+ coupling relies on the presence of a sufficient number of suitable coordinating side chains and/or structural water molecules. These ligands accomplish a shift in the relative binding energetics, which under some physiological conditions may be sufficient to provide Na+ dependence.

A second shell residue modulates a conserved ATP-binding site with radically different affinities for ATP.

BACKGROUND: Prediction of ligand binding and design of new function in enzymes is a time-consuming and expensive process. Crystallography gives the impression that proteins adopt a fixed shape, yet enzymes are functionally dynamic. Molecular dynamics offers the possibility of probing protein movement while predicting ligand binding. Accordingly, we choose the bacterial F1Fo ATP synthase ε subunit to unravel why ATP affinity by ε subunits from Bacillus subtilis and Bacillus PS3 differs ~500-fold, despite sharing identical sequences at the ATP-binding site. METHODS: We first used the Bacillus PS3 ε subunit structure to model the B. subtilis ε subunit structure and used this to explore the utility of molecular dynamics (MD) simulations to predict the influence of residues outside the ATP binding site. To verify the MD predictions, point mutants were made and ATP binding studies were employed. RESULTS: MD simulations predicted that E102 in the B. subtilis ε subunit, outside of the ATP binding site, influences ATP binding affinity. Engineering E102 to alanine or arginine revealed a ~10 or ~54 fold increase in ATP binding, respectively, confirming the MD prediction that E102 drastically influences ATP binding affinity. CONCLUSIONS: These findings reveal how MD can predict how changes in the "second shell" residues around substrate binding sites influence affinity in simple protein structures. Our results reveal why seemingly identical ε subunits in different ATP synthases have radically different ATP binding affinities. GENERAL SIGNIFICANCE: This study may lead to greater utility of molecular dynamics as a tool for protein design and exploration of protein design and function.