ALS-linked PFN1 variants exhibit loss and gain of functions in the context of formin-induced actin polymerization.

Profilin-1 (PFN1) plays important roles in modulating actin dynamics through binding both monomeric actin and proteins enriched with polyproline motifs. Mutations in PFN1 have been linked to the neurodegenerative disease amyotrophic lateral sclerosis (ALS). However, whether ALS-linked mutations affect PFN1 function has remained unclear. To address this question, we employed an unbiased proteomics analysis in mammalian cells to identify proteins that differentially interact with mutant and wild-type (WT) PFN1. These studies uncovered differential binding between two ALS-linked PFN1 variants, G118V and M114T, and select formin proteins. Furthermore, both variants augmented formin-mediated actin assembly relative to PFN1 WT. Molecular dynamics simulations revealed mutation-induced changes in the internal dynamic couplings within an alpha helix of PFN1 that directly contacts both actin and polyproline, as well as structural fluctuations within the actin- and polyproline-binding regions of PFN1. These data indicate that ALS-PFN1 variants have the potential for heightened flexibility in the context of the ternary actin-PFN1-polyproline complex during actin assembly. Conversely, PFN1 C71G was more severely destabilized than the other PFN1 variants, resulting in reduced protein expression in both transfected and ALS patient lymphoblast cell lines. Moreover, this variant exhibited loss-of-function phenotypes in the context of actin assembly. Perturbations in actin dynamics and assembly can therefore result from ALS-linked mutations in PFN1. However, ALS-PFN1 variants may dysregulate actin polymerization through different mechanisms that depend upon the solubility and stability of the mutant protein.

Friction-Limited Folding of Disulfide-Reduced Monomeric SOD1.

The folding reaction of a stable monomeric variant of Cu/Zn superoxide dismutase (mSOD1), an enzyme responsible for the conversion of superoxide free radicals into hydrogen peroxide and oxygen, is known to be among the slowest folding processes that adhere to two-state behavior. The long lifetime, ∼10 s, of the unfolded state presents ample opportunities for the polypeptide chain to transiently sample nonnative structures before the formation of the productive folding transition state. We recently observed the formation of a nonnative structure in a peptide model of the C-terminus of SOD1, a sequence that might serve as a potential source of internal chain friction-limited folding. To test for friction-limited folding, we performed a comprehensive thermodynamic and kinetic analysis of the folding mechanism of mSOD1 in the presence of the viscogens glycerol and glucose. Using a, to our knowledge, novel analysis of the folding reactions, we found the disulfide-reduced form of the protein that exposes the C-terminal sequence, but not its disulfide-oxidized counterpart that protects it, experiences internal chain friction during folding. The sensitivity of the internal friction to the disulfide bond status suggests that one or both of the cross-linked regions play a critical role in driving the friction-limited folding. We speculate that the molecular mechanisms giving rise to the internal friction of disulfide-reduced mSOD1 might play a role in the amyotrophic lateral sclerosis-linked aggregation of SOD1.

Structural Organization and Dynamics of Homodimeric Cytohesin Family Arf GTPase Exchange Factors in Solution and on Membranes.

Membrane dynamic processes require Arf GTPase activation by guanine nucleotide exchange factors (GEFs) with a Sec7 domain. Cytohesin family Arf GEFs function in signaling and cell migration through Arf GTPase activation on the plasma membrane and endosomes. In this study, the structural organization of two cytohesins (Grp1 and ARNO) was investigated in solution by size exclusion-small angle X-ray scattering and negative stain-electron microscopy and on membranes by dynamic light scattering, hydrogen-deuterium exchange-mass spectrometry and guanosine diphosphate (GDP)/guanosine triphosphate (GTP) exchange assays. The results suggest that cytohesins form elongated dimers with a central coiled coil and membrane-binding pleckstrin-homology (PH) domains at opposite ends. The dimers display significant conformational heterogeneity, with a preference for compact to intermediate conformations. Phosphoinositide-dependent membrane recruitment is mediated by one PH domain at a time and alters the conformational dynamics to prime allosteric activation by Arf-GTP. A structural model for membrane targeting and allosteric activation of full-length cytohesin dimers is discussed.

In Situ Measurements of Polymer Micellization Kinetics with Millisecond Temporal Resolution.

