Benchmarking all-atom simulations using hydrogen exchange

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.

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Worth the weight: Sub-Pocket EXplorer (SubPEx), a weighted-ensemble method to enhance binding-pocket conformational sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Computational Modeling of Molecular Structures Guided by Hydrogen-Exchange Data

Data produced by hydrogen-exchange monitoring experiments have been used in structural studies of molecules for several decades. Despite uncertainties about the structural determinants of hydrogen exchange itself, such data have successfully helped guide the structural modeling of challenging molecular systems, such as membrane proteins or large macromolecular complexes. As hydrogen-exchange monitoring provides information on the dynamics of molecules in solution, it can complement other experimental techniques in so-called integrative modeling approaches. However, hydrogen-exchange data have often only been used to qualitatively assess molecular structures produced by computational modeling tools. In this paper, we look beyond qualitative approaches and survey the various paradigms under which hydrogen-exchange data have been used to quantitatively guide the computational modeling of molecular structures. Although numerous prediction models have been proposed to link molecular structure and hydrogen exchange, none of them has been widely accepted by the structural biology community. Here, we present as many hydrogen-exchange prediction models as we could find in the literature, with the aim of providing the first exhaustive list of its kind. From purely structure-based models to so-called fractional-population models or knowledge-based models, the field is quite vast. We aspire for this paper to become a resource for practitioners to gain a broader perspective on the field and guide research toward the definition of better prediction models. This will eventually improve synergies between hydrogen-exchange monitoring and molecular modeling.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Hydrogen/Deuterium Exchange Measurements May Provide an Incomplete View of Protein Dynamics: a Case Study on Cytochrome c

Many aspects of protein function rely on conformational fluctuations. Hydrogen/deuterium exchange (HDX) mass spectrometry (MS) provides a window into these dynamics. Despite the widespread use of HDX-MS, it remains unclear whether this technique provides a truly comprehensive view of protein dynamics. HDX is mediated by H-bond-opening/closing events, implying that HDX methods provide an H-bond-centric view. This raises the question if there could be fluctuations that leave the H-bond network unaffected, thereby rendering them undetectable by HDX-MS. We explore this issue in experiments on cytochrome c (cyt c). Compared to the Fe(II) protein, Fe(III) cyt c shows enhanced deuteration on both the distal and proximal sides of the heme. Previous studies have attributed the enhanced dynamics of Fe(III) cyt c to the facile and reversible rupture of the distal M80-Fe(III) bond. Using molecular dynamics (MD) simulations, we conducted a detailed analysis of various cyt c conformers. Our MD data confirm that rupture of the M80-Fe(III) contact triggers major reorientation of the distal Ω loop. Surprisingly, this event takes place with only miniscule H-bonding alterations. In other words, the distal loop dynamics are almost "HDX-silent". Moreover, distal loop movements cannot account for enhanced dynamics on the opposite (proximal) side of the heme. Instead, enhanced deuteration of Fe(III) cyt c is attributed to sparsely populated conformers where both the distal (M80) and proximal (H18) coordination bonds have been ruptured, along with opening of numerous H-bonds on both sides of the heme. We conclude that there can be major structural fluctuations that are only weakly coupled to changes in H-bonding, making them virtually impossible to track by HDX-MS. In such cases, HDX-MS may provide an incomplete view of protein dynamics.

Protein folding - seeing is deceiving

This Perspective is intended to raise questions about the conventional interpretation of protein folding. According to the conventional interpretation, developed over many decades, a protein population can visit a vast number of conformations under unfolding conditions, but a single dominant native population emerges under folding conditions. Accordingly, folding comes with a substantial loss of conformational entropy. How is this price paid? The conventional answer is that favorable interactions between and among the side chains can compensate for entropy loss, and moreover, these interactions are responsible for the structural particulars of the native conformation. Challenging this interpretation, the Perspective introduces a proposal that high energy (i.e., unfavorable) excluding interactions winnow the accessible population substantially under physical-chemical conditions that favor folding. Both steric clash and unsatisfied hydrogen bond donors and acceptors are classified as excluding interactions, so called because conformers with such disfavored interactions will be largely excluded from the thermodynamic population. Both excluding interactions and solvent factors that induce compactness are somewhat nonspecific, yet together they promote substantial chain organization. Moreover, proteins are built on a backbone scaffold consisting of α-helices and strands of β-sheet, where the number of hydrogen bond donors and acceptors is exactly balanced. These repetitive secondary structural elements are the only two conformers that can be both completely hydrogen-bond satisfied and extended indefinitely without encountering a steric clash. Consequently, the number of fundamental folds is limited to no more than ~10,000 for a protein domain. Once excluding interactions are taken into account, the issue of "frustration" is largely eliminated and the Levinthal paradox is resolved. Putting the "bottom line" at the top: it is likely that hydrogen-bond satisfaction represents a largely under-appreciated parameter in protein folding models.

