Artificial protein cavities as specific ligand-binding templates: characterization of an engineered heterocyclic cation-binding site that preserves the evolved specificity of the parent protein.

Cavity complementation has been observed in many proteins, where an appropriate small molecule binds to a cavity-forming mutant. Here, the binding of compounds to the W191G cavity mutant of cytochrome c peroxidase is characterized by X-ray crystallography and binding thermodynamics. Unlike cavities created by removal of hydrophobic side-chains, the W191G cavity does not bind neutral or hydrophobic compounds, but displays a strong specificity for heterocyclic cations, consistent with the role of the protein to stabilize a tryptophan radical at this site. Ligand dissociation constants for the protonated cationic state ranged from 6 microM for 2-amino-5-methylthiazole to 1 mM for neutral ligands, and binding was associated with a large enthalpy-entropy compensation. X-ray structures show that each of 18 compounds with binding behavior bind specifically within the artificial cavity and not elsewhere in the protein. The compounds make multiple hydrogen bonds to the cavity walls using a subset of the interactions seen between the protein and solvent in the absence of ligand. For all ligands, every atom that is capable of making a hydrogen bond does so with either protein or solvent. The most often seen interaction is to Asp235, and most compounds bind with a specific orientation that is defined by their ability to interact with this residue. Four of the ligands do not have conventional hydrogen bonding atoms, but were nevertheless observed to orient their most polar CH bond towards Asp235. Two of the larger ligands induce disorder in a surface loop between Pro190 and Asn195 that has been identified as a mobile gate to cavity access. Despite the predominance of hydrogen bonding and electrostatic interactions, the small variation in observed binding free energies were not correlated readily with the strength, type or number of hydrogen bonds or with calculated electrostatic energies alone. Thus, as with naturally occurring binding sites, affinities to W191G are likely to be due to a subtle balance of polar, non-polar, and solvation terms. These studies demonstrate how cavity complementation and judicious choice of site can be used to produce a protein template with an unusual ligand-binding specificity.     

Binding site multiplicity with fatty acid ligands: Implications for the regulation of PKR kinase autophosphorylation with palmitate

Saturated long chain‐free fatty acids (FFAs), especially palmitate, have been implicated in apoptosis by inhibiting the activity of PKR (double‐stranded RNA‐dependent protein kinase). We recently found evidence that palmitate interacts directly with the kinase domain of PKR, subsequently inhibiting the autophosphorylation of PKR. To investigate the interactions of palmitate with PKR and its effects on PKR autophosphorylation, we performed extensive unbiased MD simulations combined with biochemical and biophysical experiments. The simulations predict multiple putative binding sites of palmitate on both the phosphorylated and unphosphorylated PKR with similar binding affinities. Ligand‐protein interactions involving a large variety of different binding modes challenge the conventional view of highly specific, single binding sites. Key interactions of palmitate involve the αC‐helix of PKR, especially near residue R307. Experimental mutation of R307 was found to affect palmitate binding and reduce its inhibitory effect. Based on this study a new allosteric mechanism is proposed where palmitate binding to the αC‐helix prevents the inactive‐to‐active transition of PKR and subsequently reduces its ability to autophosphorylate. Proteins 2014; 82:2429–2442. © 2014 Wiley Periodicals, Inc.     

Structure of a novel PTH‐related peptide hPTH′ and its interaction with the PTH receptor

We have previously shown that a recombinant human PTH fragment, Pro‐Pro‐[Arg11] hPTH (1–34)‐Pro‐Pro‐Asp (hPTH′), could be a potentially better and more cost‐effective therapeutic agent than PTH (1–34) on osteoporosis. In this report, we characterized the solution conformations of hPTH′ by NMR spectroscopy and modeled the interactions between the hPTH′ and the PTH receptor. By comparing it with PTH (1–34) structures and their respective interactions with the PTH receptor, we identified two segments of helix extending from Ile5 to Met8 and from Glu22 to Gln29 with a divided kink between the two helixes around Arg20. Mutated arginine makes hPTH′ fill the receptor cavity better as well as forms hydrogen bonds with Val193. Understanding the ligand receptor interactions will help us design small molecules to better mimic the activities of PTH. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

