Publications

Biomolecular condensates are essential for cellular functionality, yet the complex interplay among the diverse molecular interactions that mediate their formation remains poorly understood. Here, using coarse-grained molecular dynamics simulations, we address the contribution of cation-π interactions to the stability of condensates formed via liquid-liquid phase separation. We found greater stabilization of up to 80% via cation-π interactions in condensates formed from peptides with higher aromatic residue content or less charge clustering. The contribution of cation-π interactions to droplet stability increases with increasing ionic strength, suggesting a trade-off between cation-π and electrostatic interactions. Cation-π interactions, therefore, can compensate for reduced electrostatic interactions, such as occurs at higher salt concentrations and in sequences with less charged residue content or clustering. Designing condensates with desired biophysical characteristics therefore requires quantification not only of the individual interactions but also cross-talks involving charge-charge, π-π, and cation-π interactions.

In eukaryotes, many DNA/RNA-binding proteins possess intrinsically disordered regions (IDRs) with large negative charge, some of which involve a consecutive sequence of aspartate (D) or glutamate (E) residues. We refer to them as D/E repeats. The functional role of D/E repeats is not well understood, though some of them are known to cause autoinhibition through intramolecular electrostatic interaction with functional domains. In this work, we investigated the impacts of D/E repeats on the target DNA search kinetics for the high-mobility group box 1 (HMGB1) protein and the artificial protein constructs of the Antp homeodomain fused with D/E repeats of varied lengths. Our experimental data showed that D/E repeats of particular lengths can accelerate the target association in the overwhelming presence of non-functional high-Affinity ligands ('decoys'). Our coarse-grained molecular dynamics (CGMD) simulations showed that the autoinhibited proteins can bind to DNA and transition into the uninhibited complex with DNA through an electrostatically driven induced-fit process. In conjunction with the CGMD simulations, our kinetic model can explain how D/E repeats can accelerate the target association process in the presence of decoys. This study illuminates an unprecedented role of the negatively charged IDRs in the target search process.

Silk is a unique, remarkably strong biomaterial made of simple protein building blocks. To date, no synthetic method has come close to reproducing the properties of natural silk, due to the complexity and insufficient understanding of the mechanism of the silk fiber formation. Here, we use a combination of bulk analytical techniques and nanoscale analytical methods, including nano-infrared spectroscopy coupled with atomic force microscopy, to probe the structural characteristics directly, transitions, and evolution of the associated mechanical properties of silk protein species corresponding to the supramolecular phase states inside the silkworm's silk gland. We found that the key step in silk-fiber production is the formation of nanoscale compartments that guide the structural transition of proteins from their native fold into crystalline β-sheets. Remarkably, this process is reversible. Such reversibility enables the remodeling of the final mechanical characteristics of silk materials. These results open a new route for tailoring silk processing for a wide range of new material formats by controlling the structural transitions and self-assembly of the silk protein's supramolecular phases.

Intrinsically disordered regions (IDRs) of proteins are often characterized by a high fraction of charged residues, but differ in their overall net charge and in the organization of the charged residues. The function-encoding information stored via IDR charge composition and organization remains elusive. Here, we aim to decipher the sequence–function relationship in IDRs by presenting a comprehensive bioinformatic analysis of the charge properties of IDRs in the human, mouse, and yeast proteomes. About 50% of the proteins comprise at least a single IDR, which is either positively or negatively charged. Highly negatively charged IDRs are longer and possess greater net charge per residue compared with highly positively charged IDRs. A striking difference between positively and negatively charged IDRs is the characteristics of the repeated units, specifically, of consecutive Lys or Arg residues (K/R repeats) and Asp or Glu (D/E repeats) residues. D/E repeats are found to be about five times longer than K/R repeats, with the longest found containing 49 residues. Long stretches of consecutive D and E are found to be more prevalent in nucleic acid-related proteins. They are less common in prokaryotes, and in eukaryotes their abundance increases with genome size. The functional role of D/E repeats and the profound differences between them and K/R repeats are discussed.

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.

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.

