Role of Backbone−Solvent Interactions in Determining Conformational Equilibria of Intrinsically Disordered Proteins

NusG-Spt5 Transcription Factors: Universal, Dynamic Modulators of Gene Expression

The accurate and efficient biogenesis of RNA by cellular RNA polymerase (RNAP) requires accessory factors that regulate the initiation, elongation, and termination of transcription. Of the many discovered to date, the elongation regulator NusG-Spt5 is the only universally conserved transcription factor. With orthologs and paralogs found in all three domains of life, this ubiquity underscores their ancient and essential regulatory functions. NusG-Spt5 proteins evolved to maintain a similar binding interface to RNAP through contacts of the NusG N-terminal domain (NGN) that bridge the main DNA-binding cleft. We propose that varying strength of these contacts, modulated by tethering interactions, either decrease transcriptional pausing by smoothing the rugged thermodynamic landscape of transcript elongation or enhance pausing, depending on which conformation of RNAP is stabilized by NGN contacts. NusG-Spt5 contains one (in bacteria and archaea) or more (in eukaryotes) C-terminal domains that use a KOW fold to contact diverse targets, tether the NGN, and control RNA biogenesis. Recent work highlights these diverse functions in different organisms. Some bacteria contain multiple specialized NusG paralogs that regulate subsets of operons via sequence-specific targeting, controlling production of antibiotics, toxins, or capsule proteins. Despite their common origin, NusG orthologs can differ in their target selection, interacting partners, and effects on RNA synthesis. We describe the current understanding of NusG-Spt5 structure, interactions with RNAP and other regulators, and cellular functions including significant recent progress from genome-wide analyses, single-molecule visualization, and cryo-EM. The recent findings highlight the remarkable diversity of function among these structurally conserved proteins.

Solutes unmask differences in clustering versus phase separation of FET proteins

Phase separation and percolation contribute to phase transitions of multivalent macromolecules. Contributions of percolation are evident through the viscoelasticity of condensates and through the formation of heterogeneous distributions of nano- and mesoscale pre-percolation clusters in sub-saturated solutions. Here, we show that clusters formed in sub-saturated solutions of FET (FUS-EWSR1-TAF15) proteins are affected differently by glutamate versus chloride. These differences on the nanoscale, gleaned using a suite of methods deployed across a wide range of protein concentrations, are prevalent and can be unmasked even though the driving forces for phase separation remain unchanged in glutamate versus chloride. Strikingly, differences in anion-mediated interactions that drive clustering saturate on the micron-scale. Beyond this length scale the system separates into coexisting phases. Overall, we find that sequence-encoded interactions, mediated by solution components, make synergistic and distinct contributions to the formation of pre-percolation clusters in sub-saturated solutions, and to the driving forces for phase separation.

Synthetic intrinsically disordered protein fusion tags that enhance protein solubility

We report the de novo design of small (<20 kDa) and highly soluble synthetic intrinsically disordered proteins (SynIDPs) that confer solubility to a fusion partner with minimal effect on the activity of the fused protein. To identify highly soluble SynIDPs, we create a pooled gene-library utilizing a one-pot gene synthesis technology to create a large library of repetitive genes that encode SynIDPs. We identify three small (<20 kDa) and highly soluble SynIDPs from this gene library that lack secondary structure and have high solvation. Recombinant fusion of these SynIDPs to three known inclusion body forming proteins rescue their soluble expression and do not impede the activity of the fusion partner, thereby eliminating the need for removal of the SynIDP tag. These findings highlight the utility of SynIDPs as solubility tags, as they promote the soluble expression of proteins in E. coli and are small, unstructured proteins that minimally interfere with the biological activity of the fused protein.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

A molecular dynamics simulation study of glycine/serine octapeptides labeled with 2,3-diazabicyclo[2.2.2]oct-2-ene fluorophore

The compound 2,3-diazabicyclo[2.2.2]oct-2-ene (DBO) is a versatile fluorophore widely used in Förster resonance energy transfer (FRET) spectroscopy studies due to its remarkable sensitivity, enabling precise donor-acceptor distance measurements, even for short peptides. Integrating time-resolved and FRET spectroscopies with molecular dynamics simulations provides a robust approach to unravel the structure and dynamics of biopolymers in a solution. This study investigates the structural behavior of three octapeptide variants: Trp-(Gly-Ser)3-Dbo, Trp-(GlyGly)3-Dbo, and Trp-(SerSer)3-Dbo, where Dbo represents the DBO-containing modified aspartic acid, using molecular dynamics simulations. Glycine- and serine-rich amino acid fragments, common in flexible protein regions, play essential roles in functional properties. Results show excellent agreement between end-to-end distances, orientational factors from simulations, and the available experimental and theoretical data, validating the reliability of the GROMOS force field model. The end-to-end distribution, modeled using three Gaussian distributions, reveals a complex shape, confirmed by cluster analysis highlighting a limited number of significant conformations dominating the peptide landscape. All peptides predominantly adopt a disordered state in the solvent, yet exhibit a compact shape, aligning with the model of disordered polypeptide chains in poor solvents. Conformations show marginal dependence on chain composition, with Ser-only chains exhibiting slightly more elongation. This study enhances our understanding of peptide behavior, providing valuable insights into their structural dynamics in solution.

Ensemble structure of the N-terminal domain (1-267) of FUS in a biomolecular condensate

Solutions of some proteins phase separate into a condensed state of high protein concentration and a dispersed state of low concentration. Such behavior is observed in living cells for a number of RNA-binding proteins that feature intrinsically disordered domains. It is relevant for cell function via the formation of membraneless organelles and transcriptional condensates. On a basic level, the process can be studied in vitro on protein domains that are necessary and sufficient for liquid-liquid phase separation (LLPS). We have performed distance distribution measurements by electron paramagnetic resonance for 13 sections in an N-terminal domain (NTD) construct of the protein fused in sarcoma (FUS), consisting of the QGSY-rich domain and the RGG1 domain, in the denatured, dispersed, and condensed state. Using 10 distance distribution restraints for ensemble modeling and three such restraints for model validation, we have found that FUS NTD behaves as a random-coil polymer under good-solvent conditions in both the dispersed and condensed state. Conformation distribution in the biomolecular condensate is virtually indistinguishable from the one in an unrestrained ensemble, with the latter one being based on only residue-specific Ramachandran angle distributions. Over its whole length, FUS NTD is slightly more compact in the condensed than in the dispersed state, which is in line with the theory for random coils in good solvent proposed by de Gennes, Daoud, and Jannink. The estimated concentration in the condensate exceeds the overlap concentration resulting from this theory. The QGSY-rich domain is slightly more extended, slightly more hydrated, and has slightly higher propensity for LLPS than the RGG1 domain. Our results support previous suggestions that LLPS of FUS is driven by multiple transient nonspecific hydrogen bonding and π-sp2 interactions.

