Assessing the chemical accuracy of protein structures via peptide acidity

Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states

Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.

Structural dynamics of pentapeptide repeat proteins

Pentapeptide repeat proteins (PRPs) represent a large superfamily with more than 38 000 sequences in nearly 3500 species, the majority belonging to cyanobacteria but represented among all branches of life. PRPs contain at least eight consecutive pentapeptide repeats with the consensus (A/C/S/V/T/L/I)(D/N/S/K/E/I/R)(L/F)(S/T/R/E/Q/K/V/D)(G/D/E/N/R/Q/K). PRPs fold into right-handed quadrilateral β helices, also known as repeat-five-residue (Rfr)-folds, with four consecutive pentapeptide repeats comprising a single coil, the ~90° change in polypeptide direction in square-shaped coils achieved by type I, II and IV β turns, and hydrogen bonds between coils establishing β ladders on each Rfr-fold face. PRPs are broadly categorized into group 1 and 2 involved in antibiotic resistance and group 3 currently having unknown functions. Motivated by their intriguing structures, we are investigating PRP biophysical characteristics, including Rfr-fold thermal stability, β turn and β ladder hydrogen bond amide exchange rates and backbone dynamics. Here, we present analysis of 20 ns molecular dynamics (MD) simulations and all atom normal mode analysis (aaNMA) calculations for four group 1 and group 2 and four group 3 PRPs whose structures have been determined by X-ray crystallography. The MD cross-correlation matrices and aaNMA indicated strong correlated motion between adjacent coils and weak coupled motion between coils separated by one or more intervening coils. Slow anticorrelated motions were detected between adjacent coils in aaNMA modes that we hypothesize are requisite to access exchange-competent states necessary to permit solvent exchange of amide hydrogens involved in β-ladder and β-turns hydrogen bonds, which can have lifetimes on the order of months.

Calcium-Mediated Control of S100 Proteins: Allosteric Communication via an Agitator/Signal Blocking Mechanism

Allosteric proteins possess dynamically coupled residues for the propagation of input signals to distant target binding sites. The input signals usually correspond to "effector is present" or "effector is not present". Many aspects of allosteric regulation remain incompletely understood. This work focused on S100A11, a dimeric EF-hand protein with two hydrophobic target binding sites. An annexin peptide (Ax) served as the target. Target binding is allosterically controlled by Ca²⁺ over a distance of ∼26 Å. Ca²⁺ promotes formation of a [Ca₄ S100 Ax₂] complex, where the Ax peptides are accommodated between helices III/IV and III'/IV'. Without Ca²⁺ these binding sites are closed, precluding interactions with Ax. The allosteric mechanism was probed by microsecond MD simulations in explicit water, complemented by hydrogen exchange mass spectrometry (HDX/MS). Consistent with experimental data, MD runs in the absence of Ca²⁺ and Ax culminated in target binding site closure. In simulations on [Ca₄ S100] the target binding sites remained open. These results capture the essence of allosteric control, revealing how Ca²⁺ prevents binding site closure. Both HDX/MS and MD data showed that the metalation sites become more dynamic after Ca²⁺ loss. However, these enhanced dynamics do not represent the primary trigger of the allosteric cascade. Instead, a labile salt bridge acts as an incessantly active "agitator" that destabilizes the packing of adjacent residues, causing a domino chain of events that culminates in target binding site closure. This agitator represents the starting point of the allosteric signal propagation pathway. Ca²⁺ binding rigidifies elements along this pathway, thereby blocking signal transmission. This blocking mechanism does not conform to the commonly held view that allosteric communication pathways generally originate at the sites where effectors interact with the protein.

Analysing the visible conformational substates of the FK506-binding protein FKBP12

The ¹H-¹⁵N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80′s loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50′s loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~10⁵-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80′s loop. Line-broadening was also decreased for the residues in the α-helix and 50′s loop, whereas line-broadening in the 40′s loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40′s loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.

Analysing the visible conformational substates of the FK506-binding protein FKBP12

The 1H-15N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80's loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50's loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~105-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80's loop. Line-broadening was also decreased for the residues in the α-helix and 50's loop, whereas line-broadening in the 40's loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40's loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.