Interdomain flexibility and interfacial integrity of β-lactamase inhibitory protein (BLIP) modulate its binding to class A β-lactamases

β-Lactamase inhibitory protein (BLIP) consists of a tandem repeat of αβ domains conjugated by an interdomain loop and can effectively bind and inactivate class A β-lactamases, which are responsible for resistance of bacteria to β-lactam antibiotics. The varied ability of BLIP to bind different β-lactamases and the structural determinants for significant enhancement of BLIP variants with a point mutation are poorly understood. Here, we investigated the conformational dynamics of BLIP upon binding to three clinically prevalent class A β-lactamases (TEM1, SHV1, and PC1) with dissociation constants between subnanomolar and micromolar. Hydrogen deuterium exchange mass spectrometry revealed that the flexibility of the interdomain region was significantly suppressed upon strong binding to TEM1, but was not significantly changed upon weak binding to SHV1 or PC1. E73M and K74G mutations in the interdomain region improved binding affinity toward SHV1 and PC1, respectively, showing significantly increased flexibility of the interdomain region compared to the wild-type and favorable conformational changes upon binding. In contrast, more rigidity of the interfacial loop 135–145 was observed in these BLIP mutants in both free and bound states. Consistently, molecular dynamics simulations of BLIP exhibited drastic changes in the flexibility of the loop 135–145 in all complexes. Our results indicated for the first time that higher flexibility of the interdomain linker, as well as more rigidity of the interfacial loop 135–145, could be desirable determinants for enhancing inhibition of BLIP to class A β-lactamases. Together, these findings provide unique insights into the design of enhanced inhibitors.

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.

Conformational Preferences of an Intrinsically Disordered Protein Domain: A Case Study for Modern Force Fields

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/E_b, ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

Conformational switch that induces GDP release from Gi

Heterotrimeric guanine nucleotide-binding proteins (G proteins) are composed of α, β, and γ subunits. Gα switches between guanosine diphosphate (GDP)-bound inactive and guanosine triphosphate (GTP)-bound active states, and Gβγ interacts with the GDP-bound state. The GDP-binding regions are composed of two sites: the phosphate-binding and guanine-binding regions. The turnover of GDP and GTP is induced by guanine nucleotide-exchange factors (GEFs), including G protein-coupled receptors (GPCRs), Ric8A, and GIV/Girdin. However, the key structural factors for stabilizing the GDP-bound state of G proteins and the direct structural event for GDP release remain unclear. In this study, we investigated structural factors affecting GDP release by introducing point mutations in selected, conserved residues in Gαi3. We examined the effects of these mutations on the GDP/GTP turnover rate and the overall conformation of Gαi3 as well as the binding free energy between Gαi3 and GDP. We found that dynamic changes in the phosphate-binding regions are an immediate factor for the release of GDP.

Dynamics of an LPS translocon induced by substrate and an antimicrobial peptide

Lipopolysaccharide (LPS) transport to the outer membrane (OM) is a crucial step in the biogenesis of microbial surface defenses. Although many features of the translocation mechanism have been elucidated, molecular details of LPS insertion via the LPS transport (Lpt) OM protein LptDE remain elusive. Here, we integrate native MS with hydrogen-deuterium exchange MS and molecular dynamics simulations to investigate the influence of substrate and peptide binding on the conformational dynamics of LptDE. Our data reveal that LPS induces opening of the LptD β-taco domain, coupled with conformational changes on β-strands adjacent to the putative lateral exit gate. Conversely, an antimicrobial peptide, thanatin, stabilizes the β-taco, thereby preventing LPS transport. Our results illustrate that LPS insertion into the OM relies on concerted opening movements of both the β-barrel and β-taco domains of LptD, and suggest a means for developing antimicrobial therapeutics targeting this essential process in Gram-negative ESKAPE pathogens.

Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different

The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.

Integrating hydrogen-deuterium exchange mass spectrometry with molecular dynamics simulations to probe lipid-modulated conformational changes in membrane proteins