The respiratory and intestinal tracts are exposed to physical and biological hazards accompanying the intake of air and food. Likewise, the vasculature is threatened by inflammation and trauma. Mucin glycoproteins and the related von Willebrand factor guard the vulnerable cell layers in these diverse systems. Colon mucins additionally house and feed the gut microbiome. Here, we present an integrated structural analysis of the intestinal mucin MUC2. Our findings reveal the shared mechanism by which complex macromolecules responsible for blood clotting, mucociliary clearance, and the intestinal mucosal barrier form protective polymers and hydrogels. Specifically, cryo-electron microscopy and crystal structures show how disulfide-rich bridges and pH-tunable interfaces control successive assembly steps in the endoplasmic reticulum and Golgi apparatus. Remarkably, a densely O-glycosylated mucin domain performs an organizational role in MUC2. The mucin assembly mechanism and its adaptation for hemostasis provide the foundation for rational manipulation of barrier function and coagulation.

Proteins with intrinsically disordered regions have a tendency to condensate via liquid-liquid phase separation both in vitro and in vivo. Such biomolecular coacervates play various significant roles in biologically important regulatory processes. The present work explores the structural and dynamic features of coacervates formed by model polyampholytes, being intrinsically disordered proteins, that differ in terms of their charged amino acid patterns. Differences in the distribution of charged amino acids along the polyampholyte sequence lead to distinctly different structural features in the dense phase and hence to different liquid properties. Increased charge clustering raises the critical temperature for phase separation and results in each polyampholyte experiencing a larger number of inter-chain contacts with neighboring proteins in the condensate. Consequently, polyampholytes with greater charge clustering adopt a much more extended conformation, having a radius of gyration up to twice that observed in the dilute bulk phase. Translational diffusion within the droplet is pronounced, being just 4-20 times slower than in the bulk, consistently with the high conformational entropy in the dense phase and high exchange rate of the network of intermolecular interactions in the condensate. Coupled to the faster diffusion, the condensate also adopts a more elongated shape and exhibits imperfect packing, which results in cavities. This study quantifies the fundamental microscopic properties of condensates including the effect of long-range electrostatic forces and particularly how they can be modulated by the charge pattern.

Thermodynamic stability is an important property of proteins that is linked to many of the trade-offs that characterize a protein molecule and therefore its function. Designing a protein with a desired stability is a complicated task given the intrinsic trade-off between enthalpy and entropy which applies for both the folded and unfolded states. Traditionally, protein stability is manipulated by point mutations which regulate the folded state enthalpy. In some cases, the entropy of the unfolded state has also been manipulated by means that drastically restrict its conformational dynamics such as engineering disulfide bonds. In this mini-review, we survey various approaches to modify protein stability by manipulating the entropy of either the unfolded or the folded states. We show that point mutations that involve elimination of long-range contacts may have a greater destabilization effect than mutations that eliminate shorter-range contacts. Protein conjugation can also affect the entropy of the unfolded state and thus the overall stability. In addition, we show that entropy can contribute to shape the folded state and yield greater protein stabilization. Hence, we argue that the entropy component can be practically manipulated both in the folded and unfolded state to modify protein stability.

DNA sequences are often recognized by multi-domain proteins that may have higher affinity and specificity than single-domain proteins. However, the higher affinity to DNA might be coupled with slower recognition kinetics. In this study, we address this balance between stability and kinetics for multi-domain Cys(2)His(2)- (C2H2-) type zinc-finger (ZF) proteins. These proteins are the most prevalent DNA-binding domain in eukaryotes and C2H2 type zinc-finger proteins (C2H2-ZFPs) constitute nearly one-half of all known and predicted transcription factors in human. Extensive contact with DNA via tandem ZF domains confers high stability on the sequence-specific complexes. However, this can limit target search efficiency, especially for low abundance ZFPs. Earlier, we found that asymmetrical distribution of electrostatic charge among the three ZF domains of the low abundance transcription factor Egr-1 facilitates its DNA search process. Here, on a diverse set of 273 human C2H2-ZFP comprised of 3-15 tandem ZF domains, we find that, in many cases, electrostatic charge and binding specificity are asymmetrically distributed among the ZF domains so that neighbouring domains have different DNA-binding properties. For proteins containing 3-6 ZF domains, we show that the low abundance proteins possess a higher degree of non-specific asymmetry and vice versa. Our findings suggest that where the electrostatics of tandem ZF domains are similar (i.e., symmetrical), the ZFPs are more abundant to optimize their DNA search efficiency. This study reveals new insights into the fundamental determinants of recognition by C2H2-ZFPs of their DNA binding sites in the cellular landscape. The importance of electrostatic asymmetry with respect to binding site recognition by C2H2-ZFPs suggests the possibility that it may also be important in other ZFP systems and reveals a new design feature for zinc finger engineering.Author summaryOptimal recognition of proteins to DNA is governed by various factors among them the thermodynamics, kinetics and specificity of the protein-DNA complex. Multi-domain DNA-binding proteins are expected to have higher affinity and specificity due to the extensive interface they form with DNA. However, larger interface may result with higher friction when these proteins scan the DNA for the target site via the sliding mechanism. A way to overcome this drawback is to have asymmetry in the protein so that the interface with DNA is smaller. Alternatively, higher abundance can also increase the search speed. Here, using computational analysis of large data set of multi-domain zinc finger DNA-binding proteins, we report a trade-off between asymmetry and abundance.

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.

