Search

Intrinsically disordered proteins play a crucial role in cellular phase separation, yet the diverse molecular forces driving phase separation are not fully understood. It is of utmost importance to understand how peptide sequence, and particularly the balance between the peptides' short- and long-range interactions with other peptides, may affect the stability, structure, and dynamics of liquid-liquid phase separation in protein condensates. Here, using coarse-grained molecular dynamics simulations, we studied the liquid properties of the condensate in a series of polymers in which the ratio of short-range dispersion interactions to long-range electrostatic interactions varied. As the fraction of mutations that participate in short-range interactions increases at the expense of long-range electrostatic interactions, a significant decrease in the critical temperature of phase separation is observed. Nevertheless, sequences with a high fraction of short-range interactions exhibit stabilization, which suggests compensation for the loss of long-range electrostatic interactions. Decreased condensate stability is coupled with decreased translational diffusion of the polymers in the condensate, which may result in the loss of liquid characteristics in the presence of a high fraction of uncharged residues. The effect of exchanging long-range electrostatic interactions for short-range interactions can be explained by the kinetics of breaking intermolecular contacts with neighboring polymers and the kinetics of intramolecular fluctuations. While both time scales are coupled and increase as electrostatic interactions are lost, for sequences that are dominated by short-range interactions, the kinetics of intermolecular contact breakage significantly slows down. Our study supports the contention that different types of interactions can maintain protein condensates, however, long-range electrostatic interactions enhance its liquid-like behavior.

DNA mismatch repair (MMR) is an important postreplication process that eliminates mispaired or unpaired nucleotides to ensure genomic replication fidelity. In humans, Msh2-Msh6 and Msh2-Msh3 are the two mismatch repair initiation factors that recognize DNA lesions. While X-ray crystal structures exist for these proteins in complex with DNA lesions, little is known about their structures during the initial search along nonspecific double-stranded DNA, because they are short-lived and difficult to determine experimentally. In this study, various computational approaches were used to sidestep these difficulties. All-atom and coarse-grained simulations based on the crystal structures of Msh2-Msh3 and Msh2-Msh6 showed no translation along the DNA, suggesting that the initial search conformation differs from the lesion-bound crystal structure. We modeled probable search-mode structures of MSH proteins and showed, using coarse-grained molecular dynamics simulations, that they can perform rotation-coupled diffusion on DNA, which is a suitable and efficient search mechanism for their function and one predicted earlier by fluorescence resonance energy transfer and fluorescence microscopy studies. This search mechanism is implemented by electrostatic interactions among the mismatch-binding domain (MBD), the clamp domains, and the DNA backbone. During simulations, their diffusion rate did not change significantly with an increasing salt concentration, which is consistent with observations from experimental studies. When the gap between their DNA-binding clamps was increased, Msh2-Msh3 diffused mostly via the clamp domains while Msh2-Msh6 still diffused using the MBD, reproducing the experimentally measured lower diffusion coefficient of Msh2-Msh6. Interestingly, Msh2-Msh3 was capable of dissociating from the DNA, whereas Msh2-Msh6 always diffused on the DNA duplex. This is consistent with the experimental observation that Msh2-Msh3, unlike Msh2-Msh6, can overcome obstacles such as nucleosomes. Our models provide a molecular picture of the different mismatch search mechanisms undertaken by Msh2-Msh6 and Msh2-Msh3, despite the similarity of their structures.

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other's substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.

Despite the continuous discovery of host and guest proteins in membraneless organelles, complex host–guest interactions hinder the understanding of the molecular grammar governing liquid–liquid phase separation. In this study, we characterized the localization and dynamic properties of guest proteins in liquid droplets using single-molecule fluorescence microscopy. Eighteen guest proteins of different sizes, structures, and oligomeric states were examined in host p53 liquid droplets. Recruitment did not significantly depend on the structural properties of the guest proteins, but was moderately correlated with their length, total charge, and number of R and Y residues. In contrast, the diffusion of disordered guest proteins was comparable to that of host p53, whereas that of folded proteins varied widely. Molecular dynamics simulations suggest that folded proteins diffuse within the voids of the liquid droplet while interacting weakly with neighboring host proteins, whereas disordered proteins adapt their structures to form tight interactions with the host proteins. Our study provides insights into the key molecular principles of the localization and dynamics of guest proteins in liquid droplets.

Replication protein A (RPA) plays a critical role in all eukaryotic DNA processing involving single-stranded DNA (ssDNA). Contrary to the notion that RPA provides solely inert protection to transiently formed ssDNA, the RPA–ssDNA complex acts as a dynamic DNA processing unit. Here, we studied the diffusion of RPA along 60 nt ssDNA using a coarse-grained model in which the ssDNA–RPA interface was modeled by both aromatic and electrostatic interactions. Our study provides direct evidence of bulge formation during the diffusion of ssDNA along RPA. Bulges can form at a few sites along the interface and store 1–7 nt of ssDNA whose release, upon bulge dissolution, leads to propagation of ssDNA diffusion. These findings thus support the reptation mechanism, which involves bulge formation linked to the aromatic interactions, whose short range nature reduces cooperativity in ssDNA diffusion. Greater cooperativity and a larger diffusion coefficient for ssDNA diffusion along RPA are observed for RPA variants with weaker aromatic interactions and for interfaces homogenously stabilized by electrostatic interactions. ssDNA propagation in the latter instance is characterized by lower probabilities of bulge formation; thus, it may fit the sliding-without-bulge model better than the reptation model. Thus, the reptation mechanism allows ssDNA mobility despite the extensive and high affinity interface of RPA with ssDNA. The short-range aromatic interactions support bulge formation while the long-range electrostatic interactions support the release of the stored excess ssDNA in the bulge and thus the overall diffusion.