Biological membranes define the boundaries of cells and are composed primarily of phospholipids and membrane proteins. It has become increasingly evident that direct interactions of membrane proteins with their surrounding lipids play key roles in regulating both protein conformations and function. However, the exact nature and structural consequences of these interactions remain difficult to track at the molecular level. Here, we present a protocol that specifically addresses this challenge. First, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of membrane proteins incorporated into nanodiscs of controlled lipid composition is used to obtain information on the lipid species that are involved in modulating the conformational changes in the membrane protein. Then molecular dynamics (MD) simulations in lipid bilayers are used to pinpoint likely lipid-protein interactions, which can be tested experimentally using HDX-MS. By bringing together the MD predictions with the conformational readouts from HDX-MS, we have uncovered key lipid-protein interactions implicated in stabilizing important functional conformations. This protocol can be applied to virtually any integral membrane protein amenable to classic biophysical studies and for which a near-atomic-resolution structure or homology model is available. This protocol takes ~4 d to complete, excluding the time for data analysis and MD simulations, which depends on the size of the protein under investigation.

Regulation of Co-transcriptional Pre-mRNA Splicing by m6A through the Low-Complexity Protein hnRNPG

N⁶-methyladenosine (m⁶A) modification occurs co-transcriptionally and impacts pre-mRNA processing; however, the mechanism of co-transcriptional m⁶A-dependent alternative splicing regulation is still poorly understood. Heterogeneous nuclear ribonucleoprotein G (hnRNPG) is an m⁶A reader protein that binds RNA through RRM and Arg-Gly-Gly (RGG) motifs. Here, we show that hnRNPG directly binds to the phosphorylated carboxy-terminal domain (CTD) of RNA polymerase II (RNAPII) using RGG motifs in its low-complexity region. Through interactions with the phosphorylated CTD and nascent RNA, hnRNPG associates co-transcriptionally with RNAPII and regulates alternative splicing transcriptome-wide. m⁶A near splice sites in nascent pre-mRNA modulates hnRNPG binding, which influences RNAPII occupancy patterns and promotes exon inclusion. Our results reveal an integrated mechanism of co-transcriptional m⁶A-mediated splicing regulation, in which an m⁶A reader protein uses RGG motifs to co-transcriptionally interact with both RNAPII and m⁶A-modified nascent pre-mRNA to modulate RNAPII occupancy and alternative splicing.

Using Small-Angle Scattering Data and Parametric Machine Learning to Optimize Force Field Parameters for Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) play important roles in many aspects of normal cell physiology, such as signal transduction and transcription, as well as pathological states, including Alzheimer's, Parkinson's, and Huntington's disease. Unlike their globular counterparts that are defined by a few structures and free energy minima, IDP/IDR comprise a large ensemble of rapidly interconverting structures and a corresponding free energy landscape characterized by multiple minima. This aspect has precluded the use of structural biological techniques, such as X-ray crystallography and nuclear magnetic resonance (NMR) for resolving their structures. Instead, low-resolution techniques, such as small-angle X-ray or neutron scattering (SAXS/SANS), have become a mainstay in characterizing coarse features of the ensemble of structures. These are typically complemented with NMR data if possible or computational techniques, such as atomistic molecular dynamics, to further resolve the underlying ensemble of structures. However, over the past 10–15 years, it has become evident that the classical, pairwise-additive force fields that have enjoyed a high degree of success for globular proteins have been somewhat limited in modeling IDP/IDR structures that agree with experiment. There has thus been a significant effort to rehabilitate these models to obtain better agreement with experiment, typically done by optimizing parameters in a piecewise fashion. In this work, we take a different approach by optimizing a set of force field parameters simultaneously, using machine learning to adapt force field parameters to experimental SAXS scattering profiles. We demonstrate our approach in modeling three biologically IDP ensembles based on experimental SAXS profiles and show that our optimization approach significantly improve force field parameters that generate ensembles in better agreement with experiment.

Protein dynamics revealed by hydrogen/deuterium exchange mass spectrometry: Correlation between experiments and simulation

Hydrogen deuterium exchange mass spectrometry (HDX-MS) is a powerful technique for studying protein dynamics, which is an important factor governing protein functions. However, the process of hydrogen/deuterium exchange (HDX) of proteins is highly complex and the underlying mechanism has not yet been fully elucidated. Meanwhile, molecular dynamics (MD) simulation is a computational technique that can be used to elucidate HDX behaviour on proteins and facilitate interpretation of HDX-MS data. This article aims to summarize the current understandings on the mechanism of HDX and its correlation with MD simulation, to discuss the recent developments in the techniques of HDX-MS and MD simulation and to extend the perspectives of these two techniques in protein dynamics study.

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.

