Thermodynamic dissection of the binding energetics of KNI-272, a potent HIV-1 protease inhibitor

Phage display identification of high-affinity ligands for human TSG101-UEV: A structural and thermodynamic study of PTAP recognition

The ubiquitin E2 variant domain of TSG101 (TSG101-UEV) plays a pivotal role in protein sorting and virus budding by recognizing PTAP motifs within ubiquitinated proteins. Disrupting TSG101-UEV/PTAP interactions has emerged as a promising strategy for the development of novel host-oriented antivirals with a broad spectrum of action. Nonetheless, finding inhibitors with good properties as therapeutic agents remains a challenge since the key determinants of binding affinity and specificity are still poorly understood. Here we present a detailed thermodynamic, structural, and dynamic characterization viral PTAP Late domain recognition by TSG101-UEV, combining isothermal titration calorimetry, X-ray diffraction structural studies, molecular dynamics simulations, and computational analysis of intramolecular communication pathways. Our analysis highlights key contributions from conserved hydrophobic contacts and water-mediated hydrogen bonds at the PTAP binding interface. We have identified additional electrostatic hotspots adjacent to the core motif that modulate affinity. Using competitive phage display screening we have improved affinity by 1-2 orders of magnitude, producing novel peptides with low micromolar affinities that combine critical elements found in the best natural binders. Molecular dynamics simulations revealed that optimized peptides engage new pockets on the UEV domain surface. This study provides a comprehensive view of the molecular forces directing TSG101-UEV recognition of PTAP motifs, revealing that binding is governed by conserved structural elements yet tuneable through targeted optimization. These insights open new venues to design inhibitors targeting TSG101-dependent pathways with potential application as novel broad-spectrum antivirals.

Site-directed mutagenesis at the Glu78 in Ec-NhaA transporter impacting ion exchange: a biophysical study

Na+/H+ antiporters facilitate the exchange of Na+ for H+ across the cytoplasmic membrane in prokaryotic and eukaryotic cells. These transporters are crucial to maintain the homeostasis of sodium ions, consequently pH, and volume of the cells. Therefore, sodium/proton antiporters are considered promising therapeutic targets in humans. The Na+/H+ antiporter in Escherichia coli (Ec-NhaA), a prototype of cation-proton antiporter (CPA) family, transports two protons and one sodium (or Li+) in opposite direction. Previous mutagenesis experiments on Ec-NhaA have proposed Asp164, Asp163, and Asp133 amino acids with the significant implication in functional and structural integrity and create site for ion-binding. However, the mechanism and the sites for the binding of the two protons remain unknown and controversial which could be critical for pH regulation. In this study, we have explored the role of Glu78 in the regulation of pH by Ec-NhaA. Although we have created various mutants, E78C has shown a considerable effect on the stoichiometry of NhaA and presented comparable phenotypes. The ITC experiment has shown the binding of ~ 5 protons in response to the transport of one lithium ion. The phenotype analysis on selective medium showed a significant expression compared to WT Ec-NhaA. This represents the importance of Glu78 in transporting the H+ across the membrane where a single mutation with Cys amino acid alters the number of H+ significantly maintaining the activity of the protein.

Direct Determination of High-Affinity Binding Constants by Continuous Injection Isothermal Titration Calorimetry

Isothermal titration calorimetry (ITC) is a method to determine thermodynamic values (ΔG, ΔH, and ΔS) for ligand-receptor binding in biological and abiological systems. It is challenging to directly determine subnanomolar dissociation constants using a standard incremental injection approach ITC (IIA-ITC) measurement. We recently demonstrated a continuous injection approach ITC (CIA-ITC) [ J. Phys. Chem. B 2021, 125, 8075-8087]enables the estimation of thermodynamic parameters in situ. In this work, we demonstrate a label-free and surface modification-free CIA-ITC to determine the complete binding thermodynamics of a ligand with a subnanomolar dissociation constant KD. The KD for desthiobiotin (DTB)-avidin binding was determined to be 6.5 pM with respect to the ligand by CIA-ITC, a quantity unsuccessfully measured with IIA-ITC and surface plasmon resonance spectroscopy (SPR). This value compares well with literature-reported spectroscopic determination of DTB-avidin binding. Criteria with respect to the concentration of the ligand and receptor and flow rate for obtaining true equilibrium dissociation constants without displacement titration are presented.

Biophysical Approaches for the Characterization of Protein-Metabolite Interactions