Significance Six decades after DNA structure was first revealed, fundamental questions remain open. How is the entwined embrace of double-stranded nucleic acids formed or disrupted? How does the energetics underlying this process influence nucleic-acid processing machineries? By combining simulations and experiments, our work addresses these questions and reveals that asymmetric base-pair dynamics drives the stepwise separation of nucleic acid duplexes, predicts the unwinding efficiency of helicases, and intimately relates the intrinsic dynamics of base pairs to the enzymatic mechanism evolved for their opening. Taken together, our data suggest a layer of regulation of the genetic material encoded in the “unwindability” of the double helix., The opening of a Watson–Crick double helix is required for crucial cellular processes, including replication, repair, and transcription. It has long been assumed that RNA or DNA base pairs are broken by the concerted symmetric movement of complementary nucleobases. By analyzing thousands of base-pair opening and closing events from molecular simulations, here, we uncover a systematic stepwise process driven by the asymmetric flipping-out probability of paired nucleobases. We demonstrate experimentally that such asymmetry strongly biases the unwinding efficiency of DNA helicases toward substrates that bear highly dynamic nucleobases, such as pyrimidines, on the displaced strand. Duplex substrates with identical thermodynamic stability are thus shown to be more easily unwound from one side than the other, in a quantifiable and predictable manner. Our results indicate a possible layer of gene regulation coded in the direction-dependent unwindability of the double helix.

Lassa virus (LASV) is a notorious human pathogen in West Africa. Its class I trimeric spike complex displays a distinct architecture, and its cell entry mechanism involves unique attributes not shared by other related viruses. We determined the crystal structure of the GP2 fusion glycoprotein from the spike complex of LASV (GP2LASV) in its post-fusion conformation. GP2LASV adopts a canonical helical bundle configuration similarly to other viruses in its family. The core packing of GP2LASV, however, is more organized compared to GP2 from other viruses reducing the formation of internal hydrophobic cavities. We demonstrate a link between the formation of such unfavorable hydrophobic cavities and the efficiencies of membrane fusion and cell entry. Our study suggests that LASV has evolved a more efficient membrane fusogen compared to other viruses from its family by optimizing the post-fusion configuration of its GP2 module.

Highly negatively charged segments containing only aspartate or glutamate residues (“D/E repeats”) are found in many eukaryotic proteins. For example, the C-terminal 30 residues of the HMGB1 protein are entirely D/E repeats. Using nuclear magnetic resonance (NMR), fluorescence, and computational approaches, we investigated how the D/E repeats causes the autoinhibition of HMGB1 against its specific binding to cisplatin-modified DNA. By varying ionic strength in a wide range (40–900 mM), we were able to shift the conformational equilibrium between the autoinhibited and uninhibited states toward either of them to the full extent. This allowed us to determine the macroscopic and microscopic equilibrium constants for the HMGB1 autoinhibition at various ionic strengths. At a macroscopic level, a model involving the autoinhibited and uninhibited states can explain the salt concentration-dependent binding affinity data. Our data at a microscopic level show that the D/E repeats and other parts of HMGB1 undergo electrostatic fuzzy interactions, each of which is weaker than expected from the macroscopic autoinhibitory effect. This discrepancy suggests that the multivalent nature of the fuzzy interactions enables strong autoinhibition at a macroscopic level despite the relatively weak intramolecular interaction at each site. Both experimental and computational data suggest that the D/E repeats interact preferentially with other intrinsically disordered regions (IDRs) of HMGB1. We also found that mutations mimicking post-translational modifications relevant to nuclear export of HMGB1 can moderately modulate DNA-binding affinity, possibly by impacting the autoinhibition. This study illuminates a functional role of the fuzzy interactions of D/E repeats.

DNA-binding proteins rely on linear diffusion along the longitudinal DNA axis, supported by their nonspecific electrostatic affinity for DNA, to search for their target recognition sites. One may therefore expect that the ability to engage in linear diffusion along DNA is universal to all DNA-binding proteins, with the detailed biophysical characteristics of that diffusion differing between proteins depending on their structures and functions. One key question is whether the linear diffusion mechanism is defined by translation coupled with rotation, a mechanism that is often termed sliding. We conduct coarse-grained and atomistic molecular dynamics simulations to investigate the minimal requirements for protein sliding along DNA. We show that coupling, while widespread, is not universal. DNA-binding proteins that slide along DNA transition to uncoupled translation-rotation (i.e., hopping) at higher salt concentrations. Furthermore, and consistently with experimental reports, we find that the sliding mechanism is the less dominant mechanism for some DNA-binding proteins, even at low salt concentrations. In particular, the toroidal PCNA protein is shown to follow the hopping rather than the sliding mechanism.

