Exploring the early stages of the pH‐induced conformational change of influenza hemagglutinin

Hemagglutinin (HA) mediates the membrane fusion process of influenza virus through its pH‐induced conformational change. However, it remains challenging to study its structure reorganization pathways in atomic details. Here, we first applied continuous constant pH molecular dynamics approach to predict the pK_a values of titratable residues in H2 subtype HA. The calculated net‐charges in HA1 globular heads increase from 0e (pH 7.5) to +14e (pH 4.5), indicating that the charge repulsion drives the detrimerization of HA globular domains. In HA2 stem regions, critical pH sensors, such as Glu103₂, His18₁, and Glu89₁, are identified to facilitate the essential structural reorganizations in the fusing pathways, including fusion peptide release and interhelical loop transition. To probe the contribution of identified pH sensors and unveil the early steps of pH‐induced conformational change, we carried out conventional molecular dynamics simulations in explicit water with determined protonation state for each titratable residue in different environmental pH conditions. Particularly, energy barriers involving previously uncharacterized hydrogen bonds and hydrophobic interactions are identified in the fusion peptide release pathway. Nevertheless, comprehensive comparisons across HA family members indicate that different HA subtypes might employ diverse pH sensor groups along with different fusion pathways. Finally, we explored the fusion inhibition mechanism of antibody CR6261 and small molecular inhibitor TBHQ, and discovered a novel druggable pocket in H2 and H5 subtypes. Our results provide the underlying mechanism for the pH‐driven conformational changes and also novel insight for anti‐flu drug development. Proteins 2014; 82:2412–2428. © 2014 Wiley Periodicals, Inc.     

pH‐dependent molecular dynamics of vesicular stomatitis virus glycoprotein G

Vesicular stomatitis virus glycoprotein G (VSV‐G) belongs to a new class of viral fusion proteins (Class III). The structure of VSV‐G has been solved in two different conformations and fusion is known to be triggered by low pH. To investigate Class III fusion mechanisms, molecular dynamics simulations were performed on the VSV‐G prefusion structure in two different protonation states: at physiological pH (pH 7) and low pH present in the endosome (pH 5). Domain IV containing the fusion loops, which need to interact with the target membrane, exhibits the highest mobility. Energetic analyses revealed weakened interaction between Domain IV and the protein core at pH 5, which can be attributed to two pairs of structurally neighboring conserved and differentially protonated residues in the Domain IV–core interface. Energetic calculations also demonstrated that the interaction between the subunits in the core of the trimeric VSV‐G is strengthened at pH 5, mainly due to newly formed interactions between the C‐terminal loop of Domain II and the N‐terminus of the adjacent subunit. A pair of interacting residues in this interface that is affected by differential protonation was shown to be the main effectors of this phenomenon. The results of this study thus enhance the mechanistic understanding of the effects of protonation changes in VSV‐G. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

Bending of the BST‐2 coiled‐coil during viral budding

BST‐2/tetherin is a human extracellular transmembrane protein that serves as a host defense factor against HIV‐1 and other viruses by inhibiting viral spreading. Structurally, BST‐2 is a homo‐dimeric coiled‐coil that is connected to the host cell membrane by N and C terminal transmembrane anchors. The C‐terminal membrane anchor of BST‐2 is inserted into the budding virus while the N‐terminal membrane anchor remains in the host cell membrane creating a viral tether. The structural mechanism of viral budding and tethering as mediated by BST‐2 is not clear. To more fully describe the mechanism of viral tethering, we created a model of BST‐2 embedded in a membrane and used steered molecular dynamics to simulate the transition from the host cell membrane associated form to the cell‐virus membrane bridging form. We observed that BST‐2 did not transition as a rigid structure, but instead bent at positions with a reduced interface between the helices of the coiled‐coil. The simulations for the human BST‐2 were then compared with simulations on the mouse homolog, which has no apparent weak spots. We observed that the mouse homolog spread the bending across the ectodomain, rather than breaking at discrete points as observed with the human homolog. These simulations support previous biochemical and cellular work suggesting some flexibility in the coiled‐coil is necessary for viral tethering, while also highlighting how subtle changes in protein sequence can influence the dynamics and stability of proteins with overall similar structure.     

