EF-Tu binding peptides identified,dissected, and affinity optimized by phage display

The highly abundant GTP binding protein elongation factor Tu (EF-Tu) fulfills multiple roles in bacterial protein biosynthesis. Phage-displayed peptides with high affinity for EF-Tu were selected from a library of approximately 4.7 x 10(11) different peptides. The lack of sequence homology among the identified EF-Tu ligands demonstrates promiscuous peptide binding by EF-Tu. Homolog shotgun scanning of an EF-Tu ligand was used to dissect peptide molecular recognition by EF-Tu. All homolog shotgun scanning selectants bound to EF-Tu with higher affinity than the starting ligand. Thus, homolog shotgun scanning can simultaneously optimize binding affinity and rapidly provide detailed structure activity relationships for multiple side chains of a polypeptide ligand. The reported peptide ligands do not compete for binding to EF-Tu with various antibiotic EF-Tu inhibitors, and could identify an EF-Tu peptide binding site distinct from the antibiotic inhibitory sites.     

Binding and stability determinants of the PPARgamma nuclear receptor-coactivator interface as revealed by shotgun alanine scanning and in vivo selection

We modified an existing selection for protein-protein interactions based on the fragment complementation of the enzyme DHFR. Using shotgun alanine scanning in conjunction with this selection, we analyzed the interaction of the nuclear receptor PPARgamma with two peptides derived from nuclear receptor coactivators SRC1 and TRAP220. A large binding epitope stretching between and including the charge clamp residues K301 and E471 of PPARgamma was identified as necessary for PPARgamma-coactivator interaction. To decouple protein stability from the propensity to form a receptor-coactivator interface, libraries of PPARgamma variants generated by shotgun scanning were further processed using a high-throughput screen measuring their in vivo stabilities. Our findings demonstrate that many of the residues that make up the binding epitope of PPARgamma are also crucial for the stability of the PPARgamma.     

Interdomain communication between the thiolation and thioesterase domains of EntF explored by combinatorial mutagenesis and selection

Thiolation (T) domains are protein way stations in natural product assembly lines. In the enterobactin synthetase, the T domain on EntF is recognized in cis by its catalytic partners: the EntF condensation (C), adenylation (A), and thioesterase (TE) domains. To assess surface features of the EntF T domain recognized by C, A, and TE, regions of the EntF T domain were submitted to shotgun alanine scanning and Ent production selection, which revealed residues that could not be substituted by Ala. EntF mutants bearing Ala in such positions were assayed in vitro for Ent production with EntEB, and for A-T, C-T, and T-TE communications. We concluded that G1027A and M1030A are specifically defective in acyl transfer from T to TE. These residues define an interaction surface between these two in cis domains in an NRPS module.     

Thermodynamics of binding of 2-methoxy-3-isopropylpyrazine and 2-methoxy-3-isobutylpyrazine to the major urinary protein

In the present study we examine the thermodynamics of binding of two related pyrazine-derived ligands to the major urinary protein, MUP-I, using a combination of isothermal titration calorimetry (ITC), X-ray crystallography, and NMR backbone (15)N and methyl side-chain (2)H relaxation measurements. Global thermodynamics data derived from ITC indicate that binding is driven by favorable enthalpic contributions, rather than the classical entropy-driven hydrophobic effect. Unfavorable entropic contributions from the protein backbone and side-chain residues in the vicinity of the binding pocket are partially offset by favorable entropic contributions at adjacent positions, suggesting a "conformational relay" mechanism whereby increased rigidity of residues on ligand binding are accompanied by increased conformational freedom of side chains in adjacent positions. The principal driving force governing ligand affinity and specificity can be attributed to solvent-driven enthalpic effects from desolvation of the protein binding pocket.     

Expanding the realm of ultrafast protein folding: gpW, a midsize natural single-domain with alpha+beta topology that folds downhill

