Evolution of the substrate specificity of an RNA ligase ribozyme from phosphorimidazole to triphosphate activation.

The acquisition of new RNA functions through evolutionary processes was essential for the diversification of RNA-based primordial biology and its subsequent transition to modern biology. However, the mechanisms by which RNAs access new functions remain unclear. Do RNA enzymes need completely new folds to support new but related functions, or is reoptimization of the active site sufficient? What are the roles of neutral and adaptive mutations in evolutionary innovation? Here, we address these questions experimentally by focusing on the evolution of substrate specificity in RNA-catalyzed RNA assembly. We use directed in vitro evolution to show that a ligase ribozyme that uses prebiotically relevant 5'-phosphorimidazole-activated substrates can be evolved to catalyze ligation with substrates that are 5'-activated with the biologically relevant triphosphate group. Interestingly, despite catalyzing a related reaction, the new ribozyme folds into a completely new structure and exhibits promiscuity by catalyzing RNA ligation with both triphosphate and phosphorimidazole-activated substrates. Although distinct in sequence and structure, the parent phosphorimidazolide ligase and the evolved triphosphate ligase ribozymes can be connected by a series of point mutations where the intermediate sequences retain at least some ligase activity. The existence of a quasi-neutral pathway between these distinct ligase ribozymes suggests that neutral drift is sufficient to enable the acquisition of new substrate specificity, thereby providing opportunities for subsequent adaptive optimization. The transition from RNA-catalyzed RNA assembly using phosphorimidazole-activated substrates to triphosphate-activated substrates may have foreshadowed the later evolution of the protein enzymes that use monomeric triphosphates (nucleoside triphosphates, NTPs) for RNA synthesis.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center.

The mechanism by which genetic information was copied prior to the evolution of ribozymes is of great interest because of its importance to the origin of life. The most effective known process for the nonenzymatic copying of an RNA template is primer extension by a two-step pathway in which 2-aminoimidazole-activated nucleotides first react with each other to form an imidazolium-bridged intermediate that subsequently reacts with the primer. Reaction kinetics, structure-activity relationships, and X-ray crystallography have provided insight into the overall reaction mechanism, but many puzzles remain. In particular, high concentrations of Mg2+ are required for efficient primer extension, but the mechanism by which Mg2+ accelerates primer extension remains unknown. By analogy with the mechanism of DNA and RNA polymerases, a role for Mg2+ in facilitating the deprotonation of the primer 3'-hydroxyl is often assumed, but no catalytic metal ion is seen in crystal structures of the primer-extension complex. To explore the potential effects of Mg2+ binding in the reaction center, we performed atomistic molecular dynamics simulations of a series of modeled complexes in which a Mg2+ ion was placed in the reaction center with inner-sphere coordination with different sets of functional groups. Our simulations suggest that coordination of a Mg2+ ion with both O3' of the terminal primer nucleotide and the pro-Sp nonbridging oxygen of the reactive phosphate of an imidazolium-bridged dinucleotide would help to pre-organize the structure of the primer/template substrate complex to favor the primer-extension reaction. Our results suggest that the catalytic metal ion may play an important role in overcoming electrostatic repulsion between a deprotonated O3' and the reactive phosphate of the bridged dinucleotide and lead to testable predictions of the mode of Mg2+ binding that is most relevant to catalysis of primer extension.

Lipid Exchange Promotes Fusion of Model Protocells.

Vesicle fusion is an important process underlying cell division, transport, and membrane trafficking. In phospholipid systems, a range of fusogens including divalent cations and depletants have been shown to induce adhesion, hemifusion, and then full content fusion between vesicles. This work shows that these fusogens do not perform the same function for fatty acid vesicles, which are used as model protocells (primitive cells). Even when fatty acid vesicles appear adhered or hemifused to each other, the intervening barriers between vesicles do not rupture. This difference is likely because fatty acids have a single aliphatic tail, and are more dynamic than their phospholipid counterparts. To address this, it is postulated that fusion could instead occur under conditions, such as lipid exchange, that disrupt lipid packing. Using both experiments and molecular dynamics simulations, it is verified that fusion in fatty acid systems can indeed be induced by lipid exchange. These results begin to probe how membrane biophysics could constrain the evolutionary dynamics of protocells.