Energetics of the loop-to-helix transition leading to the coiled-coil structure of influenza virus hemagglutinin HA2 subunits

Fusion of influenza virus with the endosomal membrane of the host cell is mediated by the homotrimer-organized glycoprotein hemagglutinin (HA). Its fusion activity is triggered by a low pH-mediated conformational change affecting the structure of the HA1 and HA2 subunits. The HA2 subunits undergo a loop-to-helix transition leading to a coiled-coil structure, a highly conserved motif for many fusion mediating viral proteins. However, experimental studies showed that the HA2 coiled-coil structure is stable at neutral and low pH, implying that there is no direct relationship between low pH and the HA2 loop-to-helix transition. To interpret this observation, we used a computational approach based on the dielectric continuum solvent model to explore the influence of water and pH on the free energy change of the transition. The computations showed that the electrostatic interaction between HA2 fragments and water is the major driving force of the HA2 loop-to-helix transition leading to the coiled-coil structure, as long as the HA1 globular domain covering the HA2 subunits in the nonfusion competent conformation is reorganized and thereby allows water molecules to interact with the whole loop segments of the HA2 subunits. Moreover, we show that the energy released by the loop-to-helix transition may account for those energies required for driving the subsequent steps of membrane fusion. Such a water-driven process may resemble a general mechanism for the formation of the highly conserved coiled-coil motif of enveloped viruses.     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.     

A pH‐ and ionic strength‐dependent conformational change in the neck region regulates DNGR‐1 function in dendritic cells

DNGR‐1 is receptor expressed by certain dendritic cell (DC) subsets and by DC precursors in mouse. It possesses a C‐type lectin‐like domain (CTLD) followed by a poorly characterized neck region coupled to a transmembrane region and short intracellular tail. The CTLD of DNGR‐1 binds F‐actin exposed by dead cell corpses and causes the receptor to signal and potentiate cross‐presentation of dead cell‐associated antigens by DCs. Here, we describe a conformational change that occurs in the neck region of DNGR‐1 in a pH‐ and ionic strength‐dependent manner and that controls cross‐presentation of dead cell‐associated antigens. We identify residues in the neck region that, when mutated, lock DNGR‐1 in one of the two conformational states to potentiate cross‐presentation. In contrast, we show that chimeric proteins in which the neck region of DNGR‐1 is replaced by that of unrelated C‐type lectin receptors fail to promote cross‐presentation. Our results suggest that the neck region of DNGR‐1 is an integral receptor component that senses receptor progression through the endocytic pathway and has evolved to maximize extraction of antigens from cell corpses, coupling DNGR‐1 function to its cellular localization.     

pH-induced conformational changes of membrane-bound influenza hemagglutinin and its effect on target lipid bilayers.

Influenza virus hemagglutinin (HA) has served as a paradigm for both pH-dependent and -independent viral membrane fusion. Although large conformational changes were observed by X-ray crystallography when soluble fragments of HA were subjected to fusion-pH conditions, it is not clear whether the same changes occur in membrane-bound HA, what the spatial relationship is between the conformationally changed HA and the target and viral membranes, and in what way HA perturbs the target membrane at low pH. We have taken a spectroscopic approach using an array of recently developed FTIR techniques to address these questions. Difference attenuated total reflection FTIR spectroscopy was employed to reveal reversible and irreversible components of the pH-induced conformational change of the membrane-bound bromelain fragment of HA, BHA. Additional proteolytic fragments of BHA were produced which permitted a tentative assignment of the observed changes to the HA1 and HA2 subunits, respectively. The membrane-bound HA1 subunit undergoes a reversible conformational change, which most likely involves the loss of a small proportion of beta-sheet at low pH. BHA was found to undergo a partially reversible tilting motion relative to the target membrane upon exposure to pH 5, indicating a previously undescribed hinge near the anchoring point to the target membrane. Time-resolved amide H/D exchange experiments revealed a more dynamic (tertiary) structure of membrane-bound BHA and its HA2, but not its HA1, subunit. Finally BHA and, to a lesser degree, HA1 perturbed the lipid bilayer of the target membrane at the interface, as assessed by spectral changes of the lipid ester carbonyl groups. These results are discussed in the context of a complementary study of HA that was bound to viral membranes through its transmembrane peptide (Gray C, Tamm LK, 1997, Protein Sci 6:1993-2006). A distinctive role for the HA1 subunit in the conformational change of HA becomes apparent from these combined studies.     

