Mutation-induced conformational remodeling explains gain- and loss-of-function phenotypes in Euglena gracilis phytoene synthase

Technical Summary: Mutation-Induced Conformational Remodeling in Euglena gracilis Phytoene Synthase

Problem Statement
Phytoene synthase (CrtB) catalyzes the rate-limiting first step in carotenoid biosynthesis, making it a critical target for metabolic engineering. While experimental screening in Euglena gracilis has identified specific mutations (Y197A, E121H, and L349F) that result in contrasting gain-of-function (enhanced carotenoid production) and loss-of-function (reduced production) phenotypes, the structural mechanisms driving these differences remain unclear. Static structural models and experimental observations alone fail to explain how single-residue substitutions alter enzymatic activity. The study addresses the gap in understanding the dynamic structural basis linking specific mutations to altered carotenoid productivity.

Methodology
The authors employed Molecular Dynamics (MD) simulations of the apo (substrate-free) form of E. gracilis CrtB to investigate the conformational dynamics of the wild-type enzyme and three variants: Y197A, E121H, and L349F. The computational approach included:

• Energetic Analysis: MM-GBSA and MM-PBSA calculations to assess mutation-dependent energetic perturbations and stability ($\Delta\Delta G$).
• Electrostatic Characterization: PROPKA analysis to determine shifts in isoelectric points (pI) and net charge at pH 7.0.
• Global Structural Metrics: Analysis of Solvent Accessible Surface Area (SASA) and Radius of Gyration ($R_g$) to evaluate global compactness and conformational sampling.
• Local Dynamics: Root Mean Square Fluctuation (RMSF) and $\Delta$RMSF analyses to map residue-specific mobility changes.
• Active Site Architecture: Detailed examination of solvent accessibility, binding pocket volumes, and hydrogen-bond networks within the Asp-rich motifs (essential for Mg²⁺ coordination) and the hydrophobic flap.

Key Results
The simulations revealed that the mutations induce distinct remodeling of the conformational energy landscape, affecting local interaction networks and active-site architecture without causing substantial global structural changes:

1. Energetic and Electrostatic Shifts:

   • E121H exhibited the largest redistribution of conformational dynamics, with positive $\Delta\Delta G$ values indicating energetic destabilization relative to the native protein. It also showed the highest predicted pI and a slight positive net charge at neutral pH.
   • Y197A displayed near-neutral to favorable energetic changes.
   • L349F induced modest energetic perturbations but showed pronounced destabilization in PROPKA folding free energy predictions.

2. Conformational Sampling and Flexibility:

   • Gain-of-function variants (Y197A, E121H): Both retained relatively accessible hydrophobic flaps. E121H showed the broadest $R_g$ distribution, while Y197A and E121H exhibited broader SASA distributions than the native protein, suggesting enhanced conformational sampling while maintaining overall compactness.
   • Loss-of-function variant (L349F): Closely resembled the native protein in global SASA and $R_g$ distributions but displayed a dramatic, localized constriction of the hydrophobic flap pocket (approx. 35% volume reduction).

3. Local Interaction Networks and Active Site Remodeling:

   • Hydrophobic Flap: The L349F mutation caused a significant reduction in flap pocket volume and solvent accessibility, consistent with a more closed conformation. In contrast, gain-of-function variants maintained greater flap exposure.
   • Asp-rich Motifs: All mutants showed increased hydrogen-bond occupancy in the Asp-rich 1 region. Notably, L349F displayed the highest occupancy despite reduced catalytic performance, suggesting that catalytic efficiency relies on a balance between stability and flexibility rather than maximal local stabilization.
   • Rigidification in L349F: Structural analysis of the L349F mutant revealed that residue Tyr127 adopted a more rigid conformation, forming a persistent interaction network involving $\pi$-cation and $\pi$-anion interactions with Arg154 and Glu158. This contrasts with the native protein and gain-of-function variants, where Tyr127 remained conformationally flexible.

Significance and Claims
The paper claims to provide a structural framework linking mutation-induced conformational remodeling around the catalytic center to experimentally observed differences in carotenoid production. The authors assert that the contrasting phenotypes arise not merely from local structural perturbations at the mutation sites, but from the redistribution of conformational sampling around functionally important motifs, specifically the hydrophobic flap and the adjacent Asp-rich 1 motif.

The study highlights that the loss-of-function phenotype in L349F is characterized by a "pronounced constriction" of the hydrophobic flap and a "rigid interaction network," whereas gain-of-function variants retain a "relatively accessible hydrophobic flap" and greater flexibility. These findings suggest that dynamic structural organization is a critical determinant of CrtB function. The authors conclude that these insights offer a basis for the rational engineering of phytoene synthases, though they explicitly note that future simulations incorporating substrates (GGPP) and Mg²⁺ ions are required to fully elucidate how these intrinsic conformational changes influence substrate recognition and catalysis.