The current report extends the facilitated diffusion model to account for conflict between the search and recognition binding modes adopted by DNA-binding proteins (DBPs) as they search DNA and subsequently recognize and bind to their specific binding site. The speed of the search dynamics is governed by the energetic ruggedness of the protein-DNA landscape, whereas the rate for the recognition process is mostly dictated by the free energy barrier for the transition between the DBP's search and recognition binding modes. We show that these two modes are negatively coupled, such that fast 1D sliding and rapid target site recognition probabilities are unlikely to coexist. Thus, a tradeoff occurs between optimizing the timescales for finding and binding the target site. We find that these two kinetic properties can be balanced to produce a fast timescale for the total target search and recognition process by optimizing frustration. Quantification of the facilitated diffusion model by including a frustration term enables it to explain several experimental observations concerning search and recognition speeds. The extended model captures experimental estimate of the energetic ruggedness of the protein-DNA landscape and predicts how various molecular properties of protein-DNA binding affect recognition kinetics. Particularly, point mutations may change the frustration and so affect protein association with DNA, thus providing a means to modulate protein-DNA affinity by manipulating the protein's association or dissociation reactions.

The alpha-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.

Protein backbone alternation due to insertion/deletion or mutation operation often results in a change of fundamental biophysical properties of proteins. The proposed work intends to encode the protein stability changes associated with single point deletions (SPDs) of amino acids in proteins. The encoding will help in the primary screening of detrimental backbone modifications before opting for expensive in vitro experimentations. In the absence of any benchmark database documenting SPDs, we curate a data set containing SPDs that lead to both folded conformations and unfolded state. We differentiate these SPD instances with the help of simple structural and physicochemical features and eventually classify the foldability resulting out of SPDs using a Random Forest classifier and an Elliptic Envelope based outlier detector. Adhering to leave one out cross validation, the accuracy of the Random Forest classifier and the Elliptic Envelope is of 99.4% and 98.1%, respectively. The newly defined database and the delineation of SPD instances based on its resulting foldability provide a head start toward finding a solution to the given problem.

Predicting the effect of a single point mutation on protein thermodynamic stability (Delta Delta G) is an ongoing challenge with high relevance for both fundamental and applicable aspects of protein science. Drawbacks that limit the predictive power of stability prediction tools include the lack of representations for the explicit energetic terms of the unfolded state. Using coarse-grained simulations and analytical modeling analysis, we found that a mutation that involves the breaking of long-range contacts may lead to an increase in the unfolded state entropy, which can lead to an overall destabilization of the protein. A bioinformatics analysis indicates that the effect of mutation on the unfolded state is greater for hydrophobic or charged (compared with polar) residues that participate in long-range contacts through a loop length longer than 18 amino acids and whose formation probabilities are relatively high.

Many mutations that cause familial hypercholesterolemia localize to ligand-binding domain 5 (LA5) of the low-density lipoprotein receptor, motivating investigation of the folding and misfolding of this small, disulfide-rich, calcium-binding domain. LA5 folding is known to involve non-native disulfide isomers, yet these folding intermediates have not been structurally characterized. To provide insight into these intermediates, we used nuclear magnetic resonance (NMR) to follow LA5 folding in real time. We demonstrate that misfolded or partially folded disulfide intermediates are indistinguishable from the unfolded state when focusing on the backbone NMR signals, which provide information on the formation of only the final, native state. However,
¹³C labeling of cysteine side chains differentiated transient intermediates from the unfolded and native states and reported on disulfide bond formation in real time. The cysteine pairings in a dominant intermediate were identified using
¹³C-edited three-dimensional NMR, and coarse-grained molecular dynamics simulations were used to investigate the preference of this disulfide set over other non-native arrangements. The transient population of LA5 species with particular non-native cysteine connectitivies during folding supports the conclusion that cysteine pairing is not random and that there is a bias toward certain structural ensembles during the folding process, even prior to the binding of calcium.

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.

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.

