Target of rapamycin FATC domain as a general membrane anchor: The FKBP‐12 like domain of FKBP38 as a case study

Increased efforts have been undertaken to better understand the formation of signaling complexes at cellular membranes. Since the preparation of proteins containing a transmembrane domain or a prenylation motif is generally challenging an alternative membrane anchoring unit that is easy to attach, water‐soluble and binds to different membrane mimetics would find broad application. The 33‐residue long FATC domain of yeast TOR1 (y1fatc) fulfills these criteria and binds to neutral and negatively charged micelles, bicelles, and liposomes. As a case study, we fused it to the FKBP506‐binding region of the protein FKBP38 (FKBP38‐BD) and used ¹H–¹⁵N NMR spectroscopy to characterize localization of the chimeric protein to micelles, bicelles, and liposomes. Based on these and published data for y1fatc, its use as a C‐terminally attachable membrane anchor for other proteins is compatible with a wide range of buffer conditions (pH circa 6–8.5, NaCl 0 to >150 mM, presence of reducing agents, different salts such as MgCl₂ and CaCl₂). The high water‐solubility of y1fatc enables its use for titration experiments against a membrane‐localized interaction partner of the fused target protein. Results from studies with peptides corresponding to the C‐terminal 17–11 residues of the 33‐residue long domain by 1D ¹H NMR and CD spectroscopy indicate that they still can interact with membrane mimetics. Thus, they may be used as membrane anchors if the full y1fatc sequence is disturbing or if a chemically synthesized y1fatc peptide shall be attached by native chemical ligation, for example, unlabeled peptide to ¹⁵N‐labeled target protein for NMR studies.     

Cripto stabilizes GRP78 on the cell membrane

The ER resident chaperone molecule GRP78 has been shown to translocate to the cell surface where it associates with Cripto and signals cell growth, playing a still partially understood role in tumorigenesis. Consequently, a better understanding of GRP78 topology and structure at the surface of cancer cells represents an important step in the development of a new class of therapeutics. Here, we used a set of programs for creation of a complex containing GRP78 and Cripto proteins. We elucidated possible interactions of GRP78, Cripto, and their complex with the membrane. Using molecular dynamics simulations, we demonstrated that Cripto binding to GRP78 completely changes the dynamics of its behavior on the membrane, not allowing GRP78 to disconnect from it, thus enabling GRP78 tumorigenic functions.     

Sequence and conformational preferences at termini of α‐helices in membrane proteins: Role of the helix environment

α‐helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These α‐helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C‐termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze α‐helices in a high‐resolution dataset of integral α‐helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C‐termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near‐helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins. Proteins 2014; 82:3420–3436. © 2014 Wiley Periodicals, Inc.     

SIMPLE estimate of the free energy change due to aliphatic mutations: superior predictions based on first principles

The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.     

Subunit-specific backbone NMR assignments of a 64 kDa trp repressor/DNA complex: A role for N-terminal residues in tandem binding

Deuterium decoupled, triple resonance NMR spectroscopy was used to analyze complexes of ²H,¹⁵N,¹³C labelled intact and (des2–7) trp repressor (2–7 trpR) from E. coli bound in tandem to an idealized 22 basepair trp operator DNA fragment and the corepressor 5-methyltryptophan. The DNA sequence used here binds two trpR dimers in tandem resulting in chemically nonequivalent environments for the two subunits of each dimer. Sequence- and subunit-specific NMR resonance assignments were made for backbone ¹HN, ¹⁵N, ¹³C positions in both forms of the protein and for¹³ C in the intact repressor. The differences in backbone chemical shifts between the two subunits within each dimer of 2–7 trpR reflect dimer-dimer contacts involving the helix-turn-helix domains and N-terminal residues consistent with a previously determined crystal structure [Lawson and Carey (1993) Nature, 366, 178–182]. Comparison of the backbone chemical shifts of DNA-bound 2–7 trpR with those of DNA-bound intact trpR reveals significant changes for those residues involved in N-terminal-mediated interactions observed in the crystal structure. In addition, our solution NMR data contain three sets of resonances for residues 2–12 in intact trpR suggesting that the N-terminus has multiple conformations in the tandem complex. Analysis of C chemical shifts using a chemical shift index (CSI) modified for deuterium isotope effects has allowed a comparison of the secondary structure of intact and 2–7 tprR. Overall these data demonstrate that NMR backbone chemical shift data can be readily used to study specific structural details of large protein complexes.     

