Elucidating the Protonation State of the γ-Secretase Catalytic Dyad

γ-Secretase (GS) is an intramembrane aspartyl protease that participates in the sequential cleavage of C99 to generate different isoforms of the amyloid-β (Aβ) peptides that are associated with the development of Alzheimer's disease. Due to its importance in the proteolytic processing of C99 by GS, we performed pH replica exchange molecular dynamics (pH-REMD) simulations of GS in its apo and substrate-bound forms to sample the protonation states of the catalytic dyad. We found that the catalytic dyad is deprotonated at physiological pH in our apo form, but the presence of the substrate at the active site displaces its monoprotonated state toward physiological pH. Our results show that Asp257 acts as the general base and Asp385 as the general acid during the cleavage mechanism. We identified different amino acids such as Lys265, Arg269, and the PAL motif interacting with the catalytic dyad and promoting changes in its acid-base behavior. Finally, we also found a significant pK_a shift of Glu280 related to the internalization of TM6-CT in the GS-apo form. Our study provides critical mechanistic insight into the GS mechanism and the basis for future research on the genesis of Aβ peptides and the development of Alzheimer's disease.

Cosolvent Simulations with Fragment-Bound Proteins Identify Hot Spots to Direct Lead Growth

In drug design, chemical groups are sequentially added to improve a weak-binding fragment into a tight-binding lead molecule. Often, the direction to make these additions is unclear, and there are numerous chemical modifications to choose. Lead development can be guided by crystal structures of the fragment-bound protein, but this alone is unable to capture structural changes like closing or opening of the binding site and any side-chain movements. Accounting for adaptation of the site requires a dynamic approach. Here, we use molecular dynamics calculations of small organic solvents with protein-fragment pairs to reveal the nearest "hot spots". These close hot spots show the direction to make appropriate additions and suggest types of chemical modifications that could improve binding affinity. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that is well established for finding binding "hot spots" in active sites and allosteric sites of proteins. We simulated 20 fragment-bound and apo forms of key pharmaceutical targets to map out hot spots for potential lead space. Furthermore, we analyzed whether the presence of a fragment facilitates the probes' binding in the lead space, a type of binding cooperativity. To the best of our knowledge, this is the first use of cosolvent MD conducted with bound inhibitors in the simulation. Our work provides a general framework to extract molecular features of binding sites to choose chemical groups for growing lead molecules. Of the 20 systems, 17 systems were well mapped by MixMD. For the three not-mapped systems, two had lead growth out into solution away from the protein, and the third had very small modifications which indicated no nearby hot spots. Therefore, our lack of mapping in three systems was appropriate given the experimental data (true-negative cases). The simulations are run for very short time scales, making this method tractable for use in the pharmaceutical industry.

pH Effects and Cooperativity among Key Titratable Residues for Escherichia coli Glycinamide Ribonucleotide Transformylase

Human glycinamide ribonucleotide transformylase (GAR Tfase) is a regulatory enzyme in the de novo purine biosynthesis pathway that has been extensively studied as an anticancer target. To some extent, inhibition of GAR Tfase selectively targets cancer cells over normal cells and inhibits purine formation and DNA replication. In this study, we investigated E. coli GAR Tfase, which shares high sequence similarity with the human GAR Tfase, and most functional residues are conserved. Herein, we aim to predict the pH-activity curve through a computational approach. We carried out pH-replica exchange molecular dynamics (pH-REMD) simulations to investigate pH-dependent functions such as structural changes, ligand binding, and catalytic activity. To compute the pH-activity curve, we identified the catalytic residues in specific protonation states, referred to as the catalytic competent protonation states (CCPS), which maintain the structure, keep ligands bound, and facilitate catalysis. Our computed population of CCPS with respect to pH matches well with the experimental pH-activity curve. To compute the microscopic pK_a values in the catalytically active conformation, we devised a thermodynamic model that considers the coupling between protonation states of CCPS residues and conformational states. These results allow us to correctly identify the general acid and base catalysts and interpret the pH-activity curve at an atomistic level.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

Predicting the binding modes and sites of metabolism of xenobiotics

Metabolism studies are an essential integral part of ADMET profiling of drug candidates to evaluate their safety and efficacy. Cytochrome P-450 (CYP) metabolizes a wide variety of xenobiotics/drugs. The binding modes of these compounds with CYP and their intrinsic reactivities decide the metabolic products. We report here a novel computational protocol, which comprises docking of ligands to heme-containing CYPs and prediction of binding energies through a newly developed scoring function, followed by analyses of the docked structures and molecular orbitals of the ligand molecules, for predicting the sites of metabolism (SOM) of ligands. The calculated binding free energies of 121 heme-containing protein-ligand docked complexes yielded a correlation coefficient of 0.84 against experiment. Molecular orbital analyses of the resultant top three unique poses of the docked complexes provided a success rate of 87% in identifying the experimentally known sites of metabolism of the xenobiotics. The SOM prediction methodology is freely accessible at .