Microtubules (MTs) are essential components of the eukaryotic cytoskeleton that serve as “highways” for intracellular trafficking. In addition to the well-known active transport of cargo by motor proteins, many MT-binding proteins seem to adopt diffusional motility as a transportation mechanism. However, because of the limited spatial resolution of current experimental techniques, the detailed mechanism of protein diffusion has not been elucidated. In particular, the precise role of tubulin tails and tail modifications in the diffusion process is unclear. Here, using coarse-grained molecular dynamics simulations validated against atomistic simulations, we explore the molecular mechanism of protein diffusion along MTs. We found that electrostatic interactions play a central role in protein diffusion; the disordered tubulin tails enhance affinity but slow down diffusion, and diffusion occurs in discrete steps. While diffusion along wild-type MT is performed in steps of dimeric tubulin, the removal of the tails results in a step of monomeric tubulin. We found that the energy barrier for diffusion is larger when diffusion on MTs is mediated primarily by the MT tails rather than the MT body. In addition, globular proteins (EB1 and PRC1) diffuse more slowly than an intrinsically disordered protein (Tau) on MTs. Finally, we found that polyglutamylation and polyglycylation of tubulin tails lead to slower protein diffusion along MTs, although polyglycylation leads to faster diffusion across MT protofilaments. Taken together, our results explain experimentally observed data and shed light on the roles played by disordered tubulin tails and tail modifications in the molecular mechanism of protein diffusion along MTs.

Huntington's disease is a neurodegenerative disorder resulting from an expanded polyglutamine (polyQ) repeat of the Huntingtin (Htt) protein. Affected tissues often contain aggregates of the N-terminal Htt exon 1 (Htt-Ex1) fragment. The N-terminal N17 domain proximal to the polyQ tract is key to enhance aggregation and modulate Htt toxicity. Htt-Ex1 is intrinsically disordered, yet it has been postulated that under physiological conditions membranes induce the N17 to adopt an α-helical structure, which then plays a key role in regulating Htt protein aggregation. The present study leverages the recently available assignment of NMR peaks in an N17Q17 construct, in order to provide a look into the changes occurring in vitro upon exposing this fragment to various brain extract fragments as well as to synthetic bilayers. Residue-specific changes were observed by 3D HNCO NMR, whose nature was further clarified with ancillary CD and aggregation studies, as well as with molecular dynamic calculations. From this combination of measurements and computations, a unified picture emerges, whereby transient structures consisting of α-helices spanning a fraction of the N17 residues form during N17Q17-membrane interactions. These interactions are fairly dynamic, but they qualitatively mimic more rigid variants that have been discussed in the literature. The nature of these interactions and their potential influence on the aggregation process of these kinds of constructs under physiological conditions are briefly assessed.

Several DNA-binding proteins, such as topoisomerases, helicases and sliding clamps, have a toroidal (i.e. ring) shape that topologically traps DNA, with this quality being essential to their function. Many DNA-binding proteins that function, for example, as transcription factors or enzymes were shown to be able to diffuse linearly (i.e. slide) along DNA during the search for their target binding sites. The protein's sliding properties and ability to search DNA, which often also involves hopping and dissociation, are expected to be different when it encircles the DNA. In this study, we explored the linear diffusion of four ring-shaped proteins of very similar structure: three sliding clamps (PCNA, β-clamp, and the gp45) and the 9-1-1 protein, with a particular focus on PCNA. Coarse-grained molecular dynamics simulations were performed to decipher the sliding mechanism adopted by these ring-shaped proteins and to determine how the molecular properties of the inner and outer ring govern its search speed. We designed in silico variants to dissect the contributions of ring geometry and electrostatics to the sliding speed of ring-shaped proteins along DNA. We found that the toroidal proteins diffuse when they are tilted relative to the DNA axis and able to rotate during translocation, but that coupling between rotation and translocation is quite weak. Their diffusion speed is affected by the shape of the inner ring and, to a lesser extent, by its electrostatic properties. However, breaking the symmetry of the electrostatic potential can result in deviation of the DNA from the center of the ring and cause slower linear diffusion. The findings are discussed in light of earlier computational and experimental studies on the sliding of clamps.

Proteins are subject to various conflicting forces that trade-off against each other. For example, during folding, the protein achieves lower enthalpy at the cost of lower entropy. Similarly, the trade-off for increased stability may be decreased flexibility, which may abolish allosteric pathways. Accordingly, stability trades-off against function, which may also trade-off against folding kinetics and mechanism. Furthermore, attaining increased stability may reduce a protein's ability to adopt novel functions. Understanding the biophysics and function of proteins requires quantification of the driving forces involved in each of the trade-offs. Indeed, quantification of the linkages in the network of trade-offs is essential to obtaining a more complete understanding of protein structure and function.

Microtubules (MTs), an essential component of the eukaryotic cytoskeleton, are a lattice of polymerized tubulin dimers and are crucial for various cellular processes. The genetic and chemical diversity of tubulin and their disordered tails gives rise to a “tubulin code”. The functional role of tubulin post-translational modifications (PTMs), which contribute to the chemical diversity of the tubulin code, is gradually being unraveled. However, variation in the length and spatial organization of tubulin poly-modifications leads to an enormous combinatorial PTM space, which is difficult to study experimentally. Hence, the impact of the combinatorial tubulin PTM space on the biophysical properties of tubulin tails and their interactions with other proteins remains elusive. Here, we combine all-atom and coarse-grained molecular dynamics simulations to elucidate the biophysical implications of the large combinatorial tubulin PTM space in the context of an MT lattice. We find that tail–body interactions are more dominant in the tubulin dimer than in an MT lattice, and are more significant for the tails of α compared with β tubulin. In addition, polyglutamylation, but not polyglycylation, expands the dimensions of the tubulin tails. Polyglutamylation also leads to a decrease in the diffusion rate of MT-associated protein EB1 on MTs, while polyglycylation often increases diffusion rate. These observations are generally not sensitive to the organization of the polymodifications. The effect of PTMs on MT charge density and tail dynamics are also discussed. Overall, this study presents a molecular quantification of the biophysical properties of tubulin tails and their polymodifications, and provides predictions on the functional importance of tubulin PTMs.

