Stabilization of a strained protein loop conformation through protein engineering

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field

Coarse-grained (CG) level molecular dynamics simulations are routinely used to study various biomolecular processes. The Martini force field is currently the most widely adopted parameter set for such simulations. The functional form of this and several other CG force fields enforces secondary protein structure support by employing a variety of harmonic potentials or restraints that favor the protein's native conformation. We propose a straightforward method to calculate the energetic consequences of transitions between predefined conformational states in systems in which multiple factors can affect protein conformational equilibria. This method is designed for use within the Martini force field and involves imposing conformational transitions by linking a Martini-inherent elastic network to the coupling parameter λ. We demonstrate the applicability of our method using the example of five biomolecular systems that undergo experimentally characterized conformational transitions between well-defined structures (Staphylococcal nuclease, C-terminal segment of surfactant protein B, LAH4 peptide, and β2-adrenergic receptor) as well as between folded and unfolded states (GCN4 leucine zipper protein). The results show that the relative free energy changes associated with protein conformational transitions, which are affected by various factors, such as pH, mutations, solvent, and lipid membrane composition, are correctly reproduced. The proposed method may be a valuable tool for understanding how different conditions and modifications affect conformational equilibria in proteins.

Exploring the functional residues in a flavin-binding fluorescent protein using deep mutational scanning

Flavin mononucleotide (FMN)-based fluorescent proteins are versatile reporters that can monitor various cellular processes in both aerobic and anaerobic conditions. However, the understanding of the role of individual amino acid residues on the protein function has been limited and has restricted the development of better functional variants. Here we examine the functional amino acid residues of Escherichia coli flavin mononucleotide binding fluorescent protein (EcFbFP) using the application of high-throughput sequencing of functional variants, termed deep mutational scanning. The variants were classified into 329 function-retained (FR) and 259 function-loss (FL) mutations, and further the mutational enrichment in each amino acid residues was weighed to find the functionally important residues of EcFbFP. We show that the crucial amino acid residues of EcFbFP lie among the FMN-binding pocket, turns and loops of the protein where conformation changes occur, and spatially clustered residues near the E56-K97 salt bridges. In addition, the mutational sensitivity of the critical residues was confirmed by site-directed mutagenesis. The deep mutational scanning of EcFbFP has demonstrated important implications for constructing better functioning protein variants.

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.

The pressure-temperature free energy-landscape of staphylococcal nuclease monitored by (1)H NMR

The thermodynamic stability of staphylococcal nuclease was studied against the variation of both temperature and pressure by utilizing (1)H NMR spectroscopy at 750 MHz in 20 mM Mes buffer containing 99.9 % (2)H(2)O, pH 5.3. Equilibrium fractions of folded and unfolded protein species were evaluated with the proton signals of two histidine residues as monitor in the pressure range of 30-3300 bar and in the temperature range of 1.5 degrees C-35 degrees C. From the multi-parameter fit of the experimental data to the Gibbs energy equation expressed as a simultaneous function of pressure and temperature, we determined the compressibility change (Deltabeta), the volume change at 1 bar (DeltaV degrees ) and the expansivity change (Deltaalpha) upon unfolding among other thermodynamic parameters: Deltabeta=0.02(+/-0.003) ml mol(-1) bar(-1); Deltaalpha=1.33(+/-0.2) ml mol(-1) K(-1); DeltaV degrees =-41.9(+/-6. 3) ml mol(-1) (at 24 degrees C); DeltaG degrees =13.18(+/-2) kJ mol(-1) (at 24 degrees C); DeltaC(p)=13.12(+/-2) kJ mol(-1) K(-1); DeltaS degrees =0.32(+/-0.05) kJ mol(-1) K(-1 )(at 24 degrees C). The result yields a three-dimensional free energy surface, i.e. the free energy-landscape of staphylococcal nuclease on the P-T plane. The significantly positive Deltabeta and Deltaalpha values suggest that, in the pressure-denatured state, staphylococcal nuclease forms a loosely packed and fluctuating structure. The slight but statistically significant difference between the unfolding transitions of the His8 and His124 environments is considered to reflect local fluctuations in the native state, leading to pre-melting of the His124 environment prior to the cooperative unfolding of the major part of the protein.

Protein retention in ion-exchange chromatography: effect of net charge and charge distribution

The charge regulated slab model is used to evaluate the salt dependence of the retention of Staphylococcal nuclease A and its mutants in cation-exchange chromatography. An important feature of this work is that the net charge of the proteins is varied in two different ways: (a) by changing the eluent pH so that the charges are created by protonation and (b) by point mutation at position 116. Since the structure of Staphylococcal nuclease and the mutants are known, the pH dependence of retention data of the different mutants gives detailed insights into the retention mechanism. Experimental results show that the salt dependence of retention is affected more strongly by changes of the eluent pH than by point mutations. This implies that the amino acid in position 116 has only a moderately strong interaction with the stationary phase surface and that a patch on one side of the protein surface is mainly responsible for the electrostatic interaction with the surface.

Solution structures of staphylococcal nuclease from multidimensional, multinuclear NMR: nuclease-H124L and its ternary complex with Ca2+ and thymidine-3',5'-bisphosphate