Surface recognition elements of membrane protein oligomerization

Although certain membrane proteins are functional as monomeric polypeptides, others must assemble into oligomers to carry out their biological roles. High-resolution membrane protein structures provide a valuable resource for examining the sequence features that facilitate-or preclude-assembly of membrane protein monomers into multimeric structures. Here we have utilized a data set of 28 high-resolution alpha-helical membrane protein structures comprising 32 nonredundant polypeptides to address this issue. The lipid-exposed surfaces of membrane proteins that have reached their fully assembled and functional biological units have been compared with those of the individual subunits that build quaternary structures. Though the overall amino acid composition of each set of surfaces is similar, a key distinction-the distribution of small-xxx-small motifs-delineates subunits from membrane proteins that have reached a functioning oligomeric state. Quaternary structure formation may therefore be dictated by small-xxx-small motifs that are not satisfied by intrachain contacts.     

Water‐mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water‐mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins—which is a well‐accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free‐energy calculation for a thermodynamic cycle of a protein folding using a liquid‐state density functional theory. Here, the solvation‐free energy for a leucine zipper formation was examined in the coiled‐coil protein GCN4‐p1, a typical model for hydrophobic interactions, which demonstrated that water‐mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.     

Interplay between the electrostatic membrane potential and conformational changes in membrane proteins

Transmembrane electrostatic membrane potential is a major energy source of the cell. Importantly, it determines the structure as well as function of charge‐carrying membrane proteins. Here, we discuss the relationship between membrane potential and membrane proteins, in particular whether the conformation of these proteins is integrally connected to the membrane potential. Together, these concepts provide a framework for rationalizing the types of conformational changes that have been observed in membrane proteins and for better understanding the electrostatic effects of the membrane potential on both reversible as well as unidirectional dynamic processes of membrane proteins.     

Structures of designed armadillo repeat proteins binding to peptides fused to globular domains

Designed armadillo repeat proteins (dArmRP) are α‐helical solenoid repeat proteins with an extended peptide binding groove that were engineered to develop a generic modular technology for peptide recognition. In this context, the term “peptide” not only denotes a short unstructured chain of amino acids, but also an unstructured region of a protein, as they occur in termini, loops, or linkers between folded domains. Here we report two crystal structures of dArmRPs, in complex with peptides fused either to the N‐terminus of Green Fluorescent Protein or to the C‐terminus of a phage lambda protein D. These structures demonstrate that dArmRPs bind unfolded peptides in the intended conformation also when they constitute unstructured parts of folded proteins, which greatly expands possible applications of the dArmRP technology. Nonetheless, the structures do not fully reflect the binding behavior in solution, that is, some binding sites remain unoccupied in the crystal and even unexpected peptide residues appear to be bound. We show how these differences can be explained by restrictions of the crystal lattice or the composition of the crystallization solution. This illustrates that crystal structures have to be interpreted with caution when protein–peptide interactions are characterized, and should always be correlated with measurements in solution.     

Protein unfolding at interfaces: slow dynamics of alpha-helix to beta-sheet transition

A two-phase sequential dynamic change in the secondary structure of hen egg lysozyme (Lys) adsorbed on solid substrates was observed. The first phase involved fast conversion of alpha-helix to random/turns (within the first minute or at very low coverage or high substrate wettability) with no perceptible change in beta-sheet content. The second phase (1-1200 min), however, involved a relatively slow conversion from alpha-helix to beta-sheet without a noticeable change in random/turns. An important finding of this work is that the concentration of lysozyme in the adsorbed state has a substantial effect on the fractional content of secondary structures. Attenuated total reflection Fourier transform infrared (ATR/FTIR) spectroscopy, along with a newly-developed optimization algorithm for predicting the content of secondary structure motifs, was used to correlate the secondary structure and the amount of adsorbed lysozyme with the surface wettability of six different flat nanoporous substrates. Although three independent variables, surface wettability, solution concentration and time for adsorption, were used to follow the fractional structural changes of lysozyme, the results were all normalized onto a single plot with the amount adsorbed as the universal independent variable. Consequently, lateral interactions among proteins likely drive the transition process. Direct intermolecular force adhesion measurements between lysozyme and different functionalized self-assembled alkanethiol monolayers confirm that hydrophobic surfaces interact strongly with proteins. The lysozyme-unfolding pathway during early adsorption appears to be similar to that predicted by published molecular modeling results.     