Intersegmental transfer that involves direct relocation of a DNA-binding protein from one nonspecific DNA site to another was previously shown to contribute to speeding up the identification of the DNA target site. This mechanism is promoted when the protein is composed of at least two domains that have different DNA binding affinities and thus show a degree of mobility. In this study, we investigate the effect of particle crowding on the ability of a multi-domain protein to perform intersegmental transfer. We show that although crowding conditions often favor 1D diffusion of proteins along DNA over 3D diffusion, relocation of one of the tethered domains to initiate intersegmental transfer is possible even under crowding conditions. The tendency to perform intersegmental transfer by a multi-domain protein under crowding conditions is much higher for larger crowding particles than smaller ones and can be even greater than under no-crowding conditions. We report that the asymmetry of the two domains is even magnified by the crowders. The observations that crowding supports intersegmental transfer serve as another example that in vivo complexity does not necessarily slow down DNA search kinetics by proteins.

In prokaryotes, the RecA protein catalyzes the repair and strand exchange of double-stranded DNA. RecA binds to single-stranded DNA (ssDNA) and forms a presynaptic complex in which the protein polymerizes around the ssDNA to form a right-handed helical nucleoprotein filament structure. In the present work, the mechanism for the formation of the RecA-ssDNA filament structure is modeled using coarse-grained molecular dynamics simulations. Information from the X-ray structure was used to model the protein itself but not its interactions; the interactions between the protein and the ssDNA were modeled solely by electrostatic, aromatic, and repulsive energies. For the present study, the monomeric, dimeric, and trimeric units of RecA and 4, 8, and 11 NT-long ssDNA, respectively, were studied. Our results indicate that monomeric RecA is not sufficient for nucleoprotein filament formation; rather, dimeric RecA is the elementary binding unit, with higher multimeric units of RecA facilitating filament formation. Our results reveal that loop region flexibility at the primary binding site of RecA is essential for it to bind the incoming ssDNA, that the aromatic residues present in the loop region play an important role in ssDNA binding, and that ATP may play a role in guiding the ssDNA by changing the electrostatic potential of the RecA protein.

The α-helical coiled coil (CC) is a common protein motif that because of the simplicity of its sequence/structure relationship, it has been studied extensively to address fundamental questions in protein science as well as to develop strategies for designing protein with novel architectures. Nevertheless, a complete understanding of CC structures and their dynamics is still far from achieved. Particularly, spontaneous sliding at interfaces of CC proteins was observed for some systems, but its mechanism and usage as an intrinsic conformational change at CCs in protein-protein interfaces is unclear. Using coarse-grained and atomistic simulations, we study various sequences of homodimeric CC, in both parallel and antiparallel configurations. Both the strength of the hydrophobic core and the existence of salt bridges at the periphery of the interface affect sliding dynamics at the CC interface. Although the energy landscape for sliding along a CC interface is different for parallel and antiparallel configurations, both are characterized by a free energy of 1–1.5 kcal/mol, depending on the residues that constitute the CC interface. These barrier heights suggest that sliding kinetics is relatively slow in CC systems and are not expected to be of long length scale, yet they can be involved in functional motions. Our study explains the sliding that has been experimentally observed for the antiparallel CC of the dynein stalk region and the nuclear pore complex and suggests that this one-dimensional motion is an intrinsic feature in CC systems that can be involved in other CC systems.

The Engrailed Homeodomain (EnHD) transcription factor of Drosophila melanogaster was fused to the enhanced green fluorescent protein (eGFP) either at its C- or N-terminus via three- or ten-residue flexible linkers. Here, we show that EnHD undergoes destabilization upon fusing it to eGFP regardless of the linker length used and whether the tethering is to its N- or C-terminus. The destabilization is reflected in melting points that are lower by up to 9°C. Thermodynamic analysis and coarse-grained molecular dynamic simulations indicate that this destabilization is due to eGFP-promoted entropic stabilization of the denatured state ensemble of EnHD. Our results provide, therefore, an example for destabilizing interdomain allostery. They are also important given the widespread use of eGFP tagging in cell biology, as they indicate that such tagging can cause unintended protein destabilization and concomitant effects.

Glycosylation plays not only a functional role but can also modify the biophysical properties of the modified protein. Usually, natural glycosylation results in protein stabilization; however, in vitro and in silico studies showed that sometimes glycosylation results in thermodynamic destabilization. Here, we applied coarse-grained and all-atom molecular dynamics simulations to understand the mechanism underlying the loss of stability of the MM1 protein by glycosylation. We show that the origin of the destabilization is a conformational distortion of the protein caused by the interaction of the monosaccharide with the protein surface. Though glycosylation creates new short-range glycan-protein interactions that stabilize the conjugated protein, it breaks long-range protein-protein interactions. This has a destabilizing effect because the probability of long- and short-range interactions forming differs between the folded and unfolded states. The destabilization originates not from simple loss of interactions but due to a trade-off between the short- and long-range interactions.

Ubiquitination is one of the most common post-translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site-specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48-linked tetra-ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half-life through proteasomal degradation.

