A phenylnorstatine inhibitor binding to HIV-1 protease: geometry, protonation, and subsite-pocket interactions analyzed at atomic resolution

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

PyRod: Tracing Water Molecules in Molecular Dynamics Simulations

Ligands entering a protein binding pocket essentially compete with water molecules for binding to the protein. Hence, the location and thermodynamic properties of water molecules in protein structures have gained increased attention in the drug design community. Including corresponding data into 3D pharmacophore modeling is essential for efficient high throughput virtual screening. Here, we present PyRod, a free and open-source Python software that allows for visualization of pharmacophoric binding pocket characteristics, identification of hot spots for ligand binding, and subsequent generation of pharmacophore features for virtual screening. The implemented routines analyze the protein environment of water molecules in molecular dynamics (MD) simulations and can differentiate between hydrogen bonded waters as well as waters in a protein environment of hydrophobic, charged, or aromatic atom groups. The gathered information is further processed to generate dynamic molecular interaction fields (dMIFs) for visualization and pharmacophoric features for virtual screening. The described software was applied to 5 therapeutically relevant drug targets, and generated pharmacophores were evaluated using DUD-E benchmarking sets. The best performing pharmacophore was found for the HIV1 protease with an early enrichment factor of 54.6. PyRod adds a new perspective to structure-based screening campaigns by providing easy-to-interpret dMIFs and purely protein-based pharmacophores that are solely based on tracing water molecules in MD simulations. Since structural information about cocrystallized ligands is not needed, screening campaigns can be followed, for which less or no ligand information is available. PyRod is freely available at https://github.com/schallerdavid/pyrod .

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein–ligand interactions and implicit COSMO solvation

Strictly uphill – in cognate docking experiments we show that a quantum mechanical description of interaction and solvation outperforms established scoring functions in sharply distinguishing the native state from decoy poses.

Atomic resolution crystal structure of Sapp2p, a secreted aspartic protease from Candida parapsilosis

The virulence of the Candida pathogens is enhanced by the production of secreted aspartic proteases, which therefore represent possible targets for drug design. Here, the crystal structure of the secreted aspartic protease Sapp2p from Candida parapsilosis was determined. Sapp2p was isolated from its natural source and crystallized in complex with pepstatin A, a classical aspartic protease inhibitor. The atomic resolution of 0.83 Å allowed the protonation states of the active-site residues to be inferred. A detailed comparison of the structure of Sapp2p with the structure of Sapp1p, the most abundant C. parapsilosis secreted aspartic protease, was performed. The analysis, which included advanced quantum-chemical interaction-energy calculations, uncovered molecular details that allowed the experimentally observed equipotent inhibition of both isoenzymes by pepstatin A to be rationalized.

Small molecule regulation of protein conformation by binding in the Flap of HIV protease

The fragment indole-6-carboxylic acid (1F1), previously identified as a flap site binder in a fragment-based screen against HIV protease (PR), has been cocrystallized with pepstatin-inhibited PR and with apo-PR. Another fragment, 3-indolepropionic acid (1F1-N), predicted by AutoDock calculations and confirmed in a novel inhibition of nucleation crystallization assay, exploits the same interactions in the flap site in two crystal structures. Both 1F1 and 1F1-N bind to the closed form of apo-PR and to pepstatin:PR. In solution, 1F1 and 1F1-N raise the Tm of apo-PR by 3.5-5 °C as assayed by differential scanning fluorimetry (DSF) and show equivalent low-micromolar binding constants to both apo-PR and pepstatin:PR, assayed by backscattering interferometry (BSI). The observed signal intensities in BSI are greater for each fragment upon binding to apo-PR than to pepstatin-bound PR, consistent with greater conformational change in the former binding event. Together, these data indicate that fragment binding in the flap site favors a closed conformation of HIV PR.

On the reliability of the corrected semiempirical quantum chemical method (PM6-DH2) for assigning the protonation states in HIV-1 protease/inhibitor complexes

A novel computational protocol for determining the most probable protonation states in protein/ligand complexes is presented. The method consists in treating large parts of the enzyme using the corrected semiempirical quantum chemical (QM) method – PM6-D2 for optimization and PM6-DH2 for single-point energies – while the rest is calculated using molecular mechanics (MM) within a hybrid QM/MM fashion. The surrounding solvent is approximated by an implicit model. This approach is applied to two model systems, two different carboxylate pairs in one general and one unique HIV-1 protease/inhibitor complex. The effect of the size of the movable QM part is investigated in a series of several sizes, 3-, 6-, 8- and 10-Å regions surrounding the inhibitor. For the smallest region (< 450 atoms) the computationally more costly DFT QM/MM optimizations are performed as a check of the correctness. Proton transfer (PT) phenomena occur at both the PM6-D2 and DFT levels, which underlines the requirement for a QM approach. The barriers of PT are checked in model carboxylic acid pairs using the highly accurate MP2 and CCSD(T) values. An important result of this study is the fine-tuning of the protocol which can be used in further applications; its limitations are also shown, pointing to future developments. The calculations reveal which protonation variants of the active site are the most stable. In conclusion, the presented protocol can also be utilized for defining probable isomers in biomolecular systems. It can also serve as a preparatory step for further interaction-energy and binding-score calculations.

Enzymatic and structural analysis of the I47A mutation contributing to the reduced susceptibility to HIV protease inhibitor lopinavir