The food and pharmaceutical additive benzoic acid induces amyloid fibrillation of an intrinsically disordered protein

Benzoic acid (BA) is a microbe-inhibiting flavoring agent used extensively as an additive in foods, pharmaceuticals, and hygiene and cosmetic products. The level of BA in foodstuffs prescribed by world bodies and governmental agencies is assumed to be safe so as to prevent adverse health effects. The safety level of BA is however controversial, and whether different conditions of its use would be generally regarded as safe (GRAS) has been rarely determined. In the quest of how food additives affect the structure and conformation of proteins, this study evaluates the interaction of BA with an intrinsically disordered protein (IDP) at pH 4.2 that matches the pH conditions applicable for the commercial use of benzoate preservatives, and examines its structural transformation by NMR, fluorescence, and high-resolution microscopy. The interaction with BA transforms the protein to a denatured aggregated mesophase that undergoes reconfiguration to yield rigid amyloid fibrils. Significantly, fibrils are observed even with 0.1 mM BA while the recommended level of its use as a preservative is in the 0.4-8 mM range. The discussion refrains from safety comments with no projection of the BA level that could be GRAS.

Backbone interactions and secondary structures in phase separation of disordered proteins

Intrinsically disordered proteins (IDPs) are one of the major drivers behind the formation and characteristics of biomolecular condensates. Due to their inherent flexibility, the backbones of IDPs are significantly exposed, rendering them highly influential and susceptible to biomolecular phase separation. In densely packed condensates, exposed backbones have a heightened capacity to interact with neighboring protein chains, which might lead to strong coupling between the secondary structures and phase separation and further modulate the subsequent transitions of the condensates, such as aging and fibrillization. In this mini-review, we provide an overview of backbone-mediated interactions and secondary structures within biomolecular condensates to underscore the importance of protein backbones in phase separation. We further focus on recent advances in experimental techniques and molecular dynamics simulation methods for probing and exploring the roles of backbone interactions and secondary structures in biomolecular phase separation involving IDPs.

Sequence-Dependent Conformational Transitions of Disordered Proteins During Condensation

Intrinsically disordered proteins (IDPs) can form biomolecular condensates through phase separation. It is recognized that the conformation of IDPs in the dense and dilute phases as well as at the interfaces of condensates can critically impact the resulting properties associated with their functionality. However, a comprehensive understanding of the conformational transitions of IDPs during condensation remains elusive. In this study, we employ a coarse-grained polyampholyte model, comprising an equal number of oppositely charged residues-glutamic acid and lysine-whereby conformations and phase behavior can be readily tuned by altering the protein sequence. By manipulating the sequence patterns from perfectly alternating to block-like, we obtain chains with ideal-like conformations to semi-compact structures in the dilute phase, while in the dense phase, the chain conformation is approximately that of an ideal chain, irrespective of the protein sequence. By performing simulations at different concentrations, we find that the chains assemble from the dilute phase through small oligomeric clusters to the dense phase, accompanied by a gradual swelling of the individual chains. We further demonstrate that these findings are applicable to several naturally occurring proteins involved in the formation of biological condensates. Concurrently, we delve deeper into the chain conformations within the condensate, revealing that chains at the interface show a strong sequence dependence, but remain more collapsed than those in the bulk-like dense phase. This study addresses critical gaps in our knowledge of IDP conformations within condensates as a function of protein sequence.

Coexistence of two coacervate phases of polyglycine in water suggested by polymer reference interaction site model theory

Mixing Gibbs energy and phase equilibria of aqueous solutions of polyglycine were studied theoretically by means of polymer reference interaction site model integral equation theory combined with the Gibbs-Duhem method. In addition to the ordinary liquid-liquid phase separation between dilute and concentrated solutions, the theoretical calculation predicted the coexistence of two coacervate phases, namely, the lower- and higher-density coacervates. The relative thermodynamic stabilities of these two phases change with the polymerization degree of polyglycine. The higher-density coacervate phase was rapidly stabilized by increasing the polymer length, and the lower-density phase became metastable at large polymers. The hydrogen bonds between the peptide chains were strengthened, and water was thermodynamically destabilized in the higher-density coacervate. A possible relation with the formation of amyloid fibril within a liquid droplet is also discussed.

Phase Transitions of Associative Biomacromolecules

Multivalent proteins and nucleic acids, collectively referred to as multivalent associative biomacromolecules, provide the driving forces for the formation and compositional regulation of biomolecular condensates. Here, we review the key concepts of phase transitions of aqueous solutions of associative biomacromolecules, specifically proteins that include folded domains and intrinsically disordered regions. The phase transitions of these systems come under the rubric of coupled associative and segregative transitions. The concepts underlying these processes are presented, and their relevance to biomolecular condensates is discussed.

Solvent quality and solvent polarity in polypeptides