The recognition of DNA-binding proteins (DBPs) to their specific site often precedes by a search technique in which proteins slide, hop along the DNA contour or perform inter-segment transfer and 3D diffusion to dissociate and re-associate to distant DNA sites. In this study, we demonstrated that the strength and nature of the non-specific electrostatic interactions, which govern the search dynamics of DBPs, are strongly correlated with the conformation of the DNA. We tuned two structural parameters, namely curvature and the extent of helical twisting in circular DNA. These two factors are mutually independent of each other and can modulate the electrostatic potential through changing the geometry of the circular DNA conformation. The search dynamics for DBPs on circular DNA is therefore markedly different compared with linear B-DNA. Our results suggest that, for a given DBP, the rotation-coupled sliding dynamics is precluded in highly curved DNA (as well as for over-twisted DNA) because of the large electrostatic energy barrier between the inside and outside of the DNA molecule. Under such circumstances, proteins prefer to hop in order to explore interior DNA sites. The change in the balance between sliding and hopping propensities as a function of DNA curvature or twisting may result in different search efficiency and speed.

DNA recognition by DNA-binding proteins, which is a pivotal event in most gene regulatory processes, is often preceded by an extensive search for the correct site. A facilitated diffusion process, in which a DBP combines 3D diffusion in solution with 1D sliding along DNA, has been suggested to explain how proteins can locate their target sites on DNA much faster than predicted by 3D diffusion alone. One of the key mechanisms in the localization of the target by a DNA-binding protein is intersegment transfer in which the protein forms a bridged intermediate between two distant DNA regions. This jumping mechanism is more enhanced when the DNA-binding protein is asymmetric in its structure or its dynamics. We suggest that asymmetry supports the “monkey bar” mechanism, in which different domains of the protein interact with different DNA regions. In this minireview, we discuss how the molecular architectures of the proteins and DNA may modulate the efficiency of monkey bar dynamics.

Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA-ssDBP and ssRNA-ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base-aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA-protein and ssRNA-protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA-ssRBP complexes compared with the ssDNA-ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (koff) for complexes involving the former.

We focus on dimeric DNA-binding proteins from two well-studied families: orthodox type II restriction endonucleases (REs) and transcription factors (TFs). Interactions of the protein's recognition sites with the DNA and, particularly, the contribution of each of the monomers to one-dimensional (1D) sliding along nonspecific DNA were studied using computational tools. Coarse-grained molecular dynamics simulations of DNA scanning by various TFs and REs provide insights into how the symmetry of a homodimer can be broken while they nonspecifically interact with DNA. The characteristics of protein sliding along DNA, such as the average sliding length, partitioning between 1D and 3D search, and the one-dimensional diffusion coefficient D1, strongly depend on the salt concentration, which in turn affects the probability of the two monomers adopting a cooperative symmetric sliding mechanism. Indeed, we demonstrate that maximal DNA search efficiency is achieved when the protein adopts an asymmetric search mode in which one monomer slides while its partner hops. We find that proteins classified as TFs have a higher affinity for the DNA, longer sliding lengths, and an increased probability of symmetric sliding in comparison with REs. Moreover, TFs can perform their biological function over a much wider range of salt concentrations than REs. Our results demonstrate that the different biological functions of DNA-binding proteins are related to the different nonspecific DNA search mechanisms they adopt.

DNA recognition by DNA-binding proteins (DBPs), which is a pivotal event in most gene regulatory processes, is often preceded by an extensive search for the correct site. A facilitated diffusion process in which a DBP combines three-dimensional diffusion in solution with one-dimensional sliding along DNA has been suggested to explain how proteins can locate their target sites on DNA much faster than predicted by three-dimensional diffusion alone. Although experimental and theoretical studies have recently advanced understanding of the biophysical principles underlying the search mechanism, the process under in vivo cellular conditions is poorly understood. In this study, we used various computational approaches to explore how the presence of obstacle proteins on the DNA influences search efficiency. At a low obstacle occupancy (i.e., when few obstacles occupy sites on the DNA), sliding by the searching DBP may be confined, which may impair search efficiency. The obstacles, however, can be bypassed during hopping events, and the number of bypasses is larger for higher obstacle occupancies. Dynamism on the part of the obstacles may even further facilitate search kinetics. Our study shows that the nature and efficiency of the search process may be governed not only by the intrinsic properties of the DBP and the salt concentration of the medium, but also by the in vivo association of DNA with other macromolecular obstacles, their location, and occupancy.

Extensive contact with DNA via multiple zinc fingers allows highly specific DNA-binding of zinc-finger-class transcription factors, but can also slow the target search process. Here we introduce recent insights into how zinc-finger proteins can rapidly scan DNA. Potential application of the new knowledge to the zinc-finger-based technology is also discussed.

Differential signaling of the type I interferon receptor (IFNAR) has been correlated with the ability of its subunit, IFNAR1, to differentially recognize a large spectrum of different ligands, which involves intricate conformational re-arrangements of multiple interacting domains. To shed light onto the structural determinants governing ligand recognition, we compared the force-induced unfolding of the IFNAR1 ectodomain when bound to interferon and when free, using the atomic force microscope and steered molecular dynamics simulations. Unexpectedly, we find that IFNAR1 is easier to mechanically unfold when bound to interferon than when free. Analysis of the structures indicated that the origin of the reduction in unfolding forces is a conformational change in IFNAR1 induced by ligand binding.

Repeat proteins have unique elongated structures that, unlike globular proteins, are quite modular. Despite their simple one-dimensional structure, repeat proteins exhibit intricate folding behavior with a complexity similar to that of globular proteins. Therefore, repeat proteins allow one to quantify fundamental aspects of the biophysics of protein folding. One important feature of repeat proteins is the interfaces between the repeating units. In particular, the distribution of stabilities within and between the repeats was previously suggested to affect their folding characteristics. In this study, we explore how the interface affects folding kinetics and cooperativity by investigating two families of repeat proteins, namely, the Ankyrin and tetratricopeptide repeat proteins, which differ in the number of interfacial contacts that are formed between their units as well as in their folding behavior. By using simple topology-based models, we show that modulating the energetic strength of the interface relative to that of the repeat itself can drastically change the protein stability, folding rate, and cooperativity. By further dissecting the interfacial contacts into several subsets, we isolated the effects of each of these groups on folding kinetics. Our study highlights the importance of interface connectivity in determining the folding behavior.

