On the active site protonation state in aspartic proteases: implications for drug design

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.

Plasmepsin Inhibitors in Antimalarial Drug Discovery: Medicinal Chemistry and Target Validation (2000 to Present)

Plasmepsins represent novel antimalarial drug targets. However, plasmepsin-based antimalarial drug discovery efforts in the past 2 decades have generally suffered some drawbacks including lack of translatability of target inhibition to potent parasite inhibition in vitro and in vivo as well as poor selectivity over the related human aspartic proteases. Most studies reported in this period have over-relied on the use of hemoglobinase plasmepsins I-IV (particularly I and II) as targets for the new inhibitors even though these are known to be nonessential at the asexual stage of parasite development. Therefore, future antimalarial drug discovery efforts seeking to identify plasmepsin inhibitors should focus on incorporating non-hemoglobinase plasmepsins such as V, IX, and X in their screening in order to maximize chances of success. Additionally, there is need to go beyond just target enzymatic activity profiling to establishing cellular activity, physicochemical as well as drug metabolism and pharmacokinetics properties and finally in vivo proof-of-concept while ensuring selectivity over related human host proteases.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Cryo-temperature effects on membrane protein structure and dynamics

Cryo-electron structures revolutionize biology, yet cooling effects are unclear. Using a simulation protocol of hot, cold, and rapidly cooled γ-secretase we identify cryo-contraction and modes relevant to Aβ production and cryo-analysis in general.

G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket

The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.

Characterizing the structural ensemble of γ-secretase using a multiscale molecular dynamics approach

γ-Secretase is an intramembrane-cleaving aspartyl protease that plays an essential role in the processing of a variety of integral membrane proteins. Its role in the ultimate cleavage step in the processing of amyloid precursor protein to form amyloid-β (Aβ) peptide makes it an important therapeutic target in Alzheimer's disease research. Significant recent advances have been made in structural studies of this critical membrane protein complex. However, details of the mechanism of activation of the enzyme complex remain unclear. Using a multiscale computational modeling approach, combining multiple coarse-grained microsecond dynamic trajectories with all-atom models, the structure and two conformational states of the γ-secretase complex were evaluated. The transition between enzymatic state 1 and state 2 is shown to critically depend on the protonation states of the key catalytic residues Asp257 and Asp385 in the active site domain. The active site formation, related to our γ-secretase state 2, is observed to involve a concerted movement of four transmembrane helices from the catalytic subunit, resulting in the required localization of the catalytic residues. Global analysis of the structural ensemble of the enzyme complex was used to identify collective fluctuations important to the mechanism of substrate recognition and demonstrate that the corresponding fluctuations observed were uncorrelated with structural changes associated with enzyme activation. Overall, this computational study provides essential insight into the role of structure and dynamics in the activation and function of γ-secretase.

Constructing Interconsistent, Reasonable, and Predictive Models for Both the Kinetic and Thermodynamic Properties of HIV-1 Protease Inhibitors

Accumulated evidence suggests that the in vivo biological potency of a ligand is more strongly correlated with the binding/unbinding kinetics than the equilibrium thermodynamics of the protein-ligand interaction (PLI). However, the existing experimental and computational techniques are largely insufficient and limited in large-scale measurements or accurate predictions of the kinetic properties of PLI. In this work, elaborate efforts have been made to develop interconsistent, reasonable, and predictive models of the association rate constant (k_on), dissociation rate constant (k_off), and equilibrium dissociation constant (K_D) of a series of HIV protease inhibitors with different structural skeletons. The results showed that nine Volsurf descriptors derived from water (OH2) and hydrophobic (DRY) probes are key molecular determinants for the kinetic and thermodynamic properties of HIV-1 protease inhibitors. To the best of our knowledge, this is the first time that interconsistent and reasonable models with strong prediction power have been established for both the kinetic and thermodynamic properties of HIV protease inhibitors.