Solution structure of barley lipid transfer protein complexed with palmitate. Two different binding modes of palmitate in the homologous maize and barley nonspecific lipid transfer proteins.

The structure of a nonspecific lipid transfer protein from barley (ns-LTPbarley) in complex with palmitate has been determined by NMR spectroscopy. The structure has been compared to the structure of ns-LTPbarley in the absence of palmitate, to the structure of ns-LTPbarley in complex with palmitoyl coenzyme A, to the structure of ns-LTPmaize in its free form, and to the maize protein complexed with palmitate. Binding of palmitate only affects the structure of ns-LTPbarley moderately in contrast to the binding of palmitoyl coenzyme A, which leads to a considerable expansion of the protein. The modes of binding palmitate to the maize and barley protein are different. Although in neither case there are major conformational changes in the protein, the orientation of the palmitate in the two proteins is exactly opposite.     

Post-translational modification of barley LTP1b: the lipid adduct lies in the hydrophobic cavity and alters the protein dynamics

NMR techniques have been used to characterise the effects of a lipid-like post-translational modification on barley lipid transfer protein (LTP1b). NMR chemical shift data indicate that the lipid-like molecule lies in the hydrophobic cavity of LTP1b, with Tyr 79 being displaced to accommodate the ligand in the cavity. The modified protein has a reduced level of backbone amide hydrogen exchange protection, presumably reflecting increased dynamics in the protein. This may result from a loosening of the protein structure and may explain the enhanced surface properties observed for LTP1b.     

Contribution of intra- and intermolecular hydrogen bonds to the conformational stability of human lysozyme(,).

In globular proteins, there are intermolecular hydrogen bonds between protein and water molecules, and between water molecules, which are bound with the proteins, in addition to intramolecular hydrogen bonds. To estimate the contribution of these hydrogen bonds to the conformational stability of a protein, the thermodynamic parameters for denaturation and the crystal structures of five Thr to Val and five Thr to Ala mutant human lysozymes were determined. The denaturation Gibbs energy (DeltaG) of Thr to Val and Thr to Ala mutant proteins was changed from 4.0 to -5.6 kJ/mol and from 1.6 to -6.3 kJ/mol, respectively, compared with that of the wild-type protein. The contribution of hydrogen bonds to the stability (DeltaDeltaG(HB)) of the Thr and other mutant human lysozymes previously reported was extracted from the observed stability changes (DeltaDeltaG) with correction for changes in hydrophobicity and side chain conformational entropy between the wild-type and mutant structures. The estimation of the DeltaDeltaG(HB) values of all mutant proteins after removal of hydrogen bonds, including protein-water hydrogen bonds, indicates a favorable contribution of the intra- and intermolecular hydrogen bonds to the protein stability. The net contribution of an intramolecular hydrogen bond (DeltaG(HB[pp])), an intermolecular one between protein and ordered water molecules (DeltaG(HB[pw])), and an intermolecular one between ordered water molecules (DeltaG(HB[ww])) could be estimated to be 8. 5, 5.2, and 5.0 kJ/mol, respectively, for a 3 A long hydrogen bond. This result shows the different contributions to protein stability of intra- and intermolecular hydrogen bonds. The entropic cost due to the introduction of a water molecule (DeltaG(H)()2(O)) could be also estimated to be about 8 kJ/mol.     

Molecular dynamics simulations of the auxin‐binding protein 1 in complex with indole‐3‐acetic acid and naphthalen‐1‐acetic acid