Utilizing synchrotron small-angle X-ray scattering (SAXS) integrated with a microfluidic device, micellization kinetics of a diblock co-polymer, poly(ethylene glycol)-b-poly(caprolactone) (PEG-b-PCL) was measured in situ with millisecond temporal and micrometer spatial resolution. The evolutionary regimes of polymer micellization - nucleation, fusion, and insertion were directly observed. The five-inlet microfluidic device provided steady continuous mixing of the polymer solution and the antisolvent. Solvent replacement was mainly dominated by lateral diffusion across the hydrodynamically focused central layer, whose thickness could be precisely designed and manipulated from mass balance of the partitioning streams. Knowing the micellization kinetics of the polymers is essential for design and optimization of self-assembled polymeric nanostructures. The technique of integrating SAXS with microfluidic devices can be translatable to other systems for a breadth of applications.

Frustration and folding of a TIM barrel protein.

Triosephosphate isomerase (TIM) barrel proteins have not only a conserved architecture that supports a myriad of enzymatic functions, but also a conserved folding mechanism that involves on- and off-pathway intermediates. Although experiments have proven to be invaluable in defining the folding free-energy surface, they provide only a limited understanding of the structures of the partially folded states that appear during folding. Coarse-grained simulations employing native centric models are capable of sampling the entire energy landscape of TIM barrels and offer the possibility of a molecular-level understanding of the readout from sequence to structure. We have combined sequence-sensitive native centric simulations with small-angle X-ray scattering and time-resolved Förster resonance energy transfer to monitor the formation of structure in an intermediate in the Sulfolobus solfataricus indole-3-glycerol phosphate synthase TIM barrel that appears within 50 µs and must at least partially unfold to achieve productive folding. Simulations reveal the presence of a major and 2 minor folding channels not detected in experiments. Frustration in folding, i.e., backtracking in native contacts, is observed in the major channel at the initial stage of folding, as well as late in folding in a minor channel before the appearance of the native conformation. Similarities in global and pairwise dimensions of the early intermediate, the formation of structure in the central region that spreads progressively toward each terminus, and a similar rate-limiting step in the closing of the ß-barrel underscore the value of combining simulation and experiment to unravel complex folding mechanisms at the molecular level.

Nonnative structure in a peptide model of the unfolded state of superoxide dismutase 1 (SOD1): Implications for ALS-linked aggregation.

Dozens of mutations throughout the sequence of the gene encoding superoxide dismutase 1 (SOD1) have been linked to toxic protein aggregation in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). A parsimonious explanation for numerous genotypes resulting in a common phenotype would be mutation-induced perturbation of the folding free-energy surface that increases the populations of high-energy states prone to aggregation. The absence of intermediates in the folding of monomeric SOD1 suggests that the unfolded ensemble is a potential source of aggregation. To test this hypothesis, here we dissected SOD1 into a set of peptides end-labeled with FRET probes to model the local behavior of the corresponding sequences in the unfolded ensemble. Using time-resolved FRET, we observed that the peptide corresponding to the Loop VII-ß8 sequence at the SOD1 C terminus was uniquely sensitive to denaturant. Utilizing a two-dimensional form of maximum entropy modeling, we demonstrate that the sensitivity to denaturant is the surprising result of a two-state-like transition from a compact to an expanded state. Variations of the peptide sequence revealed that the compact state involves a nonnative interaction between the disordered N terminus and the hydrophobic C terminus of the peptide. This nonnative intramolecular structure could serve as a precursor for intermolecular association and result in aggregation associated with ALS. We propose that this precursor would provide a common molecular target for therapeutic intervention in the dozens of ALS-linked SOD1 mutations.

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.

Transient Kinetic Analysis of SWR1C-Catalyzed H2A.Z Deposition Unravels the Impact of Nucleosome Dynamics and the Asymmetry of Histone Exchange.

The SWR1C chromatin remodeling enzyme catalyzes ATP-dependent replacement of nucleosomal H2A with the H2A.Z variant, regulating key DNA-mediated processes such as transcription and DNA repair. Here, we investigate the transient kinetic mechanism of the histone exchange reaction, employing ensemble FRET, fluorescence correlation spectroscopy (FCS), and the steady-state kinetics of ATP hydrolysis. Our studies indicate that SWR1C modulates nucleosome dynamics on both the millisecond and microsecond timescales, poising the nucleosome for the dimer exchange reaction. The transient kinetic analysis of the remodeling reaction performed under single turnover conditions unraveled a striking asymmetry in the ATP-dependent replacement of nucleosomal dimers, promoted by localized DNA unwrapping. Taken together, our transient kinetic studies identify intermediates and provide crucial insights into the SWR1C-catalyzed dimer exchange reaction and shed light on how the mechanics of H2A.Z deposition might contribute to transcriptional regulation in vivo.

