Genome editing retraces the evolution of toxin resistance in the monarch butterfly.

Identifying the genetic mechanisms of adaptation requires the elucidation of links between the evolution of DNA sequence, phenotype, and fitness1. Convergent evolution can be used as a guide to identify candidate mutations that underlie adaptive traits2-4, and new genome editing technology is facilitating functional validation of these mutations in whole organisms1,5. We combined these approaches to study a classic case of convergence in insects from six orders, including the monarch butterfly (Danaus plexippus), that have independently evolved to colonize plants that produce cardiac glycoside toxins6-11. Many of these insects evolved parallel amino acid substitutions in the α-subunit (ATPα) of the sodium pump (Na+/K+-ATPase)7-11, the physiological target of cardiac glycosides12. Here we describe mutational paths involving three repeatedly changing amino acid sites (111, 119 and 122) in ATPα that are associated with cardiac glycoside specialization13,14. We then performed CRISPR-Cas9 base editing on the native Atpα gene in Drosophila melanogaster flies and retraced the mutational path taken across the monarch lineage11,15. We show in vivo, in vitro and in silico that the path conferred resistance and target-site insensitivity to cardiac glycosides16, culminating in triple mutant 'monarch flies' that were as insensitive to cardiac glycosides as monarch butterflies. 'Monarch flies' retained small amounts of cardiac glycosides through metamorphosis, a trait that has been optimized in monarch butterflies to deter predators17-19. The order in which the substitutions evolved was explained by amelioration of antagonistic pleiotropy through epistasis13,14,20-22. Our study illuminates how the monarch butterfly evolved resistance to a class of plant toxins, eventually becoming unpalatable, and changing the nature of species interactions within ecological communities2,6-11,15,17-19.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

Oncogenic mutations on Rac1 affect global intrinsic dynamics underlying GTP and PAK1 binding.

Rac1 is a small member of the Rho GTPase family. One of the most important downstream effectors of Rac1 is a serine/threonine kinase, p21-activated kinase 1 (PAK1). Mutational activation of PAK1 by Rac1 has oncogenic signaling effects. Here, although we focus on Rac1-PAK1 interaction by atomic-force-microscopy-based single-molecule force spectroscopy experiments, we explore the effect of active mutations on the intrinsic dynamics and binding interactions of Rac1 by Gaussian network model analysis and molecular dynamics simulations. We observe that Rac1 oncogenic mutations are at the hinges of three global modes of motion, suggesting the mechanical changes as potential markers of oncogenicity. Indeed, the dissociation of wild-type Rac1-PAK1 complex shows two distinct unbinding dynamic states that are reduced to one with constitutively active Q61L and oncogenic Y72C mutant Rac1, as revealed by single-molecule force spectroscopy experiments. Q61L and Y72C mutations change the mechanics of the Rac1-PAK1 complex by increasing the elasticity of the protein and slowing down the transition to the unbound state. On the other hand, Rac1's intrinsic dynamics reveal more flexible GTP and PAK1-binding residues on switches I and II with Q61L, Y72C, oncogenic P29S and Q61R, and negative T17N mutations. The cooperativity in the fluctuations of GTP-binding sites around the p-loop and switch I decreases in all mutants, mostly in Q61L, whereas some PAK1-binding residues display enhanced coupling with GTP-binding sites in Q61L and Y72C and within each other in P29S. The predicted binding free energies of the modeled Rac1-PAK1 complexes show that the change in the dynamic behavior likely means a more favorable PAK1 interaction. Overall, these findings suggest that the active mutations affect intrinsic functional dynamic events and alter the mechanics underlying the binding of Rac1 to GTP and upstream and downstream partners including PAK1.

History, rare, and multiple events of mechanical unfolding of repeat proteins.

Mechanical unfolding of proteins consisting of repeat domains is an excellent tool to obtain large statistics. Force spectroscopy experiments using atomic force microscopy on proteins presenting multiple domains have revealed that unfolding forces depend on the number of folded domains (history) and have reported intermediate states and rare events. However, the common use of unspecific attachment approaches to pull the protein of interest holds important limitations to study unfolding history and may lead to discarding rare and multiple probing events due to the presence of unspecific adhesion and uncertainty on the pulling site. Site-specific methods that have recently emerged minimize this uncertainty and would be excellent tools to probe unfolding history and rare events. However, detailed characterization of these approaches is required to identify their advantages and limitations. Here, we characterize a site-specific binding approach based on the ultrastable complex dockerin/cohesin III revealing its advantages and limitations to assess the unfolding history and to investigate rare and multiple events during the unfolding of repeated domains. We show that this approach is more robust, reproducible, and provides larger statistics than conventional unspecific methods. We show that the method is optimal to reveal the history of unfolding from the very first domain and to detect rare events, while being more limited to assess intermediate states. Finally, we quantify the forces required to unfold two molecules pulled in parallel, difficult when using unspecific approaches. The proposed method represents a step forward toward more reproducible measurements to probe protein unfolding history and opens the door to systematic probing of rare and multiple molecule unfolding mechanisms.

Cullin neddylation may allosterically tune polyubiquitin chain length and topology.