Using molecular dynamics and thermodynamic integration, we report on the solvation process of seven polypeptides (GLY, ALA, ILE, ASN, LYS, ARG, GLU) in water and in cyclohexane. The polypeptides are selected to cover the full hydrophobic scale while varying their chain length from tri- to undeca-homopeptides, providing indications on possible non-additivity effects as well as the role of the peptide backbone in the overall stability of the polypeptides. The use of different solvents and different polypeptides allows us to investigate the relation between solvent quality - the capacity of a given solvent to fold/unfold a given biopolymer often described on a scale ranging from "good" to "poor"; and solvent polarity - related to the specific interactions of any solvent with respect to a reference solvent. Undeca-glycine is found to be the only polypeptide to have a stable collapse in water (polar solvent), with the other hydrophobic polypeptides displaying repeated folding and unfolding events in water, with polar polypeptides presenting even more complex behavior. By contrast, all polypeptides are found to keep an extended conformation in cyclohexane, irrespective of their polarity. All considered polypeptides are also found to have favorable solvation free energy independent of the solvent polarity and their intrinsic hydrophobicity, clearly highlighting the prominent stabilizing role of the peptide backbone - with the solvation process largely enthalpically dominated in polar polypeptides and partially entropically driven for hydrophobic polypeptides. Our study thus reveals the complexity of the solvation process of polypeptides defying the common view "like dissolves like", with the solute polarity playing the most prominent role. The absence of mirror symmetry upon the inversion of polarities of both the solvent and the polypeptides is confirmed.

Effects of Conformational Constraint on Peptide Solubility Limits

Liquid-liquid phase separation of proteins preferentially involves intrinsically disordered proteins or disordered regions. Understanding the solution chemistry of these phase separations is key to learning how to quantify and manipulate systems that involve such processes. Here, we investigate the effect of cyclization on the liquid-liquid phase separation of short polyglycine peptides. We simulated separate aqueous systems of supersaturated cyclic and linear GGGGG and observed spontaneous liquid-liquid phase separation in each of the solutions. The cyclic GGGGG phase separates less robustly than linear GGGGG and has a higher aqueous solubility, even though linear GGGGG has a more favorable single molecule solvation free energy. The versatile and abundant interpeptide contacts formed by the linear GGGGG stabilize the condensed droplet phase, driving the phase separation in this system. In particular, we find that van der Waals close contact interactions are enriched in the droplet phase as opposed to electrostatic interactions. An analysis of the change in backbone conformational entropy that accompanies the phase transition revealed that cyclic peptides lose significantly less entropy in this process as expected. However, we find that the enhanced interaction enthalpy of linear GGGGG in the droplet phase is enough to compensate for a larger decrease in conformational entropy.

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

Combinatorial and antagonistic effects of tubulin glutamylation and glycylation on katanin microtubule severing

Microtubules have spatiotemporally complex posttranslational modification patterns. How cells interpret this tubulin modification code is largely unknown. We show that C. elegans katanin, a microtubule severing AAA ATPase mutated in microcephaly and critical for cell division, axonal elongation, and cilia biogenesis, responds precisely, differentially, and combinatorially to three chemically distinct tubulin modifications-glycylation, glutamylation, and tyrosination-but is insensitive to acetylation. Glutamylation and glycylation are antagonistic rheostats with glycylation protecting microtubules from severing. Katanin exhibits graded and divergent responses to glutamylation on the α- and β-tubulin tails, and these act combinatorially. The katanin hexamer central pore constrains the polyglutamate chain patterns on β-tails recognized productively. Elements distal to the katanin AAA core sense α-tubulin tyrosination, and detyrosination downregulates severing. The multivalent microtubule recognition that enables katanin to read multiple tubulin modification inputs explains in vivo observations and illustrates how effectors can integrate tubulin code signals to produce diverse functional outcomes.

Negative Thermal Expansion and Disorder-to-Order Collapse of an Intrinsically Disordered Protein under Marginally Denaturing Conditions

Recent work with intrinsically disordered proteins (IDPs) has projected a myriad of their survival instincts based mainly on the total charge content, the abundance of polar residues, and the paucity of hydrophobic amino acids. This work uses a plant IDP AtPP16-1 (Arabidopsis thaliana phloem protein class 16-1), whose solution NMR structure was determined by us recently, to show legitimate negative thermal expansion (NTE) of the native state. The thermal expansion continues to be negative even as the tertiary structure is perturbed by ultralow levels of urea up to 0.4 M. The NTE of these subdenatured states is called apparent NTE because they are prone to undergo conformational changes with temperature. Hydrodynamic shrinkage of the NTE IDP is also observed by dynamic light scattering (DLS) and NMR-measured global rotational correlation time (τ_c). The protein with denatured tertiary structure but otherwise preserved native-state secondary structure collapses to a dynamically rigid state. The data are mainly based on thermal coefficients of chemical shift and nuclear relaxation measured by heteronuclear NMR. The hydrodynamic shrinkage and collapse under marginally varying solvent compositions that may arise from unstable tertiary structure and dynamic disorder of chain segments across the backbone could be a generic property of IDPs.

The Molecular and Functional Interaction Between Membrane-Bound Organelles and Membrane-Less Condensates

A major recent advance in cell biology is the mechanistic and kinetic understanding of biogenesis of many membrane-less condensates. As membrane-less condensates and membrane-bound organelles are two major approaches used by the eukaryotic cells to organize cellular contents, it is not surprising that these membrane-less condensates interact with the membrane-bound organelles and are dynamically regulated by the cellular signaling, metabolic states, and proteostasis network. In this review, I will discuss recent progress in the biogenesis of membrane-less condensates and their connections with well-studied membrane-bound organelles. Future work will reveal the molecular and functional connectome among different condensates and membrane-bound organelles.

Deciphering how naturally occurring sequence features impact the phase behaviours of disordered prion-like domains

Prion-like low-complexity domains (PLCDs) have distinctive sequence grammars that determine their driving forces for phase separation. Here we uncover the physicochemical underpinnings of how evolutionarily conserved compositional biases influence the phase behaviour of PLCDs. We interpret our results in the context of the stickers-and-spacers model for the phase separation of associative polymers. We find that tyrosine is a stronger sticker than phenylalanine, whereas arginine is a context-dependent auxiliary sticker. In contrast, lysine weakens sticker-sticker interactions. Increasing the net charge per residue destabilizes phase separation while also weakening the strong coupling between single-chain contraction in dilute phases and multichain interactions that give rise to phase separation. Finally, glycine and serine residues act as non-equivalent spacers, and thus make the glycine versus serine contents an important determinant of the driving forces for phase separation. The totality of our results leads to a set of rules that enable comparative estimates of composition-specific driving forces for PLCD phase separation.

Conformational and Solvation Dynamics of an Amyloidogenic Intrinsically Disordered Domain of a Melanosomal Protein