Auxin‐binding protein 1 (ABP1) is suggested to be an auxin receptor which plays an important role in several processes in green plants. Maize ABP1 was simulated with the natural auxin indole‐3‐acetic acid (IAA) and the synthetic analog naphthalen‐1‐acetic acid (NAA), to elucidate the role of the KDEL sequence and the helix at the C‐terminus. The KDEL sequence weakens the intermolecular interactions between the monomers but stabilizes the C‐terminal helix. Conformational changes at the C‐terminus occur within the KDEL sequence and are influenced by the binding of the simulated ligands. This observation helps to explain experimental findings on ABP1 interactions with antibodies that are modulated by the presence of auxin, and supports the hypothesis that ABP1 acts as an auxin receptor. Stable hydrogen bonds between the monomers are formed between Glu40 and Glu62, Arg10 and Thr97, Lys39, and Glu62 in all simulations. The amino acids Ile22, Leu25, Trp44, Pro55, Ile130, and Phe149 are located in the binding pocket and are involved in hydrophobic interactions with the ring system of the ligand. Trp151 is stably involved in a face to end interaction with the ligand. The calculated free energy of binding using the linear interaction energy approach showed a higher binding affinity for NAA as compared to IAA. Our simulations confirm the asymmetric behavior of the two monomers, the stronger interaction of NAA than IAA and offers insight into the possible mechanism of ABP1 as an auxin receptor. Proteins 2014; 82:2744–2755. © 2014 Wiley Periodicals, Inc.     

Solution structure of a lipid transfer protein extracted from rice seeds. Comparison with homologous proteins.

Nuclear magnetic resonance (NMR) spectroscopy was used to determine the three dimensional structure of rice nonspecific lipid transfer protein (ns-LTP), a 91 amino acid residue protein belonging to the broad family of plant ns-LTP. Sequence specific assignment was obtained for all but three HN backbone 1H resonances and for more than 95% of the 1H side-chain resonances using a combination of 1H 2D NOESY; TOCSY and COSY experiments at 293 K. The structure was calculated on the basis of four disulfide bridge restraints, 1259 distance constraints derived from 1H-1H Overhauser effects, 72 phi angle restraints and 32 hydrogen-bond restraints. The final solution structure involves four helices (H1: Cys3-Arg18, H2: Ala25-Ala37, H3: Thr41-Ala54 and H4: Ala66-Cys73) followed by a long C-terminal tail (T) with no observable regular structure. N-capping residues (Thr2, Ser24, Thr40), whose side-chain oxygen atoms are involved in hydrogen bonds with i + 3 amide proton additionally stabilize the N termini of the first three helices. The fourth helix involving Pro residues display a mixture of alpha and 3(10) conformation. The rms deviation of 14 final structures with respect to the average structure is 1.14 +/- 0.16 A for all heavy atoms (C, N, O and S) and 0.72 +/- 0.01 A for the backbone atoms. The global fold of rice ns-LTP is close to the previously published structures of wheat, barley and maize ns-LTPs exhibiting nearly identical pattern of the numerous sequence specific interactions. As reported previously for different four-helix topology proteins, hydrophobic, hydrogen bonding and electrostatic mechanisms of fold stabilization were found for the rice ns-LTP. The sequential alignment of 36 ns-LTP primary structures strongly suggests that there is a uniform pattern of specific long-range interactions (in terms of sequence), which stabilize the fold of all plant ns-LTPs.     

A novel lipid transfer protein from the dill Anethum graveolens L.: isolation,structure, heterologous expression,and functional characteristics

A novel lipid transfer protein, designated as Ag‐LTP, was isolated from aerial parts of the dill Anethum graveolens L. Structural, antimicrobial, and lipid binding properties of the protein were studied. Complete amino acid sequence of Ag‐LTP was determined. The protein has molecular mass of 9524.4 Da, consists of 93 amino acid residues including eight cysteines forming four disulfide bonds. The recombinant Ag‐LTP was overexpressed in Escherichia coli and purified. NMR investigation shows that the Ag‐LTP spatial structure contains four α ‐helices, forming the internal hydrophobic cavity, and a long C‐terminal tail. The measured volume of the Ag‐LTP hydrophobic cavity is equal to ~800 A³, which is much larger than those of other plant LTP1s. Ag‐LTP has weak antifungal activity and unpronounced lipid binding specificity but effectively binds plant hormone jasmonic acid. Our results afford further molecular insight into biological functions of LTP in plants. Copyright © 2015 European Peptide Society and John Wiley & Sons, Ltd.     

Effect of Y220C Mutation on p53 and Its Rescue Mechanism: A Computer Chemistry Approach