Intrinsically disordered regions, terminal tails, and flexible linkers are abundant in DNA-binding proteins and play a crucial role by increasing the affinity and specificity of DNA binding. Disordered tails often undergo a disorder-to-order transition during interactions with DNA and improve both the kinetics and thermodynamics of specific DNA binding. The DNA search by proteins that interact nonspecifically with DNA can be supported by disordered tails as well. The disordered tail may increase the overall protein-DNA interface and thus increase the affinity of the protein to the DNA and its sliding propensity while slowing linear diffusion. The exact effect of the disordered tails on the sliding rate depends on the degree of positive charge clustering, as has been shown for homeodomains and p53 transcription factors. The disordered tails, which may be viewed as DNA recognizing subdomains, can facilitate intersegment transfer events that occur via a "monkey bar" mechanism in which the domains bridge two different DNA fragments simultaneously. The "monkey bar" mechanism can be facilitated by internal disordered linkers in multidomain proteins that mediate the cross-talks between the constituent domains and especially their brachiation dynamics and thus their overall capability to search DNA efficiently. The residue sequence of the disordered tails has unique characteristics that were evolutionarily selected to achieve the optimized function that is unique to each protein. Perturbation of the electrostatic characteristics of the disordered tails by post-translational modifications, such as acetylation and phosphorylation, may affect protein affinity to DNA and therefore can serve to regulate DNA recognition. Modifying the disordered protein tails or the flexibility of the inter-domain linkers of multidomain proteins may affect the cross-talk between the constituent domains so as to facilitate the search kinetics of non-specific DNA sequences and increase affinity to the specific sequences.

The p53 protein is a homotetrameric transcription factor whose monomers comprise several domains. Although its organization with and without DNA was elucidated recently, characterizing the p53-DNA complex at the atomic level remains challenging because of its many disordered regions. Here we use computational models to predict the wiring of the four chains composing p53 and study its sliding dynamics along DNA in different oligomeric states. We find that helical sliding along the major groove is the most feasible DNA search mechanism for a large range of salt concentrations. Tighter packing of the tetrameric core domain is associated with a greater nonspecific affinity for DNA and the slowest linear diffusion dynamics along DNA. C-tails facilitate linear diffusion but restrict the association of two primary dimers into a tetramer. This restriction can disappear at higher salt concentrations, which decrease the affinity of C-tails for DNA, or upon interaction of the C-tail with other DNA segments. Our results support evidence for the positive regulation of p53 function by the C-tails and suggest that posttranslational charge modifications may alter the affinity of the tails for DNA. Conversely, the N-termini have little effect on sliding characteristics. Changes in the electrostatic potentials of the core domain via missense mutations corresponding to cancer development can also affect sliding by p53. Our study provides molecular insight into the role of various p53 domains during DNA search and indicates that the complex interdomain and protein-DNA cross-talks in which p53 engages may be related to its repertoire of cellular functions.

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.

More than 70% of eukaryotic proteins are composed of multiple domains. However, most studies of the search for DNA focus on individual protein domains and do not consider potential cross talk within a multidomain transcription factor. In this study, the molecular features of the DNA search mechanism were explored for two multidomain transcription factors: human Pax6 and Oct-1. Using a simple computational model, we compared a DNA search of multidomain proteins with a search of isolated domains. Furthermore, we studied how manipulating the binding affinity of a single domain to DNA can affect the overall DNA search of the multidomain protein. Tethering the two domains via a flexible linker increases their affinity to the DNA, resulting in a higher propensity for sliding along the DNA, which is more significant for the domain with the weaker DNA-binding affinity. In this case, the domain that binds DNA more tightly anchors the multidomain protein to the DNA and, via the linker, increases the local concentration of the weak DNA-binding domain (DBD). The tethered domains directly exchange between two parallel DNA molecules via a bridged intermediate, where intersegmental transfer is promoted by the weaker DBD. We found that, in general, the relative affinity of the two domains can significantly affect the cross talk between them and thus their overall capability to search DNA efficiently. The results we obtained by examining various multidomain DNA-binding proteins support the necessity of discrepancies between the DNA-binding affinities of the constituent domains.

The search through nonspecific DNA for a specific site by proteins is known to be facilitated by sliding, hopping, and intersegment transfer between separate DNA strands, yet the driving forces of these protein dynamics from the molecular perspective are unclear. In this study, molecular features of the DNA search mechanism were explored for three homologous proteins (the HoxD9, Antp, and NK-2 homeodomains) using a simple computational model in which protein-DNA interactions are represented solely by electrostatic forces. In particular, we studied the impact that disordered N-terminal tails (N-tails), which are more common in DNA-binding proteins than in other proteins, have on the efficiency of DNA search. While the three homeodomain proteins were found to use similar binding interfaces in specific and nonspecific interactions with DNAs, their different electrostatic potentials affect the nature of their sliding dynamics. The different lengths and net charges of the N-tails of the homeodomains affect their motion along the DNA. The presence of an N-tail increases sliding propensity but slows linear diffusion along the DNA. When the search is performed in the presence of two parallel DNA molecules, a direct transfer, which is facilitated by the protein tail, from one nonspecific DNA to another occurs. The tailed proteins jump between two DNA molecules through an intermediate in which the recognition helix of the protein is adsorbed to one DNA fragment and the N-tail is adsorbed to the second, suggesting a "monkey bar" mechanism. Our study illustrates how the molecular architecture of proteins controls the efficiency of DNA scanning.

