Publications

Aromatic residues can participate in various biomolecular interactions, such as π-π, cation−π, and CH−π interactions, which are essential for protein structure and function. Here, we re-evaluate the geometry and energetics of these interactions using quantum mechanical (QM) calculations, focusing on pairwise interactions involving the aromatic amino acids Phe, Tyr, and Trp and the cationic amino acids Arg and Lys. Our findings reveal that π-π interactions, while energetically favorable, are less abundant in structured proteins than commonly assumed and are often overshadowed by previously underappreciated, yet prevalent, CH−π interactions. Cation−π interactions, particularly those involving Arg, show strong binding energies and a specific geometric preference toward stacked conformations, despite the global QM minimum, suggesting that a rather perpendicular T-shape conformation should be more favorable. Our results support a more nuanced understanding of protein stabilization via interactions involving aromatic residues. On the one hand, our results challenge the traditional emphasis on π-π interactions in structured proteins by showing that CH−π and cation−π interactions contribute significantly to their structure. On the other hand, π-π interactions appear to be key stabilizers in solvated regions and thus may be particularly important to the stabilization of intrinsically disordered proteins.

Mechanical energy, specifically in the form of ultrasound, can induce pressure variations and temperature fluctuations when applied to an aqueous media. These conditions can both positively and negatively affect protein complexes, consequently altering their stability, folding patterns, and self-assembling behavior. Despite much scientific progress, our current understanding of the effects of ultrasound on the self-assembly of amyloidogenic proteins remains limited. In the present study, we demonstrate that when the amplitude of the delivered ultrasonic energy is sufficiently low, it can induce refolding of specific motifs in protein monomers, which is sufficient for primary nucleation; this has been revealed by MD. These ultrasound-induced structural changes are initiated by pressure perturbations and are accelerated by a temperature factor. Furthermore, the prolonged action of low-amplitude ultrasound enables the elongation of amyloid protein nanofibrils directly from natively folded monomeric lysozyme protein, in a controlled manner, until it reaches a critical length. Using solution X-ray scattering, we determined that nanofibrillar assemblies, formed either under the action of sound or from natively fibrillated lysozyme, share identical structural characteristics. Thus, these results provide insights into the effects of ultrasound on fibrillar protein self-assembly and lay the foundation for the potential use of sound energy in protein chemistry.

The kinetics of proteinDNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. ProteinDNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRYDNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific proteinDNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance proteinDNA recognition despite its intrinsically slow kinetics.

The Origin Recognition Complex (ORC) seeds replication-fork formation by binding to DNA replication origins, which in budding yeast contain a 17bp DNA motif. High resolution structure of the ORC-DNA complex revealed two base-interacting elements: a disordered basic patch (Orc1-BP4) and an insertion helix (Orc4-IH). To define the ORC elements guiding its DNA binding in vivo, we mapped genomic locations of 38 designed ORC mutants, revealing that different ORC elements guide binding at different sites. At silencing-associated sites lacking the motif, ORC binding and activity were fully explained by a BAH domain. Within replication origins, we reveal two dominating motif variants showing differential binding modes and symmetry: a non-repetitive motif whose binding requires Orc1-BP4 and Orc4-IH, and a repetitive one where another basic patch, Orc1-BP3, can replace Orc4-IH. Disordered basic patches are therefore key for ORC-motif binding in vivo, and we discuss how these conserved, minor-groove interacting elements can guide specific ORC-DNA recognition.

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.

Peptide-RNA coacervates can result in the concentration and compartmentalization of simple biopolymers. Given their primordial relevance, peptide-RNA coacervates may have also been a key site of early protein evolution. However, the extent to which such coacervates might promote or suppress the exploration of novel peptide conformations is fundamentally unknown. To this end, we used electron paramagnetic resonance spectroscopy (EPR) to characterize the structure and dynamics of an ancient and ubiquitous nucleic acid binding element, the helix-hairpin-helix (HhH) motif, alone and in the presence of RNA, with which it forms coacervates. Double electron-electron resonance (DEER) spectroscopy applied to singly labeled peptides containing one HhH motif revealed the presence of dimers, even in the absence of RNA. Moreover, dimer formation is promoted upon RNA binding and was detectable within peptide-RNA coacervates. DEER measurements of spin-diluted, doubly labeled peptides in solution indicated transient α-helical character. The distance distributions between spin labels in the dimer and the signatures of α-helical folding are consistent with the symmetric (HhH)2-Fold, which is generated upon duplication and fusion of a single HhH motif and traditionally associated with dsDNA binding. These results support the hypothesis that coacervates are a unique testing ground for peptide oligomerization and that phase-separating peptides could have been a resource for the construction of complex protein structures via common evolutionary processes, such as duplication and fusion.

Despite the continuous discovery of host and guest proteins in membraneless organelles, complex hostguest interactions hinder the understanding of the molecular grammar governing liquidliquid 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 translationrotation (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.

Replication protein A (RPA) plays a critical role in all eukaryotic DNA processing involving singlestranded 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.

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.

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 (k(off)) for complexes involving the former.Author summary Quantifying bimolecular self-assembly is pivotal to understanding cellular function. In recent years, a large progress has been made in understanding the structure and biophysics of protein-protein interactions. Particularly, various computational tools are available for predicting these structures and to estimate their stability and the driving forces of their formation. The understating of the interactions between proteins and nucleic acids, however, is still limited, presumably due to the involvement of non-specific interactions as well as the high conformational plasticity that may demand an induced-fit mechanism. In particular, the interactions between proteins and single-stranded nucleic acids (i.e., single-stranded DNA and RNA) is very challenging due to their high flexibility. Furthermore, the interface between proteins and single-stranded nucleic acids is often chemically more heterogeneous than the interface between proteins and double-stranded DNA. In this study, we developed a coarse-grained computational model to predict the structure of complexes between proteins and single-stranded nucleic acids. The model was applied to estimate binding affinities and the estimated binding energies agreed well with the corresponding experimental binding affinities. The kinetics of association as well as the specificity of the complexes between proteins and ssDNA are different than those with ssRNA, mostly due to differences in their conformational flexibility.

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.

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.

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

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.

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

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.

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.

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.

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

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.