Structural-Energetic Basis for Coupling between Equilibrium Fluctuations and Phosphorylation in a Protein Native Ensemble

A finely balanced order-disorder equilibrium sculpts the folding-binding landscape of an antibiotic sequestering protein

TipA, a MerR family transcription factor from Streptomyces lividans, promotes antibiotic resistance by sequestering broad-spectrum thiopeptide-based antibiotics, thus counteracting their inhibitory effect on ribosomes. TipAS, a minimal binding motif which is expressed as an isoform of TipA, harbors a partially disordered N-terminal subdomain that folds upon binding multiple antibiotics. The extent and nature of the underlying molecular heterogeneity in TipAS that shapes its promiscuous folding-function landscape is an open question and is critical for understanding antibiotic-sequestration mechanisms. Here, combining equilibrium and time-resolved experiments, statistical modeling, and simulations, we show that the TipAS native ensemble exhibits a pre-equilibrium between binding-incompetent and binding-competent substates, with the fully folded state appearing only as an excited state under physiological conditions. The binding-competent state characterized by a partially structured N-terminal subdomain loses structure progressively in the physiological range of temperatures, swells on temperature increase, and displays slow conformational exchange across multiple conformations. Binding to the bactericidal antibiotic thiostrepton follows a combination of induced-fit and conformational-selection-like mechanisms, via partial binding and concomitant stabilization of the binding-competent substate. These ensemble features are evolutionarily conserved across orthologs from select bacteria that infect humans, underscoring the functional role of partial disorder in the native ensemble of antibiotic-sequestering proteins belonging to the MerR family.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

Phosphorylation-Competent Metastable State of Escherichia coli Toxin HipA

Phosphorylation is a key post-translational modification that alters the functional state of many proteins. The Escherichia coli toxin HipA, which phosphorylates glutamyl-tRNA synthetase and triggers bacterial persistence under stress, becomes inactivated upon autophosphorylation of Ser150. Interestingly, Ser150 is phosphorylation-incompetent in the crystal structure of HipA since it is deeply buried ("in-state"), although in the phosphorylated state it is solvent exposed ("out-state"). To be phosphorylated, a minor population of HipA must exist in the phosphorylation-competent "out-state" (solvent-exposed Ser150), not detected in the crystal structure of unphosphorylated HipA. Here we report a molten-globule-like intermediate of HipA at low urea (∼4 kcal/mol unstable than natively folded HipA). The intermediate is aggregation-prone, consistent with a solvent exposed Ser150 and its two flanking hydrophobic neighbors (Val/Ile) in the "out-state". Molecular dynamics simulations showed the HipA "in-out" pathway to contain multiple free energy minima with an increasing degree of Ser150 solvent exposure with the free energy difference between the "in-state" and the metastable exposed state(s) to be ∼2-2.5 kcal/mol, with unique sets of hydrogen bonds and salt bridges associated with the metastable loop conformations. Together, the data clearly identify the existence of a phosphorylation-competent metastable state of HipA. Our results not only suggest a mechanism of HipA autophosphorylation but also add to a number of recent reports on unrelated protein systems where the common proposed mechanism for phosphorylation of buried residues is their transient exposure even without phosphorylation.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.