Commonly used FRET fluorophores promote collapse of an otherwise disordered protein

The dimensions that unfolded proteins, including intrinsically disordered proteins (IDPs), adopt in the absence of denaturant remain controversial. We developed an analysis procedure for small-angle X-ray scattering (SAXS) profiles and used it to demonstrate that even relatively hydrophobic IDPs remain nearly as expanded in water as they are in high denaturant concentrations. In contrast, as demonstrated here, most fluorescence resonance energy transfer (FRET) measurements have indicated that relatively hydrophobic IDPs contract significantly in the absence of denaturant. We use two independent approaches to further explore this controversy. First, using SAXS we show that fluorophores employed in FRET can contribute to the observed discrepancy. Specifically, we find that addition of Alexa-488 to a normally expanded IDP causes contraction by an additional 15%, a value in reasonable accord with the contraction reported in FRET-based studies. Second, using our simulations and analysis procedure to accurately extract both the radius of gyration (R_g) and end-to-end distance (R_ee) from SAXS profiles, we tested the recent suggestion that FRET and SAXS results can be reconciled if the R_g and R_ee are "uncoupled" (i.e., no longer simply proportional), in contrast to the case for random walk homopolymers. We find, however, that even for unfolded proteins, these two measures of unfolded state dimensions remain proportional. Together, these results suggest that improved analysis procedures and a correction for significant, fluorophore-driven interactions are sufficient to reconcile prior SAXS and FRET studies, thus providing a unified picture of the nature of unfolded polypeptide chains in the absence of denaturant.

A protein interaction free energy model based on amino acid residue contributions: Assessment of point mutation stability of T4 lysozyme

Here we present a model to estimate the interaction free energy contribution of each amino acid residue of a given protein. Protein interaction energy is described in terms of per-residue interaction factors, μ. Multibody interactions are implicitly captured in μ through the combination of amino acid terms (γ) guided by local conformation indices (σ). The model enables construction of an interaction factor heat map for a protein in a given fold, allows prima facie assessment of the degree of residue-residue interaction, and facilitates a qualitative and quantitative evaluation of protein association properties. The model was used to compute thermal stability of T4 bacteriophage lysozyme mutants across seven sites. Qualitative assessment of mutational effects provides a straightforward rationale regarding whether a particular site primarily perturbs native or non-native states, or both. The presented model was found to be in good agreement with experimental mutational data (R 2 = 0.73) and suggests an approach by which to convert structure space into energy space.

Structural Basis for Teneurin Function in Circuit-Wiring: A Toxin Motif at the Synapse

Teneurins (TENs) are cell-surface adhesion proteins with critical roles in tissue development and axon guidance. Here, we report the 3.1-Å cryoelectron microscopy structure of the human TEN2 extracellular region (ECR), revealing a striking similarity to bacterial Tc-toxins. The ECR includes a large β barrel that partially encapsulates a C-terminal domain, which emerges to the solvent through an opening in the mid-barrel region. An immunoglobulin (Ig)-like domain seals the bottom of the barrel while a β propeller is attached in a perpendicular orientation. We further show that an alternatively spliced region within the β propeller acts as a switch to regulate trans-cellular adhesion of TEN2 to latrophilin (LPHN), a transmembrane receptor known to mediate critical functions in the central nervous system. One splice variant activates trans-cellular signaling in a LPHN-dependent manner, whereas the other induces inhibitory postsynaptic differentiation. These results highlight the unusual structural organization of TENs giving rise to their multifarious functions.

Trajectory-based training enables protein simulations with accurate folding and Boltzmann ensembles in cpu-hours

An ongoing challenge in protein chemistry is to identify the underlying interaction energies that capture protein dynamics. The traditional trade-off in biomolecular simulation between accuracy and computational efficiency is predicated on the assumption that detailed force fields are typically well-parameterized, obtaining a significant fraction of possible accuracy. We re-examine this trade-off in the more realistic regime in which parameterization is a greater source of error than the level of detail in the force field. To address parameterization of coarse-grained force fields, we use the contrastive divergence technique from machine learning to train from simulations of 450 proteins. In our procedure, the computational efficiency of the model enables high accuracy through the precise tuning of the Boltzmann ensemble. This method is applied to our recently developed Upside model, where the free energy for side chains is rapidly calculated at every time-step, allowing for a smooth energy landscape without steric rattling of the side chains. After this contrastive divergence training, the model is able to de novo fold proteins up to 100 residues on a single core in days. This improved Upside model provides a starting point both for investigation of folding dynamics and as an inexpensive Bayesian prior for protein physics that can be integrated with additional experimental or bioinformatic data.