Mutation causes inactivation of ‘p53’ tumor suppressor protein in almost fifty percent of cancers in humans. Outside the DNA-binding surface of p53, Y220C is the most common cancerous mutation. Previous studies have shown that a surface cavity is created by this mutation which destabilizes p53. PhiKan083, a carbazole derivative capable of binding with that cavity, and slows down its thermal denaturation rate. We investigated, theoretically, on mechanisms of structural stability loss due to Y220C mutation and mechanisms of stability restoration by PhiKan083 at the atomic level. From this study it is found that in Tp53C, Tyr220 has five electrostatic interactions with residues Val 147, Prol51, Pro153 and Pro223 located on S3/S4 loop and S7/S8 loop. The S7/S8 loop is stabilized by these electrostatic interactions. Due to the Y220C mutation all these electrostatic interactions are lost. As a result the structural fluctuation occurs at S7/S8 loop, and the loop is displaced from its original position after 6 ns MD simulation. When PhiKan083 is present (inserted) at the mutation site it provides five electrostatic interactions with Pro155, Glu221 and Thr230, and two hydrogen bonds with Leu145 and Asp228, respectively. These interactions provided by Pkikan083 stabilized the S7/S8 loop, and as a result it couldn’t be displaced. Our results showed that due to Y220C mutation p53 became destabilized through structural fluctuations surrounding the mutation site. When PhiKan083 is present at the Y220C mutation site (in 2vuk), it provides electrostatic and hydrogen bonding interactions among residue-220, its neighboring residues and PhiKan08. These interactions give additional stability to Y220C mutant p53, thus Y220C mutant p53 doesn’t destabilize.     

Molecular dynamics simulations of a winter flounder “antifreeze” polypeptide in aqueous solution

A winter flounder antifreeze polypeptide (HPLC-6) has been studied in vacuo and in aqueous solution using molecular dynamics computer simulation techniques. The helical conformation of this polypeptide was found to be stable both in vacuum and in solution. The major stabilizing interactions were found to be the main-chain hydrogen bonds, a salt-bridge interaction, and solute–solvent hydrogen bonds. A significant bending in the middle of the polypeptide chain was observed both in vacuo and in solvent at 300 K. Possible causes of the bending are discussed. From simulations of mutant polypeptide molecules in vacuo, it is concluded that the bend in the native polypeptide was caused by side chain to backbone hydrogen bond competition involving the Thr 24 side chain and facilitated by strains on the helix resulting from the Lys 18-Glu 22 salt bridge. © 1993 John Wiley & Sons, Inc.     

Molecular dynamics simulation study of AG10 and tafamidis binding to the Val122Ile transthyretin variant

Molecular dynamics (MD) simulations were used to investigate the binding of four ligands to the Val122Ile mutant of the protein transthyretin. Dissociation, misfolding, and subsequent aggregation of mutated transthyretin proteins are associated with the disease Familial Amyloidal Cardiomyopathy. The ligands investigated were the drug candidate AG10 and its decarboxy and N-methyl derivatives along with the drug tafamidis. These ligands bound to the receptor in two halogen binding pockets (HBP) designated AB and A’B’. Inter-ligand distances, solvent accessible surface areas, root mean squared deviation measurements, and extracted structures showed very little change in the AG10 ligands' conformations or locations within the HBP during the MD simulation. In addition, the AG10 ligands experienced stable, two-point interactions with the protein by forming hydrogen bonds with Ser-117 residues in both the AB and A’B’ binding pockets and Lysine-15 residues found near the surface of the receptor. Distance measurements showed these H-bonds formed simultaneously during the MD simulation. Removal of the AG10 carboxylate functional group to form decarboxy-AG10 disrupted this two-point interaction causing the ligand in the AB pocket to undergo a conformational change during the MD simulation. Likewise, addition of a methyl group to the AG10 hydrazone functional group also disrupted the two-point interaction by decreasing hydrogen bonding interactions with the receptor. Finally, MD simulations showed that the tafamidis ligands experienced fewer hydrogen bonding interactions than AG10 with the protein receptor. The tafamidis ligand in pocket A’B’ was also found to move deeper into the HBP during the MD simulation.     