The key feature explaining the rapid recognition of a DNA target site by its protein lies in the 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. Using coarse-grained molecular dynamics and Monte Carlo simulations, we show that the crowding particles can enhance facilitated diffusion and accelerate search kinetics. This effect originates from a trade-off between 3D and 1D diffusion. The 3D diffusion coefficient is lower under crowded conditions, but it has little influence because the excluded volume effect of molecular crowding restricts its use. Largely prevented from using 3D diffusion, the searching protein dramatically increases its use of the hopping search mode, which results in a higher linear diffusion coefficient. The coefficient of linear diffusion also increases under crowded conditions as a result of increased collisions between the crowding particles and the searching protein. Overall, less 3D diffusion coupled with an increase in the use of the hopping and speed of 1D diffusion :results in faster search kinetics Under crowded conditions. Our study shows, that the search kinetics and mechanism are modulated not only by the crowding occupancy but also by the properties of the crowding particles and the salt concentration.

Understanding the molecular mechanism and the fast kinetics of DNA target site recognition by a protein is essential to decipher genetic activity in the cell. The speed of searching DNA may depend on the structural complexity of the proteins and the DNA molecules as well as the cellular environment. Coarse-grained (CG) molecular dynamics simulations are powerful means to investigate the molecular details of the search performed by protein to locate the target sites. Recent studies showed how different proteins scan DNA and how the search efficiency can be enhanced and regulated by the protein properties. In this review, we discuss computational approaches to study the physical chemistry of DNA search processes using CG molecular dynamics simulations and their advantage in covering long time-scale biomolecular processes. WIREs Comput Mol Sci 2016, 6:515–531. doi: 10.1002/wcms.1262. For further resources related to this article, please visit the WIREs website.

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 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.

Many novel applications in bioelectronics rely on the interaction between biomolecules and electronically conducting substrates. However, crucial knowledge about the relation between electronic transport via peptides and their amino-acid composition is still absent. Here, we report results of electronic transport measurements via several homopeptides as a function of their structural properties and temperature. We demonstrate that the conduction through the peptide depends on its length and secondary structure as well as on the nature of the constituent amino acid and charge of its residue. We support our experimental observations with high-level electronic structure calculations and suggest off-resonance tunneling as the dominant conduction mechanism via extended peptides. Our findings indicate that both peptide composition and structure can affect the efficiency of electronic transport across peptides.

Water molecules are abundant in protein-DNA interfaces, especially in their nonspecific complexes. In this study, we investigated the organization and energetics of the interfacial water by simplifying the geometries of the proteins and the DNA to represent them as two equally and oppositely charged planar surfaces immersed in water. We found that the potential of mean force for bringing the two parallel surfaces into close proximity comprises energetic barriers whose properties strongly depend on the charge density of the surfaces. We demonstrated how the organization of the water molecules into discretized layers and the corresponding energetic barriers to dehydration can be modulated by the charge density on the surfaces, salt, and the structure of the surfaces. The 1-2 layers of ordered water are tightly bound to the charged surfaces representing the nonspecific protein-DNA complex. This suggests that water might mediate one-dimensional diffusion of proteins along DNA (sliding) by screening attractive electrostatic interactions between the positively charged molecular surface on the protein and the negatively charged DNA backbone and, in doing so, reduce intermolecular friction in a manner that smoothens the energetic landscape for sliding, and facilitates the 1D diffusion of the protein.

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.

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 D₁, 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.

Abstract The following sections are included: Introduction Methods Statistical analysis of acetylation and phosphorylation in DNA-binding proteins Simulation model for searching the DNA of modified proteins Results The occurrence of acetylation and phosphorylation at disordered regions The occurrence of acetylation and phosphorylation at disordered tails The interplay between post-translational modifications of disordered Conclusions Acknowledgments

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.