The solution structures of staphylococcal nuclease (nuclease) H124L and its ternary complex, (nuclease-H124L).pdTp.Ca2+, were determined by ab initio dynamic simulated annealing using 1925 NOE, 119 phi, 20 chi 1 and 112 hydrogen bond constraints for the free protein, and 2003 NOE, 118 phi, 20 chi 1 and 114 hydrogen bond constraints for the ternary complex. In both cases, the final structures display only small deviations from idealized covalent geometry. In structured regions, the overall root-mean-square deviations from mean atomic coordinates are 0.46 (+/- 0.05) A and 0.41 (+/- 0.05) A for the backbone heavy atoms of nuclease and its ternary complex, respectively. The backbone conformations of residues in the loop formed by Arg81-Gly86, which is adjacent to the active site, are more precisely defined in the ternary complex than in unligated nuclease. Also, the protein side chains that show NOEs and evidence for hydrogen bonds to pdTp (Arg35, Lys84, Tyr85, Arg87, Tyr113, and Tyr115) are better defined in the ternary complex. As has been observed previously in the X-ray structures of nuclease-WT, the binding of pdTp causes the backbone of Tyr113 to change from an extended to a left-handed alpha-helical conformation. The NMR structures reported here were compared with available X-ray structures: nuclease-H124L [Truckses et al. (1996) Protein Sci., 5, 1907-1916] and the ternary complex of wild-type staphylococcal nuclease [Loll and Lattman (1989) Proteins Struct. Funct. Genet., 5, 183-201]. Overall, the solution structures of nuclease-H124L are consistent with these crystal structures, but small differences were observed between the structures in the solution and crystal environments. These included differences in the conformations of certain side chains, a reduction in the extent of helix 1 in solution, and many fewer hydrogen bonds involving side chains in solution.

Coupling between trans/cis proline isomerization and protein stability in staphylococcal nuclease

The nucleases A produced by two strains of Staphylococcus aureus, which have different stabilities, differ only in the identity of the single amino acid at residue 124. The nuclease from the Foggi strain of S. aureus (by convention nuclease WT), which contains His124, is 1.9 kcal.mol-1 less stable (at pH 5.5 and 20 degrees C) than the nuclease from the V8 strain (by convention nuclease H124L), which contains Leu124. In addition, the population of the trans conformer at the Lys116-Pro117 peptide bond, as observed by NMR spectroscopy, is different for the two variants: about 15% for nuclease WT and 9% for nuclease H124L. In order to improve our understanding of the origin of these differences, we compared the properties of WT and H124L with those of the H124A and H124I variants. We discovered a correlation between effects of different residues at this position on protein stability and on stabilization of the cis configuration of the Lys116-Pro117 peptide bond. In terms of free energy, approximately 17% of the increase in protein stability manifests itself as stabilization of the cis configuration at Lys116-Pro117. This result implies that the differences in stability arise mainly from structural differences between the cis configurational isomers at Pro117 of the different variants at residue 124. We solved the X-ray structure of the cis form of the most stable variant, H124L, and compared it with the published high-resolution X-ray structure of the cis form of the most stable variant, WT (Hynes TR, Fox RO, 1991, Proteins Struct Funct Genet 10:92-105). The two structures are identical within experimental error, except for the side chain at residue 124, which is exposed in the models of both variants. Thus, the increased stability and changes in the trans/cis equilibrium of the Lys116-Pro117 peptide bond observed in H124L relative to WT are due to subtle structural changes that are not observed by current structure determination technique. Residue 124 is located in a helix. However, the stability changes are too large and follow the wrong order of stability to be explained simply by differences in helical propensity. A second site of conformational heterogeneity in native nuclease is found at the His46-Pro47 peptide bond, which is approximately 80% trans in both WT and H124L. Because proline to glycine substitutions at either residue 47 or 117 remove the structural heterogeneity at that position and increase protein stability, we determined the X-ray structures of H124L + P117G and H124L + P47G + P117G and the kinetic parameters of H124L, H124L + P47G, H124L + P117G, and H124L + P47G + P117G. The individual P117G and P47G mutations cause decreases in nuclease activity, with kcat affected more than Km, and their effects are additive. The P117G mutation in nuclease H124L leads to the same local conformational rearrangement described for the P117G mutant of WT (Hynes TR, Hodel A, Fox RO, 1994, Biochemistry 33:5021-5030). In both P117G mutants, the loop formed by residues 112-117 is located closer to the adjacent loop formed by residues 77-85, and residues 115-118 adopt a type I' beta-turn conformation with the Lys116-Gly117 peptide bond in the trans configuration, as compared with the parent protein in which these residues have a typeVIa beta-turn conformation with the Lys116-Pro117 peptide bond in the cis configuration. Addition of the P47G mutation appears not to cause any additional structural changes. However, the electron density for part of the loop containing this peptide bond was not strong enough to be interpreted.

Proline cis-trans isomerization in staphylococcal nuclease: multi-substrate free energy perturbation calculations

Staphylococcal nuclease A exists in two folded forms that differ in the isomerization state of the Lys 116-Pro 117 peptide bond. The dominant form (90% occupancy) adopts a cis peptide bond, which is observed in the crystal structure. NMR studies show that the relatively small difference in free energy between the cis and trans forms (delta Gcis-->trans approximately 1.2 kcal/mol) results from large and nearly compensating differences in enthalpy and entropy (delta Hcis-->trans approximately delta TScis-->trans approximately 10 kcal/mol). There is evidence from X-ray crystal structures that the structural differences between the cis and the trans forms of nuclease are confined to the conformation of residues 112-117, a solvated protein loop. Here, we obtain a thermodynamic and structural description of the conformational equilibrium of this protein loop through an exhaustive conformational search that identified several substates followed by free energy simulations between the substrates. By partitioning the search into conformational substates, we overcame the multiple minima problem in this particular case and obtained precise and reproducible free energy values. The protein and water environment was implicitly modeled by appropriately chosen nonbonded terms between the explicitly treated loop and the rest of the protein. These simulations correctly predicted a small free energy difference between the cis and trans forms composed of larger, compensating differences in enthalpy and entropy. The structural predictions of these simulations were qualitatively consistent with known X-ray structures of nuclease variants and yield a model of the unknown minor trans conformation.