Investigation of the induced-fit mechanism and catalytic activity of the human cytomegalovirus protease homodimer via molecular dynamics simulations

Human cytomegalovirus (HCMV) is a highly species-specific DNA virus infecting up to 80% of the general population. The viral genome contains the open reading frame UL80, which encodes the full-length 80 kDa HCMV serine protease and its substrate. Full-length HCMV protease is composed of an N-terminal 256-amino-acid proteolytic domain, called assemblin, a linker region, and a C-terminal structural domain, the assembly protein precursor. Biochemical studies have shown that dimerization activates assemblin because of an induced stabilization of the oxyanion hole (Arg166). Thus, we performed here molecular dynamics (MD) simulations on HCMV protease models to study the induced-fit mechanism of the enzyme upon the binding of substrates and peptidyl inhibitors, and structural and energetic factors that are responsible for the catalytic activity of the enzyme dimer. Long and stable trajectories were obtained for the models of the monomeric and dimeric states, free in solution and bound to a peptidyl-activated carbonyl inhibitor, with very good agreement between theoretical and experimental results. Our results suggest that HCMV protease is indeed a novel example of serine protease that operates by an induced-fit mechanism. Also, in agreement with mutagenesis studies, our MD simulations suggest that the dimeric form is necessary to activate the enzyme because of an induced stabilization of the oxyanion hole.