Computer design of microfluidic mixers for protein/RNA folding studies.

Kinetic studies of biological macromolecules increasingly use microfluidic mixers to initiate and monitor reaction progress. A motivation for using microfluidic mixers is to reduce sample consumption and decrease mixing time to microseconds. Some applications, such as small-angle x-ray scattering, also require large (>10 micron) sampling areas to ensure high signal-to-noise ratios and to minimize parasitic scattering. Chaotic to marginally turbulent mixers are well suited for these applications because this class of mixers provides a good middle ground between existing laminar and turbulent mixers. In this study, we model various chaotic to marginally turbulent mixing concepts such as flow turning, flow splitting, and vortex generation using computational fluid dynamics for optimization of mixing efficiency and observation volume. Design iterations show flow turning to be the best candidate for chaotic/marginally turbulent mixing. A qualitative experimental test is performed on the finalized design with mixing of 10 M urea and water to validate the flow turning unsteady mixing concept as a viable option for RNA and protein folding studies. A comparison of direct numerical simulations (DNS) and turbulence models suggests that the applicability of turbulence models to these flow regimes may be limited.

Structural Dynamics Control Allosteric Activation of Cytohesin Family Arf GTPase Exchange Factors.

Membrane dynamic processes including vesicle biogenesis depend on Arf guanosine triphosphatase (GTPase) activation by guanine nucleotide exchange factors (GEFs) containing a catalytic Sec7 domain and a membrane-targeting module such as a pleckstrin homology (PH) domain. The catalytic output of cytohesin family Arf GEFs is controlled by autoinhibitory interactions that impede accessibility of the exchange site in the Sec7 domain. These restraints can be relieved through activator Arf-GTP binding to an allosteric site comprising the PH domain and proximal autoinhibitory elements (Sec7-PH linker and C-terminal helix). Small-angle X-ray scattering and negative-stain electron microscopy were used to investigate the structural organization and conformational dynamics of cytohesin-3 (Grp1) in autoinhibited and active states. The results support a model in which hinge dynamics in the autoinhibited state expose the activator site for Arf-GTP binding, while subsequent C-terminal helix unlatching and repositioning unleash conformational entropy in the Sec7-PH linker to drive exposure of the exchange site.

Atomistic structural ensemble refinement reveals non-native structure stabilizes a sub-millisecond folding intermediate of CheY.

The dynamics of globular proteins can be described in terms of transitions between a folded native state and less-populated intermediates, or excited states, which can play critical roles in both protein folding and function. Excited states are by definition transient species, and therefore are difficult to characterize using current experimental techniques. Here, we report an atomistic model of the excited state ensemble of a stabilized mutant of an extensively studied flavodoxin fold protein CheY. We employed a hybrid simulation and experimental approach in which an aggregate 42 milliseconds of all-atom molecular dynamics were used as an informative prior for the structure of the excited state ensemble. This prior was then refined against small-angle X-ray scattering (SAXS) data employing an established method (EROS). The most striking feature of the resulting excited state ensemble was an unstructured N-terminus stabilized by non-native contacts in a conformation that is topologically simpler than the native state. Using these results, we then predict incisive single molecule FRET experiments as a means of model validation. This study demonstrates the paradigm of uniting simulation and experiment in a statistical model to study the structure of protein excited states and rationally design validating experiments.

Selenomethionine Quenching of Tryptophan Fluorescence Provides a Simple Probe of Protein Structure.

