The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes.

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

Vectorial channeling as a mechanism for translational control by functional prions and condensates.

Translation of messenger RNA (mRNA) is regulated through a diverse set of RNA-binding proteins. A significant fraction of RNA-binding proteins contains prion-like domains which form functional prions. This raises the question of how prions can play a role in translational control. Local control of translation in dendritic spines by prions has been invoked in the mechanism of synaptic plasticity and memory. We show how channeling through diffusion and processive translation cooperate in highly ordered mRNA/prion aggregates as well as in less ordered mRNA/protein condensates depending on their substructure. We show that the direction of translational control, whether it is repressive or activating, depends on the polarity of the mRNA distribution in mRNA/prion assemblies which determines whether vectorial channeling can enhance recycling of ribosomes. Our model also addresses the effect of changes of substrate concentration in assemblies that have been suggested previously to explain translational control by assemblies through the introduction of a potential of mean force biasing diffusion of ribosomes inside the assemblies. The results from the model are compared with the experimental data on translational control by two functional RNA-binding prions, CPEB involved in memory and Rim4 involved in gametogenesis.

Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-ß fibrils.

Amyloid fibrillization is an exceedingly complex process in which incoming peptide chains bind to the fibril while concertedly folding. The coupling between folding and binding is not fully understood. We explore the molecular pathways of association of Aß40 monomers to fibril tips by combining time-resolved in situ scanning probe microscopy with molecular modeling. The comparison between experimental and simulation results shows that a complex supported by nonnative contacts is present in the equilibrium structure of the fibril tip and impedes fibril growth in a supersaturated solution. The unraveling of this frustrated state determines the rate of fibril growth. The kinetics of growth of freshly cut fibrils, in which the bulk fibril structure persists at the tip, complemented by molecular simulations, indicate that this frustrated complex comprises three or four monomers in nonnative conformations and likely is contained on the top of a single stack of peptide chains in the fibril structure. This pathway of fibril growth strongly deviates from the common view that the conformational transformation of each captured peptide chain is templated by the previously arrived peptide. The insights into the ensemble structure of the frustrated complex may guide the search for suppressors of Aß fibrillization. The uncovered dynamics of coupled structuring and assembly during fibril growth are more complex than during the folding of most globular proteins, as they involve the collective motions of several peptide chains that are not guided by a funneled energy landscape.

Mass Spectrometry of RNA-Binding Proteins during Liquid-Liquid Phase Separation Reveals Distinct Assembly Mechanisms and Droplet Architectures.

Liquid-liquid phase separation (LLPS) of heterogeneous ribonucleoproteins (hnRNPs) drives the formation of membraneless organelles, but structural information about their assembled states is still lacking. Here, we address this challenge through a combination of protein engineering, native ion mobility mass spectrometry, and molecular dynamics simulations. We used an LLPS-compatible spider silk domain and pH changes to control the self-assembly of the hnRNPs FUS, TDP-43, and hCPEB3, which are implicated in neurodegeneration, cancer, and memory storage. By releasing the proteins inside the mass spectrometer from their native assemblies, we could monitor conformational changes associated with liquid-liquid phase separation. We find that FUS monomers undergo an unfolded-to-globular transition, whereas TDP-43 oligomerizes into partially disordered dimers and trimers. hCPEB3, on the other hand, remains fully disordered with a preference for fibrillar aggregation over LLPS. The divergent assembly mechanisms revealed by ion mobility mass spectrometry of soluble protein species that exist under LLPS conditions suggest structurally distinct complexes inside liquid droplets that may impact RNA processing and translation depending on biological context.

A generalized Flory-Stockmayer kinetic theory of connectivity percolation and rigidity percolation of cytoskeletal networks.

Actin networks are essential for living cells to move, reproduce, and sense their environments. The dynamic and rheological behavior of actin networks is modulated by actin-binding proteins such as α-actinin, Arp2/3, and myosin. There is experimental evidence that actin-binding proteins modulate the cooperation of myosin motors by connecting the actin network. In this work, we present an analytical mean field model, using the Flory-Stockmayer theory of gelation, to understand how different actin-binding proteins change the connectivity of the actin filaments as the networks are formed. We follow the kinetics of the networks and estimate the concentrations of actin-binding proteins that are needed to reach connectivity percolation as well as to reach rigidity percolation. We find that Arp2/3 increases the actomyosin connectivity in the network in a non-monotonic way. We also describe how changing the connectivity of actomyosin networks modulates the ability of motors to exert forces, leading to three possible phases of the networks with distinctive dynamical characteristics: a sol phase, a gel phase, and an active phase. Thus, changes in the concentration and activity of actin-binding proteins in cells lead to a phase transition of the actin network, allowing the cells to perform active contraction and change their rheological properties.