Crystal structure of a trimeric form of the KV7.1 (KCNQ1) A‐domain tail coiled‐coil reveals structural plasticity and context dependent changes in a putative coiled‐coil trimerization motif

Coiled-coils are widespread protein–protein interaction motifs typified by the heptad repeat (abcdefg)_n in which “a” and “d” positions are hydrophobic residues. Although identification of likely coiled-coil sequences is robust, prediction of strand order remains elusive. We present the X-ray crystal structure of a short form (residues 583–611), “Q1-short,” of the coiled-coil assembly specificity domain from the voltage-gated potassium channel Kv7.1 (KCNQ1) determined at 1.7 Å resolution. Q1-short lacks one and half heptads present in a previously studied tetrameric coiled-coil construct, Kv7.1 585–621, “Q1-long.” Surprisingly, Q1-short crystallizes as a trimer. In solution, Q1-short self-assembles more poorly than Q1-long and depends on an R-h-x-x-h-E motif common to trimeric coiled-coils. Addition of native sequences that include “a” and “d” positions C-terminal to Q1-short overrides the R-h-x-x-h-E motif influence and changes assembly state from a weakly associated trimer to a strongly associated tetramer. These data provide a striking example of a naturally occurring amino sequence that exhibits context-dependent folding into different oligomerization states, a three-stranded versus a four-stranded coiled-coil. The results emphasize the degenerate nature of coiled-coil energy landscapes in which small changes can have drastic effects on oligomerization. Discovery of these properties in an ion channel assembly domain and prevalence of the R-h-x-x-h-E motif in coiled-coil assembly domains of a number of different channels that are thought to function as tetrameric assemblies raises the possibility that such sequence features may be important for facilitating the assembly of intermediates en route to the final native state.     

Crosstalk between PKA and PKG controls pH‐dependent host cell egress of Toxoplasma gondii

Toxoplasma gondii encodes three protein kinase A catalytic (PKAc1‐3) and one regulatory (PKAr) subunits to integrate cAMP‐dependent signals. Here, we show that inactive PKAc1 is maintained at the parasite pellicle by interacting with acylated PKAr. Either a conditional knockdown of PKAr or the overexpression of PKAc1 blocks parasite division. Conversely, down‐regulation of PKAc1 or stabilisation of a dominant‐negative PKAr isoform that does not bind cAMP triggers premature parasite egress from infected cells followed by serial invasion attempts leading to host cell lysis. This untimely egress depends on host cell acidification. A phosphoproteome analysis suggested the interplay between cAMP and cGMP signalling as PKAc1 inactivation changes the phosphorylation profile of a putative cGMP‐phosphodiesterase. Concordantly, inhibition of the cGMP‐dependent protein kinase G (PKG) blocks egress induced by PKAc1 inactivation or environmental acidification, while a cGMP‐phosphodiesterase inhibitor circumvents egress repression by PKAc1 or pH neutralisation. This indicates that pH and PKAc1 act as balancing regulators of cGMP metabolism to control egress. These results reveal a crosstalk between PKA and PKG pathways to govern egress in T. gondii.     

On the pH‐optimum of activity and stability of proteins