With an estimate of hundred thousands of protein molecules per cell and the number of metabolites several orders of magnitude higher, protein-metabolite interactions are omnipresent. In vitro analyses are one of the main pillars on the way to establish a solid understanding of how these interactions contribute to maintaining cellular homeostasis. A repertoire of biophysical techniques is available by which protein-metabolite interactions can be quantitatively characterized in terms of affinity, specificity, and kinetics in a broad variety of solution environments. Several of those provide information on local or global conformational changes of the protein partner in response to ligand binding. This review chapter gives an overview of the state-of-the-art biophysical toolbox for the study of protein-metabolite interactions. It briefly introduces basic principles, highlights recent examples from the literature, and pinpoints promising future directions.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying "magic bullets" amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand-target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural-energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

Curcumin inhibited hepatitis B viral entry through NTCP binding

Hepatitis B virus (HBV) has been implicated in hepatitis and hepatocellular carcinoma. Current agents (nucleos(t)ide analogs and interferons) could only attenuate HBV infection. A combination of agents targeting different stages of viral life cycle (e.g., entry, replication, and cccDNA stability) was expected to eradicate the infection. Curcumin (CCM) was investigated for inhibitory action toward HBV attachment and internalization. Immortalized hepatocyte-like cells (imHCs), HepaRG and non-hepatic cells served as host cells for binding study with CCM. CCM decreased viral load, HBeAg, HBcAg (infectivity), intracellular HBV DNA, and cccDNA levels. The CCM-induced suppression of HBV entry was directly correlated with the density of sodium-taurocholate co-transporting polypeptide (NTCP), a known host receptor for HBV entry. The site of action of CCM was confirmed using TCA uptake assay. The affinity between CCM and NTCP was measured using isothermal titration calorimetry (ITC). These results demonstrated that CCM interrupted HBV entry and would therefore suppress HBV re-infection.

Influence of the disordered domain structure of MeCP2 on its structural stability and dsDNA interaction

Methyl-CpG binding protein 2 (MeCP2) is a transcriptional regulator and a chromatin-associated structural protein. MeCP2 deregulation results in two neurodevelopmental disorders: MeCP2 dysfunction is associated with Rett syndrome, while excess of activity is associated with MeCP2 duplication syndrome. MeCP2 is an intrinsically disordered protein (IDP) constituted by six structural domains with variable, small percentage of well-defined secondary structure. Two domains, methyl-CpG binding domain (MBD) and transcription repressor domain (TRD), are the elements responsible for dsDNA binding ability and recruitment of the gene transcription/silencing machinery, respectively. Previously we studied the influence of the completely disordered, MBD-flanking domains (N-terminal domain, NTD, and intervening domain, ID) on the structural and functional features of the MBD (Claveria-Gimeno, R. et al. Sci Rep. 2017, 7, 41,635). Here we report the biophysical study of the influence of the remaining domains (transcriptional repressor domain, TRD, and C-terminal domains, CTDα and CTDβ) on the structural stability of MBD and the dsDNA binding capabilities of MBD and ID. The influence of distant disordered domains on MBD properties makes it necessary to consider the NTD-MBD-ID variant as the minimal protein construct for studying dsDNA/chromatin binding properties, while the full-length protein should be considered for transcriptional regulation studies.

Binding Thermodynamics to Intrinsically Disordered Protein Domains

Many proteins are intrinsically disordered or contain one or more disordered domains. These domains can participate in binding interactions with other proteins or small ligands. Binding to intrinsically disordered protein domains requires the folding or structuring of those regions such that they can establish well-defined stoichiometric interactions. Since, in such a situation binding is coupled to folding, the energetics of those two events is reflected in the measured binding thermodynamics. In this protocol, we illustrate the thermodynamic differences between binding coupled to folding and binding independent of folding for the same protein. As an example, we use the HIV-1 envelope glycoprotein gp120 that contains structured as well as disordered domains. In the experiments presented, the binding of gp120 to molecules that bind to disordered regions and trigger structuring (CD4 or MAb 17b) and to molecules that bind to structured regions and do not induce conformational structuring (MAb b12) is discussed.

Identification of UQCRB as an oxymatrine recognizing protein using a T7 phage display screen

ETHNOPHARMACOLOGICAL RELEVANCE

Sophora flavescens Aiton (Radix Sophorae Flavescentis, Kushen) is used in traditional Chinese medicine to treat chronic hepatitis B (CHB), and has the ability to clear heat and dampness from the body. Oxymatrine is one of the major bioactive compounds extracted from Sophora flavescens Aiton and constitutes more than 90% of the oxymatrine injection commonly used for CHB treatment in clinics in China.

