Computational modelling of diatom silicic acid transporters predicts a conserved fold with implications for their function and evolution

Diatoms are an important group of algae that can produce intricate silicified cell walls (frustules). The complex process of silicification involves a set of enigmatic integral membrane proteins that are thought to actively transport the soluble precursor of biosilica, dissolved silicic acid. Full-length silicic acid transporters are found widely across the diatoms while homologous shorter proteins have now been identified in a range of other organisms. It has been suggested that modern silicic acid transporters arose from the union of such partial sequences. Here, we present a computational study of the silicic acid transporters and related transporter-like sequences to help understand the structure, function and evolution of this class of membrane protein. The AlphaFold software predicts that all of the protein sequences studied here share a common fold in the membrane domain which is entirely different from the predicted folds of non-homologous silicic acid transporters from plants. Substrate docking reveals how conserved polar residues could interact with silicic acid at a central solvent-accessible binding site, consistent with an alternating access mechanism of transport. The structural conservation between these proteins supports a model where modern silicon transporters evolved from smaller ancestral proteins by gene fusion.