Quantification of CH and NH/π-Stacking Interactions in Cells Using Nuclear Magnetic Resonance Spectroscopy.

π-Stacking, a type of noncovalent interactions involving aromatic residues, plays an important role in protein folding and function. In this work, an attempt has been made to measure CH/π and NH/π stacking interactions in a protein in Escherichia coli cells using a combined double-mutant cycle and nuclear magnetic resonance spectroscopy method. The results show that the CH/π and NH/π stacking interactions are generally weaker in cells than those in the buffer. The transient intermolecular noncovalent interactions between the protein and the complex cellular environment may compete with and thus weaken the stacking interactions in the protein. The weakening of stacking interactions can enhance the local conformational opening of proteins in E. coli cells. This is evident from the faster rates of amide hydrogen/deuterium exchange observed in cells than in the buffer, for residues that undergo local conformational opening. This study highlights the influence of the cellular environment on π-stacking and the conformational dynamics of proteins.

Protein Allostery Study in Cells Using NMR Spectroscopy.

Protein allostery is commonly observed in vitro. But how protein allostery behaves in cells is unknown. In this work, a protein monomer-dimer equilibrium system was built with the allosteric effect on the binding characterized using NMR spectroscopy through mutations away from the dimer interface. A chemical shift linear fitting method was developed that enabled us to accurately determine the dissociation constant. A total of 28 allosteric mutations were prepared and grouped to negative allosteric, nonallosteric, and positive allosteric modulators. ∼ 50% of mutations displayed the allosteric-state changes when moving from a buffered solution into cells. For example, there were no positive allosteric modulators in the buffered solution but eight in cells. The change in protein allostery is correlated with the interactions between the protein and the cellular environment. These interactions presumably drive the surrounding macromolecules in cells to transiently bind to the monomer and dimer mutational sites and change the free energies of the two species differently which generate new allosteric effects. These surrounding macromolecules create a new protein allostery pathway that is only present in cells.

Osmolytes Can Destabilize Proteins in Cells by Modulating Electrostatics and Quinary Interactions.

Although numerous in vitro studies have shown that osmolytes are capable of stabilizing proteins, their effect on protein folding in vivo has been less understood. In this work, we investigated the effect of osmolytes, including glycerol, sorbitol, betaine, and taurine, on the folding of a protein GB3 variant in E. coli cells using NMR spectroscopy. 400 mM osmolytes were added to E. coli cells; only glycerol stabilizes the folded protein, whereas betaine and taurine considerably destabilize the protein through modulating folding and unfolding rates. Further investigation indicates that betaine and taurine can enhance the quinary interaction between the protein and cellular environment and manifestly weaken the electrostatic attraction in protein salt bridges. The combination of the two factors causes destabilization of the protein in E. coli cells. These factors counteract the preferential exclusion mechanism that is adopted by osmolytes to stabilize proteins.

Generating Five Independent Molecular Alignments for Simultaneous Protein Structure and Dynamics Determination Using Nuclear Magnetic Resonance Spectroscopy.

Residual dipolar couplings (RDCs) are commonly used in NMR for protein structure and dynamics studies, but it is challenging to generate five independent RDC data sets (required for simultaneous structure and dynamics determination) for most protein molecules in the magnetic field. In this work, a reporter protein with a lanthanide tag is introduced to create five independent alignments. This reporter protein is then attached to target proteins where five independent sets of RDCs are also obtained for the target proteins. The fitting of RDCs provides important information about the structure and dynamics of the target proteins. The method is simple and effective and, in principle, can be used to generate complete sets of RDCs for different protein molecules.

Polyol and sugar osmolytes can shorten protein hydrogen bonds to modulate function.