AIM OF THE STUDY

We aim to analyze the protein binding target of oxymatrine in treating CHB by screening a T7 phage display cDNA library of human CHB and examine the biochemistry of protein-ligand binding between oxymatrine and its ligands.

MATERIALS AND METHODS

A T7 phage cDNA library of human CHB was biopanned by affinity selection using oxymatrine as bait. The interaction of oxymatrine with its candidate binding protein was investigated by affinity assay, molecular docking, Isothermal Titration Calorimetry (ITC) and Surface Plasmon Resonance (SPR).

RESULTS

A library of potential oxymatrine binding peptides was generated. Ubiquinol-cytochrome c reductase binding protein (UQCRB) was one of the candidate binding proteins of oxymatrine. UQCRB-displaying T7 phage binding numbers in the oxymatrine group were significantly higher than that in the control group, biotin group, and matrine group (p<0.05 or p<0.01). Three-dimensional structure modeling of the UQCRB with oxymatrine showed that their binding interfaces matched and oxymatrine inserted into a deeper pocket of UQCRB, which mainly involved amino acid residues Tyr21, Arg33, Tyr83, Glu84, Asp86, Pro88, and Glu91. The binding affinity constant (Kb) from SPR was 4.2mM. The Kb from ITC experiment was 3.9mM and stoichiometry was fixed as 1, which fit very well with the result of SPR. The binding of oxymatrine to UQCRB was driven by strong enthalpy forces such as hydrogen bonds and polar interactions as the heat released was about 157kcal/mol and ΔG was less than zero.

CONCLUSIONS

In this study, using the T7 phage display system, we have identified UQCRB as a direct binding protein of oxymatrine. Furthermore, the specificity and molecular interaction of oxymatrine with UQCRB were also determined. The binding of UQCRB to oxymatrine suggests that UQCRB is a potential target of oxymatrine in treating CHB. These results provide new understanding into the mechanism of oxymatrine and insights into the strategy on the treatment of CHB.

Enthalpy screen of drug candidates

The enthalpic and entropic contributions to the binding affinity of drug candidates have been acknowledged to be important determinants of the quality of a drug molecule. These quantities, usually summarized in the thermodynamic signature, provide a rapid assessment of the forces that drive the binding of a ligand. Having access to the thermodynamic signature in the early stages of the drug discovery process will provide critical information towards the selection of the best drug candidates for development. In this paper, the Enthalpy Screen technique is presented. The enthalpy screen allows fast and accurate determination of the binding enthalpy for hundreds of ligands. As such, it appears to be ideally suited to aid in the ranking of the hundreds of hits that are usually identified after standard high throughput screening.

The Ec-NhaA antiporter switches from antagonistic to synergistic antiport upon a single point mutation

The Na⁺, Li⁺/H⁺ antiporter of Escherichia coli (Ec-NhaA) maintains pH, Na⁺ homeostasis in enterobacteria. We used isothermal titration calorimetry to perform a detailed thermodynamic analysis of Li⁺ binding to Ec-NhaA and several of its mutants. We found that, in line with the canonical alternative access mechanistic model of secondary transporters, Li⁺/H⁺ binding to the antiporter is antagonistically coupled. Binding of Li⁺ displaces 2 H⁺ from the binding site. The process is enthalpically driven, the enthalpic gain just compensating for an entropic loss and the buffer-associated enthalpic changes dominate the overall free-energy change. Li⁺ binding, H⁺ release and antiporter activity were all affected to the same extent by mutations in the Li⁺ binding site (D163E, D163N, D164N, D164E), while D133C changed the H⁺/Li⁺ stoichiometry to 4. Most striking, however, was the mutation, A167P, which converted the Ec-NhaA antagonistic binding into synergistic binding which is only known to occur in Cl⁻/H⁺ antiporter.

Insights into Protein-Ligand Interactions: Mechanisms, Models, and Methods