Response to Comment on “Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water”

Best et al. claim that we provide no convincing basis to assert that a discrepancy remains between FRET and SAXS results on the dimensions of disordered proteins under physiological conditions. We maintain that a clear discrepancy is apparent in our and other recent publications, including results shown in the Best et al. comment. A plausible origin is fluorophore interactions in FRET experiments.

New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions

In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

Conformational Heterogeneity and FRET Data Interpretation for Dimensions of Unfolded Proteins

A mathematico-physically valid formulation is required to infer properties of disordered protein conformations from single-molecule Förster resonance energy transfer (smFRET). Conformational dimensions inferred by conventional approaches that presume a homogeneous conformational ensemble can be unphysical. When all possible-heterogeneous as well as homogeneous-conformational distributions are taken into account without prejudgment, a single value of average transfer efficiency 〈E〉 between dyes at two chain ends is generally consistent with highly diverse, multiple values of the average radius of gyration 〈R_g〉. Here we utilize unbiased conformational statistics from a coarse-grained explicit-chain model to establish a general logical framework to quantify this fundamental ambiguity in smFRET inference. As an application, we address the long-standing controversy regarding the denaturant dependence of 〈R_g〉 of unfolded proteins, focusing on Protein L as an example. Conventional smFRET inference concluded that 〈R_g〉 of unfolded Protein L is highly sensitive to [GuHCl], but data from SAXS suggested a near-constant 〈R_g〉 irrespective of [GuHCl]. Strikingly, our analysis indicates that although the reported 〈E〉 values for Protein L at [GuHCl] = 1 and 7 M are very different at 0.75 and 0.45, respectively, the Bayesian R_g² distributions consistent with these two 〈E〉 values overlap by as much as 75%. Our findings suggest, in general, that the smFRET-SAXS discrepancy regarding unfolded protein dimensions likely arise from highly heterogeneous conformational ensembles at low or zero denaturant, and that additional experimental probes are needed to ascertain the nature of this heterogeneity.

Innovative scattering analysis shows that hydrophobic disordered proteins are expanded in water

A substantial fraction of the proteome is intrinsically disordered, and even well-folded proteins adopt non-native geometries during synthesis, folding, transport, and turnover. Characterization of intrinsically disordered proteins (IDPs) is challenging, in part because of a lack of accurate physical models and the difficulty of interpreting experimental results. We have developed a general method to extract the dimensions and solvent quality (self-interactions) of IDPs from a single small-angle x-ray scattering measurement. We applied this procedure to a variety of IDPs and found that even IDPs with low net charge and high hydrophobicity remain highly expanded in water, contrary to the general expectation that protein-like sequences collapse in water. Our results suggest that the unfolded state of most foldable sequences is expanded; we conjecture that this property was selected by evolution to minimize misfolding and aggregation.

Dielectric screening effect of electronic polarization and intramolecular hydrogen bonding

Recent site-resolved hydrogen exchange measurements have uncovered significant discrepancies between simulations and experimental data during protein folding, including the excessive intramolecular hydrogen bonds in simulations. This finding indicates a possibility that intramolecular charge-charge interactions have not included sufficient dielectric screening effect of the electronic polarization. Scaling down peptide atomic charges according to the optical dielectric constant is tested in this study. As a result, the number of intramolecular hydrogen bonds is lower than using unscaled atomic charges while reaching the same levels of helical contents or β-hairpin backbone hydrogen bonds, because van der Waals interactions contribute substantially to peptide folding in water. Reducing intramolecular charge-charge interactions and hydrogen bonding increases conformational search efficiency. In particular, it reduces the equilibrium helical content in simulations using AMBER force field and the energy barrier in folding simulations using CHARMM force field.

Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field

An accurate potential energy model is crucial for biomolecular simulations. Despite many recent improvements of classical protein force fields, there are remaining key issues: much weaker temperature dependence of folding/unfolding equilibrium and overly collapsed unfolded or disordered states. For the latter problem, a new water model (TIP4P-D) has been proposed to correct the significantly underestimated water dispersion interactions. Here, using TIP4P-D, we reveal problems in current force fields through failures in folding model systems (a polyalanine peptide, Trp-cage, and the GB1 hairpin). By using residue-specific parameters to achieve better match between amino acid sequences and native structures and adding a small H-bond correction to partially compensate the missing many-body effects in α-helix formation, the new RSFF2+ force field with the TIP4P-D water model can excellently reproduce experimental melting curves of both α-helical and β-hairpin systems. The RSFF2+/TIP4P-D method also gives less collapsed unfolded structures and describes well folded proteins simultaneously.