Ligand loading at the surface of an optical biosensor and its effect upon the kinetics of protein–protein interactions

Optical biosensors are finding increasing use in the determination of kinetic and equilibrium constants for a variety of biomolecular interactions. Usually these biosensors require one biomolecule, the ligand, to be covalently attached to a hydrogel matrix which itself is bonded to the sensing surface. The ligands partner, the ligate, then binds from solution resulting in a measurable change in response which the instrument records as a function of time. Although in many cases, optical biosensors are used in order to obtain parameters that relate to interactions in solution, it is becoming clear that measurements involving the interaction of ligate with immobilized ligands on surfaces require careful experimental design. Here we report on how the density of ligand loading within the hydogel matrix affects the measured interaction kinetics. It is found that crowding of ligand within this matrix results in a significant reduction in the measured association rate constant, with a corresponding effect in the calculated overall affinity. However, measurements at low ligand loadings show association rate constants that are comparable to those measured in solution. Clearly, where this comparison is required, it is important to perform measurements under such conditions. © 1997 John Wiley & Sons, Ltd.     

Thermodynamics,the structure of integral membrane proteins,and transport

Membranes are structures whose lipid and protein components are at, or close to, equilibrium in the plane of the membrane, but are not at equilibrium across the membrane. The thermodynamic tendency of ionic and highly polar molecules to be in contact with water rather than with nonpolar media (hydrophilic interactions) is important in determining these equilibrium and nonequilibrium states. In this paper, we speculate about the structures and orientations of integral proteins in a membrane, and about how the equilibrium and nonequilibrium features of such structures and orientations might be influenced by the special mechanisms of biosynthesis, processing, and membrane insertion of these proteins. The relevance of these speculations to the mechanisms of the translocation event in membrane transport is discussed, and specific protein models of transport that have been proposed are analyzed.     

Translocons, thermodynamics, and the folding of membrane proteins

Recent three-dimensional structures of helical membrane proteins present new challenges for the prediction of structure from amino acid sequence. Membrane proteins reside stably in a thermodynamic free energy minimum after release into the membrane's bilayer fabric from the translocon complex. This means that structure prediction is primarily a problem of physical chemistry. But the folding processes within the translocon must also be considered. A distilled overview of the physical principles of membrane protein stability is presented, and extended to encompass translocon-assisted folding.     

A fast and simple method for probing the interaction of peptides and proteins with lipids and membrane-mimetics using GB1 fusion proteins and NMR spectroscopy

The expression of peptides and proteins as fusions to the B1 domain of streptococcal protein G (GB1) is very popular since GB1 often improves the solubility of the target protein and because the first purification step using IgG affinity chromatography is simple and efficient. However, the following protease digest is not always complete or can result in a digest of the target protein. In addition, a further purification step such as RP-HPLC has to be used to get rid of the GB1 tag and undigested fusion protein. Because the protease digest and the following purification step are not only time-consuming but generally also expensive, we tested if GB1 fusion proteins can directly be used for NMR interaction studies using lipids or membrane-mimetics. Based on NMR binding studies using only the GB1 part, this fusion tag does not significantly interact with different membrane-mimetics such as micelles, bicelles, or liposomes. Thus spectral changes observed using GB1-fusion proteins indicate lipid- and membrane interactions of the target protein. The method was initially established to probe membrane interactions of a large number of mutants of the FATC domain of the ser/thr kinase TOR. To demonstrate the usefulness of the approach, we show NMR binding data for the wild type protein and a leucine to alanine mutant.     

Flexible docking of peptides to proteins using CABS‐dock

Molecular docking of peptides to proteins can be a useful tool in the exploration of the possible peptide binding sites and poses. CABS‐dock is a method for protein–peptide docking that features significant conformational flexibility of both the peptide and the protein molecules during the peptide search for a binding site. The CABS‐dock has been made available as a web server and a standalone package. The web server is an easy to use tool with a simple web interface. The standalone package is a command‐line program dedicated to professional users. It offers a number of advanced features, analysis tools and support for large‐sized systems. In this article, we outline the current status of the CABS‐dock method, its recent developments, applications, and challenges ahead.     

