skip to main content


Title: Deep Analysis of Residue Constraints (DARC): identifying determinants of protein functional specificity
Abstract

Protein functional constraints are manifest as superfamily and functional-subgroup conserved residues, and as pairwise correlations. Deep Analysis of Residue Constraints (DARC) aids the visualization of these constraints, characterizes how they correlate with each other and with structure, and estimates statistical significance. This can identify determinants of protein functional specificity, as we illustrate for bacterial DNA clamp loader ATPases. These load ring-shaped sliding clamps onto DNA to keep polymerase attached during replication and contain one δ, three γ, and one δ’ AAA+ subunits semi-circularly arranged in the order δ-γ123-δ’. Only γ is active, though both γ and δ’ functionally influence an adjacent γ subunit. DARC identifies, as functionally-congruent features linking allosterically the ATP, DNA, and clamp binding sites: residues distinctive of γ and of γ/δ’ that mutually interact in trans, centered on the catalytic base; several γ/δ’-residues and six γ/δ’-covariant residue pairs within the DNA binding N-termini of helices α2 and α3; and γ/δ’-residues associated with the α2 C-terminus and the clamp-binding loop. Most notable is a trans-acting γ/δ’ hydroxyl group that 99% of other AAA+ proteins lack. Mutation of this hydroxyl to a methyl group impedes clamp binding and opening, DNA binding, and ATP hydrolysis—implying a remarkably clamp-loader-specific function.

 
more » « less
Award ID(s):
1817869
NSF-PAR ID:
10154198
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
10
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Brine shrimp (Artemia) are the only animals to thrive at sodium concentrations above 4 M. Salt excretion is powered by the Na+,K+-ATPase (NKA), a heterodimeric (αβ) pump that usually exports 3Na+in exchange for 2 K+per hydrolyzed ATP.Artemiaexpress several NKA catalytic α-subunit subtypes. High-salinity adaptation increases abundance of α2KK, an isoform that contains two lysines (Lys308 and Lys758 in transmembrane segments TM4 and TM5, respectively) at positions where canonical NKAs have asparagines (Xenopusα1’s Asn333 and Asn785). Using de novo transcriptome assembly and qPCR, we found thatArtemiaexpress two salinity-independent canonical α subunits (α1NNand α3NN), as well as two β variants, in addition to the salinity-controlled α2KK. These β subunits permitted heterologous expression of the α2KKpump and determination of its CryoEM structure in a closed, ion-free conformation, showing Lys758 residing within the ion-binding cavity. We used electrophysiology to characterize the function of α2KKpumps and compared it to that ofXenopusα1 (and its α2KK-mimicking single- and double-lysine substitutions). The double substitution N333K/N785K confers α2KK-like characteristics toXenopusα1, and mutant cycle analysis reveals energetic coupling between these two residues, illustrating how α2KK’s Lys308 helps to maintain high affinity for external K+when Lys758 occupies an ion-binding site. By measuring uptake under voltage clamp of the K+-congener86Rb+, we prove that double-lysine-substituted pumps transport 2Na+and 1 K+per catalytic cycle. Our results show how the two lysines contribute to generate a pump with reduced stoichiometry allowingArtemiato maintain steeper Na+gradients in hypersaline environments.

     
    more » « less
  2. Many viruses utilize ringed packaging ATPases to translocate double-stranded DNA into procapsids during replication. A critical step in the mechanochemical cycle of such ATPases is ATP binding, which causes a subunit within the motor to grip DNA tightly. Here, we probe the underlying molecular mechanism by which ATP binding is coupled to DNA gripping and show that a glutamate-switch residue found in AAA+ enzymes is central to this coupling in viral packaging ATPases. Using free-energy landscapes computed through molecular dynamics simulations, we determined the stable conformational state of the ATPase active site in ATP- and ADP-bound states. Our results show that the catalytic glutamate residue transitions from an active to an inactive pose upon ATP hydrolysis and that a residue assigned as the glutamate switch is necessary for regulating this transition. Furthermore, we identified via mutual information analyses the intramolecular signaling pathway mediated by the glutamate switch that is responsible for coupling ATP binding to conformational transitions of DNA-gripping motifs. We corroborated these predictions with both structural and functional experimental measurements. Specifically, we showed that the crystal structure of the ADP-bound P74-26 packaging ATPase is consistent with the structural coupling predicted from simulations, and we further showed that disrupting the predicted signaling pathway indeed decouples ATPase activity from DNA translocation activity in the φ29 DNA packaging motor. Our work thus establishes a signaling pathway that couples chemical and mechanical events in viral DNA packaging motors.

     
    more » « less
  3. Hsp104 is a large AAA+ molecular machine that can rescue proteins trapped in amorphous aggregates and stable amyloids by drawing substrate protein into its central pore. Recent cryo-EM studies image Hsp104 at high resolution. We used hydrogen exchange mass spectrometry analysis (HX MS) to resolve and characterize all of the functionally active and inactive elements of Hsp104, many not accessible to cryo-EM. At a global level, HX MS confirms the one noncanonical interprotomer interface in the Hsp104 hexamer as a marker for the spiraled conformation revealed by cryo-EM and measures its fast conformational cycling under ATP hydrolysis. Other findings enable reinterpretation of the apparent variability of the regulatory middle domain. With respect to detailed mechanism, HX MS determines the response of each Hsp104 structural element to the different bound adenosine nucleotides (ADP, ATP, AMPPNP, and ATPγS). They are distinguished most sensitively by the two Walker A nucleotide-binding segments. Binding of the ATP analog, ATPγS, tightly restructures the Walker A segments and drives the global open-to-closed/extended transition. The global transition carries part of the ATP/ATPγS-binding energy to the somewhat distant central pore. The pore constricts and the tyrosine and other pore-related loops become more tightly structured, which seems to reflect the energy-requiring directional pull that translocates the substrate protein. ATP hydrolysis to ADP allows Hsp104 to relax back to its lowest energy open state ready to restart the cycle.

     
    more » « less
  4. The F-ATP synthase, consisting of F 1 and F O motors connected by a central rotor and the stators, is the enzyme responsible for synthesizing the majority of ATP in all organisms. The F 1 (αβ) 3 ring stator contains three catalytic sites. Single-molecule F 1 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 F 1 F O , the membrane-bound F O 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 F 1 . Single-molecule studies of F 1 F O embedded in lipid bilayer nanodisks showed that the c-ring transiently stopped F 1 -ATPase-driven rotation every 36° (at each c-subunit in the c 10 -ring of E. coli F 1 F O ) 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. 
    more » « less
  5. Abstract Motivation

    The analysis of sequence conservation patterns has been widely utilized to identify functionally important (catalytic and ligand-binding) protein residues for over a half-century. Despite decades of development, on average state-of-the-art non-template-based functional residue prediction methods must predict ∼25% of a protein’s total residues to correctly identify half of the protein’s functional site residues. The overwhelming proportion of false positives results in reported ‘F-Scores’ of ∼0.3. We investigated the limits of current approaches, focusing on the so-far neglected impact of the specific choice of homologs included in multiple sequence alignments (MSAs).

    Results

    The limits of conservation-based functional residue prediction were explored by surveying the binding sites of 1023 proteins. A straightforward conservation analysis of MSAs composed of randomly selected homologs sampled from a PSI-BLAST search achieves average F-Scores of ∼0.3, a performance matching that reported by state-of-the-art methods, which often consider additional features for the prediction in a machine learning setting. Interestingly, we found that a simple combinatorial MSA sampling algorithm will in almost every case produce an MSA with an optimal set of homologs whose conservation analysis reaches average F-Scores of ∼0.6, doubling state-of-the-art performance. We also show that this is nearly at the theoretical limit of possible performance given the agreement between different binding site definitions. Additionally, we showcase the progress in this direction made by Selection of Alignment by Maximal Mutual Information (SAMMI), an information-theory-based approach to identifying biologically informative MSAs. This work highlights the importance and the unused potential of optimally composed MSAs for conservation analysis.

    Supplementary information

    Supplementary data are available at Bioinformatics online.

     
    more » « less