Molecular recognition, which is the process of biological macromolecules interacting with each other or various small molecules with a high specificity and affinity to form a specific complex, constitutes the basis of all processes in living organisms. Proteins, an important class of biological macromolecules, realize their functions through binding to themselves or other molecules. A detailed understanding of the protein–ligand interactions is therefore central to understanding biology at the molecular level. Moreover, knowledge of the mechanisms responsible for the protein-ligand recognition and binding will also facilitate the discovery, design, and development of drugs. In the present review, first, the physicochemical mechanisms underlying protein–ligand binding, including the binding kinetics, thermodynamic concepts and relationships, and binding driving forces, are introduced and rationalized. Next, three currently existing protein-ligand binding models—the “lock-and-key”, “induced fit”, and “conformational selection”—are described and their underlying thermodynamic mechanisms are discussed. Finally, the methods available for investigating protein–ligand binding affinity, including experimental and theoretical/computational approaches, are introduced, and their advantages, disadvantages, and challenges are discussed.

Synthetic, structural mimetics of the β-hairpin flap of HIV-1 protease inhibit enzyme function

Small-molecule mimetics of the β-hairpin flap of HIV-1 protease (HIV-1 PR) were designed based on a 1,4-benzodiazepine scaffold as a strategy to interfere with the flap-flap protein-protein interaction, which functions as a gated mechanism to control access to the active site. Michaelis-Menten kinetics suggested our small-molecules are competitive inhibitors, which indicates the mode of inhibition is through binding the active site or sterically blocking access to the active site and preventing flap closure, as designed. More generally, a new bioactive scaffold for HIV-1PR inhibition has been discovered, with the most potent compound inhibiting the protease with a modest K(i) of 11 μM.

Characterization of protein-protein interactions by isothermal titration calorimetry

The analysis of protein-protein interactions has attracted the attention of many researchers from both a fundamental point of view and a practical point of view. From a fundamental point of view, the development of an understanding of the signaling events triggered by the interaction of two or more proteins provides key information to elucidate the functioning of many cell processes. From a practical point of view, understanding protein-protein interactions at a quantitative level provides the foundation for the development of antagonists or agonists of those interactions. Isothermal Titration Calorimetry (ITC) is the only technique with the capability of measuring not only binding affinity but the enthalpic and entropic components that define affinity. Over the years, isothermal titration calorimeters have evolved in sensitivity and accuracy. Today, TA Instruments and MicroCal market instruments with the performance required to evaluate protein-protein interactions. In this methods paper, we describe general procedures to analyze heterodimeric (porcine pancreatic trypsin binding to soybean trypsin inhibitor) and homodimeric (bovine pancreatic α-chymotrypsin) protein associations by ITC.

Structural analysis of the regulation of the DYNLL/LC8 binding to Nek9 by phosphorylation

The NIMA family protein kinases Nek9/Nercc1, Nek6, and Nek7 constitute a signaling module activated in early mitosis involved in the control of spindle organization. DYNLL/LC8 (dynein light chain 8) was originally described as a component of the dynein complex, but the recent discovery of multiple interaction partners for LC8 has suggested that it has a general role as a dimerization hub that organizes different protein partners. Recent experiments suggested that LC8 binding to Nek9 was regulated by Nek9 autophosphorylation on Ser(944), a residue immediately located N-terminal to the LC8 conserved (K/R)xTQT binding motif, and that this was crucial for the control of signal transduction through the Nek/Nek6/7 module. In the present work, we present two crystal structures of LC8 with a peptide corresponding to the Nek9 binding region with and without a phosphorylation on Ser(944). Structural analysis of LC8 with both Nek9 peptides, together with different biophysical experiments, explains the observed diminished binding affinity of Nek9 to LC8 upon phosphorylation on Ser(944) within the Nek9 sequence, thus shedding light into a novel phosphorylation regulatory mechanism that interferes with LC8 protein · protein complex formation.

Malarial dihydrofolate reductase as a paradigm for drug development against a resistance-compromised target

Malarial dihydrofolate reductase (DHFR) is the target of antifolate antimalarial drugs such as pyrimethamine and cycloguanil, the clinical efficacy of which have been compromised by resistance arising through mutations at various sites on the enzyme. Here, we describe the use of cocrystal structures with inhibitors and substrates, along with efficacy and pharmacokinetic profiling for the design, characterization, and preclinical development of a selective, highly efficacious, and orally available antimalarial drug candidate that potently inhibits both wild-type and clinically relevant mutated forms of Plasmodium falciparum (Pf) DHFR. Important structural characteristics of P218 include pyrimidine side-chain flexibility and a carboxylate group that makes charge-mediated hydrogen bonds with conserved Arg122 (PfDHFR-TS amino acid numbering). An analogous interaction of P218 with human DHFR is disfavored because of three species-dependent amino acid substitutions in the vicinity of the conserved Arg. Thus, P218 binds to the active site of PfDHFR in a substantially different fashion from the human enzyme, which is the basis for its high selectivity. Unlike pyrimethamine, P218 binds both wild-type and mutant PfDHFR in a slow-on/slow-off tight-binding mode, which prolongs the target residence time. P218, when bound to PfDHFR-TS, resides almost entirely within the envelope mapped out by the dihydrofolate substrate, which may make it less susceptible to resistance mutations. The high in vivo efficacy in a SCID mouse model of P. falciparum malaria, good oral bioavailability, favorable enzyme selectivity, and good safety characteristics of P218 make it a potential candidate for further development.