Amino acid substitutions at protein–protein interfaces that modulate the oligomeric state

Despite similarities in their sequence and structure, there are a number of homologous proteins that adopt various oligomeric states. Comparisons of these homologous protein pairs, in terms of residue substitutions at the protein–protein interfaces, have provided fundamental characteristics that describe how proteins interact with each other. We have prepared a dataset composed of pairs of related proteins with different homo‐oligomeric states. Using the protein complexes, the interface residues were identified, and using structural alignments, the shadow‐interface residues have been defined as the surface residues that align with the interface residues. Subsequently, we investigated residue substitutions between the interfaces and the shadow interfaces. Based on the degree of the contributions to the interactions, the aligned sites of the interfaces and shadow interfaces were divided into primary and secondary sites; the primary sites are the focus of this work. The primary sites were further classified into two groups (i.e. exposed and buried) based on the degree to which the residue is buried within the shadow interfaces. Using these classifications, two simple mechanisms that mediate the oligomeric states were identified. In the primary‐exposed sites, the residues on the shadow interfaces are replaced by more hydrophobic or aromatic residues, which are physicochemically favored at protein–protein interfaces. In the primary‐buried sites, the residues on the shadow interfaces are replaced by larger residues that protrude into other proteins. These simple rules are satisfied in 23 out of 25 Structural Classification of Proteins (SCOP) families with a different‐oligomeric‐state pair, and thus represent a basic strategy for modulating protein associations and dissociations. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Identification of compounds with binding affinity to proteins via magnetization transfer from bulk water*

A powerful screening by NMR methodology (WaterLOGSY), based on transfer of magnetization from bulk water, for the identification of compounds that interact with target biomolecules (proteins, RNA and DNA fragments) is described. The method exploits efficiently the large reservoir of H₂O magnetization. The high sensitivity of the technique reduces the amount of biomolecule and ligands needed for the screening, which constitutes an important requirement for high throughput screening by NMR of large libraries of compounds. Application of the method to a compound mixture against the cyclin-dependent kinase 2 (cdk2) protein is presented.     

Using collections of structural models to predict changes of binding affinity caused by mutations in protein–protein interactions

Protein–protein interactions (PPIs) in all the molecular aspects that take place both inside and outside cells. However, determining experimentally the structure and affinity of PPIs is expensive and time consuming. Therefore, the development of computational tools, as a complement to experimental methods, is fundamental. Here, we present a computational suite: MODPIN, to model and predict the changes of binding affinity of PPIs. In this approach we use homology modeling to derive the structures of PPIs and score them using state‐of‐the‐art scoring functions. We explore the conformational space of PPIs by generating not a single structural model but a collection of structural models with different conformations based on several templates. We apply the approach to predict the changes in free energy upon mutations and splicing variants of large datasets of PPIs to statistically quantify the quality and accuracy of the predictions. As an example, we use MODPIN to study the effect of mutations in the interaction between colicin endonuclease 9 and colicin endonuclease 2 immune protein from Escherichia coli. Finally, we have compared our results with other state‐of‐art methods.     

Thermodynamics of binding by calmodulin correlates with target peptide α‐helical propensity

In this work, we have examined contributions to the thermodynamics of calmodulin (CaM) binding from the intrinsic propensity for target peptides to adopt an α‐helical conformation. CaM target sequences are thought to commonly reside in disordered regions within proteins. Using the ability of TFE to induce α‐helical structure as a proxy, the six peptides studied range from having almost no propensity to adopt α‐helical structure through to a very high propensity. This despite all six peptides having similar CaM‐binding affinities. Our data indicate there is some correlation between the deduced propensities and the thermodynamics of CaM binding. This finding implies that molecular recognition features, such as CaM target sequences, may possess a broad range of propensities to adopt local structure. Given that these peptides bind to CaM with similar affinities, the data suggest that having a higher propensity to adopt α‐helical structure does not necessarily result in tighter binding, and that the mechanism of CaM binding is very dependent on the nature of the substrate sequence. Proteins 2013. © 2012 Wiley Periodicals, Inc.