A comparative QM/MM study of the reaction mechanism of the Hepatitis C virus NS3/NS4A protease with the three main natural substrates NS5A/5B, NS4B/5A and NS4A/4B

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

The Q41R mutation in the HCV-protease enhances the reactivity towards MAVS by suppressing non-reactive pathways

Recent experimental findings pointed out a new mutation in the HCV protease, Q41R, responsible for a significant enhancement of the enzyme's reactivity towards the mitochondrial antiviral-signaling protein (MAVS). The Q41R mutation is located rather far from the active site, and its involvement in the overall reaction mechanism is thus unclear. We used classical molecular dynamics and QM/MM to study the acylation reaction of HCV NS3/4A protease variants bound to MAVS and the NS4A/4B substrate and uncovered the indirect mechanism by which the Q41R mutation plays a critical role in the efficient cleavage of the substrate. Our simulations reveal that there are two major conformations of the MAVS H1'(p) residue for the wild type protease and only one conformation for the Q41R mutant. The conformational space of H1'(p) is restricted by the Q41R mutation due to a π-π stacking between H1'(p) and R41 as well as a strong hydrogen bond between the backbone of H57 and the side chain of R41. Further QM/MM calculations indicate that the complex with the conformation ruled out by the Q41R substitution is a non-reactive species due to its higher free energy barrier for the acylation reaction. Based on our calculations, we propose a kinetic mechanism that explains experimental data showing an increase of apparent rate constants for MAVS cleavage in Q41R mutants. Our model predicts that the non-reactive conformation of the enzyme-substrate complex modulates reaction kinetics like an uncompetitive inhibitor.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

Identification of novel small molecules as inhibitors of hepatitis C virus by structure-based virtual screening

Hepatitis C virus (HCV) NS3/NS4A serine protease is essential for viral replication, which is regarded as a promising drug target for developing direct-acting anti-HCV agents. In this study, sixteen novel compounds with cell-based HCV replicon activity ranging from 3.0 to 28.2 μM (IC₅₀) were successfully identified by means of structure-based virtual screening. Compound 5 and compound 11, with an IC₅₀ of 3.0 μM and 5.1 μM, respectively, are the two most potent molecules with low cytotoxicity.

Dissection of the difference between the group I metal ions in inhibiting GSK3β: a computational study

Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase that requires two cofactor Mg(2+) ions for catalysis in regulating many important cellular signals. Experimentally, Li(+) is a competitive inhibitor of GSK3β relative to Mg(2+), while this mechanism is not experienced with other group I metal ions. Herein, we use native Mg(2)(2+)-Mg(1)(2+) GSK3β and its Mg(2)(2+)-M(1)(+) (M = Li, Na, K, and Rb) derivatives to investigate the effect of metal ion substitution on the mechanism of inhibition through two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics (MD) simulations. The results of ONIOM calculations elucidate that the interaction of Na(+), K(+), and Rb(+) with ATP is weaker compared to that of Mg(2+) and Li(+) with ATP, and the critical triphosphate moiety of ATP undergoes a large conformational change in the Na(+), K(+), and Rb(+) substituted systems. As a result, the three metal ions (Na(+), K(+), and Rb(+)) are not stable and depart from the active site, while Mg(2+) and Li(+) can stabilize in the active site, evident in MD simulations. Comparisons of Mg(2)(2+)-Mg(1)(2+) and Mg(2)(2+)-Li(1)(+) systems reveal that the inline phosphor-transfer of ATP and the two conserved hydrogen bonds between Lys85 and ATP, together with the electrostatic potential at the Li(1)(+) site, are disrupted in the Mg(2)(2+)-Li(1)(+) system. These computational results highlight the possible mechanism why Li(+) inhibits GSK3β.