Extreme entropy-enthalpy compensation in a drug-resistant variant of HIV-1 protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

The application of bioisosteres in drug design for novel drug discovery: focusing on acid protease inhibitors

INTRODUCTION

A bioisostere is a powerful concept for medicinal chemistry. It allows the improvement of the stability; oral absorption; membrane permeability; and absorption, distribution, metabolism and excretion (ADME) of drug candidate, while retaining their biological properties. The term 'bioisostere' is derived from 'isostere', whose physical and chemical properties, such as steric size, hydrophobicity, and electronegativity, are similar to those of a functional or atomic group, and is considered to possess biological properties. Here, the authors highlight the recent applications of bioisosteres in drug design, mainly based on our drug discovery studies.

AREAS COVERED

This review discusses the application of bioisosteres for novel drug discovery with focus on the authors' drug discovery studies such as renin, HIV-protease, and β-secretase inhibitors. The authors highlight that some bioisosteres can form the scaffolding for drug candidates, namely substrate transition state, amide/ester, and carboxylic acid bioisosteres. Moreover, the authors propose the new terms 'electron-donor bioisostere' and 'conformational bioisostere' for drug discovery.

EXPERT OPINION

The authors discuss the importance of bioisostere's design concept based on specific interaction with the corresponding biomolecule. In addition, some strategies for drug discovery based on the bioisostere concept are introduced. Many bioisosteres, which are recognized by corresponding target biomolecules as exhibiting similar biological properties, have been reported to date; most of the recently developed bioisosteres were designed by cheminformatics approaches. Some molecular design softwares and databases are introduced.

Placevent: an algorithm for prediction of explicit solvent atom distribution-application to HIV-1 protease and F-ATP synthase

We have created a simple algorithm for automatically predicting the explicit solvent atom distribution of biomolecules. The explicit distribution is coerced from the three-dimensional (3D) continuous distribution resulting from a 3D reference interaction site model (3D-RISM) calculation. This procedure predicts optimal location of solvent molecules and ions given a rigid biomolecular structure and the solvent composition. We show examples of predicting water molecules near the KNI-272 bound form of HIV-1 protease and predicting both sodium ions and water molecules near the rotor ring of F-adenosine triphosphate (ATP) synthase. Our results give excellent agreement with experimental structure with an average prediction error of 0.39-0.65 Å. Further, unlike experimental methods, this method does not suffer from the partial occupancy limit. Our method can be performed directly on 3D-RISM output within minutes. It is extremely useful for examining multiple specific solvent-solute interactions, as a convenient method for generating initial solvent structures for molecular dynamics calculations, and may assist in refinement of experimental structures. © 2012 Wiley Periodicals, Inc.

Synergistic substrate binding determines the stoichiometry of transport of a prokaryotic H(+)/Cl(-) exchanger

Active exchangers dissipate the gradient of one substrate to accumulate nutrients, export xenobiotics and maintain cellular homeostasis. Mechanistic studies suggested that all exchangers share two fundamental properties: substrate binding is antagonistic and coupling is maintained by preventing shuttling of the empty transporter. The CLC Cl⁻: H⁺ exchangers control the homeostasis of cellular compartments in most living organisms but their transport mechanism remains unclear. We show that substrate binding to CLC-ec1 is synergistic rather than antagonistic: chloride binding induces protonation of a critical glutamate. The simultaneous binding of H⁺ and Cl⁻ gives rise to a fully-loaded state incompatible with conventional mechanisms. Mutations in the Cl⁻ transport pathway identically alter the stoichiometries of Cl⁻: H⁺ exchange and binding. We propose that the thermodynamics of synergistic substrate binding determine the stoichiometry of transport rather than the kinetics of conformational changes and ion binding.

Finding a better path to drug selectivity