Biological macromolecules evolved to perform their function in specific cellular environment (subcellular compartments or tissues); therefore, they should be adapted to the biophysical characteristics of the corresponding environment, one of them being the characteristic pH. Many macromolecular properties are pH dependent, such as activity and stability. However, only activity is biologically important, while stability may not be crucial for the corresponding reaction. Here, we show that the pH‐optimum of activity (the pH of maximal activity) is correlated with the pH‐optimum of stability (the pH of maximal stability) on a set of 310 proteins with available experimental data. We speculate that such a correlation is needed to allow the corresponding macromolecules to tolerate small pH fluctuations that are inevitable with cellular function. Our findings rationalize the efforts of correlating the pH of maximal stability and the characteristic pH of subcellular compartments, as only pH of activity is subject of evolutionary pressure. In addition, our analysis confirmed the previous observation that pH‐optimum of activity and stability are not correlated with the isoelectric point, pI, or with the optimal temperature. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Cell‐free production of trimeric influenza hemagglutinin head domain proteins as vaccine antigens

In order to effectively combat pandemic influenza threats, there is a need for more rapid and robust vaccine production methods. In this article, we demonstrate E. coli‐based cell‐free protein synthesis (CFPS) as a method to rapidly produce domains from the protein hemagglutinin (HA), which is present on the surface of the influenza virus. The portion of the HA coding sequence for the “head” domain from the 2009 pandemic H1N1 strain was first optimized for E. coli expression. The protein domain was then produced in CFPS reactions and purified in soluble form first as a monomer and then as a trimer by a C‐terminal addition of the T4 bacteriophage foldon domain. Production of soluble trimeric HA head domain was enhanced by introducing stabilizing amino acid mutations to the construct in order to avoid aggregation. Trimerization was verified using size exclusion HPLC, and the stabilized HA head domain trimer was more effectively recognized by antibodies from pandemic H1N1 influenza vaccine recipients than was the monomer and also bound to sialic acids more strongly, indicating that the trimers are correctly formed and could be potentially effective as vaccines. Biotechnol. Bioeng. 2012; 109: 2962–2969. © 2012 Wiley Periodicals, Inc.     

Structure of HIV‐1 reverse transcriptase/d4TTP complex: Novel DNA cross‐linking site and pH‐dependent conformational changes

Stavudine (d4T, 2′,3′‐didehydro‐2′,3′‐dideoxythymidine) was one of the first chain‐terminating nucleoside analogs used to treat HIV infection. We present the first structure of the active, triphosphate form of d4T (d4TTP) bound to a catalytic complex of HIV‐1 RT/dsDNA template‐primer. We also present a new strategy for disulfide (S–S) chemical cross‐linking between N⁶ of a modified adenine at the second overhang base to I63C in the fingers subdomain of RT. The cross‐link site is upstream of the duplex‐binding region of RT, however, the structure is very similar to published RT structures with cross‐linking to Q258C in the thumb, which suggests that cross‐linking at either site does not appreciably perturb the RT/DNA structures. RT has a catalytic maximum at pH 7.5. We determined the X‐ray structures of the I63C‐RT/dsDNA/d4TTP cross‐linked complexes at pH 7, 7.5, 8, 8.5, 9, and 9.5. We found small (~0.5 Å), pH‐dependent motions of the fingers subdomain that folds in to form the dNTP‐binding pocket. We propose that the pH‐activity profile of RT relates to this motion of the fingers. Due to side effects of neuropathy and lipodystrophy, use of d4T has been stopped in most countries, however, chemical modification of d4T might lead to the development of a new class of nucleoside analogs targeting RNA and DNA polymerases.     

Conformational analysis of the full‐length M2 protein of the influenza A virus using solid‐state NMR

