Significance: The molecular basis of morphological and physiological adaptations in snakes is largely unknown. Here, we study these phenotypes using the genome of the Burmese python ( Python molurus bivittatus ), a model for extreme phenotypic plasticity and metabolic adaptation. We discovered massive rapid changes in gene expression that coordinate major changes in organ size and function after feeding. Many significantly responsive genes are associated with metabolism, development, and mammalian diseases. A striking number of genes experienced positive selection in ancestral snakes. Such genes were related to metabolism, development, lungs, eyes, heart, kidney, and skeletal structure—all highly modified features in snakes. Snake phenotypic novelty seems to be driven by the system-wide coordination of protein adaptation, gene expression, and changes in genome structure.

Boiga dendrophila (mangrove catsnake) is a colubrid snake that lives in Southeast Asian lowland rainforests and mangrove swamps and that preys primarily on birds. We have isolated, purified, and sequenced a novel toxin from its venom, which we named denmotoxin. It is a monomeric polypeptide of 77 amino acid residues with five disulfide bridges. In organ bath experiments, it displayed potent postsynaptic neuromuscular activity and irreversibly inhibited indirectly stimulated twitches in chick biventer cervicis nerve-muscle preparations. In contrast, it induced much smaller and readily reversible inhibition of electrically induced twitches in mouse hemidiaphragm nerve-muscle preparations. More precisely, the chick muscle alpha(1)betagammadelta-nicotinic acetylcholine receptor was 100-fold more susceptible compared with the mouse receptor. These data indicate that denmotoxin has a bird-specific postsynaptic activity. We chemically synthesized denmotoxin, crystallized it, and solved its crystal structure at 1.9 A by the molecular replacement method. The toxin structure adopts a non-conventional three-finger fold with an additional (fifth) disulfide bond in the first loop and seven additional residues at its N terminus, which is blocked by a pyroglutamic acid residue. This is the first crystal structure of a three-finger toxin from colubrid snake venom and the first fully characterized bird-specific toxin. Denmotoxin illustrates the relationship between toxin specificity and the primary prey type that constitutes the snake's diet.

Snakebite envenoming remains a devastating and neglected tropical disease, claiming over 100,000 lives annually and causing severe complications and long-lasting disabilities for many more^1,2. Three-finger toxins (3FTx) are highly toxic components of elapid snake venoms that can cause diverse pathologies, including severe tissue damage³ and inhibition of nicotinic acetylcholine receptors, resulting in life-threatening neurotoxicity⁴. At present, the only available treatments for snakebites consist of polyclonal antibodies derived from the plasma of immunized animals, which have high cost and limited efficacy against 3FTxs^5,6,7. Here we used deep learning methods to de novo design proteins to bind short-chain and long-chain α-neurotoxins and cytotoxins from the 3FTx family. With limited experimental screening, we obtained protein designs with remarkable thermal stability, high binding affinity and near-atomic-level agreement with the computational models. The designed proteins effectively neutralized all three 3FTx subfamilies in vitro and protected mice from a lethal neurotoxin challenge. Such potent, stable and readily manufacturable toxin-neutralizing proteins could provide the basis for safer, cost-effective and widely accessible next-generation antivenom therapeutics. Beyond snakebite, our results highlight how computational design could help democratize therapeutic discovery, particularly in resource-limited settings, by substantially reducing costs and resource requirements for the development of therapies for neglected tropical diseases.