skip to main content


Title: In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non‐native interactions
Abstract

We examine the influence of cellular interactions in all‐atom models of a section of theHomo sapienscytoplasm on the early folding events of the three‐helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30‐microsecond simulations. We did however capture the formation of all three helices in a compact native‐like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta‐sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild‐type sequence. This study demonstrates that even fast‐folding proteins can get stuck in non‐native states in the cell, making them useful models for protein–chaperone interactions and early stages of aggregate formation relevant to cellular disease.

 
more » « less
Award ID(s):
2205665
NSF-PAR ID:
10480207
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Protein Science
Volume:
32
Issue:
11
ISSN:
0961-8368
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Proteins in the cellular milieu reside in environments crowded by macromolecules and other solutes. Although crowding can significantly impact the protein folded state stability, most experiments are conducted in dilute buffered solutions. To resolve the effect of crowding on protein stability, we use19F nuclear magnetic resonance spectroscopy to follow the reversible, two‐state unfolding thermodynamics of the N‐terminal Src homology 3 domain of theDrosophilasignal transduction protein drk in the presence of polyethylene glycols (PEGs) of various molecular weights and concentrations. Contrary to most current theories of crowding that emphasize steric protein–crowder interactions as the main driving force for entropically favored stabilization, our experiments show that PEG stabilization is accompanied by significant heat release, and entropy disfavors folding. Using our newly developed model, we find that stabilization by ethylene glycol and small PEGs is driven by favorable binding to the folded state. In contrast, for larger PEGs, chemical or soft PEG–protein interactions do not play a significant role. Instead, folding is favored by excluded volume PEG–protein interactions and an exothermic nonideal mixing contribution from release of confined PEG and water upon folding. Our results indicate that crowding acts through molecular interactions subtler than previously assumed and that interactions between solution components with both the folded and unfolded states must be carefully considered.

     
    more » « less
  2. Studies of proteins from one organism in another organism’s cells have shown that such exogenous proteins stick more, pointing toward coevolution of the cytoplasm and protein surface to minimize stickiness. Here we flip this question around by asking whether exogenous proteins can assemble efficiently into their target complexes in a non-native cytoplasm. We use as our model system the assembly of BtubA and BtubB from Prosthecobacter hosted in human U-2 OS cells. BtubA and B evolved from eukaryotic tubulins after horizontal gene transfer, but they have low surface sequence identity with the homologous human tubulins and do not respond to tubulin drugs such as nocodazole. In U-2 OS cells, BtubA and B assemble efficiently into dimers compared to in vitro, and the wild-type BtubA and B proteins subsequently are able to form microtubules as well. We find that generic crowding effects (Ficoll 70 in vitro) contribute significantly to efficient dimer assembly when compared to sticking interactions (U-2 OS cell lysate in vitro), consistent with the notion that a generic mechanism such as crowding can be effective at driving assembly of exogenous proteins, even when protein-cytoplasm quinary structure and sticking have been modified in a non-native cytoplasm. A simple Monte Carlo model of in vitro and in-cell interactions, treating BtubA and B as sticky dipoles in a matrix of sticky or nonsticky crowders, rationalizes all the experimental trends with two adjustable parameters and reveals nucleation as the likely mechanism for the time-scale separation between dimer- and tubule formation in-cell and in vitro. 
    more » « less
  3. Schneidman-Duhovny, Dina (Ed.)
    We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology. 
    more » « less
  4. Abstract

    Like many Gram‐negative pathogens,Shigellarely on a type three secretion system (T3SS) for injection of effector proteins directly into eukaryotic host cells to initiate and sustain infection. Protein secretion through the needle‐like type three secretion apparatus (T3SA) requires ATP hydrolysis by the T3SS ATPase Spa47, making it a likely target for in vivo regulation of T3SS activity and an attractive target for small molecule therapeutics against shigellosis. Here, we developed a model of an activated Spa47 homo‐hexamer, identifying two distinct regions at each protomer interface that we hypothesized to provide intermolecular interactions supporting Spa47 oligomerization and enzymatic activation. Mutational analysis and a series of high‐resolution crystal structures confirm the importance of these residues, as many of the engineered mutants are unable to form oligomers and efficiently hydrolyze ATP in vitro. Furthermore, in vivo evaluation ofShigellavirulence phenotype uncovered a strong correlation between T3SS effector protein secretion, host cell membrane disruption, and cellular invasion by the tested mutant strains, suggesting that perturbation of the identified interfacial residues/interactions influences Spa47 activity through preventing oligomer formation, which in turn regulatesShigellavirulence. The most impactful mutations are observed within the conserved Site 2 interface where the native residues support oligomerization and likely contribute to a complex hydrogen bonding network that organizes the active site and supports catalysis. The critical reliance on these conserved residues suggests that aspects of T3SS regulation may also be conserved, providing promise for the development of a cross‐species therapeutic that broadly targets T3SS ATPase oligomerization and activation.

     
    more » « less
  5. Abstract

    The liquid–liquid phase separation (LLPS) of Tau has been postulated to play a role in modulating the aggregation property of Tau, a process known to be critically associated with the pathology of a broad range of neurodegenerative diseases including Alzheimer's Disease. Taucan undergo LLPS by homotypic interaction through self‐coacervation (SC) or by heterotypic association through complex‐coacervation (CC) between Tau and binding partners such as RNA. What is unclear is in what way the formation mechanisms for self and complex coacervation of Tau are similar or different, and the addition of a binding partner to Tau alters the properties of LLPS and Tau. A combination ofin vitroexperimental and computational study reveals that the primary driving force for both Tau CC and SC is electrostatic interactions between Tau‐RNA or Tau‐Tau macromolecules. The liquid condensates formed by the complex coacervation of Tau and RNA have distinctly higher micro‐viscosity and greater thermal stability than that formed by the SC of Tau. Our study shows that subtle changes in solution conditions, including molecular crowding and the presence of binding partners, can lead to the formation of different types of Tau condensates with distinct micro‐viscosity that can coexist as persistent and immiscible entities in solution. We speculate that the formation, rheological properties and stability of Tau droplets can be readily tuned by cellular factors, and that liquid condensation of Tau can alter the conformational equilibrium of Tau.

     
    more » « less