Stress-Triggered Phase Separation Is an Adaptive, Evolutionarily Tuned Response

In eukaryotic cells, diverse stresses trigger coalescence of RNA-binding proteins into stress granules. In vitro, stress-granule-associated proteins can demix to form liquids, hydrogels, and other assemblies lacking fixed stoichiometry. Observing these phenomena has generally required conditions far removed from physiological stresses. We show that poly(A)-binding protein (Pab1 in yeast), a defining marker of stress granules, phase separates and forms hydrogels in vitro upon exposure to physiological stress conditions. Other RNA-binding proteins depend upon low-complexity regions (LCRs) or RNA for phase separation, whereas Pab1's LCR is not required for demixing, and RNA inhibits it. Based on unique evolutionary patterns, we create LCR mutations, which systematically tune its biophysical properties and Pab1 phase separation in vitro and in vivo. Mutations that impede phase separation reduce organism fitness during prolonged stress. Poly(A)-binding protein thus acts as a physiological stress sensor, exploiting phase separation to precisely mark stress onset, a broadly generalizable mechanism.

Preferential binding effects on protein structure and dynamics revealed by coarse-grained Monte Carlo simulation

The effect of preferential binding of solute molecules within an aqueous solution on the structure and dynamics of the histone H3.1 protein is examined by a coarse-grained Monte Carlo simulation. The knowledge-based residue-residue and hydropathy-index-based residue-solvent interactions are used as input to analyze a number of local and global physical quantities as a function of the residue-solvent interaction strength (f). Results from simulations that treat the aqueous solution as a homogeneous effective solvent medium are compared to when positional fluctuations of the solute molecules are explicitly considered. While the radius of gyration (Rg) of the protein exhibits a non-monotonic dependence on solvent interaction over a wide range of f within an effective medium, an abrupt collapse in Rg occurs in a narrow range of f when solute molecules rapidly bind to a preferential set of sites on the protein. The structure factor S(q) of the protein with wave vector (q) becomes oscillatory in the collapsed state, which reflects segmental correlations caused by spatial fluctuations in solute-protein binding. Spatial fluctuations in solute binding also modify the effective dimension (D) of the protein in fibrous (D ∼ 1.3), random-coil (D ∼ 1.75), and globular (D ∼ 3) conformational ensembles as the interaction strength increases, which differ from an effective medium with respect to the magnitude of D and the length scale.

Optimization of the GBMV2 implicit solvent force field for accurate simulation of protein conformational equilibria

Accurate treatment of solvent environment is critical for reliable simulations of protein conformational equilibria. Implicit treatment of solvation, such as using the generalized Born (GB) class of models arguably provides an optimal balance between computational efficiency and physical accuracy. Yet, GB models are frequently plagued by a tendency to generate overly compact structures. The physical origins of this drawback are relatively well understood, and the key to a balanced implicit solvent protein force field is careful optimization of physical parameters to achieve a sufficient level of cancellation of errors. The latter has been hampered by the difficulty of generating converged conformational ensembles of non-trivial model proteins using the popular replica exchange sampling technique. Here, we leverage improved sampling efficiency of a newly developed multi-scale enhanced sampling technique to re-optimize the generalized-Born with molecular volume (GBMV2) implicit solvent model with the CHARMM36 protein force field. Recursive optimization of key GBMV2 parameters (such as input radii) and protein torsion profiles (via the CMAP torsion cross terms) has led to a more balanced GBMV2 protein force field that recapitulates the structures and stabilities of both helical and β-hairpin model peptides. Importantly, this force field appears to be free of the over-compaction bias, and can generate structural ensembles of several intrinsically disordered proteins of various lengths that seem highly consistent with available experimental data. © 2017 Wiley Periodicals, Inc.

Building a More Predictive Protein Force Field: A Systematic and Reproducible Route to AMBER-FB15

The increasing availability of high-quality experimental data and first-principles calculations creates opportunities for developing more accurate empirical force fields for simulation of proteins. We developed the AMBER-FB15 protein force field by building a high-quality quantum chemical data set consisting of comprehensive potential energy scans and employing the ForceBalance software package for parameter optimization. The optimized potential surface allows for more significant thermodynamic fluctuations away from local minima. In validation studies where simulation results are compared to experimental measurements, AMBER-FB15 in combination with the updated TIP3P-FB water model predicts equilibrium properties with equivalent accuracy, and temperature dependent properties with significantly improved accuracy, in comparison with published models. We also discuss the effect of changing the protein force field and water model on the simulation results.

