Search

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.

International audience; 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 flycasting mechanism based on the capillarity approximation and polymer chain statistics. The model shows that flycasting is most effective when the protein unfolding barrier is small and the part of the chain which extends towards the target is relatively rigid. These features are often seen in known examples of flycasting 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.

In the cell, protein complexes form by relying on specific interactions between their monomers. Excluded volume effects due to molecular crowding would lead to correlations between molecules even without specific interactions. What is the interplay of these effects in the crowded cellular environment? We study dimerization of a model homodimer when the mondimers are free and when they are tethered to each other. We consider a structured environment: Two monomers first diffuse into a cavity of size L and then fold and bind within the cavity. The folding and binding are simulated by using molecular dynamics based on a simplified topology based model. The confinement in the cell is described by an effective molecular concentration C ∼ L −3 . A two-state coupled folding and binding behavior is found. We show the maximal rate of dimerization occurred at an effective molecular concentration C op ≃ 1 mM, which is a relevant cellular concentration. In contrast, for tethered chains the rate keeps at a plateau when C < C op but then decreases sharply when C > C op . For both the free and tethered cases, the simulated variation of the rate of dimerization and thermodynamic stability with effective molecular concentration agrees well with experimental observations. In addition, a theoretical argument for the effects of confinement on dimerization is also made.

Energy landscape theory requires that the protein-folding mechanism is generally globally directed or funneled toward the native state. The collective nature of transition state ensembles further suggests that sufficient averaging of the native interactions can occur so that the knowledge of the native topology may suffice for predicting the mechanism. Nevertheless, while simple homogeneously weighted native topology-based models predict the folding mechanisms for many proteins, for other proteins knowledge of the native topology, by itself, seems not to suffice in determining the folding mechanism. Simulations of proteins with differing topologies reveal that the failure of homogeneously weighted topology-based models can, however, be completely understood within the framework of a funneled energy landscape and can be quantified by comparing the fluctuation of entropy cost for forming contacts to the expected fluctuations in contact energy. To be precise, we find the transition state ensembles of proteins with all-α topologies, which are more uniform in the specific entropy cost of contact formation, have transition state ensembles that are more readily perturbed by differences in energetic weights than are the transition state ensembles of proteins with significant amounts of β-structure, where the specific entropy costs of contact formation are more widely distributed. This behavior is consistent with a random-field Ising model analogy that follows from the free energy functional approach to folding.

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 protein's 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.

The same energy landscape principles associated with the folding of proteins into their monomeric conformations should also describe how these proteins oligomerize into domain-swapped conformations. We tested this hypothesis by using a simplified model for the epidermal growth factor receptor pathway substrate 8 src homology 3 domain protein, both of whose monomeric and domain-swapped structures have been solved. The model, which we call the symmetrized Gō-type model, incorporates only information regarding the monomeric conformation in an energy function for the dimer to predict the domain-swapped conformation. A striking preference for the correct domain-swapped structure was observed, indicating that overall monomer topology is a main determinant of the structure of domain-swapped dimers. Furthermore, we explore the free energy surface for domain swapping by using our model to characterize the mechanism of oligomerization.

Protein recognition and binding, which result in either transient or long-lived complexes, play a fundamental role in many biological functions, but sometimes also result in pathologic aggregates. We use a simplified simulation model to survey a range of systems where two highly flexible protein chains form a homodimer. In all cases, this model, which corresponds to a perfectly funneled energy landscape for folding and binding, reproduces the macroscopic experimental observations on whether folding and binding are coupled in one step or whether intermediates occur. Owing to the minimal frustration principle, we find that, as in the case of protein folding, the native topology is the major factor that governs the choice of binding mechanism. Even when the monomer is stable on its own, binding sometimes occurs fastest through unfolded intermediates, thus showing the speedup envisioned in the fly-casting scenario for molecular recognition.

We address the stability of multicharged finite systems driven by Coulomb forces beyond the Rayleigh instability limit. Our exploration of the nuclear dynamics of heavily charged Morse clusters enabled us to vary the range of the pair potential and of the fissibility parameter, which results in distinct fragmentation patterns and in the angular distributions of the fragments. The Rayleigh instability limit separates between nearly binary (or tertiary) spatially unisotropic fission and spatially isotropic Coulomb explosion into a large number of small, ionic fragments. Implications are addressed for a broad spectrum of dynamics in chemical physics, radiation physics of ultracold gases, and biophysics, involving the fission of clusters and droplets, the realization of Coulomb explosion of molecular clusters, the isotropic expansion of optical molasses, and the Coulomb instability of “isolated” proteins.

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala 12 ) 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 Ala 12 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 Ala 12 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.

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.

It is widely appreciated that water molecules play an invaluable role in governing the structure, stability, dynamics, and function of biomolecules. The hydration forces are responsible for packing and stabilization of the protein structure. Particularly, water participates in many hydrogen bond networks and screening electrostatic interactions. However, the exact range of processes mediated by water is far from being understood, and it is only in the recent years that water has been quantitatively treated as an integral component of biomolecular systems. In this issue of PNAS, Papoian et al. (1) report a significant improvement in protein structure prediction by adding a water knowledge-based potential to an established Hamiltonian for protein structure prediction. “Wetting” the Hamiltonian improves the predicted structures, especially for large proteins, when long-range interactions between polar or charged groups are mediated by water molecules. There are a variety of experimental and theoretical studies acknowledging the active role of solvent in protein stability and dynamics. Experimentally, x-ray, neutron diffraction (2), NMR (3, 4), and femtosecond fluorescence (5) measurements reveal the binding sites, structure, and dynamics of water. Theoretically, molecular dynamics simulations offer a detailed atomic description of both the biomolecule and the solvent as well as the time dependency of their dynamics (6–9). These all-atom explicit–solvent simulations of proteins, however, do not appear to provide sufficient conformational coverage to tackle many equilibrium and long-time-scale kinetic properties. Thus, a complementary approach is to adopt simplified models that trade high structural resolution of the water molecule (10, 11) (e.g., generalized born model for the water) or of both the polypeptide and the solvent (e.g., structure-based models) for enhancing conformational sampling (12, 13). Desolvation during folding processes was studied by using a structure-based (Go) model for the Src homology 3 (SH3) …