Protein ubiquitination is central to the regulation of various pathways in eukaryotes. The process of ubiquitination and its cellular outcome were investigated in hundreds of proteins to date. Despite this, the evolution of this regulatory mechanism has not yet been addressed comprehensively. Here, we quantify the rates of evolutionary changes of ubiquitination and SUMOylation (Small Ubiquitin-like MOdifier) sites. We estimate the time at which they first appeared, and compare them to acetylation and phosphorylation sites and to unmodified residues. We observe that the various modification sites studied exhibit similar rates. Mammalian ubiquitination sites are weakly more conserved than unmodified lysine residues, and a higher degree of relative conservation is observed when analyzing bona fide ubiquitination sites. Various reasons can be proposed for the limited level of excess conservation of ubiquitination, including shifts in locations of the sites, the presence of alternative sites, and changes in the regulatory pathways. We observe that disappearance of sites may be compensated by the presence of a lysine residue in close proximity, which is significant when compared to evolutionary patterns of unmodified lysine residues, especially in disordered regions. This emphasizes the importance of analyzing a window in the vicinity of functional residues, as well as the capability of the ubiquitination machinery to ubiquitinate residues in a certain region. Using prokaryotic orthologs of ubiquitinated proteins, we study how ubiquitination sites were formed, and observe that while sometimes sequence additions and rearrangements are involved, in many cases the ubiquitination machinery utilizes an already existing sequence without significantly changing it. Finally, we examine the evolution of ubiquitination, which is linked with other modifications, to infer how these complex regulatory modules have evolved. Our study gives initial insights into the formation of ubiquitination sites, their degree of conservation in various species, and their co-evolution with other posttranslational modifications.

DNA-binding proteins execute many different genetic tasks, such as transcription, repression, replication, and the repair of damaged DNA. Given the size of the genome, the large number of molecular species in vivo, and the crowded environment in which these processes take place, the remarkable efficiency and specificity of protein–DNA recognition present a major theoretical puzzle. In particular, two fundamental issues of great interest in protein– DNA recognition are: how proteins achieve a high degree of selectivity among a huge background of competing nonspecific DNA sequences and the aspect of the rapid recognition rate of the DNA target sequence by the protein.

Conjugating flexible polymers (such as oligosaccharides) to proteins or confining a protein in a restricted volume often increases protein thermal stability. In this communication, we investigate the interplay between conjugation and confinement which is not trivial as the magnitude and the mechanism of stabilization are different in each instance. Using coarse-grained computational approach the folding biophysics is studied when the protein is placed in a sphere of variable radius and is conjugated to 0-6 mono- or penta-saccharides. We observe a synergistic effect on thermal stability when short oligosaccharides are attached and the modified protein is confined in a small cage. However, when large oligosaccharides are added, a conflict between confinement and glycosylation arises as the stabilizing effect of the cage is dramatically reduced and it is almost impossible to further stabilize the protein beyond the mild stabilization induced by the sugars.

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.

Asparagine glycosylation is one of the most common and important post-translational modifications of proteins in eukaryotic cells. N-Glycosylation occurs when a triantennary glycan precursor is transferred en bloc to a nascent polypeptide (harboring the N-X-T/S sequon) as the peptide is cotranslationally translocated into the endoplasmic reticulum (ER). In addition to facilitating binding interactions with components of the ER proteostasis network, N-glycans can also have intrinsic effects on protein folding by directly altering the folding energy landscape. Previous work from our laboratories (Hanson et al. Proc. Natl. Acad. Sci. U.S.A. 2009, 109, 3131-3136; Shental-Bechor, D.; Levy, Y. Proc. Natl. Acad. Sci. U.S.A. 2008, 105, 8256-8261) suggested that the three sugar residues closest to the protein are sufficient for accelerating protein folding and stabilizing the resulting structure in vitro; even a monosaccharide can have a dramatic effect. The highly conserved nature of these three proximal sugars in N-glycans led us to speculate that introducing an N-glycosylation site into a protein that is not normally glycosylated would stabilize the protein and increase its folding rate in a manner that does not depend on the presence of specific stabilizing protein-saccharide interactions. Here, we test this hypothesis experimentally and computationally by incorporating an N-linked GlcNAc residue at various positions within the Pin WW domain, a small β-sheet-rich protein. The results show that an increased folding rate and enhanced thermodynamic stability are not general, context-independent consequences of N-glycosylation. Comparison between computational predictions and experimental observations suggests that generic glycan-based excluded volume effects are responsible for the destabilizing effect of glycosylation at highly structured positions. However, this reasoning does not adequately explain the observed destabilizing effect of glycosylation within flexible loops. Our data are consistent with the hypothesis that specific, evolved protein-glycan contacts must also play an important role in mediating the beneficial energetic effects on protein folding that glycosylation can confer.