Template-based docking of a prolactin receptor proline-rich motif octapeptide to FKBP12: implications for cytokine receptor signaling.

A conserved proline-rich motif (PRM) in the cytoplasmic domain of cytokine receptors has been suggested to be a signaling switch regulated by the action of the FK506 binding protein (FKBP) family of peptidylprolyl isomerases (O'Neal KD, Yu-Lee LY, Shearer WT, 1995, Ann NY Acad Sci 766:282-284). We have docked the prolactin receptor PRM (Ile1-Phe2-Pro3-Pro4-Val5-Pro6-Gly7-Pro8) to the ligand binding site of FKBP12. The procedure involved conformational search restricted by NMR restraints (O'Neal KD et al., 1996, Biochem J 315:833-844), energy minimization of the octapeptide conformers so obtained, template-based docking of a selected conformer to FKBP12, and energy refinement of the resulting complex. The template used was the crystal structure of a cyclic FK506-peptide hybrid bound to FKBP12. Val5-Pro6 of the PRM was taken to be the biologically relevant Xaa-Pro bond. The docked conformer is stabilized by two intramolecular hydrogen bonds, N7H7-->O4 and N2H2-->O8, and two intermolecular ones, Ile56; N-H-->O = C:Pro6 and Tyr82:O-H-->O = C:Gly7. This conformer features a Type I beta-turn and has extensive hydrophobic contacts with the FKBP12 binding surface. The observed interactions support the hypothesis that FKBP12 catalyzes cis-trans isomerization in the PRM when it is part of the longer cytoplasmic domain of a cytokine receptor, and suggest a significant role for the PRM in signal transduction.     

Crystal structure and molecular dynamics simulation of ubiquitin-like domain of murine parkin

Parkin is the gene product identified as the major cause of autosomal recessive juvenile Parkinsonism (AR-JP). Parkin, a ubiquitin ligase E3, contains a unique ubiquitin-like domain in its N-terminus designated Uld which is assumed to be a interaction domain with the Rpn 10 subunit of 26S proteasome. To elucidate the structural and functional role of Uld in parkin at the atomic level, the X-ray crystal structure of murine Uld was determined and a molecular dynamics simulation of wild Uld and its five mutants (K27N, R33Q, R42P, K48A and V56E) identified from AR-JP patients was performed. Murine Uld consists of two alpha helices [Ile23-Arg33 (alpha1) and Val56-Gln57 (alpha2)] and five beta strands [Met1-Phe7 (beta1), Tyr11-Asp18 (beta2), Leu41-Phe45 (beta3), Lys48-Pro51 (beta4) and Ser65-Arg72 (beta5)] and its overall structure is essentially the same as that of human ubiquitin with a 1.22 A rmsd for the backbone atoms of residues 1-76; however, the sequential identity and similarity between both molecules are 32% and 63%, respectively. This close resemblance is due to the core structure built by same hydrogen bond formations between and within the backbone chains of alpha1 and beta1/2/5 secondary structure elements and by nearly the same hydrophobic interactions formed between the nonpolar amino acids of their secondary structures. The side chain NetaH of Lys27 on the alpha1 helix was crucial to the stabilization of the spatial orientations of beta3 and beta4 strands, possible binding region with Rpn 10 subunit, through three hydrogen bonds. The MD simulations showed the K27N and R33Q mutations increase the structural fluctuation of these beta strands including the alpha1 helix. Reversely, the V56E mutant restricted the spatial flexibility at the periphery of the short alpha2 helix by the interactions between the polar atoms of Glu56 and Ser19 residues. However, a large fluctuation of beta4 strand with respect to beta5 strand was induced in the R42P mutant, because of the impossibility of forming paired hydrogen bonds of Pro for Arg42 in wild Uld. The X-ray structure showed that the side chains of Asp39, Gln40 and Arg42 at the N-terminal periphery of beta3 strand protrude from the molecular surface of Uld and participate in hydrogen bonds with the polar residues of neighboring Ulds. Thus, the MD simulation suggests that the mutation substitution of Pro for Arg42 not only causes the large fluctuation of beta3 strand in the Uld but also leads to the loss of the ability of Uld to trap the Rpn 10 subunit. In contrast, the MD simulation of K48A mutant showed little influence on the beta3-beta4 loop structure, but a large fluctuation of Lys48 side chain, suggesting the importance of flexibility of this side chain for the interaction with the Rpn 10 subunit. The present results would be important in elucidating the impaired proteasomal binding mechanism of parkin in AR-JP.     