Biophysics in drug discovery: impact, challenges and opportunities

Over the past 25 years, biophysical technologies such as X-ray crystallography, nuclear magnetic resonance spectroscopy, surface plasmon resonance spectroscopy and isothermal titration calorimetry have become key components of drug discovery platforms in many pharmaceutical companies and academic laboratories. There have been great improvements in the speed, sensitivity and range of possible measurements, providing high-resolution mechanistic, kinetic, thermodynamic and structural information on compound-target interactions. This Review provides a framework to understand this evolution by describing the key biophysical methods, the information they can provide and the ways in which they can be applied at different stages of the drug discovery process. We also discuss the challenges for current technologies and future opportunities to use biophysical methods to solve drug discovery problems.

Unraveling HIV protease flaps dynamics by Constant pH Molecular Dynamics simulations

The active site of HIV protease (HIV-PR) is covered by two flaps. These flaps are known to be essential for the catalytic activity of the HIV-PR, but their exact conformations at the different stages of the enzymatic pathway remain subject to debate. Understanding the correct functional dynamics of the flaps might aid the development of new HIV-PR inhibitors. It is known that, the HIV-PR catalytic efficiency is pH-dependent, likely due to the influence of processes such as charge transfer and protonation/deprotonation of ionizable residues. Several Molecular Dynamics (MD) simulations have reported information about the HIV-PR flaps. However, in MD simulations the protonation of a residue is fixed and thus it is not possible to study the correlation between conformation and protonation state. To address this shortcoming, this work attempts to capture, through Constant pH Molecular Dynamics (CpHMD), the conformations of the apo, substrate-bound and inhibitor-bound HIV-PR, which differ drastically in their flap arrangements. The results show that the HIV-PR flaps conformations are defined by the protonation of the catalytic residues Asp25/Asp25' and that these residues are sensitive to pH changes. This study suggests that the catalytic aspartates can modulate the opening of the active site and substrate binding.

Influence of surface charge, binding site residues and glycosylation on Thielavia terrestris cutinase biochemical characteristics

Cutinases are esterases of industrial importance for applications in recycling and surface modification of polyesters. The cutinase from Thielavia terrestris (TtC) is distinct in terms of its ability to retain its stability and activity in acidic pH. Stability and activity in acidic pHs are desirable for esterases as the pH of the reaction tends to go down with the generation of acid. The pH stability and activity are governed by the charged state of the residues involved in catalysis or in substrate binding. In this study, we performed the detailed structural and biochemical characterization of TtC coupled with surface charge analysis to understand its acidic tolerance. The stability of TtC in acidic pH was rationalized by evaluating the contribution of charge interactions to the Gibbs free energy of unfolding at varying pHs. The activity of TtC was found to be limited by substrate binding affinity, which is a function of the surface charge. Additionally, the presence of glycosylation affects the biochemical characteristics of TtC owing to steric interactions with residues involved in substrate binding.

Conformational Dynamics and Binding Free Energies of Inhibitors of BACE-1: From the Perspective of Protonation Equilibria

BACE-1 is the β-secretase responsible for the initial amyloidogenesis in Alzheimer’s disease, catalyzing hydrolytic cleavage of substrate in a pH-sensitive manner. The catalytic mechanism of BACE-1 requires water-mediated proton transfer from aspartyl dyad to the substrate, as well as structural flexibility in the flap region. Thus, the coupling of protonation and conformational equilibria is essential to a full in silico characterization of BACE-1. In this work, we perform constant pH replica exchange molecular dynamics simulations on both apo BACE-1 and five BACE-1-inhibitor complexes to examine the effect of pH on dynamics and inhibitor binding properties of BACE-1. In our simulations, we find that solution pH controls the conformational flexibility of apo BACE-1, whereas bound inhibitors largely limit the motions of the holo enzyme at all levels of pH. The microscopic pK_a values of titratable residues in BACE-1 including its aspartyl dyad are computed and compared between apo and inhibitor-bound states. Changes in protonation between the apo and holo forms suggest a thermodynamic linkage between binding of inhibitors and protons localized at the dyad. Utilizing our recently developed computational protocol applying the binding polynomial formalism to the constant pH molecular dynamics (CpHMD) framework, we are able to obtain the pH-dependent binding free energy profiles for various BACE-1-inhibitor complexes. Our results highlight the importance of correctly addressing the binding-induced protonation changes in protein-ligand systems where binding accompanies a net proton transfer. This work comprises the first application of our CpHMD-based free energy computational method to protein-ligand complexes and illustrates the value of CpHMD as an all-purpose tool for obtaining pH-dependent dynamics and binding free energies of biological systems.