All ultrafast folding proteins known to date are either very small in size (less than 45 residues), have an alpha-helix bundle topology, or have been artificially engineered. In fact, many of them share two or even all three features. Here we show that gpW, a natural 62-residue alpha+beta protein expected to fold slowly in a two-state fashion, folds in microseconds (i.e., from tau = 33 micros at 310 K to tau = 1.7 micros at 355 K). Thermodynamic analyses of gpW reveal probe dependent thermal denaturation, complex coupling between two denaturing agents, and differential scanning calorimetry (DSC) thermogram characteristic of folding over a negligible thermodynamic folding barrier. The free energy surface analysis of gpW folding kinetics also produces a marginal folding barrier of about thermal energy ( RT) at the denaturation midpoint. From these results we conclude that gpW folds in the downhill regime and is close to the global downhill limit. This protein seems to be poised toward downhill folding by a loosely packed hydrophobic core with low aromatic content, large stabilizing contributions from local interactions, and abundance of positive charges on the native surface. These special features, together with a complex functional role in bacteriophage lambda assembly, suggest that gpW has been engineered to fold downhill by natural selection.     

Critical residue that promotes protein dimerization: a story of partially exposed Phe25 in 14-3-3σ

Many proteins exist and function as oligomers. While hydrophobic interactions have been recognized as the major driving force for oligomerization, detailed molecular mechanisms for the assembly are unknown. Here, we used 14-3-3σ as a model protein and investigated the role of hydrophobic residues at the dimeric interface using MD simulations and coimmunoprecipitations. We found that a half-exposed and half-buried residue in the interface, Phe(25), plays a more important role in promoting homodimerization than the hydrophobic core residues by organizing both favorable hydrophobic and hydrophilic interactions. Phe(25) is critical in packing and stabilizing hydrophobic core residues. We conclude that the structural stability of hydrophobic cores is critical for a stable homodimer complex and this stable property can be bestowed by residues outside of hydrophobic core. The important organizing activity of Phe(25) for homodimerization of 14-3-3σ originates from its unique physical location, rigidity, size, and hydrophobicity. Thus, hydrophobic residues that are not deeply buried at the oligomeric interface may play important but different roles from the buried core residues and they may promote oligomerization by organizing co-operativity of core and other residues for favorable hydrophobic and electrostatic interactions.     

Conformational dynamics and stability of HP35 studied with 2D IR vibrational echoes

Two-dimensional infrared (2D IR) vibrational echo spectroscopy was used to measure the fast dynamics of two variants of chicken villin headpiece 35 (HP35). The CN of cyanophenylalanine residues inserted in the hydrophobic core were used as a vibrational probe. Experiments were performed on both singly (HP35-P) and doubly CN-labeled peptide (HP35-P(2)) within the wild-type sequence, as well as on HP-35 containing a singly labeled cyanophenylalanine and two norleucine mutations (HP35-P NleNle). There is a remarkable similarity between the dynamics measured in singly and doubly CN-labeled HP35, demonstrating that the presence of an additional CN vibrational probe does not significantly alter the dynamics of the small peptide. The substitution of two lysine residues by norleucines markedly improves the stability of HP35 by replacing charged with nonpolar residues, stabilizing the hydrophobic core. The results of the 2D IR experiments reveal that the dynamics of HP35-P are significantly faster than those of HP35-P NleNle. These observations suggest that the slower structural fluctuations in the hydrophobic core, indicating a more tightly structured core, may be an important contributing factor to HP35-P NleNle's increased stability.     

The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.     

Modulating protein structure with fluorous amino acids: increased stability and native-like structure conferred on a 4-helix bundle protein by hexafluoroleucine