Polyol and sugar osmolytes are commonly used in therapeutic protein formulations. How they may affect protein structure and function is an important question. In this work, through NMR measurements, we show that glycerol and sorbitol (polyols), as well as glucose (sugar), can shorten protein backbone hydrogen bonds. The hydrogen bond shortening is also captured by molecular dynamics simulations, which suggest a hydrogen bond competition mechanism. Specifically, osmolytes weaken hydrogen bonds between the protein and solvent to strengthen those within the protein. Although the hydrogen bond change is small, with the average experimental cross hydrogen bond 3hJNC' coupling of two proteins GB3 and TTHA increased by ~ 0.01 Hz by the three osmolytes (160 g/L), its effect on protein function should not be overlooked. This is exemplified by the PDZ3-peptide binding where several intermolecular hydrogen bonds are formed and osmolytes shift the equilibrium towards the bound state.

Arginine side chain stacking with peptide plane stabilizes the protein helix conformation in a cooperative way.

A combined experimental and computational study is performed for arginine side chain stacking with the protein α-helix. Theremostability measurements of Aristaless homeodomain, a helical protein, suggest that mutating the arginine residue R106, R137 or R141, which has the guanidino side chain stacking with the peptide plane, to alanine, destabilizes the protein. The R-PP stacking has an energy of ∼0.2-0.4 kcal/mol. This stacking interaction mainly comes from dispersion and electrostatics, based on MP2 calculations with the energy decomposition analysis. The calculations also suggest that the stacking stabilizes 2 backbone-backbone h-bonds (i→i-4 and i-3→i-7) in a cooperative way. Desolvation and electrostatic polarization are responsible for cooperativity with the i→i-4 and i-3→i-7 h-bonds, respectively. This cooperativity is supported by a protein α-helices h-bond survey in the pdb databank where stacking shortens the corresponding h-bond distances.

Direct Observation of CH/CH van der Waals Interactions in Proteins by NMR.

van der Waals interactions are important to protein stability and function. These interactions are usually identified empirically based on protein 3D structures. In this work, we performed a solution nuclear magnetic resonance (NMR) spectroscopy study of van der Waals interactions by detecting the through-space vdw JCC-coupling between protein aliphatic side chain groups. Specifically, vdw JCC-coupling values up to ∼0.5 Hz were obtained between the methyl and nearby aliphatic groups in protein GB3, providing direct experimental evidence for the van der Waals interactions. Quantum mechanical calculations suggest that the J-coupling is correlated with the exchange-repulsion term of van der Waals interaction. NMR detection of vdw JCC-coupling offers a new tool to characterize such interactions in proteins.

Entropy Drives the Formation of Salt Bridges in the Protein GB3.

Salt bridges are very common in proteins. But what drives the formation of protein salt bridges is not clear. In this work, we determined the strength of four salt bridges in the protein GB3 by measuring the ΔpKa values of the basic residues that constitute the salt bridges with a highly accurate NMR titration method at different temperatures. The results show that the ΔpKa values increase with temperature, thus indicating that the salt bridges are stronger at higher temperatures. Fitting of ΔpKa values to the van't Hoff equation yields positive ΔH and ΔS values, thus indicating that entropy drives salt-bridge formation. Molecular dynamics simulations show that the protein and solvent make opposite contributions to ΔH and ΔS. Specifically, the enthalpic gain contributed from the protein is more than offset by the enthalpic loss contributed from the solvent, whereas the entropic gain originates from the desolvation effect.

Quinary Interactions Weaken the Electric Field Generated by Protein Side-Chain Charges in the Cell-like Environment.

The intramolecular electric field (e-field) generated by protein GB3 side-chain charges K/E10, K/E19, and D/K40 was measured in the absence or presence of macromolecular crowding. The e-field responds differently to different crowding agents-dextran, Ficoll, BSA, and E. coli cell lysate. Dextran and Ficoll have no effect on the e-field. The lysate generally weakens the e-field but the amplitude of weakening varies greatly. For example, the e-field by K19 is reduced by 67% in the presence of 90 g/L lysate, corresponding to a charge change from 0.9 to 0.3 e for K19, whereas the e-fields by D/K40 are weakened only by ∼7% under the same lysate concentration. The extent of the e-field weakening by BSA is in between that by Ficoll (dextran) and lysate. Further investigations suggest that the e-field weakening mechanism by lysate is similar to that by NaCl. That is, the e-field generated by a protein surface charge affects the distribution of lysate which creates a reaction field and weakens the protein e-field. Our study indicates that the protein electrostatic property can be changed significantly due to quinary interaction with the cell environment.