Computational and theoretical advances in studies of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are increasingly realized to play diverse biological roles, ranging from molecular signaling to the formation of membraneless organelles. Their high degree of disorder makes them more challenging to study using the techniques of conventional structural biology, because any observable will be averaged over a heterogeneous ensemble of structures. Molecular simulations and theory are therefore a natural complement to experiment for studying the structure, dynamics and function of IDPs. The diverse time and length scales relevant to the roles played by IDPs require flexibility in the techniques applied. Here, I summarize some of the developments in simulation and theory in recent years, which have been driven by the desire to better capture IDP properties at different time- and length-scales. I also provide an outlook for how methods can be improved in the future and emerging problems which may be addressed by theory and simulation.

Evidence for the principle of minimal frustration in the evolution of protein folding landscapes

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Refined Parameterization of Nonbonded Interactions Improves Conformational Sampling and Kinetics of Protein Folding Simulations

Recent advances in computational technology have enabled brute-force molecular dynamics (MD) simulations of protein folding using physics-based molecular force fields. The extensive sampling of protein conformations afforded by such simulations revealed, however, considerable compaction of the protein conformations in the unfolded state, which is inconsistent with experiment. Here, we show that a set of surgical corrections to nonbonded interactions between amine nitrogen-carboxylate oxygen and aliphatic carbon-carbon atom pairs can considerably improve the realism of protein folding simulations. Specifically, we show that employing our corrections in ∼500 μs all-atom replica-exchange MD simulations of the WW domain and villin head piece proteins increases the size of the denatured proteins' conformations and does not destabilize the native conformations of the proteins. In addition to making the folded conformations a global minimum of the respective free energy landscapes at room temperature, our corrections also make the free energy landscape smoother, considerably accelerating the folding kinetics and, hence, reducing the computational expense of a protein folding simulation.

Theoretical Insights into the Biophysics of Protein Bi-stability and Evolutionary Switches

Deciphering the effects of nonsynonymous mutations on protein structure is central to many areas of biomedical research and is of fundamental importance to the study of molecular evolution. Much of the investigation of protein evolution has focused on mutations that leave a protein’s folded structure essentially unchanged. However, to evolve novel folds of proteins, mutations that lead to large conformational modifications have to be involved. Unraveling the basic biophysics of such mutations is a challenge to theory, especially when only one or two amino acid substitutions cause a large-scale conformational switch. Among the few such mutational switches identified experimentally, the one between the G_A all-α and G_B α+β folds is extensively characterized; but all-atom simulations using fully transferrable potentials have not been able to account for this striking switching behavior. Here we introduce an explicit-chain model that combines structure-based native biases for multiple alternative structures with a general physical atomic force field, and apply this construct to twelve mutants spanning the sequence variation between G_A and G_B. In agreement with experiment, we observe conformational switching from G_A to G_B upon a single L45Y substitution in the GA98 mutant. In line with the latent evolutionary potential concept, our model shows a gradual sequence-dependent change in fold preference in the mutants before this switch. Our analysis also indicates that a sharp G_A/G_B switch may arise from the orientation dependence of aromatic π-interactions. These findings provide physical insights toward rationalizing, predicting and designing evolutionary conformational switches.

The biological functions of globular proteins are intimately related to their folded structures and their associated conformational fluctuations. Evolution of new structures is an important avenue to new functions. Although many mutations do not change the folded state, experiments indicate that a single amino acid substitution can lead to a drastic change in the folded structure. The physics of this switch-like behavior remains to be elucidated. Here we develop a computational model for the relevant physical forces, showing that mutations can lead to new folds by passing through intermediate sequences where the old and new folds occur with varying probabilities. Our approach helps provide a general physical account of conformational switching in evolution and mutational effects on conformational dynamics.

Can Local Probes Go Global? A Joint Experiment-Simulation Analysis of λ(6-85) Folding

The process of protein folding is known to involve global motions in a cooperative affair; the structure of most of the protein sequences is gained or lost over a narrow range of temperature, denaturant, or pressure perturbations. At the same time, recent simulations and experiments reveal a complex structural landscape with a rich set of local motions and conformational changes. We couple experimental kinetic and thermodynamic measurements with specifically tailored analysis of simulation data to isolate local versus global folding probes. We find that local probes exhibit lower melting temperatures, smaller surface area changes, and faster kinetics compared to global ones. We also see that certain local probes of folding match the global behavior more closely than others. Our work highlights the importance of using multiple probes to fully characterize protein folding dynamics by theory and experiment.

