Yaakov Levy

Effect of glycosylation on protein folding: A close look at thermodynamic stabilization

Glycosylation is one of the most common posttranslational modifications to occur in protein biosynthesis, yet its effect on the thermodynamics and kinetics of proteins is poorly understood. A minimalist model based on the native protein topology, in which each amino acid and sugar ring was represented by a single bead, was used to study the effect of glycosylation on protein folding. We studied in silico the folding of 63 engineered SH3 domain variants that had been glycosylated with different numbers of conjugated polysaccharide chains at different sites on the proteins surface. Thermal stabilization of the protein by the polysaccharide chains was observed in proportion to the number of attached chains. Consistent with recent experimental data, the degree of thermal stabilization depended on the position of the glycosylation sites, but only very weakly on the size of the glycans. A thermodynamic analysis showed that the origin of the enhanced protein stabilization by glycosylation is destabilization of the unfolded state rather than stabilization of the folded state. The higher free energy of the unfolded state is enthalpic in origin because the bulky polysaccharide chains force the unfolded ensemble to adopt more extended conformations by prohibiting formation of a residual structure. The thermodynamic stabilization induced by glycosylation is coupled with kinetic stabilization. The effects introduced by the glycans on the biophysical properties of proteins are likely to be relevant to other protein polymeric conjugate systems that regularly occur in the cell as posttranslational modifications or for biotechnological purposes.

Protein sliding along DNA: dynamics and structural characterization.

Efficient search of DNA by proteins is fundamental to the control of cellular regulatory processes. It is currently believed that protein sliding, hopping, and transfer between adjacent DNA segments, during which the protein nonspecifically interacts with DNA, are central to the speed of their specific recognition. In this study, we focused on the structural and dynamic features of proteins when they scan the DNA. Using a simple computational model that represents protein-DNA interactions by electrostatic forces, we identified that the protein makes use of identical binding interfaces for both nonspecific and specific DNA interactions. Accordingly, in its one-dimensional diffusion along the DNA, the protein is bound at the major groove and performs a helical motion, which is stochastic and driven by thermal diffusion. A microscopic structural insight into sliding from our model, which is governed by electrostatic forces, corroborates previous experimental studies suggesting that the active site of some regulatory proteins continually faces the interior of the DNA groove while sliding along sugar-phosphate rails. The diffusion coefficient of spiral motion along the major groove of the DNA is not affected by salt concentration, but the efficiency of the search can be significantly enhanced by increasing salt concentration due to a larger number of hopping events. We found that the most efficient search comprises approximately 20% sliding along the DNA and approximately 80% hopping and three-dimensional diffusion. The presented model that captures various experimental features of facilitated diffusion has the potency to address other questions regarding the nature of DNA search, such as the sliding characteristics of oligomeric and multidomain DNA-binding proteins that are ubiquitous in the cell.

Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

Assessment of CASP8 structure predictions for template free targets

The biennial CASP experiment is a crucial way to evaluate, in an unbiased way, the progress in predicting novel 3D protein structures. In this article, we assess the quality of prediction of template free models, that is, ab initio prediction of 3D structures of proteins based solely on the amino acid sequences, that is, proteins that did not have significant sequence identity to any protein in the Protein Data Bank. There were 13 targets in this category and 102 groups submitted predictions. Analysis was based on the GDT_TS analysis, which has been used in previous CASP experiments, together with a newly developed method, the OK_Rank, as well as by visual inspection. There is no doubt that in recent years many obstacles have been removed on the long and elusive way to deciphering the protein‐folding problem. Out of the 13 targets, six were predicted well by a number of groups. On the other hand, it must be stressed that for four targets, none of the models were judged to be satisfactory. Thus, for template free model prediction, as evaluated in this CASP, successes have been achieved for most targets; however, a great deal of research is still required, both in improving the existing methods and in development of new approaches. Proteins 2009.

DNA search efficiency is modulated by charge composition and distribution in the intrinsically disordered tail

Intrinsically disordered tails are common in DNA-binding proteins and can affect their search efficiency on nonspecific DNA by promoting the brachiation dynamics of intersegment transfer. During brachiation, the protein jumps between distant DNA regions via an intermediate state in which the tail and globular moieties are bound to different DNA segments. While the disordered tail must be long and positively charged to facilitate DNA search, the effect of its residue sequence on brachiation is unknown. We explored this issue using the NK-2 and Antp homeodomain transcription factors. We designed 566 NK-2 tail-variants and 55 Antp tail-variants having different net charges and positive charge distributions and studied their dynamics and DNA search efficiencies using coarse-grained molecular dynamics simulations. More intersegment transfers occur when the tail is moderately positively charged and the positive charges are clustered together in the middle of the tail or towards its N terminus. The presence of a negatively charged residue does not significantly affect protein brachiation, although it is likely that the presence of many negatively charged residues will complicate the DNA search mechanism. A bioinformatic analysis of 1,384 wild-type homeodomains illustrates that the charge composition and distribution in their N-tail sequences are consistent with an optimal charge pattern to promote intersegment transfer. Our study thus indicates that the residue sequence of the disordered tails of DNA-binding proteins has unique characteristics that were evolutionarily selected to achieve optimized function and suggests that the sequence-structure-function paradigm known for structured proteins is valid for intrinsically disordered proteins as well.