Thr729 in human topoisomerase I modulates anti-cancer drug resistance by altering protein domain communications as suggested by molecular dynamics simulations

The role of Thr729 in modulating the enzymatic function of human topoisomerase I has been characterized by molecular dynamics (MD) simulation. In detail, the structural–dynamical behaviour of the Thr729Lys and the Thr729Pro mutants have been characterized because of their in vivo and in vitro functional properties evidenced in the accompanying paper. Both mutants can bind to the DNA substrate and are enzymatically active, but while Thr729Lys is resistant even at high concentration of the camptothecin (CPT) anti-cancer drug, Thr729Pro shows only a mild reduction in drug sensitivity and in DNA binding. MD simulations show that the Thr729Lys mutation provokes a structural perturbation of the CPT-binding pocket. On the other hand, the Thr729Pro mutant maintains the wild-type structural scaffold, only increasing its rigidity. The simulations also show the complete abolishment, in the Thr729Lys mutant, of the protein communications between the C-terminal domain (where the active Tyr723 is located) and the linker domain, that plays an essential role in the control of the DNA rotation, thus explaining the distributive mode of action displayed by this mutant.     

Effect of the N2 residue on the stability of the α-helix for all 20 amino acids

N2 is the second position in the alpha-helix. All 20 amino acids were placed in the N2 position of a synthetic helical peptide (CH(3)CO-[AXAAAAKAAAAKAAGY]-NH(2)) and the helix content was measured by circular dichroism spectroscopy at 273K. The dependence of peptide helicity on N2 residue identity has been used to determine a free-energy scale by analysis with a modified Lifson-Roig helix-coil theory that includes a parameter for the N2 energy (n2). The rank order of DeltaDeltaG((relative to Ala)) is Glu(-), Asp(-) > Ala > Glu(0), Leu, Val, Gln, Thr, Ile, Ser, Met, Asp(0), His(0), Arg, Cys, Lys, Phe > Asn, > Gly, His(+), Pro, Tyr. The results correlate very well with N2 propensities in proteins, moderately well with N1 and helix interior preferences, and not at all with N-cap preferences. The strongest energetic effects result from interactions with the helix dipole, which favors negative charges at the helix N terminus. Hydrogen bonds to side chains at N2, such as Gln, Ser, and Thr, are weak, despite occurring frequently in protein crystal structures, in contrast to the N-cap position. This is because N-cap hydrogen bonds are close to linear, whereas N2 hydrogen bonds have poor geometry. These results can be used to modify protein stability rationally, help design helices, and improve prediction of helix location and stability.     

Binding of ATP at the active site of human pancreatic glucokinase--nucleotide-induced conformational changes with possible implications for its kinetic cooperativity

Glucokinase (GK) is the central player in glucose-stimulated insulin release from pancreatic β-cells, and catalytic activation by α-D-glucose binding has a key regulatory function. Whereas the mechanism of this activation is well understood, on the basis of crystal structures of human GK, there are no similar structural data on ATP binding to the ligand-free enzyme and how it affects its conformation. Here, we report on a conformational change induced by the binding of adenine nucleotides to human pancreatic GK, as determined by intrinsic tryptophan fluorescence, using the catalytically inactive mutant form T228M to correct for the inner filter effect. Adenosine-5'-(β,γ-imido)triphosphate and ATP bind to the wild-type enzyme with apparent [L](0.5) (ligand concentration at half-maximal effect) values of 0.27±0.02 mm and 0.78±0.14 mm, respectively. The change in protein conformation was further supported by ATP inhibition of the binding of the fluorescent probe 8-anilino-1-naphthalenesulfonate and limited proteolysis by trypsin, and by molecular dynamic simulations. The simulations provide a first insight into the dynamics of the binary complex with ATP, including motion of the flexible surface/active site loop and partial closure of the active site cleft. In the complex, the adenosine moiety is packed between two α-helices and stabilized by hydrogen bonds (with Thr228, Thr332, and Ser336) and hydrophobic interactions (with Val412 and Leu415). Combined with enzyme kinetic analyses, our data indicate that the ATP-induced changes in protein conformation may have implications for the kinetic cooperativity of the enzyme.     