There has recently been much interest in exploiting the unusual properties associated with fluorocarbons to modulate the physicochemical properties of proteins. Here we present a detailed investigation into the effect on structure and stability of systematically repacking the hydrophobic core of a model protein with the extensively fluorinated (fluorous) amino acid l-5,5,5,5',5',5'-hexafluoroleucine (hFLeu). The starting point was a 27-residue peptide, alpha(4)-H, that adopts an antiparallel 4-alpha-helix bundle structure, and in which the hydrophobic core comprises six layers of leucine residues introduced at the "a" and "d" positions of the canonical heptad repeat. A series of peptides were synthesized in which the central two (alpha(4)-F(2))(,) four (alpha(4)-F(4)), or all six layers (alpha(4)-F(6)) of the core were substituted hFLeu. The free energy of unfolding increases by 0.3 (kcal/mol)/hFLeu on repacking the central two layers and by an additional 0.12 (kcal/mol)/hFLeu on repacking additional layers, so that alpha(4)-F(6) is approximately 25% more stable than the nonfluorinated protein alpha(4)-H. One-dimensional proton, two-dimensional (1)H-(15)N HSQC, and (19)F NMR spectroscopies were used to examine the effect of fluorination on the conformational dynamics of the peptide. Unexpectedly, increasing the degree of fluorination also appears to result in peptides that possess a more structured backbone and less fluid hydrophobic core. The latter only occurs in alpha(4)-F(4) and alpha(4)-F(6), suggesting that crowding of the hFLeu residues may restrict the amplitude and/or time scales for rotation of the side chains.     

Thermodynamical properties of reaction intermediates during apoplastocyanin folding in time domain

Two intermediates observed for the folding process of apoplastocyanin (apoPC) were investigated by using a photoinduced triggering system combined with the transient grating and transient lens methods. The thermodynamic quantities, enthalpy, heat capacity, partial volume, and thermal expansion volume changes during the protein folding reaction were measured in time domain for the first time. An interesting observation is the positive enthalpy changes during the folding process. This positive enthalpy change must be compensated by positive entropy changes, which could be originated from the dehydration effect of hydrophobic residues and/or the translational entropy gain of bulk water molecules. Observed negative heat capacity change was explained by the dehydration effect of hydrophilic residues and/or motional confinement of amino acid side chains and water molecules in apoPC. The signs of the volume change and thermal expansion volume were different for two processes and these changes were interpreted in terms of the different relative contributions of the hydration and the dehydration of the hydrophilic residues. These results indicated two-step hydrophobic collapses in the early stage of the apoPC folding, but the nature of the dynamics was different.     

Role of water in mediating the assembly of Alzheimer amyloid-beta Abeta16-22 protofilaments

The role of water in promoting the formation of protofilaments (the basic building blocks of amyloid fibrils) is investigated using fully atomic molecular dynamics simulations. Our model protofilament consists of two parallel beta-sheets of Alzheimer Amyloid-beta 16-22 peptides (Ac-K(16)-L(17)-V(18)-F(19)-F(20)-A(21)-E(22)-NH2). Each sheet presents a distinct hydrophobic and hydrophilic face and together self-assemble to a stable protofilament with a core consisting of purely hydrophobic residues (L(17), F(19), A(21)), with the two charged residues (K(16), E(22)) pointing to the solvent. Our simulations reveal a subtle interplay between a water mediated assembly and one driven by favorable energetic interactions between specific residues forming the interior of the protofilament. A dewetting transition, in which water expulsion precedes hydrophobic collapse, is observed for some, but not all molecular dynamics trajectories. In the trajectories in which no dewetting is observed, water expulsion and hydrophobic collapse occur simultaneously, with protofilament assembly driven by direct interactions between the hydrophobic side chains of the peptides (particularly between F-F residues). For those same trajectories, a small increase in the temperature of the simulation (on the order of 20 K) or a modest reduction in the peptide-water van der Waals attraction (on the order of 10%) is sufficient to induce a dewetting transition, suggesting that the existence of a dewetting transition in simulation might be sensitive to the details of the force field parametrization.     

Hydrophobic core repacking in a coiled-coil dimer via phage display: insights into plasticity and specificity at a protein-protein interface