Extremely high affinity and selectivity are two of the most sought-after properties of drug molecules. Selectivity has been difficult to achieve, especially for targets that belong to large families of structurally and functionally related proteins. There are essentially two ways by which selectivity can be improved during lead optimization: a chemical modification of the lead compound that improves the affinity towards the target to a higher extent than to off-target molecules; and a chemical modification that lowers the affinity of the lead compound towards off-target molecules. Maximal selectivity is achieved when both mechanisms can be combined synergistically. As we discuss here, analysis of several protease inhibitors that vary in a single functionality indicates that nonpolar functionalities preferentially follow the first mechanism, whereas polar functionalities follow the second, and that those features are imprinted in their thermodynamic signatures.

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.

Accurate ensemble molecular dynamics binding free energy ranking of multidrug-resistant HIV-1 proteases

Accurate calculation of important thermodynamic properties, such as macromolecular binding free energies, is one of the principal goals of molecular dynamics simulations. However, single long simulation frequently produces incorrectly converged quantitative results due to inadequate sampling of conformational space in a feasible wall-clock time. Multiple short (ensemble) simulations have been shown to explore conformational space more effectively than single long simulations, but the two methods have not yet been thermodynamically compared. Here we show that, for end-state binding free energy determination methods, ensemble simulations exhibit significantly enhanced thermodynamic sampling over single long simulations and result in accurate and converged relative binding free energies that are reproducible to within 0.5 kcal/mol. Completely correct ranking is obtained for six HIV-1 protease variants bound to lopinavir with a correlation coefficient of 0.89 and a mean relative deviation from experiment of 0.9 kcal/mol. Multidrug resistance to lopinavir is enthalpically driven and increases through a decrease in the protein-ligand van der Waals interaction, principally due to the V82A/I84V mutation, and an increase in net electrostatic repulsion due to water-mediated disruption of protein-ligand interactions in the catalytic region. Furthermore, we correctly rank, to within 1 kcal/mol of experiment, the substantially increased chemical potency of lopinavir binding to the wild-type protease compared to saquinavir and show that lopinavir takes advantage of a decreased net electrostatic repulsion to confer enhanced binding. Our approach is dependent on the combined use of petascale computing resources and on an automated simulation workflow to attain the required level of sampling and turn around time to obtain the results, which can be as little as three days. This level of performance promotes integration of such methodology with clinical decision support systems for the optimization of patient-specific therapy.

Computational study of the resistance shown by the subtype B/HIV-1 protease to currently known inhibitors

Human immunodeficiency virus type 1 protease (HIV-1 PR) is an essential enzyme in the HIV-1 life cycle. As such, this protein represents a major drug target in AIDS therapy, but emerging resistance to antiretroviral inhibitor cocktails, caused by high viral mutation rates, represents a significant challenge in AIDS treatment. Many mutations are not located within the active site or binding pocket, nor they do significantly modify the three-dimensional structural organization of the enzyme; hence, the mechanism(s) by which they alter inhibitor affinity for the protease remains uncertain. In this article, we present an all-atom computational analysis of the dynamic residue-residue coordination between the active site residues and the rest of the protein and of the energetic properties of different HIV-1 PR complexes. We analyze both the wild-type form and mutated forms that induce drug resistance. In particular, the results show differences between the wild type and the mutants in their mechanism of dynamic coordination, in the signal propagation between the active site residues and the rest of the protein, and in the energy networks responsible for the stabilization of the bound inhibitor conformation. Finally, we propose a dynamic and energetic explanation for HIV-1 protease drug resistance, and, through this model, we identify a possible new site that could be helpful in the design of a new family of HIV-1 PR allosteric inhibitors.

How much binding affinity can be gained by filling a cavity?

Binding affinity optimization is critical during drug development. Here, we evaluate the thermodynamic consequences of filling a binding cavity with functionalities of increasing van der Waals radii (-H, -F, -Cl, and CH(3)) that improve the geometric fit without participating in hydrogen bonding or other specific interactions. We observe a binding affinity increase of two orders of magnitude. There appears to be three phases in the process. The first phase is associated with the formation of stable van der Waals interactions. This phase is characterized by a gain in binding enthalpy and a loss in binding entropy, attributed to a loss of conformational degrees of freedom. For the specific case presented in this article, the enthalpy gain amounts to -1.5 kcal/mol while the entropic losses amount to +0.9 kcal/mol resulting in a net 3.5-fold affinity gain. The second phase is characterized by simultaneous enthalpic and entropic gains. This phase improves the binding affinity 25-fold. The third phase represents the collapse of the trend and is triggered by the introduction of chemical functionalities larger than the binding cavity itself [CH(CH(3))(2)]. It is characterized by large enthalpy and affinity losses. The thermodynamic signatures associated with each phase provide guidelines for lead optimization.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TRbeta) vs. TRalpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TRbeta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331beta) in the TRbeta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3'-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TRbeta selectivity. TR x-ray structures reveal better fit of ligand with the TRalpha LBC. The TRbeta LBC, however, expands relative to TRalpha in the presence of Triac (549 A(3) vs. 461 A(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TRbeta and permits greater flexibility of the Triac carboxylate group in TRbeta than in TRalpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TRbeta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Automated site preparation in physics-based rescoring of receptor ligand complexes