Protein flexibility is thought to play key roles in numerous biological processes, including antibody affinity maturation, signal transduction, and enzyme catalysis, yet only limited information is available regarding the molecular details linking protein dynamics with function. A single point mutation at the distal site of the endogenous tissue inhibitor of metalloproteinase 1 (TIMP-1) enables this clinical target protein to tightly bind and inhibit membrane type 1 matrix metalloproteinase (MT1-MMP) by increasing only the association constant. The high-resolution X-ray structure of this complex determined at 2 could not explain the mechanism of enhanced binding and pointed to a role for protein conformational dynamics. Molecular dynamics (MD) simulations reveal that the high-affinity TIMP-1 mutants exhibit significantly reduced binding interface flexibility and more stable hydrogen bond networks. This was accompanied by a redistribution of the ensemble of substrates to favorable binding conformations that fit the enzyme catalytic site. Apparently, the decrease in backbone flexibility led to a lower entropy cost upon formation of the complex. This work quantifies the effect of a single point mutation on the protein conformational dynamics and function of TIMP-1. Here we argue that controlling the intrinsic protein dynamics of MMP endogenous inhibitors may be utilized for rationalizing the design of selective novel protein inhibitors for this class of enzymes.

Protein ubiquitination controls the cellular fate of numerous eukaryotic proteins. Despite its importance, many fundamental questions remain regarding its mechanism. One such question is how ubiquitination alters the biophysical properties of the modified protein and whether these alterations are significant in the cellular context. In this study, we investigate the effects of ubiquitination on the folding thermodynamics and mechanism of various substrates using computational tools and find that ubiquitination changes the thermal stability of modified proteins in a manner relevant to cellular processes. These changes depend on the substrate modification site and on the type of ubiquitination. Ubiquitination of the substrate Ubc7 at the residues that are modified in vivo prior to proteasomal degradation uniquely results in significant thermal destabilization and a local unwinding near the modification site, which indicates that ubiquitination possibly facilitates the unfolding process and improves substrate degradation efficiency. With respect to the substrate p19^4inkd, our results support a synergetic effect of ubiquitination and phosphorylation on the degradation process via enhanced thermal destabilization. Our study implies that, in addition to its known role as a recognition signal, the ubiquitin attachment may be directly involved in the cellular process it regulates by changing the biophysical properties of the substrate.

Glycosylation is among the most common post-translational modifications that proteins undergo that may affect many of their activities. It may also modify the underlying energy landscape of glycoproteins in a way that their altered biophysical characteristics are linked to their bioactivity. Yet, the capability of glycosylation to modify thermodynamic and kinetic properties varies greatly between glycoproteins. Deciphering the 'glycosylation code' that dictates the interplay between the nature of the carbohydrates or the proteins and the biophysical properties of the glycosylated proteins is essential. In this article, we discuss how the size, number, and structure of glycans, as well as the attachment sites, may modulate the folding of glycoproteins. Understanding the cross-talks between the protein and the attached glycans at the molecular level may assist in tailoring the biophysical properties of proteins in general.

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 p-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. (C) 2009 Elsevier Ltd. All rights reserved.

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 ∼ 20% sliding along the DNA and ∼ 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.

The complexity of the mechanisms by which proteins fold has been shown by many studies to be governed by their native-state topologies. This was manifested in the ability of the native topology-based model to capture folding mechanisms and the success of folding rate predictions based on various topological measures, such as the contact order. However, while the finer details of topological complexity have been thoroughly examined and related to folding kinetics, simpler characteristics of the protein, such as its overall shape, have been largely disregarded. In this study, we investigated the folding of proteins with an unusual elongated geometry that differs substantially from the common globular structure. To study the effect of the elongation degree on the folding kinetics, we used repeat proteins, which become more elongated as they include more repeating units. Some of these have apparently anomalous experimental folding kinetics, with rates that are often less than expected on the basis of rates for globular proteins possessing similar topological complexity. Using experimental folding rates and a larger set of rates obtained from simulations, we have shown that as the protein becomes increasingly elongated, its folding kinetics becomes slower and deviates more from the rate expected on the basis of topology measures fitted for globular proteins. The observed slow kinetics is a result of a more complex pathway in which stable intermediates composed of several consecutive repeats can appear. We thus propose a novel measure, an elongation-sensitive contact order, that takes into account both the extent of elongation and the topological complexity of the protein. This new measure resolves the apparent discrimination between the folding of globular and elongated repeat proteins. Our study extends the current capabilities of folding-rate predictions by unifying the kinetics of repeat and globular proteins.

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 Φ_Func, and demonstrate that local unfolding may be due, in part, to competing intra-protein interactions.