Computational design and biochemical characterization of maize nonspecific lipid transfer protein variants for biosensor applications

Lipid transfer proteins (LTPs) are a family of proteins that bind and transfer lipids. Utilizing the maize LTP, we have successfully engineered fluorescent reagentless biosensors for the natural ligand of LTPs; this was achieved by using computational protein design to remove a disulfide bridge and attaching a thio-reactive fluorophore. Conformational change induced by ligand titration is thought to affect the fluorescence of the fluorophore, allowing detection of ligand binding. Fluorescence measurements show that our LTP variants have affinity to palmitate that is consistent with wild-type LTP. These molecules have the potential to be utilized as scaffolds to design hydrophobic ligand biosensors or to serve as drug carriers.     

Contribution of Thr29 to the thermodynamic stability of goat alpha-lactalbumin as determined by experimental and theoretical approaches

The Thr29 residue in the hydrophobic core of goat alpha-lactalbumin (alpha-LA) was substituted with Val (Thr29Val) and Ile (Thr29Ile) to investigate the contribution of Thr29 to the thermodynamic stability of the protein. We carried out protein stability measurements, X-ray crystallographic analyses, and free energy calculations based on molecular dynamics simulation. The equilibrium unfolding transitions induced by guanidine hydrochloride demonstrated that the Thr29Val and Thr29Ile mutants were, respectively, 1.9 and 3.2 kcal/mol more stable than the wild-type protein (WT). The overall structures of the mutants were almost identical to that of WT, in spite of the disruption of the hydrogen bonding between the side-chain O-H group of Thr29 and the main-chain C=O group of Glu25. To analyze the stabilization mechanism of the mutants, we performed free energy calculations. The calculated free energy differences were in good agreement with the experimental values. The stabilization of the mutants was mainly caused by solvation loss in the denatured state. Furthermore, the O-H group of Thr29 favorably interacts with the C=O group of Glu25 to form hydrogen bonds and, simultaneously, unfavorably interacts electrostatically with the main-chain C=O group of Thr29. The difference in the free energy profile of the unfolding path between WT and the Thr29Ile mutant is discussed in light of our experimental and theoretical results.     

Distal pocket polarity in ligand binding to myoglobin: structural and functional characterization of a threonine68(E11) mutant

Site-directed mutagenesis studies have confirmed that the distal histidine in myoglobin stabilizes bound O2 by hydrogen bonding and have suggested that it is the polar character of the imidazole side chain rather than its size that limits the rate of ligand entry into the protein. We constructed an isosteric Val68 to Thr replacement in pig myoglobin (i) to investigate whether the O2 affinity could be increased by the introduction of a second hydrogen-bonding group into the distal heme pocket and (ii) to examine the influence of polarity on the ligand binding rates more rigorously. The 1.9-A crystal structure of Thr68 aquometmyoglobin confirms that the mutant and wild-type proteins are essentially isostructural and reveals that the beta-OH group of Thr68 is in a position to form hydrogen-bonding interactions both with the coordinated water molecule and with the main chain greater than C=O of residue 64. The rate of azide binding to the ferric form of the Thr68 mutant was 60-fold lower than that for the wild-type protein, consistent with the proposed stabilization of the coordinated water molecule. However, bound O2 is destabilized in the ferrous form of the mutant protein. The observed 17-fold lowering of the O2 affinity may be a consequence of the hydrogen-bonding interaction made between the Thr68 beta-OH group and the carbonyl oxygen of residue 64. Overall association rate constants for O2, NO, and alkyl isocyanide binding to ferrous pig myoglobin were 3-10-fold lower for the mutant compared to the wild-type protein, whereas that for CO binding was little affected.(ABSTRACT TRUNCATED AT 250 WORDS)