Hydrogen atoms are not typically observable in X-ray crystal structures, but inferring their locations is often important in structure-based drug design. In addition, protonation states of the protein can change in response to ligand binding, as can the orientations of OH groups, a subtle form of "induced fit." We implement and evaluate an automated procedure for optimizing polar hydrogens in protein-binding sites in complex with ligands. Specifically, we apply the previously described Independent Cluster Decomposition Algorithm (ICDA) algorithm (Li et al., Proteins 2007;66:824-837), which assigns the ionization states of titratable residues, the amide orientations of Asn/Gln side chains, the imidazole ring orientation in His, and the orientations of OH/SH groups, in a unified algorithm. We test the utility of this method for identifying nativelike ligand poses using 247 protein-ligand complexes from an established database of docked decoys. Pose selection is performed with a physics-based scoring function based on a molecular mechanics energy function and a Generalized Born implicit solvent model. The use of the ICDA receptor preparation protocol, implemented with no knowledge of the native ligand pose, increases the accuracy of pose selection significantly, with the average RMSD over all complexes decreasing from 2.7 to 1.5 A when applying ICDA. Large improvements are seen for specific classes of binding sites with titratable groups, such as aspartyl proteases.

Structure of HIV-1 protease in complex with potent inhibitor KNI-272 determined by high-resolution X-ray and neutron crystallography

HIV-1 protease is a dimeric aspartic protease that plays an essential role in viral replication. To further understand the catalytic mechanism and inhibitor recognition of HIV-1 protease, we need to determine the locations of key hydrogen atoms in the catalytic aspartates Asp-25 and Asp-125. The structure of HIV-1 protease in complex with transition-state analog KNI-272 was determined by combined neutron crystallography at 1.9-A resolution and X-ray crystallography at 1.4-A resolution. The resulting structural data show that the catalytic residue Asp-25 is protonated and that Asp-125 (the catalytic residue from the corresponding diad-related molecule) is deprotonated. The proton on Asp-25 makes a hydrogen bond with the carbonyl group of the allophenylnorstatine (Apns) group in KNI-272. The deprotonated Asp-125 bonds to the hydroxyl proton of Apns. The results provide direct experimental evidence for proposed aspects of the catalytic mechanism of HIV-1 protease and can therefore contribute substantially to the development of specific inhibitors for therapeutic application.

Do enthalpy and entropy distinguish first in class from best in class?

A drug molecule should bind to its target with high affinity and selectivity. Because the binding affinity is a combined function of the binding enthalpy and the binding entropy, extremely high affinity requires that both terms contribute favorably to binding. The binding enthalpy, however, is notoriously more difficult to optimize than the binding entropy, a fact that has resulted in thermodynamically unbalanced molecules that do not achieve optimal potency. In fact, with current technologies, the enthalpic optimization of drug candidates may take years and only appear in second-generation products. Within that context, it is not surprising that structure/activity relationships (SARs) that explicitly incorporate the interplay between enthalpy and entropy and accelerate the optimization process are being developed and gaining popularity.

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.

Peptide inhibition of HIV-1: current status and future potential

More than 2 decades of intensive research has focused on defining replication mechanisms of HIV type 1 (HIV-1), the etiologic agent of AIDS. The delineation of strategies for combating this viral infection has yielded many innovative approaches toward this end. HIV-1 is a lentivirus in the family retroviridae that is relatively small with regard to both structure and genome size, having a diploid RNA genome of approximately 9 kb, with only three major genes and several gene products resulting from alternate splicing and translational frameshifting. Most marketed drugs for treating AIDS are inhibitors of HIV-1 reverse transcriptase or protease enzymes, but new targets include the integrase enzyme, cell surface interactions that facilitate viral entry, and also virus particle maturation and assembly. The emergence of drug-resistant variants of HIV-1 has been the main impediment to successful treatment of AIDS. Thus, there is a pressing need to develop novel treatment strategies targeting multiple stages of the virus life-cycle. Research efforts aimed at developing successful means for combating HIV-1 infection have included development of peptide inhibitors of HIV-1. This article summarizes past and current endeavors in the development of peptides that inhibit replication of HIV-1 and the role of peptide inhibitors in the search for new anti-HIV drugs.

