Search for: All records

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. 

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. 

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. 

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. 

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. 

The phase problem in X-ray crystallography arises from the fact that only the intensities, and not the phases, of the diffracting electromagnetic waves are measured directly. Molecular replacement can often estimate the relative phases of reflections starting with those derived from a template structure, which is usually a previously solved structure of a similar protein. The key factor in the success of molecular replacement is finding a good template structure. When no good solved template exists, predicted structures based partially on templates can sometimes be used to generate models for molecular replacement, thereby extending the lower bound of structural and sequence similarity required for successful structure determination. Here, the effectiveness is examined of structures predicted by a state-of-the-art prediction algorithm, the Associative memory, Water-mediated, Structure and Energy Model Suite ( AWSEM-Suite ), which has been shown to perform well in predicting protein structures in CASP13 when there is no significant sequence similarity to a solved protein or only very low sequence similarity to known templates. The performance of AWSEM-Suite structures in molecular replacement is discussed and the results show that AWSEM-Suite performs well in providing useful phase information, often performing better than I-TASSER-MR and the previous algorithm AWSEM-Template .