Enhanced structural insights into the folding energy landscape of the N-terminal dimerization domain of Escherichia coli tryptophan repressor, [2-66]2 TR, were obtained from a combined experimental and theoretical analysis of its equilibrium folding reaction. Previous studies have shown that the three intertwined helices in [2-66]2 TR are sufficient to drive the formation of a stable dimer for the full-length protein, [2-107]2 TR. The monomeric and dimeric folding intermediates that appear during the folding reactions of [2-66]2 TR have counterparts in the folding mechanism of the full-length protein. The equilibrium unfolding energy surface on which the folding and dimerization reactions occur for [2-66]2 TR was examined with a combination of native-state hydrogen exchange analysis, pepsin digestion and matrix-assisted laser/desorption mass spectrometry performed at several concentrations of protein and denaturant. Peptides corresponding to all three helices in [2-66]2 TR show multi-layered protection patterns consistent with the relative stabilities of the dimeric and monomeric folding intermediates. The observation of protection exceeding that offered by the dimeric intermediate in segments from all three helices implies that a segment-swapping mechanism may be operative in the monomeric intermediate. Protection greater than that expected from the global stability for a single amide hydrogen in a peptide from the C-helix possibly and another from the A-helix may reflect non-random structure, possibly a precursor for segment swapping, in the urea-denatured state. Native topology-based model simulations that correspond to a funnel energy landscape capture both the monomeric and dimeric intermediates suggested by the HX MS data and provide a rationale for the progressive acquisition of secondary structure in their conformational ensembles.

We studied the energetics and fragmentation patterns of multicharged (A⁺)_n Morse clusters (n = 55-321), with a total cluster charge Z = n. The Morse pair-potential parameters were characterized by the dissociation energy D = 1-10 eV, range parameter α = 1-3 Å^-1, and interatomic equilibrium separation R_e = 1-3 Å. The potential energies ε (per particle) of these multicharged Morse clusters at their equilibrium configuration (with bond length r ₀) were analyzed in terms of the liquid drop model. This resulted in the relation ε =(ā_C⁰/r₀)n ^2/3+(ā_v⁰D/αr₀) +[ā_s⁰D/(αr₀)^3/2]n ^-1/3, where the reduced parameters ā_C⁰ (for the Coulomb energy), ā_v⁰ (for the interior energy) and ā_s⁰ (for the surface energy) are independent of the Morse pair-potential parameters. The Rayleigh fissibility parameter X = E(Coulomb)/2E(surface), which determines the fragmentation pattern (i.e., X 1 for Coulomb explosion), was expressed in the form X=(Z²/n)[(2ā_s ⁰/ā_C⁰)(D/α^3/2r ₀^1/2)]^-1. The application of this result to the Coulomb instability of multicharged globular proteins reveals that X

Minding your p's and q's has become as important to protein-folding theorists as it is for those being instructed in the rules of etiquette. To assess the quality of structural reaction coordinates in predicting the transition-state ensemble (TSE) of protein folding, we benchmarked the accuracy of four structural reaction coordinates against the kinetic measure P _fold, whose value of 0.50 defines the stochastic separatrix for a two-state folding mechanism. For two proteins that fold by a simple two-state mechanism, c-src SH3 and Cl-2, the Φ-values of the TSEs predicted by native topology-based reaction coordinates (including Q, the fraction of native contacts) are almost identical to those of the TSE based on P_fold, with correlation coefficients of >0.90. For proteins with complex folding mechanisms that have especially broad, asymmetrical free energy barriers such as the designed 3-ankyrin repeating protein (3ANK) or proteins with distinct intermediates such as cyanovirin-N (CV-N), we show that the ensemble of structures with P_fold = 0.50 generally does not include the chemically relevant transition states. This weakness of P_fold limits its usefulness in protein folding studies. For such systems, elucidating the essential features of folding mechanisms requires using multiple reaction coordinates, although the number is still rather small. At the same time, for simple folding mechanisms, there is no indication of superiority for P _fold over structurally chosen and thermodynamically relevant reaction coordinates that correctly measure the degree of nativeness.