The coiled-coil, which consists of two or more interwoven amphiphilic alpha-helices, is formed by sequences that have a characteristic heptad repeat (abcdefg) where a and d are hydrophobic residues. Most efforts to elucidate the origins of coiled-coil pairing selectivity have focused on electrostatic interactions among side chains that flank the core (positions e and g) and on polar side chains that occur occasionally at core positions. We have used phage display to explore another source of coiled-coil specificity: steric matching among nonpolar side chains in the core. We introduced a destabilizing Leu-->Ala mutation into the core of one helix in a known heterodimer and then screened a phage-based library of potential partner helices in search of compensating mutations. We identified a new heterodimer pair (30 residues/helix) that is comparable in stability to the GCN4-p1 homodimer (33 residues/helix). Furthermore, the Leu-->Ala mutant shows specificity for its phage-derived partner over the original partner despite their similar sequences. These results show that a phage-based approach can provide unique insights on coiled-coil pairing preferences that should facilitate both the analysis of natural sequences and the development of specific dimerization motifs that are orthogonal to one another.     

Computational insight into small molecule inhibition of cyclophilins

Cyclophilins (Cyp) are a family of cellular enzymes possessing peptidyl-prolyl isomerase activity, which catalyze the cis-trans interconversion of proline-containing peptide bonds. The two most abundant family members, CypA and CypB, have been identified as valid drug targets for a wide range of diseases, including HCV, HIV, and multiple cancers. However, the development of small molecule inhibitors that possess nM potency and high specificity for a particular Cyp is difficult given the complete conservation of all active site residues between the enzymes. Monte Carlo statistical sampling coupled to free energy perturbation theory (MC/FEP) calculations have been carried out to elucidate the origin of the experimentally observed nM inhibition of CypA by acylurea-based derivatives and the >200-fold in vitro selectivity between CypA and CypB from aryl 1-indanylketone-based μM inhibitors. The computed free-energies of binding were in close accord with those derived from experiments. Binding affinity values for the inhibitors were determined to be dependent upon the stabilization strength of the nonbonded interactions provided toward two catalytic residues: Arg55 and Asn102 in CypA and the analogous Arg63 and Asn110 residues in CypB. Fine-tuning of the hydrophobic interactions allowed for enhanced potency among derivatives. The aryl 1-indanylketones are predicted to differentiate between the cyclophilins by using distinct binding motifs that exploit subtle differences in the active site arrangements. Ideas for the development of new selective compounds with the potential for advancement to low-nanomolar inhibition are presented.     

Fluoroalkyl and alkyl chains have similar hydrophobicities in binding to the "hydrophobic wall" of carbonic anhydrase

The hydrophobic effect, the free-energetically favorable association of nonpolar solutes in water, makes a dominant contribution to binding of many systems of ligands and proteins. The objective of this study was to examine the hydrophobic effect in biomolecular recognition using two chemically different but structurally similar hydrophobic groups, aliphatic hydrocarbons and aliphatic fluorocarbons, and to determine whether the hydrophobicity of the two groups could be distinguished by thermodynamic and biostructural analysis. This paper uses isothermal titration calorimetry (ITC) to examine the thermodynamics of binding of benzenesulfonamides substituted in the para position with alkyl and fluoroalkyl chains (H(2)NSO(2)C(6)H(4)-CONHCH(2)(CX(2))(n)CX(3), n = 0-4, X = H, F) to human carbonic anhydrase II (HCA II). Both alkyl and fluoroalkyl substituents contribute favorably to the enthalpy and the entropy of binding; these contributions increase as the length of chain of the hydrophobic substituent increases. Crystallography of the protein-ligand complexes indicates that the benzenesulfonamide groups of all ligands examined bind with similar geometry, that the tail groups associate with the hydrophobic wall of HCA II (which is made up of the side chains of residues Phe131, Val135, Pro202, and Leu204), and that the structure of the protein is indistinguishable for all but one of the complexes (the longest member of the fluoroalkyl series). Analysis of the thermodynamics of binding as a function of structure is compatible with the hypothesis that hydrophobic binding of both alkyl and fluoroalkyl chains to hydrophobic surface of carbonic anhydrase is due primarily to the release of nonoptimally hydrogen-bonded water molecules that hydrate the binding cavity (including the hydrophobic wall) of HCA II and to the release of water molecules that surround the hydrophobic chain of the ligands. This study defines the balance of enthalpic and entropic contributions to the hydrophobic effect in this representative system of protein and ligand: hydrophobic interactions, here, seem to comprise approximately equal contributions from enthalpy (plausibly from strengthening networks of hydrogen bonds among molecules of water) and entropy (from release of water from configurationally restricted positions).     

