Unravelling the Graded Millisecond Allosteric Activation Mechanism of Imidazole Glycerol Phosphate Synthase

Deciphering the molecular mechanisms of enzymatic allosteric regulation requires the structural characterization of key functional states and also their time evolution toward the formation of the allosterically activated ternary complex. The transient nature and usually slow millisecond timescale interconversion between these functional states hamper their detailed experimental and computational characterization. Here, we design a computational strategy tailored to reconstruct millisecond timescale events to describe the graded allosteric activation of imidazole glycerol phosphate synthase (IGPS) in the ternary complex. IGPS is a heterodimeric bienzyme complex responsible for the hydrolysis of glutamine to glutamate in the HisH subunit and delivering ammonia for the cyclase activity in HisF. Despite significant advances in understanding the underlying allosteric mechanism, essential molecular details of the long-range millisecond allosteric activation pathway of wild-type IGPS remain hidden. Without using a priori information of the active state, our simulations uncover how IGPS, with the allosteric effector bound in HisF, spontaneously captures glutamine in a catalytically inactive HisH conformation, subsequently attains a closed HisF:HisH interface, and finally forms the oxyanion hole in HisH for efficient glutamine hydrolysis. We show that effector binding in HisF dramatically decreases the conformational barrier associated with the oxyanion hole formation in HisH, in line with the experimentally observed 4500-fold activity increase in glutamine production. The formation of the allosterically active state is controlled by time-evolving dynamic communication networks connecting the effector and substrate binding sites. This computational strategy can be generalized to study other unrelated enzymes undergoing millisecond timescale allosteric transitions.