Formation of insoluble amyloid plaques in the vascular and hippocampal areas of the brain characterizes Alzheimer’s disease, a devastating neurodegenerative disorder causing dementia. Site-specific hydrolytic catalysis of β-secretase, or BACE-1, is responsible for production of oligomerative amyloid β-peptide. As the catalytic activity of BACE-1 is pH-dependent and its structural dynamics are intrinsic to the catalysis, we examine the dependence of dynamics of BACE-1 on solution pH and its implications on the catalytic mechanism of BACE-1. Also, we highlight the importance of accurate description of protonation states of the titratable groups in computer-aided drug discovery targeting BACE-1. We hope the understanding of pH dependence of the dynamics and inhibitor binding properties of BACE-1 will aid the structure-based inhibitor design efforts against Alzheimer’s disease.

Context-dependent GluN2B-selective inhibitors of NMDA receptor function are neuroprotective with minimal side effects

Stroke remains a significant problem despite decades of work on neuroprotective strategies. NMDA receptor (NMDAR) antagonists are neuroprotective in preclinical models, but have been clinically unsuccessful, in part due to side effects. Here we describe a prototypical GluN2B-selective antagonist with an IC50 value that is 10-fold more potent at acidic pH 6.9 associated with ischemic tissue compared to pH 7.6, a value close to the pH in healthy brain tissue. This should maximize neuroprotection in ischemic tissue while minimizing on-target side effects associated with NMDAR blockade in noninjured brain regions. We have determined the mechanism underlying pH-dependent inhibition and demonstrate the utility of this approach in vivo. We also identify dicarboxylate dimers as a novel proton sensor in proteins. These results provide insight into the molecular basis of pH-dependent neuroprotective NMDAR block, which could be beneficial in a wide range of neurological insults associated with tissue acidification.

Computer-aided structure-based design of multitarget leads for Alzheimer's disease

Alzheimer's disease is a neurodegenerative pathology with unmet clinical needs. A highly desirable approach to this syndrome would be to find a single lead that could bind to some or all of the selected biomolecules that participate in the amyloid cascade, the most accepted route for Alzheimer disease genesis. In order to circumvent the challenge posed by the sizable differences in the binding sites of the molecular targets, we propose a computer-assisted protocol based on a pharmacophore and a set of required interactions with the targets that allows for the automated screening of candidates. We used a combination of docking and molecular dynamics protocols in order to discard nonbinders, optimize the best candidates, and provide a rationale for their potential as inhibitors. To provide a proof of concept, we proceeded to screen the literature and databases, a task that allowed us to identify a set of carbazole-containing compounds that initially showed affinity only for the cholinergic targets in our experimental assays. Two cycles of design based on our protocol led to a new set of analogues that were synthesized and assayed. The assay results revealed that the designed inhibitors had improved affinities for BACE-1 by more than 3 orders of magnitude and also displayed amyloid aggregation inhibition and affinity for AChE and BuChE, a result that led us to a group of multitarget amyloid cascade inhibitors that also could have a positive effect at the cholinergic level.

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (φ, ψ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C(β) atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (φ, ψ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (φ, ψ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C(β) atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.