Multimeric proteins are ubiquitous in many cellular processes that require high levels of regulation. Eukaryotic gene expression is often regulated by a mechanism of combinatorial control that involves the binding of dimeric transcription factors to DNA together with the coordinated activity of additional proteins. In this study, we investigated the dimerization of the Arc-repressor on DNA with the aim of achieving microscopic insight into the possible advantages of interacting with DNA as a complex rather than as a monomeric single-domain protein. We used a computational coarse-grained model in which the protein dynamics was governed by native interactions and protein-DNA interactions were dictated by electrostatic forces. Inspired by previous experimental work that showed an enhanced refolding rate for the Arc-repressor in the presence of DNA and other polyanions, we focused on the mechanism and kinetics of the assembly of Arc monomers in the presence of single- (ssDNA) and double-stranded DNA (dsDNA) molecules in a low-salt concentration environment. The electrostatic interactions that attract the protein to the dsDNA were shown to be fundamental in colocalizing the unfolded Arc chains and in accelerating refolding. Arc monomers bind the dsDNA efficiently and nonspecifically, and search for each other via one-dimensional diffusion. The fastest folding of Arc is observed for DNA of 30 bp. Longer DNA is significantly less efficient in accelerating the Arc refolding rate, since the two subunits search distinct regions of the one-dimensional DNA and are therefore much less colocalized. The probability that the two unfolded chains will meet on 200 bp DNA is similar to that in the bulk. The colocalization of Arc subunits on ssDNA results in much faster folding compared to that obtained on dsDNA of the same length. Differences in the rate of Arc refolding, cooperativity, and the structure of its transition state ensemble introduced by ssDNA and dsDNA molecules demonstrate the important role of colocalization in biological self-assembly processes.

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.

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.

Conformational heterogeneity in proteins is known to often be the key to their function. We present a coarse grained model to explore the interplay between protein structure, folding and function which is applicable to allosteric or non-allosteric proteins. We employ the model to study the detailed mechanism of the reversible conformational transition of Adenylate Kinase (AKE) between the open to the closed conformation, a reaction that is crucial to the protein's catalytic function. We directly observe high strain energy which appears to be correlated with localized unfolding during the functional transition. This work also demonstrates that competing native interactions from the open and closed form can account for the large conformational transitions in AKE. We further characterize the conformational transitions with a new measure Phi(Func), and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.

Although engineering of transcription factors and DNA-modifying enzymes has drawn substantial attention for artificial gene regulation and genome editing, most efforts focus on affinity and specificity of the DNA-binding proteins, typically overlooking the kinetic properties of these proteins. However, a simplistic pursuit of high affinity can lead to kinetically deficient proteins that spend too much time at nonspecific sites before reaching their targets on DNA. We demonstrate that structural dynamic knowledge of the DNA-scanning process allows for kinetically and thermodynamically balanced engineering of DNA-binding proteins. Our current study of the zinc-finger protein Egr-1 (also known as Zif268) and its nuclease derivatives reveals kinetic and thermodynamic roles of the dynamic conformational equilibrium between two modes during the DNA-scanning process: one mode suitable for search and the other for recognition. By mutagenesis, we were able to shift this equilibrium, as confirmed by NMR spectroscopy. Using fluorescence and biochemical assays as well as computational simulations, we analyzed how the shifts of the conformational equilibrium influence binding affinity, target search kinetics, and efficiency in displacing other proteins from the target sites. A shift toward the recognition mode caused an increase in affinity for DNA and a decrease in search efficiency. In contrast, a shift toward the search mode caused a decrease in affinity and an increase in search efficiency. This accelerated site-specific DNA cleavage by the zinc-finger nuclease, without enhancing off-target cleavage. Our study shows that appropriate modulation of the dynamic conformational ensemble can greatly improve zinc-finger technology, which has used Egr-1 (Zif268) as a major scaffold for engineering.

ssDNA binding proteins (SSBs) protect ssDNA from chemical and enzymatic assault that can derail DNA processing machinery. Complexes between SSBs and ssDNA are often highly stable, but predicting their structures is challenging, mostly because of the inherent flexibility of ssDNA and the geometric and energetic complexity of the interfaces that it forms. Here, we report a newly developed coarse-grained model to predict the structure of SSB–ssDNA complexes. The model is successfully applied to predict the binding modes of six SSBs with ssDNA strands of lengths of 6–65 nt. In addition to charge–charge interactions (which are often central to governing protein interactions with nucleic acids by means of electrostatic complementarity), an essential energetic term to predict SSB–ssDNA complexes is the interactions between aromatic residues and DNA bases. For some systems, flexibility is required from not only the ssDNA but also, the SSB to allow it to undergo conformational changes and the penetration of the ssDNA into its binding pocket. The association mechanisms can be quite varied, and in several cases, they involve the ssDNA sliding along the protein surface. The binding mechanism suggests that coarse-grained models are appropriate to study the motion of SSBs along ssDNA, which is expected to be central to the function carried out by the SSBs.