Fluorescence spectroscopy, relying on intrinsic protein fluorophores, is one of the most widely used methods for studying protein folding, protein-ligand interactions, and protein dynamics. Tryptophan is usually the fluorophore of choice, given its sensitivity to its environment and having the highest quantum yield of the natural amino acids; however, changes in tryptophan fluorescence can be difficult to interpret in terms of specific structural changes. The introduction of quenchers of tryptophan fluorescence can provide information about specific structures, particularly if quenching is short-range; however, the most commonly employed quencher is histidine, and it is effective only when the imidazole side chain is protonated, thus limiting the pH range over which this approach can be employed. In addition, histidine is not always a conservative substitution and is likely to be destabilizing if inserted into the hydrophobic core of proteins. Here we illustrate the use of a Trp-selenomethionine (MSe) pair as a specific probe of protein structure. MSe requires a close approach to Trp to quench its fluorescence, and this effect can be exploited to design specific probes of α-helix and ß-sheet formation. The approach is illustrated using equilibrium and time-resolved fluorescence measurements of designed peptides and globular proteins. MSe is easily incorporated into proteins and provides a conservative replacement for hydrophobic side chains, and MSe quenching of Trp fluorescence is pH-independent. The oxidized form of MSe, selenomethionine selenoxide, is also an efficient quencher of Trp fluorescence.

Dye-Sensitized Core/Active Shell Upconversion Nanoparticles for Optogenetics and Bioimaging Applications.

Near-infrared (NIR) dye-sensitized upconversion nanoparticles (UCNPs) can broaden the absorption range and boost upconversion efficiency of UCNPs. Here, we achieved significantly enhanced upconversion luminescence in dye-sensitized core/active shell UCNPs via the doping of ytterbium ions (Yb(3+)) in the UCNP shell, which bridged the energy transfer from the dye to the UCNP core. As a result, we synergized the two most practical upconversion booster effectors (dye-sensitizing and core/shell enhancement) to amplify upconversion efficiency. We demonstrated two biomedical applications using these UCNPs. By using dye-sensitized core/active shell UCNP embedded poly(methyl methacrylate) polymer implantable systems, we successfully shifted the optogenetic neuron excitation window to a biocompatible and deep tissue penetrable 800 nm wavelength. Furthermore, UCNPs were water-solubilized with Pluronic F127 with high upconversion efficiency and can be imaged in a mouse model.

Sequential Closure of Loop Structures Forms the Folding Nucleus during the Refolding Transition of the Escherichia coli Adenylate Kinase Molecule.

The ensemble of conformers of globular protein molecules immediately following transfer from unfolding to folding conditions is assumed to be collapsed though still disordered, as the first steps of the folding pathway are initiated. In order to test the hypothesis that long loop closure transitions are part of the initiation of the folding pathway, our groups are studying the initiation of the folding transition of a model protein by time-resolved excitation energy transfer (trFRET) detected fast kinetics experiments. Site-specific double labeling is used to study the timing of conformational transitions of individual loop forming chain segments at the microsecond time regime. Previously, it was shown that at least three long loops in the Escherichia coli adenylate kinase (AK) molecule close within the first 5 ms of folding of AK, while the main global folding transition occurs in a time regime of seconds. In order to enhance the time resolution of the kinetics experiments to the microsecond time regime and determine the rate of closure of the two N terminal loops (loop I residues 1-26 and loop II residues 29-72), we applied a continuous flow based double kinetics experiment. These measurements enabled us to obtain a microsecond series of transient time dependent distributions of distances between the ends of the labeled loops. Analysis of the trFRET experiments show that the N terminal loop (loop I) is closed within less than 60 µs after the initiation of refolding. Loop II is also mostly closed within that time step but shows an additional small reduction of the mean end-to-end distance in a second phase at a rate of 0.005 µs(-1). This second phase can either reflect tightening of a loosely closed loop in the ensemble of conformers or may reflect two subpopulations in the ensemble, which differ in the rate of closure of loop II, but not in the rate of closure of loop I. This study shows the very fast closure of long loops in the otherwise disordered backbone and fine details of the very early hidden pretransition state steps that are essential for the fast and efficient folding of the protein molecule.

Selenomethionine, p-cyanophenylalanine pairs provide a convenient, sensitive, non-perturbing fluorescent probe of local helical structure.

The use of selenomethionine (MSe)-p-cyanophenylalanine (FCN) pairs to probe protein structure is demonstrated. MSe quenches FCN fluorescence via electron transfer. Both residues can be incorporated recombinantly or by peptide synthesis. Time-resolved and steady-state fluorescence measurements demonstrate that MSe-FCN pairs provide specific local probes of helical structure.

Tailoring dye-sensitized upconversion nanoparticle excitation bands towards excitation wavelength selective imaging.

