Bertrand García-Moreno E

Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrCr). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pKa values of the active site residues suggest an increase in the pKa of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009.

Molecular mechanisms of pH-driven conformational transitions of proteins: Insights from continuum electrostatics calculations of acid unfolding

The acid unfolding of staphylococcal nuclease (SNase) is very cooperative (Whitten and García‐Moreno, Biochemistry 2000;39:14292–14304). As many as seven hydrogen ions (H+) are bound preferentially by the acid‐unfolded state relative to the native (N) state in the pH range 3.2–3.9. To investigate the mechanism of acid unfolding, structure‐based pKa calculations were performed with a variety of continuum electrostatic methods. The calculations reproduced successfully the H+ binding properties of the N state between pH 5 and 9, but they systematically overestimated the number of H+ bound upon acid unfolding. The calculated pKa values of all carboxylic residues in the N state were more depressed than they should be. The discrepancy between the observed and the calculated H+ uptake upon acid unfolding was not improved by using high protein dielectric constants, structures relaxed with molecular dynamics, or other empirical modifications implemented previously by others to maximize agreement between measured and calculated pKa values. This suggests an important role for conformational fluctuations of the backbone as important determinants of pKa values of carboxylic groups. Because no global or subglobal conformational changes have been observed previously for SNase under acidic conditions above the acid‐unfolding region, these fluctuations must be local. The acid unfolding of SNase does not seem to involve the disruption of the N state by accruement of intramolecular repulsive interactions, nor the protonation of key ion paired carboxylic residues. It is more consistent with modest contributions from many H+ binding groups, with an important role for local conformational fluctuations in the coupling between H+ binding and the global structural transition. Proteins 2006.

Electrostatic effects in unfolded staphylococcal nuclease

Structure‐based calculations of pK a values and electrostatic free energies of proteins assume that electrostatic effects in the unfolded state are negligible. In light of experimental evidence showing that this assumption is invalid for many proteins, and with increasing awareness that the unfolded state is more structured and compact than previously thought, a detailed examination of electrostatic effects in unfolded proteins is warranted. Here we address this issue with structure‐based calculations of electrostatic interactions in unfolded staphylococcal nuclease. The approach involves the generation of ensembles of structures representing the unfolded state, and calculation of Coulomb energies to Boltzmann weight the unfolded state ensembles. Four different structural models of the unfolded state were tested. Experimental proton binding data measured with a variant of nuclease that is unfolded under native conditions were used to establish the validity of the calculations. These calculations suggest that weak Coulomb interactions are an unavoidable property of unfolded proteins. At neutral pH, the interactions are too weak to organize the unfolded state; however, at extreme pH values, where the protein has a significant net charge, the combined action of a large number of weak repulsive interactions can lead to the expansion of the unfolded state. The calculated pK a values of ionizable groups in the unfolded state are similar but not identical to the values in small peptides in water. These studies suggest that the accuracy of structure‐based calculations of electrostatic contributions to stability cannot be improved unless electrostatic effects in the unfolded state are calculated explicitly.

Human monoclonal antibody stability and activity at vaginal pH

Antibodies can be delivered topically to the vagina to protect against pregnancy and sexually transmitted infections, but the acidity of vaginal secretions (pH 3.5-4.5) might inactivate them. To address this question, both experimental and computational methods were used to evaluate the effects of pH on human monoclonal antibody (MAb) stability and activity. To determine the acid-sensitivity of their antigen binding sites, human MAbs against human sperm (H6-3C4) and gp120 of HIV (1511) were tested by ELISA for binding to human sperm and recombinant gp120, respectively, at pH 3.0-7.0, after storing them for 1 or 20 h at the same pH. Binding was unaltered by acidic pH> or =4 even after 20 h, and at pH 3.5 both MAbs retained > or =40% antigen binding activity. A humanized MAb against HSV-2 glycoprotein B expressed both in Chinese hamster ovary (CHO) cells and in soybean cells was incubated for 1 or 24 h at pH 3.5-7.6, brought to neutral pH, and tested for ability to block HSV-2 infection of foreskin fibroblast cells. Loss in blocking activity occurred only when antibodies were incubated at pH 3.5 for 24 h and was independent of the expression cell type. Using empirical structure-based methods, net charge, Z, and electrostatic contributions to free energy, DeltaDeltaG(el), were calculated as a function of pH for 1 human and 8 murine F(ab)s. The calculations indicate that Z changes slowly between pH 5.0 and 9.0 and that DeltaDeltaG(el) is nearly constant between pH 4.0 and 10 for all the F(ab)s and, therefore, human antibodies should remain stable in this pH range. Taken together, our data and empirical calculations suggest that vaginally applied human MAbs are likely to remain stable and active throughout the duration they are likely to reside in the vagina.

Thermodynamic principles for the engineering of pH-driven conformational switches and acid insensitive proteins.

The general thermodynamic principles behind pH driven conformational transitions of biological macromolecules are well understood. What is less obvious is how they can be used to engineer pH switches in proteins. The acid unfolding of staphylococcal nuclease (SNase) was used to illustrate different factors that can affect pH-driven conformational transitions. Acid unfolding is a structural transition driven by preferential H(+) binding to the acid unfolded state (U) over the native (N) state of a protein. It is the result of carboxylic groups that titrate with more normal pK(a) values in the U state than in the N state. Acid unfolding profiles of proteins reflect a balance between electrostatic and non-electrostatic contributions to stability. Several strategies were used in attempts to turn SNase into an acid insensitive protein: (1) enhancing global stability of the protein with mutagenesis or with osmolytes, (2) use of high salt concentrations to screen Coulomb interactions, (3) stabilizing the N state through specific anion effects, (4) removing Asp or Glu residues that titrate with depressed pK(a) values in the N state, and (5) removing basic residues that might have strong repulsive interactions in the N state at low pH. The only effective way to engineer acid resistance in SNase is not through modulation of pK(a) values of Asp/Glu but by enhancing the global stability of the protein. Modulation of pH-driven conformational transitions by selective manipulation of the electrostatic component of the switch is an extremely difficult undertaking.

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pKa values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pKa values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pKa in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pKa values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. 15N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pKa calculations: (1) The timescale of the backbones response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pKa values of ionizable residues in proteins. Proteins 2014; 82:3132–3143.

Charges in Hydrophobic Environments: A Strategy for Identifying Alternative States in Proteins

In the V23E variant of staphylococcal nuclease, Glu-23 has a pKa of 7.5. At low pH, Glu-23 is neutral and buried in the hydrophobic interior of the protein. Crystal structures and NMR spectroscopy experiments show that when Glu-23 becomes charged, the protein switches into an open state in which strands β1 and β2 separate from the β-barrel; the remaining structure is unaffected. In the open state the hydrophobic interior of the protein is exposed to bulk water, allowing Glu-23 to become hydrated. This illustrates several key aspects of protein electrostatics: (1) The apparent pKa of an internal ionizable group can reflect the average of the very different pKa values (open ≈4.5, closed ≫7.5) sampled in the different conformational states. (2) The high apparent dielectric constant reported by the pKa value of internal ionizable group reflects conformational reorganization. (3) The apparent pKa of internal groups can be governed by large conformational changes. (4) A single charge buried in the hydrophobic interior of a protein is sufficient to convert what might have been a transient, partially unfolded state into the dominant state in solution. This suggests a general strategy for examining inaccessible regions of the folding landscape and for engineering conformational switches driven by small changes in pH. These data also constitute a benchmark for stringent testing of the ability of computational algorithms to predict pKa values of internal residues and to reproduce pH-driven conformational transitions of proteins.