Fast and Reliable Thermodynamic Approach for Determining the Protonation State of the Asp Dyad

Dynamical Basis of Allosteric Activation for the Plasmodium falciparum Protein Kinase G

The Plasmodium falciparum cGMP-dependent protein kinase (PfPKG) is required for the progression of the Plasmodium's life cycle and is therefore a promising malaria drug target. PfPKG includes four cGMP-binding domains (CBD-A to CBD-D). CBD-D plays a crucial role in PfPKG regulation as it is the primary determinant for the inhibition and cGMP-dependent activation of the catalytic domain. Hence, it is critical to understand how CBD-D is allosterically regulated by cGMP. Although the apo versus holo conformational changes of CBD-D have been reported, information on the intermediates of the activation pathway is currently lacking. Here, we employed molecular dynamics simulations to model four key states along the thermodynamic cycle for the cGMP-dependent activation of the PfPKG CBD-D domain. The simulations were compared to NMR data, and they revealed that the PfPKG CBD-D activation pathway samples a compact intermediate in which the N- and C-terminal helices approach the central β-barrel. In addition, by comparing the cGMP-bound active and inactive states, the essential binding interactions that differentiate these states were identified. The identification of structural and dynamical features unique to the cGMP-bound inactive state provides a promising basis to design PfPKG-selective allosteric inhibitors as a viable treatment for malaria.

Mechanism of allosteric inhibition in the Plasmodium falciparum cGMP-dependent protein kinase

Most malaria deaths are caused by the protozoan parasite Plasmodium falciparum Its life cycle is regulated by a cGMP-dependent protein kinase (PfPKG), whose inhibition is a promising antimalaria strategy. Allosteric kinase inhibitors, such as cGMP analogs, offer enhanced selectivity relative to competitive kinase inhibitors. However, the mechanisms underlying allosteric PfPKG inhibition are incompletely understood. Here, we show that 8-NBD-cGMP is an effective PfPKG antagonist. Using comparative NMR analyses of a key regulatory domain, PfD, in its apo, cGMP-bound, and cGMP analog-bound states, we elucidated its inhibition mechanism of action. Using NMR chemical shift analyses, molecular dynamics simulations, and site-directed mutagenesis, we show that 8-NBD-cGMP inhibits PfPKG not simply by reverting a two-state active versus inactive equilibrium, but by sampling also a distinct inactive "mixed" intermediate. Surface plasmon resonance indicates that the ability to stabilize a mixed intermediate provides a means to effectively inhibit PfPKG, without losing affinity for the cGMP analog. Our proposed model may facilitate the rational design of PfPKG-selective inhibitors for improved management of malaria.

Proton transfer and drug binding details revealed in neutron diffraction studies of wild-type and drug resistant HIV-1 protease

HIV-1 protease is an essential therapeutic target for the design and development of antiviral inhibitors to treat AIDS. We used room temperature neutron crystallography to accurately determine hydrogen atom positions in several protease complexes with clinical drugs, amprenavir and darunavir. Hydrogen bonding interactions were carefully mapped to provide an unprecedented picture of drug binding to the protease target. We demonstrate that hydrogen atom positions within the enzyme catalytic site can be altered by introducing drug resistant mutations and by protonating surface residues that trigger proton transfer reactions between the catalytic Asp residues and the hydroxyl group of darunavir. When protein perdeuteration is not feasible, we validate the use of initial H/D exchange with unfolded protein and partial deuteration in pure D2O with hydrogenous glycerol to maximize deuterium incorporation into the protein, with no detrimental effects on the growth of quality crystals suitable for neutron diffraction experiments.

Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling

Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR.

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.

Consistent Prediction of Mutation Effect on Drug Binding in HIV-1 Protease Using Alchemical Calculations

Despite a large number of antiretroviral drugs targeting HIV-1 protease for inhibition, mutations in this protein during the course of patient treatment can render them inefficient. This emerging resistance inspired numerous computational studies of the HIV-1 protease aimed at predicting the effect of mutations on drug binding in terms of free binding energy Δ G, as well as in mechanistic terms. In this study, we analyze ten different protease-inhibitor complexes carrying major resistance-associated mutations (RAMs) G48V, I50V, and L90M using molecular dynamics simulations. We demonstrate that alchemical free energy calculations can consistently predict the effect of mutations on drug binding. By explicitly probing different protonation states of the catalytic aspartic dyad, we reveal the importance of the correct choice of protonation state for the accuracy of the result. We also provide insight into how different mutations affect drug binding in their specific ways, with the unifying theme of how all of them affect the crucial drug binding regions of the protease.

Determination of the Bridging Ligand in the Active Site of Tyrosinase

Tyrosinase is a type-3 copper enzyme that is widely distributed in plants, fungi, insects, and mammals. Developing high potent inhibitors against tyrosinase is of great interest in diverse fields including tobacco curing, food processing, bio-insecticides development, cosmetic development, and human healthcare-related research. In the crystal structure of Agaricus bisporus mushroom tyrosinase, there is an oxygen atom bridging the two copper ions in the active site. It is unclear whether the identity of this bridging oxygen is a water molecule or a hydroxide anion. In the present study, we theoretically determine the identity of this critical bridging oxygen by performing first-principles hybrid quantum mechanics/molecular mechanics/Poisson-Boltzmann-surface area (QM/MM-PBSA) calculations along with a thermodynamic cycle that aim to improve the accuracy. Our results show that the binding with water molecule is energy favored and the QM/MM-optimized structure is very close to the crystal structure, whereas the binding with hydroxide anions causes the increase of energy and significant structural changes of the active site, indicating that the identity of the bridging oxygen must be a water molecule rather than a hydroxide anion. The different binding behavior between water and hydroxide anions may explain why molecules with a carboxyl group or too many negative charges have lower inhibitory activity. In light of this, the design of high potent active inhibitors against tyrosinase should satisfy both the affinity to the copper ions and the charge neutrality of the entire molecule.