Drug Target Identification Using Affinity Core‐Shell Magnetic Nanoparticles and Mass Spectrometry

Affinity core‐shell magnetic nanoparticles (MNPs) were prepared for identifying the target proteins of drugs in the cell lysate when used in combination with nano‐high‐performance liquid chromatography tandem mass spectrometry (HPLC‐MS/MS)‐based shotgun proteomic analysis. A number of new potential targets of cyclosporine A (CsA) could be identified, owing to the high efficacy of the affinity MNPs in drug target identification. To the best of our knowledge, this is the first time to reveal such an abundant target spectrum of CsA.     

基于定点突变的抗克百威单链抗体的亲和力成熟

为获得亲和力更高的抗克百威(CBF)单链抗体(scFv), 从抗CBF scFv氨基酸序列出发, 通过同源模建获得抗体模型, 找出抗体中的活性口袋区域, 进而将小分子药物与抗体进行分子对接, 发现疏水作用和氢键对于抗体亲和力具有重要作用. 进一步对口袋内亲水氨基酸残基HArg40和LHis38进行模拟替换, 再进行分子对接分析, 发现当以亮氨酸为突变氨基酸时, 对接评分最高. 在此基础上, 通过构建突变scFv基因及可溶性表达, 采用ELISA法对进化后的单链抗体(evoscFv)进行了鉴定. 结果表明, evoscFv对CBF的IC₅₀值为18.11 μg/L, 低于野生型抗体的27.25 μg/L, 亲和解离常数K_d为4.06×10^-8 mol/L, 相对亲和力比野生型scFv提高了2.23倍, 说明通过分子对接分析及对抗体活性口袋中氨基酸残基进行替换, 获得了一个亲和力更高的突变体抗体.     

Interactions of peptide mimics of hyaluronic acid with the receptor for hyaluronan mediated motility (RHAMM)

Summary Using the hyaluronic acid (HA) binding region of the receptor for hyaluronan-mediated motility (RHAMM) as a model, a molecular perspective for peptide mimicry of the natural ligand was established by comparing the interaction sites of HA and unnatural peptide–ligands to RHAMM. This was accomplished by obtaining a series of octapeptide–ligands through screening experiments that bound to the HA binding domains of RHAMM (amino acids 517–576) and could be displaced by HA. These molecules were computationally docked onto a three-dimensional NMR based model of RHAMM. The NMR model showed that RHAMM(517–576) was a set of three helices, two of which contained the HA binding domains (HABDs) flanking a central groove. The structure was stabilized by hydrophobic interactions from four pairs of Val and Ile side chains extending into the groove. The presence of solvent exposed, positively charged side chains spaced 11Å apart matched the spacing of negative charges on HA. Docking experiments using flexible natural and artificial ligands demonstrated that HA and peptide–mimetics preferentially bound to the second helix that contains HABD-2. Three salt bridges between HA carboxylates and Lys548, Lys553 and Lys560 and two hydrophobic interactions involving Val538 and Val559 were predicted to stabilize the RHAMM-HA complex. The high affinity peptides and HA utilized the same charged residues, with additional contacts to other basic residues. However, hydrophobic contacts do not contribute to affinity for peptide ligand-RHAMM complexes. These results offer insight into how selectivity is achieved in the binding of HA to RHAMM, and how peptide competitors may compete for binding with HA on a single hyaladherin.     

Dispersion interactions govern the strong thermal stability of a protein