The conformational plasticity of intrinsically disordered proteins (IDPs) allows them to adopt a range of conformational states that can be important for their biological functions. The driving force for the conformational preference of an IDP emanates from an intricate interplay between chain-chain and chain-solvent interactions. Using ultrafast femtosecond and picosecond time-resolved fluorescence measurements, we characterized the conformational and solvation dynamics around the N- and C-terminal segments of a disordered repeat domain of a melanosomal protein Pmel17 that forms functional amyloid responsible for melanin biosynthesis. Our time-resolved fluorescence anisotropy results revealed slight compaction and slower rotational dynamics around the amyloidogenic C-terminal segment when compared to the proline-rich N-terminal segment of the repeat domain. The compaction of the C-terminal region was also associated with the restrained mobility of hydration water as indicated by our solvation dynamics measurements. Our findings indicate that sequence-dependent chain-solvent interactions govern both the conformational and solvation dynamics that are crucial in directing the conversion of a highly dynamic IDP into an ordered amyloid assembly. Such an interplay of amino acid composition-dependent conformational and solvation dynamics might have important physicochemical consequences in specific water-protein, ion-protein, and protein-protein interactions involved in amyloid formation and phase transitions.

Glycine-Rich Peptides from FUS Have an Intrinsic Ability to Self-Assemble into Fibers and Networked Fibrils

Glycine-rich regions feature prominently in intrinsically disordered regions (IDRs) of proteins that drive phase separation and the regulated formation of membraneless biomolecular condensates. Interestingly, the Gly-rich IDRs seldom feature poly-Gly tracts. The protein fused in sarcoma (FUS) is an exception. This protein includes two 10-residue poly-Gly tracts within the prion-like domain (PLD) and at the interface between the PLD and the RNA binding domain. Poly-Gly tracts are known to be highly insoluble, being potent drivers of self-assembly into solid-like fibrils. Given that the internal concentrations of FUS and FUS-like molecules cross the high micromolar and even millimolar range within condensates, we reasoned that the intrinsic insolubility of poly-Gly tracts might be germane to emergent fluid-to-solid transitions within condensates. To assess this possibility, we characterized the concentration-dependent self-assembly for three non-overlapping 25-residue Gly-rich peptides derived from FUS. Two of the three peptides feature 10-residue poly-Gly tracts. These peptides form either long fibrils based on twisted ribbon-like structures or self-supporting gels based on physical cross-links of fibrils. Conversely, the peptide with similar Gly contents but lacking a poly-Gly tract does not form fibrils or gels. Instead, it remains soluble across a wide range of concentrations. Our findings highlight the ability of poly-Gly tracts within IDRs that drive phase separation to undergo self-assembly. We propose that these tracts are likely to contribute to nucleation of fibrillar solids within dense condensates formed by FUS.

Thermodynamic Compensation in Peptides Following Liquid-Liquid Phase Separation

Liquid-liquid phase separation of proteins often incorporates intrinsically disordered proteins or those with disordered regions. Examining these processes via the entropy change is desirable for establishing a quantitative foundation with which to probe and understand these phase transitions. Of interest is the effect of residue sequence on the entropy of the peptide backbone. In this work we model these systems via all atom simulations of liquid-liquid phase separation of peptides. Systems of supersaturated pentapeptides separate into a peptide-dense liquid droplet phase as well as a dilute (saturated) aqueous phase. An analysis of the change in backbone conformational entropy associated with the phase transition was performed. We examined systems of four different pentapeptides (GGGGG, GGQGG, GGNGG, and GGVGG) in order to explore the effect of sequence variation on the conformational entropy, as well as the effect of side chain variation on the physical characteristics of the droplet phases. We find that the loss of conformational entropy that accompanies aqueous → droplet transitions is more than compensated by a decrease in interaction enthalpy as contributions to the free energy change for the process.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Applications of KBFF20

Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native β-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20

A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.

A Data-Driven Hydrophobicity Scale for Predicting Liquid-Liquid Phase Separation of Proteins

An accurate model for macroscale disordered assemblies of biological macromolecules such as those formed in so-called membraneless organelles would greatly assist in studying their structure, function, and dynamics. Recent evidence has suggested that liquid-liquid phase separation (LLPS) underlies the formation of membraneless organelles. While the general mechanism of exchange of macromolecule/water for macromolecule/macromolecule interactions is known to be the driving force for LLPS, the specific interactions involved are not well understood. One way that protein-water and protein-protein interactions have been understood historically is via hydrophobicity scales. However, these scales are typically optimized for describing these relative interactions in certain cases, such as protein folding or insertion of proteins into membranes. To better describe the relative interactions of proteins that undergo LLPS, we have developed a new, data-driven hydrophobicity scale. To determine the new scale, we used coarse-grained molecular dynamics simulations using the hydrophobicity scale coarse-grained model, which relates the interactions between amino acids to their hydrophobicity. Hydrophobicity values were determined via the force-balance method on a library of proteins that includes unfolded, intrinsically disordered, and phase-separating proteins (PSP). The resulting hydrophobicity scale can better predict whether a given protein will undergo LLPS at physiological conditions by using coarse-grained molecular dynamics simulations than existing hydrophobicity scales. This new scale confirms the importance of π-π interactions between amino acids as important drivers of LLPS. This new hydrophobicity scale provides a convenient and compact description of protein-protein interactions for proteins that undergo LLPS and could be used to develop new models to describe interactions between PSP and other components, such as nucleic acids.

Modulating Microtubules: A Molecular Perspective on the Effects of Tail Modifications

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

Therapeutic peptides for the treatment of cystic fibrosis: Challenges and perspectives