Carboxyl-peptide plane stacking is important for stabilization of buried E305 of Trichoderma reesei Cel5A.

Hydrogen bonds or salt bridges are usually formed to stabilize the buried ionizable residues. However, such interactions do not exist for two buried residues D271 and E305 of Trichoderma reesei Cel5A, an endoglucanase. Mutating D271 to alanine or leucine improves the enzyme thermostability quantified by the temperature T50 due to the elimination of the desolvation penalty of the aspartic acid. However, the same mutations for E305 decrease the enzyme thermostability. Free energy calculations based on the molecular dynamics simulation predict the thermostability of D271A, D271L, and E305A (compared to WT) in line with the experimental observation but overestimate the thermostability of E305L. Quantum mechanical calculations suggest that the carboxyl-peptide plane stacking interactions occurring to E305 but not D271 are important for the carboxyl group stabilization. For the protonated carboxyl group, the interaction energy can be as much as about -4 kcal/mol for parallel stacking and about -7 kcal/mol for T-shaped stacking. For the deprotonated carboxyl group, the largest interaction energies for parallel stacking and T-shaped stacking are comparable, about -7 kcal/mol. The solvation effect generally weakens the interaction, especially for the charged system. A search of the carboxyl-peptide plane stacking in the PDB databank indicates that parallel stacking but not T-shaped stacking is quite common, and the most probable distance between the two stacking fragments is close to the value predicted by the QM calculations. This work highlights the potential role of carboxyl amide π-π stacking in the stabilization of aspartic acid and glutamic acid in proteins.

The slowdown of the endoglucanase Trichoderma reesei Cel5A-catalyzed cellulose hydrolysis is related to its initial activity.

One important feature of hydrolysis of cellulose by cellulases is that the reaction slows down quickly after it starts. In this work, we investigate the slowdown mechanism at the early stage of the reaction using endoglucanase Tr. Cel5A-catalyzed phosphate acid-swollen cellulose (PASC) hydrolysis as a model system. Specifically, we focus on the effect of enzyme adsorption on the reaction slowdown. Nineteen single mutations are introduced (with the assistance of molecular dynamics simulations) to perturb the enzyme PASC interaction, yielding the adsorption partitioning coefficient Kr that ranged from 0.12 to 0.39 L/g, compared to that of the wild type (0.26 L/g). Several residues, including T18, K26, Y26, H229, and T300, are demonstrated to be important for adsorption of the enzyme to PASC. The kinetic measurements show that the slowdown of the hydrolysis is not correlated with the adsorption quantified by the partitioning coefficient Kr but is anticorrelated with the initial activity. This result suggests that the mutants with higher activity are more prone to being trapped or deplete the most reactive substrate faster and the adsorption plays no apparent role in the reaction slowdown. The initial activity of Cel5A against PASC is correlated with the enzyme specific activity against a soluble substrate p-nitrophenyl cellobioside.

Protein apparent dielectric constant and its temperature dependence from remote chemical shift effects.

A NMR protocol is introduced that permits accurate measurement of minute, remote chemical shift perturbations (CSPs), caused by a mutation-induced change in the electric field. Using protein GB3 as a model system, (1)H(N) CSPs in K19A and K19E mutants can be fitted to small changes in the electric field at distal sites in the protein using the Buckingham equation, yielding an apparent dielectric constant εa of 8.6 ± 0.8 at 298 K. These CSPs, and their derived εa value, scale strongly with temperature. For example, CSPs at 313 K are about ∼30% smaller than those at 278 K, corresponding to an effective εa value of about 7.3 at 278 K and 10.5 at 313 K. Molecular dynamics simulations in explicit solvent indicate that solvent water makes a significant contribution to εa.