Rubredoxin from the hyperthermophile Pyrococcus furiosus (Pf Rd) is an extremely thermostable protein, which makes it an attractive subject of protein folding and stability studies. A fundamental question arises as to what the reason for such extreme stability is and how it can be elucidated from a complex set of interatomic interactions. We addressed this issue first theoretically through a computational analysis of the hydrophobic core of the protein and its mutants, including the interactions taking place inside the core. Here we show that a single mutation of one of phenylalanine's residues inside the protein's hydrophobic core results in a dramatic decrease in its thermal stability. The calculated unfolding Gibbs energy as well as the stabilization energy differences between a few core residues follows the same trend as the melting temperature of protein variants determined experimentally by microcalorimetry measurements. NMR spectroscopy experiments have shown that the only part of the protein affected by mutation is the reasonably rearranged hydrophobic core. It is hence concluded that stabilization energies, which are dominated by London dispersion, represent the main source of stability of this protein.     

A structural study of model peptides derived from HIV-1 integrase central domain

The HIV-1 integrase (IN) catalyzes the integration of viral DNA in the human genome. In vitro the enzyme displays an equilibrium of monomers, dimers, tetramers and larger oligomers. However, its functional oligomeric form in vivo is not known. We report a study of the auto-associative properties of three peptides denoted K156, E156 and E159. These derive from the alpha4 helix of the IN catalytic core. The alpha4 helix is an amphipatic helix exposed at the surface of the protein and could be involved in the oligomerization process through its hydrophobic face. The peptides were obtained from the replacement of several amino acid residues by more helicogenic ones in the alpha4 helix peptide. K156 carries the basic residues Lys156 and Lys159, which have been shown important for the binding of IN to viral DNA. In E156 and E159 they are replaced with the acidic residue Glu. A fourth peptide K(E)156 obtained from the replacement of hydrophobic residues with Glu in K156 in order to abolish the auto-associative properties is used as a negative control. The capacity shown by peptides for alpha-helical formation is demonstrated by circular dichroism (CD) analysis performed in aqueous solution and in aqueous trifluoroethanol (TFE) mixtures. Both electrospray ionization mass spectrometry (ESI-MS) and glutaraldehyde chemical cross-linking show that peptides adopt different solvent-dependent equilibriums of monomers, dimers, trimers and tetramers. Oligomerization of peptides in aqueous solution is related to their ability to form helical structures. Addition of a small amount of TFE (<10%) stimulates helix stabilization and the interhelical hydrophobic contacts. Higher amounts of TFE alter the hydrophobic contacts and disrupt the oligomeric species. In addition to hydrophobic interactions, the patterns indicate that the biologically important Lys156 and Lys159 residues also participate in helix association. K(E)156 despite its ability to adopt a helical structure is unable to associate into oligomers, demonstrating the importance of hydrophobic contacts for oligomerization. Thus, the designed peptides provide us information on the functional properties of the alpha4 IN that seems to hold a dual role in DNA recognition and protein oligomerization.     

Structure-based drug design: the discovery of novel nonpeptide orally active inhibitors of human renin

BACKGROUND: The aspartic proteinase renin plays an important physiological role in the regulation of blood pressure. It catalyses the first step in the conversion of angiotensinogen to the hormone angiotensin II. In the past, potent peptide inhibitors of renin have been developed, but none of these compounds has made it to the end of clinical trials. Our primary aim was to develop novel nonpeptide inhibitors. Based on the available structural information concerning renin-substrate interactions, we synthesized inhibitors in which the peptide portion was replaced by lipophilic moieties that interact with the large hydrophobic S1/S3-binding pocket in renin. RESULTS: Crystal structure analysis of renin-inhibitor complexes combined with computational methods were employed in the medicinal-chemistry optimisation process. Structure analysis revealed that the newly designed inhibitors bind as predicted to the S1/S3 pocket. In addition, however, these compounds interact with a hitherto unrecognised large, distinct, sub-pocket of the enzyme that extends from the S3-binding site towards the hydrophobic core of the enzyme. Binding to this S3(sp) sub-pocket was essential for high binding affinity. This unprecedented binding mode guided the drug-design process in which the mostly hydrophobic interactions within subsite S3(sp) were optimised. CONCLUSIONS: Our design approach led to compounds with high in vitro affinity and specificity for renin, favourable bioavailability and excellent oral efficacy in lowering blood pressure in primates. These renin inhibitors are therefore potential therapeutic agents for the treatment of hypertension and related cardiovascular diseases.