A structural dynamics model for how CPEB3 binding to SUMO2 can regulate translational control in dendritic spines.

A prion-like RNA-binding protein, CPEB3, can regulate local translation in dendritic spines. CPEB3 monomers repress translation, whereas CPEB3 aggregates activate translation of its target mRNAs. However, the CPEB3 aggregates, as long-lasting prions, may raise the problem of unregulated translational activation. Here, we propose a computational model of the complex structure between CPEB3 RNA-binding domain (CPEB3-RBD) and small ubiquitin-like modifier protein 2 (SUMO2). Free energy calculations suggest that the allosteric effect of CPEB3-RBD/SUMO2 interaction can amplify the RNA-binding affinity of CPEB3. Combining with previous experimental observations on the SUMOylation mode of CPEB3, this model suggests an equilibrium shift of mRNA from binding to deSUMOylated CPEB3 aggregates to binding to SUMOylated CPEB3 monomers in basal synapses. This work shows how a burst of local translation in synapses can be silenced following a stimulation pulse, and explores the CPEB3/SUMO2 interplay underlying the structural change of synapses and the formation of long-term memories.

Exploring the interplay between fibrillization and amorphous aggregation channels on the energy landscapes of tau repeat isoforms.

Filaments made up of different isoforms of tau protein are associated with a variety of neurodegenerative diseases. Filaments made up of the 4R-tau isoform, which has four repeat regions (R1 to R4), are found in patients suffering from Alzheimer's disease, while filaments made of the 3R-tau isoform, which contains only three repeat units (R1, R3, and R4), are found in patients with Pick's disease (frontotemporal dementia). In this work, a predictive coarse-grained protein force field, the associative memory water-mediated structure and energy model (AWSEM), is used to study the energy landscapes of nucleation of the two different fibrils derived from patients with Pick's and Alzheimer's diseases. The landscapes for nucleating both fibril types contain amorphous oligomers leading to branched structures as well as prefibrillar oligomers. These two classes of oligomers differ in their structural details: The prefibrillar oligomers have more parallel in-register ß-strands, which ultimately lead to amyloid fibrils, while the amorphous oligomers are characterized by a near random ß-strand stacking, leading to a distinct amorphous phase. The landscape topography suggests that there must be significant structural reordering, or "backtracking," to transit from the amorphous aggregation channel to the fibrillization channel. Statistical mechanical perturbation theory allows us to evaluate the effects of changing concentration on the aggregation free-energy landscapes and to predict the effects of phosphorylation, which is known to facilitate the aggregation of tau repeats.

Braiding topology and the energy landscape of chromosome organization proteins.

Assemblies of structural maintenance of chromosomes (SMC) proteins and kleisin subunits are essential to chromosome organization and segregation across all kingdoms of life. While structural data exist for parts of the SMC-kleisin complexes, complete structures of the entire complexes have yet to be determined, making mechanistic studies difficult. Using an integrative approach that combines crystallographic structural information about the globular subdomains, along with coevolutionary information and an energy landscape optimized force field (AWSEM), we predict atomic-scale structures for several tripartite SMC-kleisin complexes, including prokaryotic condensin, eukaryotic cohesin, and eukaryotic condensin. The molecular dynamics simulations of the SMC-kleisin protein complexes suggest that these complexes exist as a broad conformational ensemble that is made up of different topological isomers. The simulations suggest a critical role for the SMC coiled-coil regions, where the coils intertwine with various linking numbers. The twist and writhe of these braided coils are coupled with the motion of the SMC head domains, suggesting that the complexes may function as topological motors. Opening, closing, and translation along the DNA of the SMC-kleisin protein complexes would allow these motors to couple to the topology of DNA when DNA is entwined with the braided coils.

Exploring the F-actin/CPEB3 interaction and its possible role in the molecular mechanism of long-term memory.