Cystic fibrosis (CF) is the most common amongst rare genetic diseases, affecting more than 70.000 people worldwide. CF is characterized by a dysfunctional chloride channel, termed cystic fibrosis conductance regulator (CFTR), which leads to the production of a thick and viscous mucus layer that clogs the lungs of CF patients and traps pathogens, leading to chronic infections and inflammation and, ultimately, lung damage. In recent years, the use of peptides for the treatment of respiratory diseases, including CF, has gained growing interest. Therapeutic peptides for CF include antimicrobial peptides, inhibitors of proteases, and modulators of ion channels, among others. Peptides display unique features that make them appealing candidates for clinical translation, like specificity of action, high efficacy, and low toxicity. Nevertheless, the intrinsic properties of peptides, together with the need of delivering these compounds locally, e.g. by inhalation, raise a number of concerns in the development of peptide therapeutics for CF lung disease. In this review, we discuss the challenges related to the use of peptides for the treatment of CF lung disease through inhalation, which include retention within mucus, proteolysis, immunogenicity and aggregation. Strategies for overcoming major shortcomings of peptide therapeutics will be presented, together with recent developments in peptide design and optimization, including computational analysis and high-throughput screening.

Evolutionary history of cotranscriptional editing in the paramyxoviral phosphoprotein gene

The phosphoprotein gene of the paramyxoviruses encodes multiple protein products. The P, V, and W proteins are generated by transcriptional slippage. This process results in the insertion of non-templated guanosine nucleosides into the mRNA at a conserved edit site. The P protein is an essential component of the viral RNA polymerase and is encoded by a faithful copy of the gene in the majority of paramyxoviruses. However, in some cases, the non-essential V protein is encoded by default and guanosines must be inserted into the mRNA in order to encode P. The number of guanosines inserted into the P gene can be described by a probability distribution, which varies between viruses. In this article, we review the nature of these distributions, which can be inferred from mRNA sequencing data, and reconstruct the evolutionary history of cotranscriptional editing in the paramyxovirus family. Our model suggests that, throughout known history of the family, the system has switched from a P default to a V default mode four times; complete loss of the editing system has occurred twice, the canonical zinc finger domain of the V protein has been deleted or heavily mutated a further two times, and the W protein has independently evolved a novel function three times. Finally, we review the physical mechanisms of cotranscriptional editing via slippage of the viral RNA polymerase.

C-Terminal Tail Polyglycylation and Polyglutamylation Alter Microtubule Mechanical Properties

Microtubules are biopolymers that perform diverse cellular functions. Microtubule behavior regulation occurs in part through post-translational modification of both the α- and β-subunits of tubulin. One class of modifications is the heterogeneous addition of glycine and/or glutamate residues to the disordered C-terminal tails (CTTs) of tubulin. Because of their prevalence in stable, high-stress cellular structures such as cilia, we sought to determine if these modifications alter microtubules' intrinsic stiffness. Here, we describe the purification and characterization of differentially modified pools of tubulin from Tetrahymena thermophila. We found that post-translational modifications do affect microtubule stiffness but do not affect the number of protofilaments incorporated into microtubules. We measured the spin dynamics of nuclei in the CTT backbone by NMR spectroscopy to explore the mechanism of this change. Our results show that the α-tubulin CTT does not protrude out from the microtubule surface, as is commonly depicted in models, but instead interacts with the dimer's surface. This suggests that the interactions of the α-tubulin CTT with the tubulin body contributes to the stiffness of the assembled microtubule, thus providing insight into the mechanism by which polyglycylation and polyglutamylation can alter microtubule mechanical properties.

The Dynamism of Intrinsically Disordered Proteins: Binding-Induced Folding, Amyloid Formation, and Phase Separation

Intrinsically disordered proteins (IDPs) or natively unfolded proteins do not undergo autonomous folding into a well-defined 3-D structure and challenge the conventional structure-function paradigm. They are involved in a multitude of critical physiological functions by adopting various structural states via order-to-disorder transitions or by maintaining their disordered characteristics in functional complexes. In recent times, there has been a burgeoning interest in the investigation of intriguing behavior of IDPs using highly multidisciplinary and complementary approaches due to the pivotal role of this unique class of protein chameleons in physiology and disease. Over the past decade or so, our laboratory has been actively investigating the unique physicochemical properties of this class of highly dynamic, flexible, rapidly interconverting proteins. We have utilized a diverse array of existing and emerging tools involving steady-state and time-resolved fluorescence, Raman spectroscopy, circular dichroism, light scattering, fluorescence microscopy, and atomic force microscopy coupled with site-directed mutagenesis and other biochemical and biophysical tools to study a variety of interesting and important aspects of IDPs. In this Feature Article, I describe our work on the conformational characteristics, solvation dynamics, binding-induced folding, amyloid formation, and liquid-liquid phase separation of a number of amyloidogenic IDPs. A series of these studies described here captures the role of conformational plasticity and dynamics in directing binding, folding, assembly, aggregation, and phase transitions implicated in physiology and pathology.

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

Physical Chemistry of the Protein Backbone: Enabling the Mechanisms of Intrinsic Protein Disorder

Over the last two decades it has become clear that well-defined structure is not a requisite for proteins to properly function. Rather, spectra of functionally competent, structurally disordered states have been uncovered requiring canonical paradigms in molecular biology to be revisited or reimagined. It is enticing and oftentimes practical to divide the proteome into structured and unstructured, or disordered, proteins. While function, composition, and structural properties largely differ, these two classes of protein are built upon the same scaffold, namely, the protein backbone. The versatile physicochemical properties of the protein backbone must accommodate structural disorder, order, and transitions between these states. In this review, we survey these properties through the conceptual lenses of solubility and conformational populations and in the context of protein-disorder mediated phenomena (e.g., phase separation, order-disorder transitions, allostery). Particular attention is paid to the results of computational studies, which, through thermodynamic decomposition and dissection of molecular interactions, can provide valuable mechanistic insight and testable hypotheses to guide further solution experiments. Lastly, we discuss changes in the dynamics of side chains and order-disorder transitions of the protein backbone as two modes or realizations of "entropic reservoirs" capable of tuning coupled thermodynamic processes.

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.

Dynamic elements and kinetics: Most favorable conformations of peptides in solution with measurements and simulations