The influenza A M2 protein forms a proton channel for virus infection and mediates virus assembly and budding. While extensive structural information is known about the transmembrane helix and an adjacent amphipathic helix, the conformation of the N‐terminal ectodomain and the C‐terminal cytoplasmic tail remains largely unknown. Using two‐dimensional (2D) magic‐angle‐spinning solid‐state NMR, we have investigated the secondary structure and dynamics of full‐length M2 (M2FL) and found them to depend on the membrane composition. In 2D ¹³C DARR correlation spectra, 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC)‐bound M2FL exhibits several peaks at β‐sheet chemical shifts, which result from water‐exposed extramembrane residues. In contrast, M2FL bound to cholesterol‐containing membranes gives predominantly α‐helical chemical shifts. Two‐dimensional J‐INADEQUATE spectra and variable‐temperature ¹³C spectra indicate that DMPC‐bound M2FL is highly dynamic while the cholesterol‐containing membranes significantly immobilize the protein at physiological temperature. Chemical‐shift prediction for various secondary‐structure models suggests that the β‐strand is located at the N‐terminus of the DMPC‐bound protein, while the cytoplasmic domain is unstructured. This prediction is confirmed by the 2D DARR spectrum of the ectodomain‐truncated M2(21–97), which no longer exhibits β‐sheet chemical shifts in the DMPC‐bound state. We propose that the M2 conformational change results from the influence of cholesterol, and the increased helicity of M2FL in cholesterol‐rich membranes may be relevant for M2 interaction with the matrix protein M1 during virus assembly and budding. The successful determination of the β‐strand location suggests that chemical‐shift prediction is a promising approach for obtaining structural information of disordered proteins before resonance assignment.     

Insertion of single‐chain variable fragment (scFv) peptide linker improves surface display of influenza hemagglutinin (HA1) on non‐recombinant Lactococcus lactis

Heterologous protein displayed on the surface of Lactococcus lactis using the binding domain of N‐acetylmuramidase (AcmA) has a potential application in vaccine delivery. In this study, we developed a non‐recombinant L. lactis surface displaying the influenza A (H1N1) 2009 hemagglutinin (HA1). Three recombinant proteins, HA1/L/AcmA, HA1/AcmA, and HA1 were overexpressed in Escherichia coli, and purified. In the binding study using flow cytometry, the HA1/L/AcmA, which contained the single‐chain variable fragment (scFv) peptide linker showed significantly higher percentage of binding counts and mean fluorescence binding intensity (MFI) (51.7 ± 1.4% and 3,594.0 ± 675.9, respectively) in comparison to the HA1/AcmA without the scFv peptide linker (41.1 ± 1.5% and 1,652.0 ± 34.1, respectively). Higher amount of HA1/L/AcmA (～2.9 × 10⁴ molecules per cell) was displayed on L. lactis when compared to HA1/AcmA (～1.1 × 10⁴ molecules per cell) in the immunoblotting analysis. The HA1/L/AcmA completely agglutinated RBCs at comparable amount of protein to that of HA1/AcmA and HA1. Computational modeling of protein structures suggested that scFv peptide linker in HA1/L/AcmA kept the HA1 and the AcmA domain separated at a much longer distance in comparison to HA1/AcmA. These findings suggest that insertion of the scFv peptide linker between HA1 and AcmA improved binding of recombinant proteins to L. lactis. Hence, insertion of scFv peptide linker can be further investigated as a potential approach for improvement of heterologous proteins displayed on the surface of L. lactis using the AcmA binding domain. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 33:154–162, 2017     

Metal‐dependent assembly of a protein nano‐cage

Short, alpha‐helical coiled coils provide a simple, modular method to direct the assembly of proteins into higher order structures. We previously demonstrated that by genetically fusing de novo–designed coiled coils of the appropriate oligomerization state to a natural trimeric protein, we could direct the assembly of this protein into various geometrical cages. Here, we have extended this approach by appending a coiled coil designed to trimerize in response to binding divalent transition metal ions and thereby achieve metal ion‐dependent assembly of a tetrahedral protein cage. Ni²⁺, Co²⁺, Cu²⁺, and Zn²⁺ ions were evaluated, with Ni²⁺ proving the most effective at mediating protein assembly. Characterization of the assembled protein indicated that the metal ion–protein complex formed discrete globular structures of the diameter expected for a complex containing 12 copies of the protein monomer. Protein assembly could be reversed by removing metal ions with ethylenediaminetetraacetic acid or under mildly acidic conditions.     