Frustration in protein–DNA binding influences conformational switching and target search kinetics

Rapid recognition of DNA target sites involves facilitated diffusion through which alternative sites are searched on genomic DNA. A key mechanism facilitating the localization of the target by a DNA-binding protein (DBP) is one-dimensional diffusion (sliding) in which electrostatic forces attract the protein to the DNA. As the protein reaches its target DNA site, it switches from purely electrostatic binding to a specific set of interactions with the DNA bases that also involves hydrogen bonding and van der Waals forces. High overlap between the DBP patches used for nonspecific and specific interactions with DNA may enable an immediate transition between the two binding modes following target site localization. By contrast, an imperfect overlap may result in greater frustration between the two potentially competing binding modes and consequently slower switching between them. A structural analysis of 125 DBPs indicates frustration between the two binding modes that results in a large difference between the orientations of the protein to the DNA when it slides compared to when it specifically interacts with DNA. Coarse-grained molecular dynamics simulations of in silico designed peptides comprising the full range of frustrations between the two interfaces show slower transition from nonspecific to specific DNA binding as the overlap between the patches involved in the two binding modes decreases. The complex search kinetics may regulate the search by eliminating trapping of the protein in semispecific sites while sliding.

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the beta-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.

Asymmetrical roles of zinc fingers in dynamic DNA-scanning process by the inducible transcription factor Egr-1

Egr-1 is an inducible transcription factor that recognizes 9-bp target DNA sites via three zinc finger domains and activates genes in response to cellular stimuli such as synaptic signals and vascular stresses. Using spectroscopic and computational approaches, we have studied structural, dynamic, and kinetic aspects of the DNA-scanning process in which Egr-1 is nonspecifically bound to DNA and perpetually changes its location on DNA. Our NMR data indicate that Egr-1 undergoes highly dynamic domain motions when scanning DNA. In particular, the zinc finger 1 (ZF1) of Egr-1 in the nonspecific complex is mainly dissociated from DNA and undergoes collective motions on a nanosecond timescale, whereas zinc fingers 2 and 3 (ZF2 and ZF3, respectively) are bound to DNA. This was totally unexpected because the previous crystallographic studies of the specific complex indicated that all of Egr-1’s three zinc fingers are equally involved in binding to a target DNA site. Mutations that are expected to enhance ZF1’s interactions with DNA and with ZF2 were found to reduce ZF1’s domain motions in the nonspecific complex suggesting that these interactions dictate the dynamic behavior of ZF1. By experiment and computation, we have also investigated kinetics of Egr-1’s translocation between two nonspecific DNA duplexes. Our data on the wild type and mutant proteins suggest that the domain dynamics facilitate Egr-1’s intersegment transfer that involves transient bridging of two DNA sites. These results shed light on asymmetrical roles of the zinc finger domains for Egr-1 to scan DNA efficiently in the nucleus.

Capillarity theory for the fly-casting mechanism.

Biomolecular folding and function are often coupled. During molecular recognition events, one of the binding partners may transiently or partially unfold, allowing more rapid access to a binding site. We describe a simple model for this fly-casting mechanism based on the capillarity approximation and polymer chain statistics. The model shows that fly casting is most effective when the protein unfolding barrier is small and the part of the chain which extends toward the target is relatively rigid. These features are often seen in known examples of fly casting in protein–DNA binding. Simulations of protein–DNA binding based on well-funneled native-topology models with electrostatic forces confirm the trends of the analytical theory.

Mutations as trapdoors to two competing native conformations of the Rop-dimer

Conformational transitions play a central role in regulating protein function. Structure-based models with multiple basins have been used to understand the mechanisms governing these transitions. A model able to accommodate multiple folding basins is proposed to explore the mutational effects in the folding of the Rop-dimer (Rop). In experiments, Rop mutants show unusually strong increases in folding rates with marginal effects on stability. We investigate the possibility of two competing conformations representing a parallel (P) and the wild-type antiparallel (AP) arrangement of the monomers as possible native conformations. We observe occupation of both distinct states and characterize the transition pathways. An interesting observation from the simulations is that, for equivalent energetic bias, the transition to the P basin (non-wild-type basin) shows a lower free-energy barrier. Thus, the rapid kinetics observed in experiments appear to be the result of two competing states with different kinetic behavior, triggered upon mutation by the opening of a trapdoor arising from Rops symmetric structure. The general concept of having competing conformations for the native state goes beyond explaining Rops mutational behaviors and can be applied to other systems. A switch between competing native structures might be triggered by external factors to allow, for example, allosteric control or signaling.