Resilience to resistance of HIV-1 protease inhibitors: profile of darunavir

The current effectiveness of HAART in the management of HIV infection is compromised by the emergence of extensively cross-resistant strains of HIV-1, requiring a significant need for new therapeutic agents. Due to its crucial role in viral maturation and therefore HIV-1 replication and infectivity, the HIV-1 protease continues to be a major development target for antiretroviral therapy. However, new protease inhibitors must have higher thresholds to the development of resistance and cross-resistance. Research has demonstrated that the binding characteristics between a protease inhibitor and the active site of the HIV-1 protease are key factors in the development of resistance. More specifically, the way in which a protease inhibitor fits within the substrate consensus volume, or "substrate envelope", appears to be critical. The currently available inhibitors are not only smaller than the native substrates, but also have a different shape. This difference in shape underlies observed patterns of resistance because primary drug-resistant mutations often arise at positions in the protease where the inhibitors protrude beyond the substrate envelope but are still in contact with the enzyme. Since all currently available protease inhibitors occupy a similar space (in spite of their structural differences) in the active site of the enzyme, the specific positions where the inhibitors protrude and contact the enzyme correspond to the locations where most mutations occur that give rise to multidrug-resistant HIV-1 strains. Detailed investigation of the structure, thermodynamics, and dynamics of the active site of the protease enzyme is enabling the identification of new protease inhibitors that more closely fit within the substrate envelope and therefore decrease the risk of drug resistance developing. The features of darunavir, the latest FDA-approved protease inhibitor, include its high binding affinity (Kd = 4.5 x 10-12 M) for the protease active site, the presence of hydrogen bonds with the backbone, and its ability to fit closely within the substrate envelope (or consensus volume). Darunavir is potent against both wild-type and protease inhibitor-resistant viruses in vitro, including a broad range of over 4,000 clinical isolates. Additionally, in vitro selection studies with wild-type HIV-1 strains have shown that resistance to darunavir develops much more slowly and is more difficult to generate than for existing protease inhibitors. Clinical studies have shown that darunavir administered with low-dose ritonavir (darunavir/ritonavir) provides highly potent viral suppression (including significant decreases in HIV viral load in patients with documented protease inhibitor resistance) together with favorable tolerability. In conclusion, as a result of its high binding affinity for and overall fit within the active site of HIV-1 protease, darunavir has a higher genetic barrier to the development of resistance and better clinical efficacy against multidrug-resistant HIV relative to current protease inhibitors. The observed efficacy, safety and tolerability of darunavir in highly treatment-experienced patients makes darunavir an important new therapeutic option for HIV-infected patients.

Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.

Compensating Enthalpic and Entropic Changes Hinder Binding Affinity Optimization

A common strategy to improve the potency of drug candidates is to introduce chemical functionalities, like hydrogen bond donors or acceptors, at positions where they are able to establish strong interactions with the target. However, it is often observed that the added functionalities do not necessarily improve potency even if they form strong hydrogen bonds. Here, we explore the thermodynamic and structural basis for those observations. KNI-10033 is a potent experimental HIV-1 protease inhibitor with picomolar affinity against the wild-type enzyme (K(d) = 13 pm). The potency of the inhibitor is the result of favorable enthalpic (DeltaH = -8.2 kcal/mol) and entropic (-TDeltaS = -6.7 kcal/mol) interactions. The replacement of the thioether group in KNI-10033 by a sulfonyl group (KNI-10075) results in a strong hydrogen bond with the amide of Asp 30B of the HIV-1 protease. This additional hydrogen bond improves the binding enthalpy by 3.9 kcal/mol; however, the enthalpy gain is completely compensated by an entropy loss, resulting in no affinity change. Crystallographic and thermodynamic analysis of the inhibitor/protease complexes indicates that the entropy losses are due to a combination of conformational and solvation effects. These results provide a set of practical guidelines aimed at overcoming enthalpy/entropy compensation and improve binding potency.