Small peptides in solution adopt a specific morphology as they function. It is of fundamental interest to examine the structural properties of these small biomolecules in solution and observe how they transition from one conformation to another and form functional structures. In this study, we have examined the structural properties of a simple dipeptide and a five-residue peptide with the application of far-UV circular dichroism (CD) spectroscopy as a function of temperature, fluorescence anisotropy, and all-atom molecular dynamics simulation. Analysis of the temperature dependent CD spectra shows that the simplest dipeptide N-acetyl-tryptophan-amide (NATA) adopts helical, beta sheet, and random coil conformations. At room temperature, NATA is found to have 5% alpha-helical, 37% beta sheet, and 58% random coil conformations. To our knowledge, this type of structural content in a simplest dipeptide has not been observed earlier. The pentapeptide (WK5) is found to have four major secondary structural elements with 8% 3₁₀ helix, 14% poly-L-proline II, 8% beta sheet, and 14% turns. A 56% unordered structural population is also present for WK5. The presence of a significant population of 3₁₀ helix in a simple pentapeptide is rarely observed. Fluorescence anisotropy decay (FAD) measurements yielded reorientation times of 45 ps for NATA and 120 ps for WK5. The fluorescence anisotropy decay measurements reveal the size differences between the two peptides, NATA and WK5, with possible contributions from differences in shape, interactions with the environment, and conformational dynamics. All-atom molecular dynamics simulations were used to model the structures and motions of these two systems in solution. The predicted structures sampled by both peptides qualitatively agree with the experimental findings. Kinetic modeling with optimal dimensionality reduction suggests that the slowest dynamic processes in the dipeptide involve sidechain transitions occurring on a 1 ns timescale. The kinetics in the pentapeptide monitors the formation of a distorted helical structure from an extended conformation on a timescale of 10 ns. Modeling of the fluorescence anisotropy decay is found to be in good agreement with the measured data and correlates with the main contributions of the measured reorientation times to individual conformers, which we define as dynamic elements. In NATA, the FAD can be well represented as a sum of contributions from representative conformers. This is not the case in WK5, where our analysis suggests the existence of coupling between conformational dynamics and global tumbling. The current study involving detailed experimental measurements and atomically detailed modeling reveals the existence of specific secondary structural elements and novel dynamical features even in the simplest peptide systems.

Dissecting the Energetics of Intrinsically Disordered Proteins via a Hybrid Experimental and Computational Approach

Intrinsically disordered proteins (IDPs) play important roles in biology, but little is known about the energetics of their inter-residue interactions. Methods that have been successfully applied to analyze the energetics of globular proteins are not applicable to the fluctuating partially ordered ensembles populated by IDPs. A combined computational experimental strategy is introduced for analyzing the energetic role of individual residues in the free state of IDPs. The approach combines experimental measurements of the binding of wild-type and mutant IDPs to their partners with alchemical free energy calculations of the structured complexes. These data allow quantitative information to be deduced about the free state via a thermodynamic cycle. The approach is validated by the analysis of the effects of mutations upon the binding free energy of the ovomucoid inhibitor third binding domain to its partners and is applied to the C-terminal domain of the measles virus nucleoprotein, a 125-residue IDP involved in the RNA transcription and replication of measles virus. The analysis reveals significant inter-residue interactions in the unbound IDP and suggests a biological role for them. The work demonstrates that advances in force fields and computational hardware have now led to the point where it is possible to develop methods, which integrate experimental and computational techniques to reveal insights that cannot be studied using either technique alone.

Thermodynamics of Conformational Transitions in a Disordered Protein Backbone Model

Conformational entropy is expected to contribute significantly to the thermodynamics of structural transitions in intrinsically disordered proteins or regions in response to protein/ligand binding, posttranslational modifications, and environmental changes. We calculated the backbone (dihedral) conformational entropy of oligoglycine (Gly_N), a protein backbone mimic and model intrinsically disordered region, as a function of chain length (N=3, 4, 5, 10, and 15) from simulations using three different approaches. The backbone conformational entropy scales linearly with chain length with a slope consistent with the entropy of folding of well-structured proteins. The entropic contributions of second-order dihedral correlations are predominantly through intraresidue ϕ-ψ pairs, suggesting that oligoglycine may be thermodynamically modeled as a system of independent glycine residues. We find the backbone conformational entropy to be largely independent of global structural parameters, like the end-to-end distance and radius of gyration. We introduce a framework referred to herein as "ensemble confinement" to estimate the loss (gain) of conformational free energy and its entropic component when individual residues are constrained to (released from) particular regions of the ϕ-ψ map. Quantitatively, we show that our protein backbone model resists ordering/folding with a significant, unfavorable ensemble confinement free energy because of the loss of a substantial portion of the absolute backbone entropy. Proteins can couple this free-energy reservoir to distal binding events as a regulatory mechanism to promote or suppress binding.

Pearl-Necklace-Like Local Ordering Drives Polypeptide Collapse

The collapse of the polypeptide backbone is an integral part of protein folding. Using polyglycine as a probe, we explore the nonequilibrium pathways of protein collapse in water. We find that the collapse depends on the competition between hydration effects and intrapeptide interactions. Once intrapeptide van der Waal interactions dominate, the chain collapses along a nonequilibrium pathway characterized by formation of pearl-necklace-like local clusters as intermediates that eventually coagulate into a single globule. By describing this coarsening through the contact probability as a function of distance along the chain, we extract a time-dependent length scale that grows in a linear fashion. The collapse dynamics is characterized by a dynamical critical exponent z ≈ 0.5 that is much smaller than the values of z = 1–2 reported for nonbiological polymers. This difference in the exponents is explained by the instantaneous formation of intrachain hydrogen bonds and local ordering that may be correlated with the observed fast folding times of proteins.

Protein plasticity driven by disorder and collapse governs the heterogeneous binding of CytR to DNA

The amplitude of thermodynamic fluctuations in biological macromolecules determines their conformational behavior, dimensions, nature of phase transitions and effectively their specificity and affinity, thus contributing to fine-tuned molecular recognition. Unique among large-scale conformational changes in proteins are temperature-induced collapse transitions in intrinsically disordered proteins (IDPs). Here, we show that CytR DNA-binding domain, an IDP that folds on binding DNA, undergoes a coil-to-globule transition with temperature in the absence of DNA while exhibiting energetically decoupled local and global structural rearrangements, and maximal thermodynamic fluctuations at the optimal bacterial growth temperature. The collapse is shown to be a continuous transition through a combination of statistical-mechanical modeling and all-atom implicit solvent simulations. Surprisingly, CytR binds single-site cognate DNA with negative cooperativity, described by Hill coefficients less than one, resulting in a graded binding response. We show that heterogeneity arising from varying binding-competent CytR conformations or orientations at the single-molecular level contributes to negative binding cooperativity at the level of bulk measurements due to the conflicting requirements of collapse transition, large fluctuations and folding-upon-binding. Our work reports strong evidence for functionally driven thermodynamic fluctuations in determining the extent of collapse and disorder with implications in protein search efficiency of target DNA sites and regulation.