Interpreting the Evolutionary Echoes of a Protein Complex Essential for Inner-Ear Mechanosensation.

The sensory epithelium of the inner ear, found in all extant lineages of vertebrates, has been subjected to over 500 million years of evolution, resulting in the complex inner ear of modern vertebrates. Inner-ear adaptations are as diverse as the species in which they are found, and such unique anatomical variations have been well studied. However, the evolutionary details of the molecular machinery that is required for hearing are less well known. Two molecules that are essential for hearing in vertebrates are cadherin-23 and protocadherin-15, proteins whose interaction with one another acts as the focal point of force transmission when converting sound waves into electrical signals that the brain can interpret. This "tip-link" interaction exists in every lineage of vertebrates, but little is known about the structure or mechanical properties of these proteins in most non-mammalian lineages. Here, we use various techniques to characterize the evolution of this protein interaction. Results show how evolutionary sequence changes in this complex affect its biophysical properties both in simulations and experiments, with variations in interaction strength and dynamics among extant vertebrate lineages. Evolutionary simulations also characterize how the biophysical properties of the complex in turn constrain its evolution and provide a possible explanation for the increase in deafness-causing mutants observed in cadherin-23 relative to protocadherin-15. Together, these results suggest a general picture of tip-link evolution in which selection acted to modify the tip-link interface, although subsequent neutral evolution combined with varying degrees of purifying selection drove additional diversification in modern tetrapods.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations.

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

Elastic versus brittle mechanical responses predicted for dimeric cadherin complexes.

Cadherins are a superfamily of adhesion proteins involved in a variety of biological processes that include the formation of intercellular contacts, the maintenance of tissue integrity, and the development of neuronal circuits. These transmembrane proteins are characterized by ectodomains composed of a variable number of extracellular cadherin (EC) repeats that are similar but not identical in sequence and fold. E-cadherin, along with desmoglein and desmocollin proteins, are three classical-type cadherins that have slightly curved ectodomains and engage in homophilic and heterophilic interactions through an exchange of conserved tryptophan residues in their N-terminal EC1 repeat. In contrast, clustered protocadherins are straighter than classical cadherins and interact through an antiparallel homophilic binding interface that involves overlapped EC1 to EC4 repeats. Here we present molecular dynamics simulations that model the adhesive domains of these cadherins using available crystal structures, with systems encompassing up to 2.8 million atoms. Simulations of complete classical cadherin ectodomain dimers predict a two-phased elastic response to force in which these complexes first softly unbend and then stiffen to unbind without unfolding. Simulated α, ß, and γ clustered protocadherin homodimers lack a two-phased elastic response, are brittle and stiffer than classical cadherins and exhibit complex unbinding pathways that in some cases involve transient intermediates. We propose that these distinct mechanical responses are important for function, with classical cadherin ectodomains acting as molecular shock absorbers and with stiffer clustered protocadherin ectodomains facilitating overlap that favors binding specificity over mechanical resilience. Overall, our simulations provide insights into the molecular mechanics of single cadherin dimers relevant in the formation of cellular junctions essential for tissue function.

Structure of the ancient TRPY1 channel from Saccharomyces cerevisiae reveals mechanisms of modulation by lipids and calcium.

Transient receptor potential (TRP) channels emerged in fungi as mechanosensitive osmoregulators. The Saccharomyces cerevisiae vacuolar TRP yeast 1 (TRPY1) is the most studied TRP channel from fungi, but the structure and details of channel modulation remain elusive. Here, we describe the full-length cryoelectron microscopy structure of TRPY1 at 3.1 Å resolution in a closed state. The structure, despite containing an evolutionarily conserved and archetypical transmembrane domain, reveals distinctive structural folds for the cytosolic N and C termini, compared with other eukaryotic TRP channels. We identify an inhibitory phosphatidylinositol 3-phosphate (PI(3)P) lipid-binding site, along with two Ca2+-binding sites: a cytosolic site, implicated in channel activation and a vacuolar lumen site, implicated in inhibition. These findings, together with data from microsecond-long molecular dynamics simulations and a model of a TRPY1 open state, provide insights into the basis of TRPY1 channel modulation by lipids and Ca2+, and the molecular evolution of TRP channels.