The prevalence of domain-swapping in nature is a manifestation of the principle of minimal frustration in that the interactions designed by evolution to stabilize the protein are also involved in this mode of binding. We previously demonstrated that the Symmetrized-Go potential accurately predicts the experimentally observed domain-swapped structure of Eps8 based solely on the structure of the monomer. There can be, however, multiple modes of domain-swapping, reflecting a higher level of frustration, which is a consequence of symmetry. The human prion and cyanovirin-N are too frustrated to form unique domain-swapped structures on the basis of the Symmetrized-Go potential. However, supplementing the completely symmetric model with intermolecular and intramolecular disulfide bonds in the prion and cyanovirin-N proteins, respectively, yielded unique domain-swapped structures with a remarkable similarity to the experimentally observed ones. These results suggest that the disulfide bonds may sometimes be critical in overcoming the intrinsic frustration of the symmetrized energy landscapes for domain-swapping. We also discuss the implications of intermolecular disulfide bonds in the formation of mammalian prion aggregates.

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 (PrP^C) and partially unfolded isoforms (PrP^PU) of the three prion protein analogs. It was found that the three energy landscapes differ in three respects: (i) the relative stability of the PrP^C and the PrP^PU states, (ii) the transition pathways from Pr^PC to PrP^PU, and (iii) the relative stability of the three helices in the PrP^C 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 PrP^C and of PrP^PU 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 PrP^C to the pathogenic isoform PrP^SC.

Abstract A set of 34 molecular dynamic (MD) simulations totaling 305 ns of simulation time of the prion protein-derived peptide PrP106?126 was performed with both explicit and implicit solvent models. The objective of these simulations is to investigate the relative stability of the α-helical conformation of the peptide and the mechanism for conversion from the helix to a random-coil structure. At neutral pH, the wild-type peptide was found to lose its initial helical structure very fast, within a few nanoseconds (ns) from the beginning of the simulations. The helix breaks up in the middle and then unwinds to the termini. The spontaneous transition into the random coil structure is governed by the hydrophobic interaction between His111 and Val122. The A117V mutation, which is linked to GSS disease, was found to destabilize the helix conformation of the peptide significantly, leading to a complete loss of helicity approximately 1 ns faster than in the wild-type. Furthermore, the A117V mutant exhibits a different mechanism for helix-coil conversion, wherein the helix begins to break up at the C-terminus and then gradually to unwind towards the N-terminus. In most simulations, the mutation was found to speed up the conversion through an additional hydrophobic interaction between Met112 and the mutated residue Val117, an interaction that did not exist in the wild-type peptide. Finally, the ?-sheet conformation of the wild-type peptide was found to be less stable at acidic pH due to a destabilization of the His111?Val122, since at acidic pH this histidine is protonated and is unlikely to participate in hydrophobic interaction. Proteins 2001;45:382?396. ? 2001 Wiley-Liss, Inc.

Conformation constraints are known to affect the flexibility and bioactivity of peptides. In this study we analyzed the effect of conformation constraints on the topography of the energy landscapes of three analogous hexapeptides. The three analogs vary in the degree of constraint imposed on their conformational motion: linear alanine hexapeptide with neutral terminals (Ala6), linear alanine hexapeptide with charged terminals (chrg-Ala6), and cyclic alanine hexapeptide (cyc-Ala6). It was found that significantly different energy landscapes characterize each of the three peptides, leading to different folding behaviors. Since all three analogs would be encoded by the same gene, these results suggest that nongenomic post-translational modifications may play an important role in determining the properties of proteins as well as of their folding pathways. In addition, the present study indicates that the complexity of those energy landscapes that are dominated by funnel topography can be captured by one or two reaction coordinates, such as conformational similarity to the native state. However, for more complex landscapes characterized by multiple basins such a description is insufficient. This study also shows that similar views of the landscape topography were obtained by principal component analysis (based only on local minima) and by topological mapping analysis (based on minima and barrier information). Both methods were able to resolve the complex landscape topographies for all three peptides. (C) 2001 American Institute of Physics.