Femtosecond Hydration Map of Intrinsically Disordered α-Synuclein

Protein hydration water plays a fundamentally important role in protein folding, binding, assembly, and function. Little is known about the hydration water in intrinsically disordered proteins that challenge the conventional sequence-structure-function paradigm. Here, by combining experiments and simulations, we show the existence of dynamical heterogeneity of hydration water in an intrinsically disordered presynaptic protein, namely α-synuclein, implicated in Parkinson's disease. We took advantage of nonoccurrence of cysteine in the sequence and incorporated a number of cysteine residues at the N-terminal segment, the central amyloidogenic nonamyloid-β component (NAC) domain, and the C-terminal end of α-synuclein. We then labeled these cysteine variants using environment-sensitive thiol-active fluorophore and monitored the solvation dynamics using femtosecond time-resolved fluorescence. The site-specific femtosecond time-resolved experiments allowed us to construct the hydration map of α-synuclein. Our results show the presence of three dynamically distinct types of water: bulk, hydration, and confined water. The amyloidogenic NAC domain contains dynamically restrained water molecules that are strikingly different from the water molecules present in the other two domains. Atomistic molecular dynamics simulations revealed longer residence times for water molecules near the NAC domain and supported our experimental observations. Additionally, our simulations allowed us to decipher the molecular origin of the dynamical heterogeneity of water in α-synuclein. These simulations captured the quasi-bound water molecules within the NAC domain originating from a complex interplay between the local chain compaction and the sequence composition. Our findings from this synergistic experimental simulation approach suggest longer trapping of interfacial water molecules near the amyloidogenic hotspot that triggers the pathological conversion into amyloids via chain sequestration, chain desolvation, and entropic liberation of ordered water molecules.

Unfolded states under folding conditions accommodate sequence-specific conformational preferences with random coil-like dimensions

Proteins are marginally stable molecules that fluctuate between folded and unfolded states. Here, we provide a high-resolution description of unfolded states under refolding conditions for the N-terminal domain of the L9 protein (NTL9). We use a combination of time-resolved Förster resonance energy transfer (FRET) based on multiple pairs of minimally perturbing labels, time-resolved small-angle X-ray scattering (SAXS), all-atom simulations, and polymer theory. Upon dilution from high denaturant, the unfolded state undergoes rapid contraction. Although this contraction occurs before the folding transition, the unfolded state remains considerably more expanded than the folded state and accommodates a range of local and nonlocal contacts, including secondary structures and native and nonnative interactions. Paradoxically, despite discernible sequence-specific conformational preferences, the ensemble-averaged properties of unfolded states are consistent with those of canonical random coils, namely polymers in indifferent (theta) solvents. These findings are concordant with theoretical predictions based on coarse-grained models and inferences drawn from single-molecule experiments regarding the sequence-specific scaling behavior of unfolded proteins under folding conditions.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Improvements to the ABSINTH Force Field for Proteins Based on Experimentally Derived Amino Acid Specific Backbone Conformational Statistics

We present an improved version of the ABSINTH implicit solvation model and force field paradigm (termed ABSINTH-C) by incorporating a grid-based term that bootstraps against experimentally derived and computationally optimized conformational statistics for blocked amino acids. These statistics provide high-resolution descriptions of the intrinsic backbone dihedral angle preferences for all 20 amino acids. The original ABSINTH model generates Ramachandran plots that are too shallow in terms of the basin structures and too permissive in terms of dihedral angle preferences. We bootstrap against the reference optimized landscapes and incorporate CMAP-like residue-specific terms that help us reproduce the intrinsic dihedral angle preferences of individual amino acids. These corrections that lead to ABSINTH-C are achieved by balancing the incorporation of the new residue-specific terms with the accuracies inherent to the original ABSINTH model. We demonstrate the robustness of ABSINTH-C through a series of examples to highlight the preservation of accuracies as well as examples that demonstrate the improvements. Our efforts show how the recent experimentally derived and computationally optimized coil-library landscapes can be used as a touchstone for quantifying errors and making improvements to molecular mechanics force fields.

Force Fields for Small Molecules

Molecular dynamics (MD) simulations have been widely applied to computer-aided drug design (CADD). While MD has been used in a variety of applications such as free energy perturbation and long-time simulations, the accuracy of the results from those methods depends strongly on the force field used. Force fields for small molecules are crucial, as they not only serve as building blocks for developing force fields for larger biomolecules but also act as model compounds that will be transferred to ligands used in CADD. Currently, a wide range of small molecule force fields based on additive or nonpolarizable models have been developed. While these nonpolarizable force fields can produce reasonable estimations of physical properties and have shown success in a variety of systems, there is still room for improvements due to inherent limitations in these models including the lack of an electronic polarization response. For this reason, incorporating polarization effects into the energy function underlying a force field is believed to be an important step forward, giving rise to the development of polarizable force fields. Recent simulations of biological systems have indicated that polarizable force fields are able to provide a better physical representation of intermolecular interactions and, in many cases, better agreement with experimental properties than nonpolarizable, additive force fields. Therefore, this chapter focuses on the development of small molecule force fields with emphasis on polarizable models. It begins with a brief introduction on the importance of small molecule force fields and their evolution from additive to polarizable force fields. Emphasis is placed on the additive CHARMM General Force Field and the polarizable force field based on the classical Drude oscillator. The theory for the Drude polarizable force field and results for small molecules are presented showing their improvements over the additive model. The potential importance of polarization for their application in a wide range of biological systems including CADD is then discussed.

Folding of poly-amino acids and intrinsically disordered proteins in overcrowded milieu induced by pH change

