Griselda Hernandez

A Billion-fold Range in Acidity for the Solvent-Exposed Amides of Pyrococcus furiosus Rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.

Conformational Dynamics in FKBP Domains: Relevance to Molecular Signaling and Drug Design

Among the 22 FKBP domains in the human genome, FKBP12.6 and the first FKBP domains (FK1) of FKBP51 and FKBP52 are evolutionarily and structurally most similar to the archetypical FKBP12. As such, the development of inhibitors with selectivity among these four FKBP domains poses a significant challenge for structure-based design. The pleiotropic effects of these FKBP domains in a range of signaling processes such as the regulation of ryanodine receptor calcium channels by FKBP12 and FKBP12.6 and steroid receptor regulation by the FK1 domains of FKBP51 and FKBP52 amply justify the efforts to develop selective therapies. In contrast to their close structural similarities, these four FKBP domains exhibit a substantial diversity in their conformational flexibility. A number of distinct conformational transitions have been characterized for FKBP12 spanning timeframes from 20 s to 10 ns and in each case these dynamics have been shown to markedly differ from the conformational behavior for one or more of the other three FKBP domains. Protein flexibilitybased inhibitor design could draw upon the transitions that are significantly populated in only one of the targeted proteins. Both the similarities and differences among these four proteins valuably inform the understanding of how dynamical effects propagate across the FKBP domains as well as potentially how such intramolecular transitions might couple to the larger scale transitions that are central to the signaling complexes in which these FKBP domains function.

Prediction of Bond Vector Autocorrelation Functions from Larmor Frequency-Selective Order Parameter Analysis of NMR Relaxation Data

Protein molecular dynamics interpretation of the standard R1, R2, and heteronuclear NOE relaxation measurements has typically been limited to a single S2 order parameter which is often insufficient to characterize the rich content of these NMR experiments. In the absence of exchange linebroadening, an optimized reduced spectral density analysis of these measurements can yield spectral density values at three distinct frequencies. Surprisingly, these three discrete spectral density values have proven to be sufficient for a Larmor frequency-selective order parameter analysis of the 223 methine and methylene H-C bonds of the B3 domain of Protein G (GB3) to accurately back-calculate the entire curve of the corresponding bond vector autocorrelation functions upon which the NMR relaxation behavior depends. The 13C relaxation values calculated from 2 μs of CHARMM36 simulation trajectories yielded the corresponding autocorrelation functions to an average rmsd of 0.44% with only three bond vectors having rmsd errors slightly greater than 1.0%. Similar quality predictions were obtained using the CHARMM22/CMAP, AMBER ff99SB, and AMBER ff99SB-ILDN force fields. Analogous predictions for the backbone 15N relaxation values were 3-fold more accurate. Excluding seven residues for which either experimental data is lacking or previous MD studies have indicated markedly divergent dynamics predictions, the CHARMM36-derived and experimentally derived 15N relaxation values for the remaining 48 amides of GB3 agree to an average of 0.016, 0.010, and 0.020 for the fast limit (Sf2) and Larmor frequency-selective (SH2 and SN2) order parameters, respectively. In contrast, for a substantial fraction of side chain positions, the statistical uncertainties obtained in the relaxation value predictions from each force field were appreciably less than the much larger differences predicted among these force fields, indicating a significant opportunity for experimental NMR relaxation measurements to provide structurally interpretable guidance for further optimizing the prediction of protein dynamics.

Nuclear Magnetic Relaxation Dispersion Measurements of Ethylenediamine-Tetraacetatoiron(III) Complexes at High pH and Ionic Strength

Abstract Nuclear magnetic relaxation dispersion measurements are reported for aqueous solutions of iron(III)-EDTA complexes at high pH and high ionic strength where the dominant species is the μ-oxo dimer. Formation of the dimeric species leaves only the outer coordination sphere contribution to the paramagnetic relaxation rate of the water protons in solution. The relaxation efficiency is further decreased by the decrease in the effective magnetic moment due to antiferromagnetic coupling. This result is in distinction to dimeric species that form with Gd(HEDTA) complexes at high pH in which antiferromagnetic coupling has not been observed to decrease the effective magnetic moment.