Cooperative folding near the downhill limit determined with amino acid resolution by hydrogen exchange

The relationship between folding cooperativity and downhill, or barrier-free, folding of proteins under highly stabilizing conditions remains an unresolved topic, especially for proteins such as λ-repressor that fold on the microsecond timescale. Under aqueous conditions where downhill folding is most likely to occur, we measure the stability of multiple H bonds, using hydrogen exchange (HX) in a λYA variant that is suggested to be an incipient downhill folder having an extrapolated folding rate constant of 2 × 10(5) s(-1) and a stability of 7.4 kcal·mol(-1) at 298 K. At least one H bond on each of the three largest helices (α1, α3, and α4) breaks during a common unfolding event that reflects global denaturation. The use of HX enables us to both examine folding under highly stabilizing, native-like conditions and probe the pretransition state region for stable species without the need to initiate the folding reaction. The equivalence of the stability determined at zero and high denaturant indicates that any residual denatured state structure minimally affects the stability even under native conditions. Using our ψ analysis method along with mutational ϕ analysis, we find that the three aforementioned helices are all present in the folding transition state. Hence, the free energy surface has a sufficiently high barrier separating the denatured and native states that folding appears cooperative even under extremely stable and fast folding conditions.

Isotope-enriched protein standards for computational amide I spectroscopy

We present a systematic isotope labeling study of the protein G mutant NuG2b as a step toward the production of reliable, structurally stable, experimental standards for amide I infrared spectroscopic simulations. By introducing isotope enriched amino acids into a minimal growth medium during bacterial expression, we induce uniform labeling of the amide bonds following specific amino acids, avoiding the need for chemical peptide synthesis. We use experimental data to test several common amide I frequency maps and explore the influence of various factors on map performance. Comparison of the predicted absorption frequencies for the four maps tested with empirical assignments to our experimental spectra yields a root-mean-square error of 6-12 cm(-1), with outliers of at least 12 cm(-1) in all models. This means that the predictions may be useful for predicting general trends such as changes in hydrogen bonding configuration; however, for finer structural constraints or absolute frequency assignments, the models are unreliable. The results indicate the need for careful testing of existing literature maps and shed light on possible next steps for the development of quantitative spectral maps.

Quantitative interpretation of FRET experiments via molecular simulation: force field and validation

Molecular simulation is a valuable and complementary tool that may assist with the interpretation of single-molecule Förster resonance energy transfer (FRET) experiments, if the energy function is of sufficiently high quality. Here we present force-field parameters for one of the most common pairs of chromophores used in experiments, AlexaFluor 488 and 594. From microsecond molecular-dynamics simulations, we are able to recover both experimentally determined equilibrium constants and association/dissociation rates of the chromophores with free tryptophan, as well as the decay of fluorescence anisotropy of a labeled protein. We find that it is particularly important to obtain a correct balance of solute-water interactions in the simulations in order to faithfully capture the experimental anisotropy decays, which provide a sensitive benchmark for fluorophore mobility. Lastly, by a combination of experiment and simulation, we address a potential complication in the interpretation of experiments on polyproline, used as a molecular ruler for FRET experiments, namely the potential association of one of the chromophores with the polyproline helix. Under conditions where simulations accurately capture the fluorescence anisotropy decay, we find at most a modest, transient population of conformations in which the chromophores associate with the polyproline. Explicit calculation of FRET transfer efficiencies for short polyprolines yields results in good agreement with experiment. These results illustrate the potential power of a combination of molecular simulation and experiment in quantifying biomolecular dynamics.

Where the complex things are: single molecule and ensemble spectroscopic investigations of protein folding dynamics

Progress in our understanding of the simple folding dynamics of small proteins and the complex dynamics of large proteins is reviewed. Recent characterizations of the folding transition path of small proteins revealed a simple dynamics explainable by the native centric model. In contrast, the accumulated data showed the substates containing residual structures in the unfolded state and partially populated intermediates, causing complexity in the early folding dynamics of small proteins. The size of the unfolded proteins in the absence of denaturants is likely expanded but still controversial. The steady progress in the observation of folding of large proteins has clarified the rapid formation of long-range contacts that seem inconsistent with the native centric model, suggesting that the folding strategy of large proteins is distinct from that of small proteins.