pH-induced structural changes of the synthetic homopolypeptides poly-E, poly-K, poly-R, and intrinsically disordered proteins (IDPs) prothymosin α (ProTα) and linker histone H1, in concentrated PEG solutions simulating macromolecular crowding conditions within the membrane-less organelles, were characterized. The conformational transitions of the studied poly-amino acids in the concentrated PEG solutions depend on the polymerization degree of these homopolypeptides, the size of their side chains, the charge distribution of the side chains, and the crowding agent concentration. The results obtained for poly-amino acids are valid for IDPs having a significant total charge. The overcrowded conditions promote a significant increase in the cooperativity of the pH-induced coil-α-helix transition of ProTα and provoke histone H1 aggregation. The most favorable conditions for the pH-induced structural transitions in concentrated PEG solutions are realized when the charged residues are grouped in blocks, and when the distance between the end of the side group carrying charge and the backbone is small. Therefore, the block-wise distribution of charged residues within the IDPs not only plays an important role in the liquid-liquid phase transitions, but may also define the expressivity of structural transitions of these proteins in the overcrowded conditions of the membrane-less organelles.

Building, Characterization, and Applications of Cuvette-FCS in Denaturant-Induced Expansion of Globular and Disordered Proteins

Fluorescence correlation spectroscopy (FCS) is a single-molecule sensitive technique with widespread applications in biophysics. However, conventional microscope-based FCS setups have limitations in performing certain experiments such as those requiring agitations such as stirring or heating, and those involving measurements in solvents with the mismatch of refractive indices. We have recently developed an FCS setup that is suitable for performing measurements inside regular cuvettes. The cuvette-FCS is suitable for performing single-molecule measurements in experiments that are regularly performed in spectrofluorometers but are generally avoided in microscope-based FCS. Here we describe building and characterization of the performance of the cuvette-FCS setup in detail. Finally, we have used a natively folded protein and an intrinsically disordered protein to demonstrate and describe how cuvette-FCS can be applied conveniently to measure urea-dependent expansion of hydrodynamic size of proteins.

Conformational preferences and phase behavior of intrinsically disordered low complexity sequences: insights from multiscale simulations

While many proteins and protein regions utilize a complex repertoire of amino acids to achieve their biological function, a subset of protein sequences are enriched in a reduced set of amino acids. These so-called low complexity (LC) sequences, specifically intrinsically disordered variants of LC sequences, have been the focus of recent investigations owing to their roles in a range of biological functions, specifically phase separation. Computational studies of LC sequences have provided rich insights into their behavior both as individual proteins in dilute solutions and as the drivers and modulators of higher-order assemblies. Here, we review how simulations performed across distinct resolutions have provided different types of insights into the biological role of LC sequences.

Quantitative Binding Behavior of Intrinsically Disordered Proteins to Nanoparticle Surfaces at Individual Residue Level

The quantitative and predictive understanding how intrinsically disordered proteins (IDPs) interact with engineered nanoparticles has potentially important implications for new therapeutics as well as nanotoxicology. Based on a recently developed solution ¹⁵ N NMR relaxation approach, the interactions between four representative IDPs with silica nanoparticles are reported at atomic detail. Each IDP possesses distinct binding modes, which can be quantitatively explained by the local amino-acid residue composition using a "free residue interaction model". The model was parameterized using the binding affinities of free proteinogenic amino acids along with long-range effects, derived by site-specific mutagenesis, that exponentially scale with distance along the primary sequence. The model, which is accessible through a web server, can be applied to predict the residue-specific binding affinities of a large number of IDPs.

Deprotonation of a Single Amino Acid Residue Induces Significant Stability in an α-Helical Heteropeptide

Stability of secondary structural elements is an integral component of a structurally stable protein. Presence of protons in the residue sequence and their immediate environment play a significant role in conformational stability. In this study, we show that removing a proton from a single amino acid residue significantly increases the stability of an α-helical heteropeptide in comparison with the unprotonated form. Far-UV circular dichroism spectroscopy, fluorescence spectroscopy, fluorescence energy transfer measurements, and over 10 μs of all-atom molecular dynamics simulations are used to provide an atomically detailed characterization of this event. There is a single histidine residue in the studied α-helical peptide sequence toward the N-terminal that interacts with a tryptophan located four residues away and quenches the fluorescence when protonated. Removing a proton from this histidine residue dequenches the tryptophan fluorescence and contributes to a significant increase in the helix stability. Atomically detailed analysis of individual residue conformations shows that the protonated histidine tends to be in closer proximity to the tryptophan, which correlates with higher helix content in the N and C termini and lower helix content in the central region of the peptide. In the presence of a neutral histidine, when tryptophan fluorescence is no longer quenched and histidine moves further away from tryptophan, the helix content remains mostly unchanged in the N-and-C termini and significantly increases in the central region. Our results strongly suggest that interactions of the tryptophan with a protonated histidine downregulate the helix population in the central segment of the helical structure compared to a neutral histidine residue. Upregulation of helix population of the central segment of this α-helical heteropeptide in the presence of a neutral histidine residue significantly increases the peptide structural stability.

Paradoxical Effect of Trehalose on the Aggregation of α-Synuclein: Expedites Onset of Aggregation yet Reduces Fibril Load

Aggregation of α-synuclein is closely connected to the pathology of Parkinson's disease. The phenomenon involves multiple steps, commenced by partial misfolding and eventually leading to mature amyloid fibril formation. Trehalose, a widely accepted osmolyte, has been shown previously to inhibit aggregation of various globular proteins owing to its ability to prevent the initial unfolding of protein. In this study, we have examined if it behaves in a similar fashion with intrinsically disordered protein α-synuclein and possesses the potential to act as therapeutic agent against Parkinson's disease. It was observed experimentally that samples coincubated with trehalose fibrillate faster compared to the case in its absence. Molecular dynamics simulations suggested that this initial acceleration is manifestation of trehalose's tendency to perturb the conformational transitions between different conformers of monomeric protein. It stabilizes the aggregation prone "extended" conformer of α-synuclein, by binding to its exposed acidic residues of the C terminus. It also favors the β-rich oligomers once formed. Interestingly, the total fibrils formed are still promisingly less since it accelerates the competing pathway toward formation of amorphous aggregates.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.