One of the key roadblocks in UCNP development is its extremely limited choices of excitation wavelengths. We report a generic design to program UCNPs to possess highly tunable dye characteristic excitation bands. Using such distinctive properties, we were able to develop a new excitation wavelength selective security imaging. This work unleashed the greater freedom of the excitation wavelengths of the upconversion nanoparticles and we believe it is a game-changer in the field and this method will enable numerous applications that are currently limited by existing UCNPs.

Methods for analysis of size-exclusion chromatography-small-angle X-ray scattering and reconstruction of protein scattering.

Size-exclusion chromatography in line with small-angle X-ray scattering (SEC-SAXS) has emerged as an important method for investigation of heterogeneous and self-associating systems, but presents specific challenges for data processing including buffer subtraction and analysis of overlapping peaks. This paper presents novel methods based on singular value decomposition (SVD) and Guinier-optimized linear combination (LC) to facilitate analysis of SEC-SAXS data sets and high-quality reconstruction of protein scattering directly from peak regions. It is shown that Guinier-optimized buffer subtraction can reduce common subtraction artifacts and that Guinier-optimized linear combination of significant SVD basis components improves signal-to-noise and allows reconstruction of protein scattering, even in the absence of matching buffer regions. In test cases with conventional SAXS data sets for cytochrome c and SEC-SAXS data sets for the small GTPase Arf6 and the Arf GTPase exchange factors Grp1 and cytohesin-1, SVD-LC consistently provided higher quality reconstruction of protein scattering than either direct or Guinier-optimized buffer subtraction. These methods have been implemented in the context of a Python-extensible Mac OS X application known as Data Evaluation and Likelihood Analysis (DELA), which provides convenient tools for data-set selection, beam intensity normalization, SVD, and other relevant processing and analytical procedures, as well as automated Python scripts for common SAXS analyses and Guinier-optimized reconstruction of protein scattering.

Non-native structure appears in microseconds during the folding of E. coli RNase H.

The folding pathway of Escherichia coli RNase H is one of the best experimentally characterized for any protein. In spite of this, spectroscopic studies have never captured the earliest events. Using continuous-flow microfluidic mixing, we have now observed the first several milliseconds of folding by monitoring the tryptophan fluorescence lifetime (60 µs dead time). Two folding intermediates are observed, the second of which is the previously characterized I(core) millisecond intermediate. The new earlier intermediate is likely on-pathway and appears to have long-range non-native structure, providing a rare example of such non-native structure formation in a folding pathway. The tryptophan fluorescence lifetimes also suggest a deviation from native packing in the second intermediate, I(core). Similar results from a fragment of RNase H demonstrate that only half of the protein is significantly involved in this early structure formation. These studies give us a view of the formation of tertiary structure on the folding pathway, which complements previous hydrogen-exchange studies that monitored only secondary structure and observed sequential native structure formation. Our results provide detailed folding information on both a timescale and a size-scale accessible to all-atom molecular dynamics simulations of protein folding.

Modulation of frustration in folding by sequence permutation.

Folding of globular proteins can be envisioned as the contraction of a random coil unfolded state toward the native state on an energy surface rough with local minima trapping frustrated species. These substructures impede productive folding and can serve as nucleation sites for aggregation reactions. However, little is known about the relationship between frustration and its underlying sequence determinants. Chemotaxis response regulator Y (CheY), a 129-amino acid bacterial protein, has been shown previously to populate an off-pathway kinetic trap in the microsecond time range. The frustration has been ascribed to premature docking of the N- and C-terminal subdomains or, alternatively, to the formation of an unproductive local-in-sequence cluster of branched aliphatic side chains, isoleucine, leucine, and valine (ILV). The roles of the subdomains and ILV clusters in frustration were tested by altering the sequence connectivity using circular permutations. Surprisingly, the stability and buried surface area of the intermediate could be increased or decreased depending on the location of the termini. Comparison with the results of small-angle X-ray-scattering experiments and simulations points to the accelerated formation of a more compact, on-pathway species for the more stable intermediate. The effect of chain connectivity in modulating the structures and stabilities of the early kinetic traps in CheY is better understood in terms of the ILV cluster model. However, the subdomain model captures the requirement for an intact N-terminal domain to access the native conformation. Chain entropy and aliphatic-rich sequences play crucial roles in biasing the early events leading to frustration in the folding of CheY.