Dendritic spines are tiny membranous protrusions on the dendrites of neurons. Dendritic spines change shape in response to input signals, thereby strengthening the connections between neurons. The growth and stabilization of dendritic spines is thought to be essential for maintaining long-term memory. Actin cytoskeleton remodeling in spines is a key element of their formation and growth. More speculatively, the aggregation of CPEB3, a functional prion that binds RNA, has been reported to be involved in the maintenance of long-term memory. Here we study the interaction between actin and CPEB3 and propose a molecular model for the complex structure of CPEB3 and an actin filament (F-actin). The results of our computational modeling, including both energetic and structural analyses, are compared with novel data from peptide array experiments. Our model of the CPEB3/F-actin interaction suggests that F-actin potentially triggers the aggregation-prone structural transition of a short CPEB3 sequence by zipping it into a beta-hairpin form. We also propose that the CPEB3/F-actin interaction might be regulated by the SUMOylation of CPEB3, based on bioinformatic searches for potential SUMOylation sites as well as SUMO interacting motifs in CPEB3. On the basis of these results and the existing literature, we put forward a possible molecular mechanism underlying long-term memory that involves CPEB3's binding to actin, its aggregation, and its regulation by SUMOylation.

The role of the Arp2/3 complex in shaping the dynamics and structures of branched actomyosin networks.

Actomyosin networks give cells the ability to move and divide. These networks contract and expand while being driven by active energy-consuming processes such as motor protein walking and actin polymerization. Actin dynamics is also regulated by actin-binding proteins, such as the actin-related protein 2/3 (Arp2/3) complex. This complex generates branched filaments, thereby changing the overall organization of the network. In this work, the spatiotemporal patterns of dynamical actin assembly accompanying the branching-induced reorganization caused by Arp2/3 were studied using a computational model (mechanochemical dynamics of active networks [MEDYAN]); this model simulates actomyosin network dynamics as a result of chemical reactions whose rates are modulated by rapid mechanical equilibration. We show that branched actomyosin networks relax significantly more slowly than do unbranched networks. Also, branched networks undergo rare convulsive movements, "avalanches," that release strain in the network. These avalanches are associated with the more heterogeneous distribution of mechanically linked filaments displayed by branched networks. These far-from-equilibrium events arising from the marginal stability of growing actomyosin networks provide a possible mechanism of the "cytoquakes" recently seen in experiments.

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations.

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.

AWSEM-Suite: a protein structure prediction server based on template-guided, coevolutionary-enhanced optimized folding landscapes.

The accurate and reliable prediction of the 3D structures of proteins and their assemblies remains difficult even though the number of solved structures soars and prediction techniques improve. In this study, a free and open access web server, AWSEM-Suite, whose goal is to predict monomeric protein tertiary structures from sequence is described. The model underlying the server's predictions is a coarse-grained protein force field which has its roots in neural network ideas that has been optimized using energy landscape theory. Employing physically motivated potentials and knowledge-based local structure biasing terms, the addition of homologous template and co-evolutionary restraints to AWSEM-Suite greatly improves the predictive power of pure AWSEM structure prediction. From the independent evaluation metrics released in the CASP13 experiment, AWSEM-Suite proves to be a reasonably accurate algorithm for free modeling, standing at the eighth position in the free modeling category of CASP13. The AWSEM-Suite server also features a front end with a user-friendly interface. The AWSEM-Suite server is a powerful tool for predicting monomeric protein tertiary structures that is most useful when a suitable structure template is not available. The AWSEM-Suite server is freely available at: https://awsem.rice.edu.

Assemblies of calcium/calmodulin-dependent kinase II with actin and their dynamic regulation by calmodulin in dendritic spines.

Calcium/calmodulin-dependent kinase II (CaMKII) plays a key role in the plasticity of dendritic spines. Calcium signals cause calcium-calmodulin to activate CaMKII, which leads to remodeling of the actin filament (F-actin) network in the spine. We elucidate the mechanism of the remodeling by combining computer simulations with protein array experiments and electron microscopic imaging, to arrive at a structural model for the dodecameric complex of CaMKII with F-actin. The binding interface involves multiple domains of CaMKII. This structure explains the architecture of the micrometer-scale CaMKII/F-actin bundles arising from the multivalence of CaMKII. We also show that the regulatory domain of CaMKII may bind either calmodulin or F-actin, but not both. This frustration, along with the multipartite nature of the binding interface, allows calmodulin transiently to strip CaMKII from actin assemblies so that they can reorganize. This observation therefore provides a simple mechanism by which the structural dynamics of CaMKII establishes the link between calcium signaling and the morphological plasticity of dendritic spines.

Forging tools for refining predicted protein structures.

