The main protease (M pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M pro , a cysteine protease, have been determined, facilitating structure-based drug design. M pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N δ (HD) and N ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.
more »
« less
Conformational dynamics of cathepsin D and binding to a small‐molecule BACE1 inhibitor
BACE1 is a major therapeutic target for prevention and treatment of Alzheimer's disease. Developing inhibitors that can selectively target BACE1 in favor of other proteases, especially cathepsin D (CatD), has presented significant challenges. Here, we investigate the conformational dynamics and protonation states of BACE1 and CatD using continuous constant pH molecular dynamics with pH replica‐exchange sampling protocol. Despite similar structure, BACE1 and CatD exhibit markedly different active site dynamics. BACE1 displays pH‐dependent flap dynamics that controls substrate accessibility, while the CatD flap is relatively rigid and remains open in the pH range 2.5–6. Interestingly, although each protease hydrolyzes peptide bonds, the protonation states of the catalytic dyads are different within the active pH range. The acidic and basic components of the BACE1 catalytic dyad are clear, while either aspartic acid of the CatD catalytic dyad could play the role of acid or base. Finally, we investigate binding of the inhibitor LY2811376 developed by Eli Lilly to BACE1 and CatD. Surprisingly, in the enzyme active pH range, LY2811376 forms a stronger salt bridge with the catalytic dyad in CatD than in BACE1, which might explain the retinal toxicity of the inhibitor related to off‐target inhibition of CatD. This work highlights the complexity and challenge in structure‐based drug design where receptor‐ligand binding induces protonation state change in both the protein and the inhibitor. © 2017 Wiley Periodicals, Inc.
more » « less- NSF-PAR ID:
- 10035366
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 38
- Issue:
- 15
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- p. 1260-1269
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of Mpro, a cysteine protease, have been determined, facilitating structure-based drug design. Mpro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, Mpro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV Mpro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 Mpro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of Mpro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored Nδ (HD) and Nϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 Mpro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.more » « less
-
more » « less
Gram-negative bacteria expressing class A β-lactamases pose a serious health threat due to their ability to inactivate all β-lactam antibiotics. The acyl–enzyme intermediate is a central milestone in the hydrolysis reaction catalyzed by these enzymes. However, the protonation states of the catalytic residues in this complex have never been fully analyzed experimentally due to inherent difficulties. To help unravel the ambiguity surrounding class A β-lactamase catalysis, we have used ultrahigh-resolution X-ray crystallography and the recently approved β-lactamase inhibitor avibactam to trap the acyl–enzyme complex of class A β-lactamase CTX-M-14 at varying pHs. A 0.83-Å-resolution CTX-M-14 complex structure at pH 7.9 revealed a neutral state for both Lys73 and Glu166. Furthermore, the avibactam hydroxylamine-
O -sulfonate group conformation varied according to pH, and this conformational switch appeared to correspond to a change in the Lys73 protonation state at low pH. In conjunction with computational analyses, our structures suggest that Lys73 has a perturbed acid dissociation constant (pKa) compared with acyl–enzyme complexes with β-lactams, hindering its function to deprotonate Glu166 and the initiation of the deacylation reaction. Further NMR analysis demonstrated Lys73 pKato be ∼5.2 to 5.6. Together with previous ultrahigh-resolution crystal structures, these findings enable us to follow the proton transfer process of the entire acylation reaction and reveal the critical role of Lys73. They also shed light on the stability and reversibility of the avibactam carbamoyl acyl–enzyme complex, highlighting the effect of substrate functional groups in influencing the protonation states of catalytic residues and subsequently the progression of the reaction. -
William DeGrado (Ed.)more » « less
Producing novel enzymes that are catalytically active in vitro and biologically functional in vivo is a key goal of synthetic biology. Previously, we reported Syn-F4, the first de novo protein that meets both criteria. Syn-F4 hydrolyzed the siderophore ferric enterobactin, and expression of Syn-F4 allowed an inviable strain of
Escherichia coli (Δfes ) to grow in iron-limited medium. Here, we describe the crystal structure of Syn-F4. Syn-F4 forms a dimeric 4-helix bundle. Each monomer comprises two long α-helices, and the loops of the Syn-F4 dimer are on the same end of the bundle (syn topology). Interestingly, there is a penetrated hole in the central region of the Syn-F4 structure. Extensive mutagenesis experiments in a previous study showed that five residues (Glu26, His74, Arg77, Lys78, and Arg85) were essential for enzymatic activity in vivo. All these residues are located around the hole in the central region of the Syn-F4 structure, suggesting a putative active site with a catalytic dyad (Glu26–His74). The complete inactivity of purified proteins with mutations at the five residues supports the putative active site and reaction mechanism. Molecular dynamics and docking simulations of the ferric enterobactin siderophore binding to the Syn-F4 structure demonstrate the dynamic property of the putative active site. The structure and active site of Syn-F4 are completely different from native enterobactin esterase enzymes, thereby demonstrating that proteins designed de novo can provide life-sustaining catalytic activities using structures and mechanisms dramatically different from those that arose in nature. -
Abstract Ab-initio molecular dynamics enables following the dynamics of biological systems from the first principles, describing the electronic structure and offering the opportunity to “watch” the evolution of biochemical processes with unique resolution, beyond the capabilities of state-of-the-art experimental techniques. This article reports the role of first-principles ( ab-initio ) molecular dynamics (MD) in the CRISPR-Cas9 genome editing revolution, achieving a profound understanding of the enzymatic function and offering valuable insights for enzyme engineering. We introduce the methodologies and explain the use of ab-initio MD simulations to establish the two-metal dependent mechanism of DNA cleavage in the RuvC domain of the Cas9 enzyme, and how a second catalytic domain, HNH, cleaves the target DNA with the aid of a single metal ion. A detailed description of how ab-initio MD is combined with free-energy methods—i.e., thermodynamic integration and metadynamics—to break and form chemical bonds is given, explaining the use of these methods to determine the chemical landscape and establish the catalytic mechanism in CRISPR-Cas9. The critical role of classical methods is also discussed, explaining theory and application of constant pH MD simulations, used to accurately predict the catalytic residues’ protonation states. Overall, first-principles methods are shown to unravel the electronic structure and reveal the catalytic mechanism of the Cas9 enzyme, providing valuable insights that can serve for the design of genome editing tools with improved catalytic efficiency or controllable activity.more » « less