Human transthyretin (TTR) is an amyloidogenic protein whose aggregation is associated with several types of amyloid diseases. The following mechanism of TTR amyloid formation has been proposed. TTR tetramer at first dissociates into native monomers, which is the rate-limiting step in fibril formation. The monomeric species then partially unfold to form amyloidogenic intermediates that subsequently undergo a downhill self-assembly process. The amyloid deposit can be facilitated by disease-associated point mutations. However, only subtle structural differences were observed between the crystal structures of the wild type and the disease-associated variants. To investigate how single-point mutations influence the effective energy landscapes of TTR monomers, molecular dynamics (MD) simulations were performed on wild-type TTR and two pathogenic variants. Principal coordinate analysis on MD-generated ensembles has revealed multiple unfolding pathways for each protein. Amyloidogenic intermediates with the dislocated C strand-loop-D strand motif were observed only on the unfolding pathways of V30M and L55P variants and not for wild-type TTR. Our study suggests that the sequence-dependent unfolding pathway plays a crucial role in the amyloidogenicity of TTR. Analyses of side chain concerted motions indicate that pathogenic mutations on "edge strands" disrupt the delicate side chain correlated motions, which in turn may alter the sequence of unfolding events.

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.

Conformational transitions are thought to be the prime mechanism of prion diseases. In this study, the energy landscapes of a wild-type prion protein (PrP) and the D178N and E200K mutant proteins were mapped, enabling the characterization of the normal isoforms (PrPC) and partially unfolded isoforms (PrPPU) of the three prion protein analogs. It was found that the three energy landscapes differ in three respects: (i) the relative stability of the PrPC and the PrPPU states, (ii) the transition pathways from PrPC to PrPPU, and (iii) the relative stability of the three helices in the PrPC state. In particular, it was found that although helix 1 (residues 144-156) is the most stable helix in wild-type PrP, its stability is dramatically reduced by both mutations. This destabilization is due to changes in the charge distribution that affects the internal salt bridges responsible for the greater stability of this helix in wild-type PrP. Although both mutations result in similar destabilization of helix 1, they a have different effect on the overall stability of PrPC and of PrPPU isoforms and on structural properties. The destabilization of helix 1 by mutations provides additional evidences to the role of this helix in the pathogenic transition from the PrPC to the pathogenic isoform PrPSC. Proteins 2002;47:458–468. © 2002 Wiley-Liss, Inc.

Multidomain transcription factors, which are especially abundant in eukaryotic genomes, are advantageous to accelerate the search kinetics for target site because they can follow the intersegment transfer via the monkey-bar mechanism in which the protein forms a bridged intermediate between two distant DNA regions. Monkey-bar dynamics highly depends on the properties of the multidomain protein (the affinity of each of the constituent domains to the DNA and the length of the linker) and the DNA molecules (their inter-distance and inter-angle). In this study, we investigate using coarse-grained molecular dynamics simulations how the local conformation of the DNA may affect the DNA search performed by a multidomain protein Pax6 in comparison to that of the isolated domains. Our results suggest that in addition to the common rotation-coupled translation along the DNA major groove, for curved DNA the tethered domains may slide in a rotation-decoupled sliding mode. Furthermore, the multidomain proteins move by longer jumps on curved DNA compared with those performed by the single domain protein. The long jumps originate from the DNA curvature bringing two sequentially distant DNA sites into close proximity with each other and they suggest that multidomain proteins may move on highly curved DNA faster than linear DNA.

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

Rapid recognition by a protein of its DNA target site is achieved through a combination of one- and three-dimensional (1D and 3D) diffusion, which allows efficient scanning of the many alternative sites. This facilitated diffusion mechanism is expected to be affected by cellular conditions, particularly crowding, given that up to 40% of the total cellular volume may by occupied by macromolecules. Both experimental and theoretical studies showed that crowding particles can enhance facilitated diffusion and accelerate search kinetics. This effect may originate from crowding forcing a trade-off between 3D and 1D diffusion. In this study, using coarse-grained molecular dynamic simulations, we investigate how the molecular properties of the crowders may modulate the effect exerted by crowding on a searcher protein. We show that crowders with an affinity to the DNA are less effective search facilitators than particles whose contribution is solely entropic. Crowders that have affinity to DNA may occupy DNA sites and thereby function as obstacles or roadblocks that slow down the searcher protein, and they may also produce a smaller excluded volume effect and so reduce usage of the hopping searching mode in favor of less-effective 3D diffusion in the bulk. We discuss how strong repulsive interactions between the crowding particles themselves may affect the overall dynamics of the crowders and their excluded volume effect. Our study shows that search kinetics and its mechanism are modulated not only by salt concentration and crowding occupancy, but also by the properties of the crowding particles.

We present an integrated experimental and computational study of the molecular mechanisms by which myristoylation affects protein folding and function, which has been little characterized to date. Myristoylation, the covalent linkage of a hydrophobic C14 fatty acyl chain to the N-terminal glycine in a protein, is a common modification that plays a critical role in vital regulated cellular processes by undergoing reversible energetic and conformational switching. Coarse-grained folding simulations for the model pH-dependent actin- and membrane-binding protein hisactophilin reveal that nonnative hydrophobic interactions of the myristoyl with the protein as well as nonnative electrostatic interactions have a pronounced effect on folding rates and thermodynamic stability. Folding measurements for hydrophobic residue mutations of hisactophilin and atomistic simulations indicate that the nonnative interactions of the myristoyl group in the folding transition state are nonspecific and robust, and so smooth the energy landscape for folding. In contrast, myristoyl interactions in the native state are highly specific and tuned for sensitive control of switching functionality. Simulations and amide hydrogen exchange measurements provide evidence for increases as well as decreases in stability localized on one side of the myristoyl binding pocket in the protein, implicating strain and altered dynamics in switching. The effects of folding and function arising from myristoylation are profoundly different from the effects of other post-translational modifications.

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.