Conjugation of Nedd8 (neddylation) to Cullins (Cul) in Cul-RING E3 ligases (CRLs) stimulates ubiquitination and polyubiquitination of protein substrates. CRL is made up of two Cul-flanked arms: one consists of the substrate-binding and adaptor proteins and the other consists of E2 and Ring-box protein (Rbx). Polyubiquitin chain length and topology determine the substrate fate. Here, we ask how polyubiquitin chains are accommodated in the limited space available between the two arms and what determines the polyubiquitin linkage topology. We focus on Cul5 and Rbx1 in three states: before Cul5 neddylation (closed state), after neddylation (open state), and after deneddylation, exploiting molecular dynamics simulations and the Gaussian Network Model. We observe that regulation of substrate ubiquitination and polyubiquitination takes place through Rbx1 rotations, which are controlled by Nedd8-Rbx1 allosteric communication. Allosteric propagation proceeds from Nedd8 via Cul5 dynamic hinges and hydrogen bonds between the C-terminal domain of Cul5 (Cul5CTD) and Rbx1 (Cul5CTD residues R538/R569 and Rbx1 residue E67, or Cul5CTD E474/E478/N491 and Rbx1 K105). Importantly, at each ubiquitination step (homogeneous or heterogeneous, linear or branched), the polyubiquitin linkages fit into the distances between the two arms, and these match the inherent CRL conformational tendencies. Hinge sites may constitute drug targets.

Allosteric Dynamic Control of Binding.

Proteins have a highly dynamic nature and there is a complex interrelation between their structural dynamics and binding behavior. By assuming various conformational ensembles, they perform both local and global fluctuations to interact with other proteins in a dynamic infrastructure adapted to functional motion. Here, we show that there is a significant association between allosteric mutations, which lead to high-binding-affinity changes, and the hinge positions of global modes, as revealed by a large-scale statistical analysis of data in the Structural Kinetic and Energetic Database of Mutant Protein Interactions (SKEMPI). We further examined the mechanism of allosteric dynamics by conducting studies on human growth hormone (hGH) and pyrin domain (PYD), and the results show how mutations at the hinge regions could allosterically affect the binding-site dynamics or induce alternative binding modes by modifying the ensemble of accessible conformations. The long-range dissemination of perturbations in local chemistry or physical interactions through an impact on global dynamics can restore the allosteric dynamics. Our findings suggest a mechanism for the coupling of structural dynamics to the modulation of protein interactions, which remains a critical phenomenon in understanding the effect of mutations that lead to functional changes in proteins.

Single-Molecule Force Spectroscopy: Experiments, Analysis, and Simulations.

The mechanical properties of cells and of subcellular components are important to obtain a mechanistic molecular understanding of biological processes. The quantification of mechanical resistance of cells and biomolecules using biophysical methods matured thanks to the development of nanotechnologies such as optical and magnetic tweezers, the biomembrane force probe, and atomic force microscopy (AFM). The quantitative nature of force spectroscopy measurements has converted AFM into a valuable tool in biophysics. Force spectroscopy allows the determination of the forces required to unfold protein domains and to disrupt individual receptor/ligand bonds. Molecular simulations as a computational microscope allow investigation of similar biological processes with an atomistic detail. In this chapter, we first provide a step-by-step protocol of force spectroscopy experiments using AFM, including sample preparation, measurements, and analysis and interpretation of the resulting dynamic force spectrum in terms of available theories. Next, we present the background for molecular dynamics (MD) simulations focusing on steered molecular dynamics (SMD) and the importance of bridging computational tools with experimental techniques.

High-speed force spectroscopy: microsecond force measurements using ultrashort cantilevers.

Complete understanding of the role of mechanical forces in biological processes requires knowledge of the mechanical properties of individual proteins and living cells. Moreover, the dynamic response of biological systems at the nano- and microscales span over several orders of magnitude in time, from sub-microseconds to several minutes. Thus, access to force measurements over a wide range of length and time scales is required. High-speed atomic force microscopy (HS-AFM) using ultrashort cantilevers has emerged as a tool to study the dynamics of biomolecules and cells at video rates. The adaptation of HS-AFM to perform high-speed force spectroscopy (HS-FS) allows probing protein unfolding and receptor/ligand unbinding up to the velocity of molecular dynamics (MD) simulations with sub-microsecond time resolution. Moreover, application of HS-FS on living cells allows probing the viscoelastic response at short time scales providing deep understanding of cytoskeleton dynamics. In this mini-review, we assess the principles and recent developments and applications of HS-FS using ultrashort cantilevers to probe molecular and cellular mechanics.

High-Speed Force Spectroscopy for Single Protein Unfolding.

Single-molecule force spectroscopy (SMFS) measurements allow for quantification of the molecular forces required to unfold individual protein domains. Atomic force microscopy (AFM) is one of the long-established techniques for force spectroscopy (FS). Although FS at conventional AFM pulling rates provides valuable information on protein unfolding, in order to get a more complete picture of the mechanism, explore new regimes, and combine and compare experiments with simulations, we need higher pulling rates and µs-time resolution, now accessible via high-speed force spectroscopy (HS-FS). In this chapter, we provide a step-by-step protocol of HS-FS including sample preparation, measurements and analysis of the acquired data using HS-AFM with an illustrative example on unfolding of a well-studied concatamer made of eight repeats of the titin I91 domain.

Structural and dynamics aspects of ASC speck assembly.

Activation of the inflammasome is accompanied by rapid formation of a micrometer-sized perinuclear structure called the ASC speck, a platform for caspase-1 activity. The ASC speck is often referred to as an aggregate and shares certain features with aggresomes. It is thus an open question whether the ASC speck formation takes place via nonspecific aggregation of hydrophobic patches or specific interactions of its domains; PYD and CARD, which belong to the death fold superfamily. Bringing together structure and dynamics studies using the Gaussian network model of PYD and CARD, and molecular dynamics simulations of the wild-type and in silico mutated PYD, with the mutational analysis on the ASC structure and its separate domains in human cells, we show that the ASC speck is an organized structure with at least two levels of distinct compaction mechanisms based on the specific interactions of PYD and CARD.