Prediction of the structures of the plant-specific regions of vascular plant cellulose synthases and correlated functional analysis

Seed plants express cellulose synthase (CESA) protein isoforms with non-redundant functions, but how the isoforms function differently is unknown. Compared to bacterial cellulose synthases, CESAs have two insertions in the large cytosolic loop: the relatively well-conserved Plant Conserved Region (P-CR) and a Class Specific Region (CSR) that varies between CESAs. Absent any atomic structure of a plant CESA, we used ab initio protein structure prediction and molecular modeling to explore how these plant-specific regions may modulate CESA function. We modeled P-CR and CSR peptides from Arabidopsis thaliana CESAs representing the six clades of seed plant CESAs. As expected, the predicted wild type P-CR structures were similar. Modeling of the mutant P-CR of Atcesa8 (R362K) (fra6) suggested that changes in local structural stability and surface electrostatics may cause the mutant phenotype. Among CSRs within CESAs required for primary wall cellulose synthesis, the amino sequence and the modeled arrangement of helices was most similar in AtCESA1 and AtCESA3. Genetic complementation of known Arabidopsis mutants showed that the CSRs of AtCESA1 and AtCESA3 can function interchangeably in vivo. Analysis of protein surface electrostatics led to ideas about how the surface charges on CSRs may mediate protein-protein interactions. Refined modeling of the P-CR and CSR regions of GhCESA1 from cotton modified their tertiary structures, spatial relationships to the catalytic domain, and preliminary predictions about CESA oligomer formation. Cumulatively, the results provide structural clues about the function of plant-specific regions of CESA.