Differential backbone dynamics of companion helices in the extended helical coiled‐coil domain of a bacterial chemoreceptor

Cytoplasmic domains of transmembrane bacterial chemoreceptors are largely extended four‐helix coiled coils. Previous observations suggested the domain was structurally dynamic. We probed directly backbone dynamics of this domain of the transmembrane chemoreceptor Tar from Escherichia coli using site‐directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. Spin labels were positioned on solvent‐exposed helical faces because EPR spectra for such positions reflect primarily polypeptide backbone movements. We acquired spectra for spin‐labeled, intact receptor homodimers solubilized in detergent or inserted into native E. coli lipid bilayers in Nanodiscs, characterizing 16 positions distributed throughout the cytoplasmic domain and on both helices of its helical hairpins, one amino terminal to the membrane‐distal tight turn (N‐helix), and the other carboxyl terminal (C‐helix). Detergent solubilization increased backbone dynamics for much of the domain, suggesting that loss of receptor activities upon solubilization reflects wide‐spread destabilization. For receptors in either condition, we observed an unanticipated difference between the N‐ and C‐helices. For bilayer‐inserted receptors, EPR spectra from sites in the membrane‐distal protein‐interaction region and throughout the C‐helix were typical of well‐structured helices. In contrast, for approximately two‐thirds of the N‐helix, from its origin as the AS‐2 helix of the membrane‐proximal HAMP domain to the beginning of the membrane‐distal protein‐interaction region, spectra had a significantly mobile component, estimated by spectral deconvolution to average approximately 15%. Differential helical dynamics suggests a four‐helix bundle organization with a pair of core scaffold helices and two more dynamic partner helices. This newly observed feature of chemoreceptor structure could be involved in receptor function.     

Optimization of norbornyl‐based carbocyclic nucleoside analogs as cyclin‐dependent kinase 2 inhibitors

We report on the discovery of norbornyl moiety as a novel structural motif for cyclin‐dependent kinase 2 (CDK2) inhibitors which was identified by screening a carbocyclic nucleoside analogue library. Three micromolar hits were expanded by the use of medicinal chemistry methods into a series of 16 novel compounds. They had prevailingly micromolar activities against CDK2 and the best compound of the series attained IC₅₀ of 190 nM. The binding modes were explored in molecular details by modeling and docking. Quantum mechanics‐based scoring was used to rationalize the affinities. In conclusion, the discovered 9‐hydroxymethylnorbornyl moiety was shown by joint experimental‐theoretical efforts to be able to serve as a novel substituent for CDK2 inhibitors. This finding opens door to the exploration of chemical space towards more effective derivatives targeting this important class of protein kinases.     

Comparison of the structures and stabilities of coiled‐coil proteins containing hexafluoroleucine and t‐butylalanine provides insight into the stabilizing effects of highly fluorinated amino acid side‐chains

Highly fluorinated analogs of hydrophobic amino acids are well known to increase the stability of proteins toward thermal unfolding and chemical denaturation, but there is very little data on the structural consequences of fluorination. We have determined the structures and folding energies of three variants of a de novo designed 4‐helix bundle protein whose hydrophobic cores contain either hexafluoroleucine (hFLeu) or t‐butylalanine (tBAla). Although the buried hydrophobic surface area is the same for all three proteins, the incorporation of tBAla causes a rearrangement of the core packing, resulting in the formation of a destabilizing hydrophobic cavity at the center of the protein. In contrast, incorporation of hFLeu, causes no changes in core packing with respect to the structure of the nonfluorinated parent protein which contains only leucine in the core. These results support the idea that fluorinated residues are especially effective at stabilizing proteins because they closely mimic the shape of the natural residues they replace while increasing buried hydrophobic surface area.