Refining predicted protein structures with all-atom molecular dynamics simulations is one route to producing, entirely by computational means, structural models of proteins that rival in quality those that are determined by X-ray diffraction experiments. Slow rearrangements within the compact folded state, however, make routine refinement of predicted structures by unrestrained simulations infeasible. In this work, we draw inspiration from the fields of metallurgy and blacksmithing, where practitioners have worked out practical means of controlling equilibration by mechanically deforming their samples. We describe a two-step refinement procedure that involves identifying collective variables for mechanical deformations using a coarse-grained model and then sampling along these deformation modes in all-atom simulations. Identifying those low-frequency collective modes that change the contact map the most proves to be an effective strategy for choosing which deformations to use for sampling. The method is tested on 20 refinement targets from the CASP12 competition and is found to induce large structural rearrangements that drive the structures closer to the experimentally determined structures during relatively short all-atom simulations of 50 ns. By examining the accuracy of side-chain rotamer states in subensembles of structures that have varying degrees of similarity to the experimental structure, we identified the reorientation of aromatic side chains as a step that remains slow even when encouraging global mechanical deformations in the all-atom simulations. Reducing the side-chain rotamer isomerization barriers in the all-atom force field is found to further speed up refinement.

Frustration and Direct-Coupling Analyses to Predict Formation and Function of Adeno-Associated Virus.

Adeno-associated virus (AAV) is a promising gene therapy vector because of its efficient gene delivery and relatively mild immunogenicity. To improve delivery target specificity, researchers use combinatorial and rational library design strategies to generate novel AAV capsid variants. These approaches frequently propose high proportions of nonforming or noninfective capsid protein sequences that reduce the effective depth of synthesized vector DNA libraries, thereby raising the discovery cost of novel vectors. We evaluated two computational techniques for their ability to estimate the impact of residue mutations on AAV capsid protein-protein interactions and thus predict changes in vector fitness, reasoning that these approaches might inform the design of functionally enriched AAV libraries and accelerate therapeutic candidate identification. The Frustratometer computes an energy function derived from the energy landscape theory of protein folding. Direct-coupling analysis (DCA) is a statistical framework that captures residue coevolution within proteins. We applied the Frustratometer to select candidate protein residues predicted to favor assembled or disassembled capsid states, then predicted mutation effects at these sites using the Frustratometer and DCA. Capsid mutants were experimentally assessed for changes in virus formation, stability, and transduction ability. The Frustratometer-based metric showed a counterintuitive correlation with viral stability, whereas a DCA-derived metric was highly correlated with virus transduction ability in the small population of residues studied. Our results suggest that coevolutionary models may be able to elucidate complex capsid residue-residue interaction networks essential for viral function, but further study is needed to understand the relationship between protein energy simulations and viral capsid metastability.

Cys-labeling kinetics of membrane protein GlpG: a role for specific SDS binding and micelle changes?

Empirically, α-helical membrane protein folding stability in surfactant micelles can be tuned by varying the mole fraction MFSDS of anionic (sodium dodecyl sulfate (SDS)) relative to nonionic (e.g., dodecyl maltoside (DDM)) surfactant, but we lack a satisfying physical explanation of this phenomenon. Cysteine labeling (CL) has thus far only been used to study the topology of membrane proteins, not their stability or folding behavior. Here, we use CL to investigate membrane protein folding in mixed DDM-SDS micelles. Labeling kinetics of the intramembrane protease GlpG are consistent with simple two-state unfolding-and-exchange rates for seven single-Cys GlpG variants over most of the explored MFSDS range, along with exchange from the native state at low MFSDS (which inconveniently precludes measurement of unfolding kinetics under native conditions). However, for two mutants, labeling rates decline with MFSDS at 0-0.2 MFSDS (i.e., native conditions). Thus, an increase in MFSDS seems to be a protective factor for these two positions, but not for the five others. We propose different scenarios to explain this and find the most plausible ones to involve preferential binding of SDS monomers to the site of CL (based on computational simulations) along with changes in size and shape of the mixed micelle with changing MFSDS (based on SAXS studies). These nonlinear impacts on protein stability highlights a multifaceted role for SDS in membrane protein denaturation, involving both direct interactions of monomeric SDS and changes in micelle size and shape along with the general effects on protein stability of changes in micelle composition.

Oleuropein derivatives from olive fruit extracts reduce α-synuclein fibrillation and oligomer toxicity.