Lopinavir (LPV) is a second-generation HIV protease inhibitor (PI) designed to overcome resistance development in patients undergoing long-term antiviral therapy. The mutation of isoleucine at position 47 of the HIV protease (PR) to alanine is associated with a high level of resistance to LPV. In this study, we show that recombinant PR containing a single I47A substitution has the inhibition constant (K(i) ) value for lopinavir by two orders of magnitude higher than for the wild-type PR. The addition of the I47A substitution to the background of a multiply mutated PR species from an AIDS patient showed a three-order-of-magnitude increase in K(i) in vitro relative to the patient PR without the I47A mutation. The crystal structure of I47A PR in complex with LPV showed the loss of van der Waals interactions in the S2/S2' subsites. This is caused by the loss of three side-chain methyl groups due to the I47A substitution and by structural changes in the A47 main chain that lead to structural changes in the flap antiparallel beta-strand. Furthermore, we analyzed possible interaction of the I47A mutation with secondary mutations V32I and I54V. We show that both mutations in combination with I47A synergistically increase the relative resistance to LPV in vitro. The crystal structure of the I47A/I54V PR double mutant in complex with LPV shows that the I54V mutation leads to a compaction of the flap, and molecular modeling suggests that the introduction of the I54V mutation indirectly affects the strain of the bound inhibitor in the PR binding cleft.

Computational analysis of plasmepsin IV bound to an allophenylnorstatine inhibitor

The plasmepsin proteases from the malaria parasite Plasmodium falciparum are attracting attention as putative drug targets. A recently published crystal structure of Plasmodium malariae plasmepsin IV bound to an allophenylnorstatine inhibitor [Clemente, J.C. et al. (2006) Acta Crystallogr. D 62, 246-252] provides the first structural insights regarding interactions of this family of inhibitors with plasmepsins. The compounds in this class are potent inhibitors of HIV-1 protease, but also show nM binding affinities towards plasmepsin IV. Here, we utilize automated docking, molecular dynamics and binding free energy calculations with the linear interaction energy LIE method to investigate the binding of allophenylnorstatine inhibitors to plasmepsin IV from two different species. The calculations yield excellent agreement with experimental binding data and provide new information regarding protonation states of active site residues as well as conformational properties of the inhibitor complexes.

Kinetic, stability, and structural changes in high-resolution crystal structures of HIV-1 protease with drug-resistant mutations L24I, I50V, and G73S

The crystal structures, dimer stabilities, and kinetics have been analyzed for wild-type human immunodeficiency virus type 1 (HIV-1) protease (PR) and resistant mutants PR(L24I), PR(I50V), and PR(G73S) to gain insight into the molecular basis of drug resistance. The mutations lie in different structural regions. Mutation I50V alters a residue in the flexible flap that interacts with the inhibitor, L24I alters a residue adjacent to the catalytic Asp25, and G73S lies at the protein surface far from the inhibitor-binding site. PR(L24I) and PR(I50V), showed a 4% and 18% lower k(cat)/K(m), respectively, relative to PR. The relative k(cat)/K(m) of PR(G73S) varied from 14% to 400% when assayed using different substrates. Inhibition constants (K(i)) of the antiviral drug indinavir for the reaction catalyzed by the mutant enzymes were about threefold and 50-fold higher for PR(L24I) and PR(I50V), respectively, relative to PR and PR(G73S). The dimer dissociation constant (K(d)) was estimated to be approximately 20 nM for both PR(L24I) and PR(I50V), and below 5 nM for PR(G73S) and PR. Crystal structures of the mutants PR(L24I), PR(I50V) and PR(G73S) were determined in complexes with indinavir, or the p2/NC substrate analog at resolutions of 1.10-1.50 Angstrom. Each mutant revealed distinct structural changes relative to PR. The mutated residues in PR(L24I) and PR(I50V) had reduced intersubunit contacts, consistent with the increased K(d) for dimer dissociation. Relative to PR, PR(I50V) had fewer interactions of Val50 with inhibitors, in agreement with the dramatically increased K(i). The distal mutation G73S introduced new hydrogen bond interactions that can transmit changes to the substrate-binding site and alter catalytic activity. Therefore, the structural alterations observed for drug-resistant mutations were in agreement with kinetic and stability changes.

Efficiency of a second-generation HIV-1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations

A subnanomolar inhibitor of human immunodeficiency virus type 1 (HIV-1) protease, designated QF34, potently inhibits the wild-type and drug-resistant enzyme. To explain its broad activity, the binding of QF34 to the wild-type HIV-1 protease is investigated by molecular dynamics simulations and compared to the binding of two inhibitors that are used clinically, saquinavir (SQV) and indinavir (IDV). Analysis of the flexibility of protease residues and inhibitor segments in the complex reveals that segments of QF34 were more mobile during the dynamics studies than the segments of SQV and IDV. The dynamics of hydrogen bonding show that QF34 forms a larger number of stable hydrogen bonds than the two inhibitors that are used clinically. Absolute binding free energies were calculated with molecular mechanics-generalized Born surface area (MM-GBSA) methodology using three protocols. The most consistent results were obtained using the single-trajectory approach, due to cancellation of errors and inadequate sampling in the separate-trajectory protocols. For all three inhibitors, energy components in favor of binding include van der Waals and electrostatic terms, whereas polar solvation and entropy terms oppose binding. Decomposition of binding energies reveals that more protease residues contribute significantly to the binding of QF34 than to the binding of SQV and IDV. Moreover, contributions from protease main chains and side chains are balanced in the case of QF34 (52:48 ratio, respectively), whereas side chain contributions prevail in both SQV and IDV (main-chain:side-chain ratios of 41:59 and 45:55, respectively). The presented results help explain the ability of QF34 to inhibit multiple resistant mutants and should be considered in the design of broad-specificity second-generation HIV-1 protease inhibitors.