Aggregation of α-synuclein (αSN) is implicated in neuronal degeneration in Parkinson's disease and has prompted searches for natural compounds inhibiting αSN aggregation and reducing its tendency to form toxic oligomers. Oil from the olive tree (Olea europaea L.) represents the main source of fat in the Mediterranean diet and contains variable levels of phenolic compounds, many structurally related to the compound oleuropein. Here, using αSN aggregation, fibrillation, size-exclusion chromatography-multiangle light scattering (SEC-MALS)-based assays, and toxicity assays, we systematically screened the fruit extracts of 15 different olive varieties to identify compounds that can inhibit αSN aggregation and oligomer toxicity and also have antioxidant activity. Polyphenol composition differed markedly among varieties. The variety with the most effective antioxidant and aggregation activities, Koroneiki, combined strong inhibition of αSN fibril nucleation and elongation with strong disaggregation activity on preformed fibrils and prevented the formation of toxic αSN oligomers. Fractionation of the Koroneiki extract identified oleuropein aglycone, hydroxyl oleuropein aglycone, and oleuropein as key compounds responsible for the differences in inhibition across the extracts. These phenolic compounds inhibited αSN amyloidogenesis by directing αSN monomers into small αSN oligomers with lower toxicity, thereby suppressing the subsequent fibril growth phase. Our results highlight the molecular consequences of differences in the level of effective phenolic compounds in different olive varieties, insights that have implications for long-term human health.

Topological constraints and modular structure in the folding and functional motions of GlpG, an intramembrane protease.

We investigate the folding of GlpG, an intramembrane protease, using perfectly funneled structure-based models that implicitly account for the absence or presence of the membrane. These two models are used to describe, respectively, folding in detergent micelles and folding within a bilayer, which effectively constrains GlpG's topology in unfolded and partially folded states. Structural free-energy landscape analysis shows that although the presence of multiple folding pathways is an intrinsic property of GlpG's modular functional architecture, the large entropic cost of organizing helical bundles in the absence of the constraining bilayer leads to pathways that backtrack (i.e., local unfolding of previously folded substructures is required when moving from the unfolded to the folded state along the minimum free-energy pathway). This backtracking explains the experimental observation of thermodynamically destabilizing mutations that accelerate GlpG's folding in detergent micelles. In contrast, backtracking is absent from the model when folding is constrained within a bilayer, the environment in which GlpG has evolved to fold. We also characterize a near-native state with a highly mobile transmembrane helix 5 (TM5) that is significantly populated under folding conditions when GlpG is embedded in a bilayer. Unbinding of TM5 from the rest of the structure exposes GlpG's active site, consistent with studies of the catalytic mechanism of GlpG that suggest that TM5 serves as a substrate gate to the active site.

Protein Frustratometer 2: a tool to localize energetic frustration in protein molecules, now with electrostatics.

The protein frustratometer is an energy landscape theory-inspired algorithm that aims at localizing and quantifying the energetic frustration present in protein molecules. Frustration is a useful concept for analyzing proteins' biological behavior. It compares the energy distributions of the native state with respect to structural decoys. The network of minimally frustrated interactions encompasses the folding core of the molecule. Sites of high local frustration often correlate with functional regions such as binding sites and regions involved in allosteric transitions. We present here an upgraded version of a webserver that measures local frustration. The new implementation that allows the inclusion of electrostatic energy terms, important to the interactions with nucleic acids, is significantly faster than the previous version enabling the analysis of large macromolecular complexes within a user-friendly interface. The webserver is freely available at URL: http://frustratometer.qb.fcen.uba.ar.

Cooperative folding of a polytopic α-helical membrane protein involves a compact N-terminal nucleus and nonnative loops.

Despite the ubiquity of helical membrane proteins in nature and their pharmacological importance, the mechanisms guiding their folding remain unclear. We performed kinetic folding and unfolding experiments on 69 mutants (engineered every 2-3 residues throughout the 178-residue transmembrane domain) of GlpG, a membrane-embedded rhomboid protease from Escherichia coli. The only clustering of significantly positive ϕ-values occurs at the cytosolic termini of transmembrane helices 1 and 2, which we identify as a compact nucleus. The three loops flanking these helices show a preponderance of negative ϕ-values, which are sometimes taken to be indicative of nonnative interactions in the transition state. Mutations in transmembrane helices 3-6 yielded predominantly ϕ-values near zero, indicating that this part of the protein has denatured-state-level structure in the transition state. We propose that loops 1-3 undergo conformational rearrangements to position the folding nucleus correctly, which then drives folding of the rest of the domain. A compact N-terminal nucleus is consistent with the vectorial nature of cotranslational membrane insertion found in vivo. The origin of the interactions in the transition state that lead to a large number of negative ϕ-values remains to be elucidated.