Multipolar interactions in the D pocket of thrombin: large differences between tricyclic imide and lactam inhibitors

Two series of tricyclic inhibitors of the serine protease thrombin, imides (+/-)-1-(+/-)-8 and lactams (+/-)-9-(+/-)-13, were analysed to evaluate contributions of orthogonal multipolar interactions with the backbone C=O moiety of Asn98 to the free enthalpy of protein-ligand complexation. The lactam derivatives are much more potent and more selective inhibitors (K(i) values between 0.065 and 0.005 microM, selectivity for thrombin over trypsin between 361- and 1609-fold) than the imide compounds (Ki values between 0.057 and 23.7 microM, selectivity for thrombin over trypsin between 3- and 67-fold). The increase in potency and selectivity is explained by the favorable occupancy of the P-pocket of thrombin by the additional isopropyl substituent in the lactam derivatives. The nature of the substituent on the benzyl ring filling the D pocket strongly influences binding potency in the imide series, with Ki values increasing in the sequence: F < OCH2O < Cl < H < OMe < OH < N(pyr)< Br. This sequence can be explained by both steric fit and the occurrence of orthogonal multipolar interactions with the backbone C[double bond, length as m-dash]O moiety of Asn98. In contrast, the substituent on the benzyl ring hardly affects the ligand potency in the lactam series. This discrepancy was clarified by the comparison of X-ray structures solved for co-crystals of thrombin with imide and lactam ligands. Whereas the benzyl substituents in the imide inhibitors are sufficiently close (< or =3.5 Angstroms) to the C=O group of Asn98 to allow for attractive orthogonal multipolar interactions, the distances in the lactam series are too large (> or =4 Angstroms) for attractive dipolar contacts to be effective.     

The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure

Summary A new simple empirical function has been developed that estimates the free energy of binding for a given protein-ligand complex of known 3D structure. The function takes into account hydrogen bonds, ionic interactions, the lipophilic protein-ligand contact surface and the number of rotatable bonds in the ligand. The dataset for the calibration of the function consists of 45 protein-ligand complexes. The new energy function reproduces the binding constants (ranging from 2.5·10^-2 to 4·10^-14 M, corresponding to binding energies between -9 and -76 kJ/mol) of the dataset with a standard deviation of 7.9 kJ/mol, corresponding to 1.4 orders of magnitude in binding affinity. The individual contributions to protein-ligand binding obtained from the scoring function are: ideal neutral hydrogen bond: -4.7 kJ/mol; ideal ionic interaction: -8.3 kJ/mol; lipophilic contact: -0.17 kJ/mol Ų; one rotatable bond in the ligand: +1.4 kJ/mol. The function also contains a constant contribution (+5.4 kJ/mol) which may be rationalized as loss of translational and rotational entropy. The function can be evaluated very fast and is therefore also suitable for application in a 3D database search or de novo ligand design program such as LUDI.     

Rapid protein-ligand costructures using chemical shift perturbations

Structure-based drug discovery requires the iterative determination of protein-ligand costructures in order to improve the binding affinity and selectivity of potential drug candidates. In general, X-ray and NMR structure determination methods are time consuming and are typically the limiting factor in the drug discovery process. The application of molecular docking simulations to filter and evaluate drug candidates has become a common method to improve the throughput and efficiency of structure-based drug design. Unfortunately, molecular docking methods suffer from common problems that include ambiguous ligand conformers or failure to predict the correct docked structure. A rapid approach to determine accurate protein-ligand costructures is described based on NMR chemical shift perturbation (CSP) data routinely obtained using 2D 1H-15N HSQC spectra in high-throughput ligand affinity screens. The CSP data is used to both guide and filter AutoDock calculations using our AutoDockFilter program. This method is demonstrated for 19 distinct protein-ligand complexes where the docked conformers exhibited an average rmsd of 1.17 +/- 0.74 A relative to the original X-ray structures for the protein-ligand complexes.     

Molecular design of potent inhibitor specific for cathepsin B based on the tertiary structure prediction.

To design a potent inhibitor specific for cathepsin B (rat liver), the tertiary structure was predicted based on the crystal structure of the papain complexed with (+)-(2S,3S)-3-(1-[N-(3-methylbutyl)amino]leucylcarbonyl)oxirane-2- carbolylic acid (E-64-c), a thiol protease inhibitor. Taking advantage of the structural characteristics of the predicted active site, seventeen inhibitors were chemically synthesized by molecular modeling, and one of them, N-(L-3-trans-propylcarbamoyloxirane-2-carbonyl)-L-isoleucyl-L-p rol ine (CA-074) was shown to be the first potent inhibitor specific for cathepsin B. The relationship between the structure and inhibitory activity is discussed based on the model structure of the cathepsin B-inhibitor complex.     

Alkali metal ion catalysis and inhibition in nucleophilic displacement reactions at phosphorus centers: ethyl and methyl paraoxon and ethyl and methyl parathion

We report on the ethanolysis of the P=O and P=S compounds ethyl and methyl paraoxon (1a and 1b) and ethyl and methyl parathion (2a and 2b). Plots of spectrophotometrically measured rate constants, kobsd versus [MOEt], the alkali ethoxide concentration, show distinct upward and downward curvatures, pointing to the importance of ion-pairing phenomena and a differential reactivity of free ions and ion pairs. Three types of reactivity and selectivity patterns have been discerned: (1) For the P=O compounds 1a and 1b, LiOEt > NaOEt > KOEt > EtO-; (2) for the P=S compound 2a, KOEt > EtO- > NaOEt > LiOEt; (3) for P=S, 2b, 18C6-crown-complexed KOEt > KOEt = EtO(-) > NaOEt > LiOEt. These selectivity patterns are characteristic of both catalysis and inhibition by alkali-metal cations depending on the nature of the electrophilic center, P=O vs P=S, and the metal cation. Ground-state (GS) vs transition-state (TS) stabilization energies shed light on the catalytic and inhibitory tendencies. The unprecedented catalytic behavior of crowned-K(+) for the reaction of 2b is noteworthy. Modeling reveals an extreme steric interaction for the reaction of 2a with crowned-K(+), which is responsible for the absence of catalysis in this system. Overall, P=O exhibits greater reactivity than P=S, increasing from 50- to 60-fold with free EtO(-) and up to 2000-fold with LiOEt, reflecting an intrinsic P=O vs P=S reactivity difference (thio effect). The origin of reactivity and selectivity differences in these systems is discussed on the basis of competing electrostatic effects and solvational requirements as function of anionic electric field strength and cation size (Eisenman's theory).     

Shape selectivity versus functional group pre-organization in molecularly imprinted polymers

A systematic survey of related molecular probes differing in shape or functional group orientation was used to compare the effects of shape selectivity versus pre-organization of functional monomers on imprinting and rebinding performance of molecularly imprinted polymers (MIPs). These studies revealed that templates with two functional group interactions with the MIPs are influenced to a larger degree by shape selective interactions than templates with three functional group interactions. For example, with two functional group interactions, increasing side chain size of compounds 1-5 increased selectivity 5-fold; while the same size change for compounds with three functional group interactions leads to a maximum 2-fold increase. Thus, the effects of shape selectivity and pre-organization of functional groups do not appear to work in concert with each other during the imprinting process or in the rebinding behavior. Furthermore, greater selectivity is generally found for templates with two functional groups, where the dominant mode of molecular recognition is shape selectivity. For example, the α value for the MIP elicited toward template compound 5 with two hydrogen bonding groups was 5-12-fold higher than equivalently shaped compounds 6-8 that have three non-covalent binding interactions (Table 3). On the other hand, pre-organization of functional groups dominated the performance of MIPs elicited toward templates with three template-functional group interactions. This is observed in Tables 6-8, where compounds with identical positioning of three functional groups all show less than an order of magnitude change in α values despite changes in shape.     

Copper(II), zinc(II), and cadmium(II) coordination polymers with 2-methylglutarato and 4,4′-dipyridyldisulfide

Three new coordination polymers of copper(II), zinc(II) and cadmium(II), Cu(H₂O)(Dpds)(2-MGA) (I), [Zn(Dpds)(2-MGA)] · 1.25H₂O (II) and [Cd(H₂O)(Dpds)(2-MGA)] · 0.25H₂O (III) (Dpds = 4,4′-dipyridyldisulfide, H₂MGA = (RS)-2-methyl glutaric acid), have been synthesized and characteried by X-ray single crystal structure determination. The Cu atoms in I are alternately bridged by Dpds ligands and 2-methylglutarato ligands to generate 1D chain. The resulted chains are assembled via S...S weak interactions into 2D layers, which are through twofold 2D parallel/2D parallel mode inclined interpenetration to induce 3D supramolecular architecture. In II, the ZnN₂O₂ tetrahedras are bridged by 2-MGA anion and Dpds ligands to form 2D (4,4) networks, which are assembled via hydrogen bonds to 3D supramolecular architecture. The centrosymmetric binuclear units Cd₂(2-MGA)₂ in III are bridged by Dpds ligands to form 1D repeated rhomboids chains, which are interlinked via S...S weak interactions into 2D layer, and the resulting 2D sheets are inclined parallel into 3D network.     

Target Self-Enhanced Selectivity in Metal-Specific DNAzymes

Highly selective recognition of metal ions by rational ligand design is challenging, and simple metal binding by biological ligands is often obscured by nonspecific interactions. In this work, binding-triggered catalysis is used and metal selectivity is greatly increased by increasing the number of metal ions involved, as exemplified in a series of in vitro selected RNA-cleaving DNAzymes. The cleavage junction is modified with a glycyl–histidine-functionalized tertiary amine moiety to provide multiple potential metal coordination sites. DNAzymes that bind 1, 2, and 3 Zn²⁺ ions, increased their selectivity for Zn²⁺ over Co²⁺ ions from approximately 20-, 1000-, to 5000-fold, respectively. This study offers important insights into metal recognition by combining rational ligand design and combinatorial selection, and it provides a set of new DNAzymes with excellent selectivity for Zn²⁺ ions.     

Structure-Based Evolution of Subtype-Selective Neurotensin Receptor Ligands

Subtype-selective agonists of the neurotensin receptor NTS2 represent a promising option for the treatment of neuropathic pain, as NTS2 is involved in the mediation of μ-opioid-independent anti-nociceptive effects. Based on the crystal structure of the subtype NTS1 and previous structure–activity relationships (SARs) indicating a potential role for the sub-pocket around Tyr11 of NT(8–13) in subtype-specific ligand recognition, we have developed new NTS2-selective ligands. Starting from NT(8–13), we replaced the tyrosine unit by β²-amino acids (type 1), by heterocyclic tyrosine bioisosteres (type 2) and peptoid analogues (type 3). We were able to evolve an asymmetric synthesis of a 5-substituted azaindolylalanine and its application as a bioisostere of tyrosine capable of enhancing NTS2 selectivity. The S-configured test compound 2 a, [(S)-3-(pyrazolo[1,5-a]pyridine-5-yl)-propionyl¹¹]NT(8–13), exhibits substantial NTS2 affinity (4.8 nm) and has a nearly 30-fold NTS2 selectivity over NTS1. The (R)-epimer 2 b showed lower NTS2 affinity but more than 600-fold selectivity over NTS1.     

Heterometallic cubane single-molecule magnets

Two new heterometallic cubane molecules have been synthesized. High-frequency electron paramagnetic resonance and magnetization measurements indicate that [Mn(3)Ni(hmp)(3)O(N(3))(3)(C(7)H(5)O(2))(3)] (1) displays a well-isolated S = 5 ground state (DeltaE > 120 K), with g = 2.0, D = -0.23 cm(-1), and ferromagnetic Mn-Mn exchange interactions competing with antiferromagnetic Ni-Mn interactions. [Mn(3)Zn(hmp)(3)O(N(3))(3)(C(3)H(5)O(2))(3)] (2) possesses a S = 6 ground state (DeltaE > 105 K), with g = 2.0, D = -0.14 cm(-1), and ferromagnetic Mn-Mn exchange interactions. Magnetization vs magnetic field data for oriented single crystals of 1 and 2 indicate that both complexes are single-molecule magnets.     

High-affinity mu opioid receptor ligands discovered by the screening of an exhaustively stereodiversified library of 1,5-enediols

In this communication, we report the synthesis of an exhaustively stereodiversified library of 16 1,5-enediols (2) and the screening of these compounds for mu opioid receptor (MOR) binding. The stereochemical configuration of 2 strongly impacted the binding affinity, and (S,S,S,R)-2 exhibited a Ki of 8.8 nM for MOR, comparable to that of endomorphin-2 (Ki = 1.2 nM). Moreover, compounds 2 exhibited 5-86-fold selectivity for MOR over delta opioid receptor (DOR) and 16-150-fold selectivity for MOR over kappa opioid receptor (KOR). Additionally, analogues of 2 were synthesized which showed the importance of the trans olefin for receptor binding but that modifications of the C-terminal amino acid were well tolerated. Ligand 11 is noteworthy because it retains only one of the amide bonds present in 1, but binds MOR with an affinity of 10 nM and 110- and 600-fold selectivity for MOR over DOR and KOR. These results demonstrate the utility of stereochemical diversity in the discovery of bioactive small molecules.     

Hepatitis C virus NS5B polymerase: QM/MM calculations show the important role of the internal energy in ligand binding

The inter- and intramolecular interactions that determine the experimentally observed binding mode of the ligand (2Z)-2-(benzoylamino)-3-[4-(2-bromophenoxy)phenyl]-2-propenoate in complex with hepatitis C virus NS5B polymerase have been studied using QM/MM calculations. DFT-based QM/MM optimizations were performed on a number of ligand conformers in the protein-ligand complex. Using these initial poses, our aim is 2-fold. First, we identify the minimum energy pose. Second, we dissect the energetic contributions to this pose using QM/MM methods. The study reveals the critical importance of internal energy for the proper energy ranking of the docked poses. Using this protocol, we successfully identified three poses that have low RMSD with respect to the crystallographic structure from among the top 20 initially docked poses. We show that the most important energetic component contributing to binding for this particular protein-ligand system is the conformational (i.e., QM internal) energy.     

Nonpeptide inhibitors of cathepsin G: optimization of a novel beta-ketophosphonic acid lead by structure-based drug design

The serine protease cathepsin G (EC 3.4.21.20; Cat G), which is stored in the azurophilic granules of neutrophils (polymorphonuclear leukocytes) and released on degranulation, has been implicated in various pathological conditions associated with inflammation. By employing high-throughput screening, we identified beta-ketophosphonic acid 1 as a moderate inhibitor of Cat G (IC(50) = 4.1 microM). We were fortunate to obtain a cocrystal of 1 with Cat G and solve its structure by X-ray crystallography (3.5 A). Structural details from the X-ray analysis of 1.Cat G served as a platform for optimization of this lead compound by structure-based drug design. With the aid of molecular modeling, substituents were attached to the 3-position of the 2-naphthyl ring of 1, which occupies the S1 pocket of Cat G, to provide an extension into the hydrophobic S3 region. Thus, we arrived at analogue 7 with an 80-fold potency improvement over 1 (IC(50) = 53 nM). From these results, it is evident that the beta-ketophosphonic acid unit can form the basis for a novel class of serine protease inhibitors.     

Development of a light-up fluorescent probe for HRAS G-quadruplex DNA

The development of small-molecule G-quadruplex DNA probes has attracted significant attention in recent years. However, G-quadruplexes can display a wide variety of topologies, which process different structures and functions. Therefore, selective discrimination one G-quadruplex structure over another is promising. Herein, we reported the design, synthesis and biological evaluation of a long-chain fatty amine functionalized triphenylamine-quinolinium conjugate 1b. Significant enhancement of the fluorescence intensity (over 180 fold) was observed when 1b bound with HRAS G-quadruplex DNA, while much weaker enhancements were presented in the presence of other G-quadruplexes (45–90-fold) and single/double-stranded DNAs (less than 20-fold), indicating 1b had an excellent selectivity to HRAS. The details of the interactions were investigated by UV–Vis, FID and CD analysis. The results show 1b could interact and stabilize HRAS structure mainly by π-π stacking binding mode. The introduced amine chain of the structure core was found to be better in the terms of inducing selectivity toward G-quadruplex structure. In addition, the application of 1b as a fluorescent agent for living cell imaging was also demonstrated.     

Strategy in inhibition of cathepsin B, a target in tumor invasion and metastasis

Cathepsin B, a cysteine protease, is an important target in fighting cancer. This enzyme has been implicated in enhancing tumor invasiveness and metastasis, therefore inhibitors for cathepsin B are highly sought as potential anticancer and antimetastatic agents. A structure-based design effort was pursued in arriving at a template for inhibition of cathepsin B. Focused compound libraries were synthesized based on this template, which were screened for cathepsin B inhibitory properties. Compound 2, 1-(2(R)-[1(S)-acetoxy-2-[2(S)-(2,4-difluoro-benzoylamino)-3-phenyl-propionylaminooxy]-2-oxo-ethyl]-pentanoyl)-pyrrolidine-2(S)-carboxylic acid benzyl ester, is the prototype of this novel class of cysteine protease inhibitor that emerged from the search. The molecule modifies the active site of cathepsin B covalently, irreversibly, and efficiently, a process for which the kinetic parameters were evaluated. A set of three judiciously altered variants of compound 2 was also synthesized to explore the details of the proposed mechanism of action by this inhibitor. Compound 2 and its analogues may prove useful tools in reversing the deleterious effect of cathepsin B in fighting cancer.     

Two Cadmium(II) Coordination Polymers based on Pamoic Acid and Different Polydentate N-donor Ligands: Syntheses,Crystal Structures,and Fluorescence Properties

Two new fluorescent coordination polymers based on pamoic acid and different polydentate N-donor ligands, namely {[Cd(PA)(TPTZ)(H₂O)](DMF)₂}_n ( 1 ) and [Cd(PA)(BIB)]_n ( 2 ) [H₂PA = pamoic acid, TPTZ = 2,4,6-tri(2-pyridyl)-1,3,5-triazine, BIB = 1,4-bis(1-imidazolyl)benzene], were synthesized and characterized. Complex 1 showed a 1D zigzag chain structure with intramolecular hydrogen bonds. The 2D supramolecular structure in 1 was formed through π–π stacking interactions and intermolecular hydrogen bonds. Complex 2 displayed a 2D network structure. Intramolecular hydrogen bonds and π–π stacking interactions were observed in 2 . By studying the fluorescence sensing performance of two coordination polymers, complex 1 exhibited high selectivity for tracking Al³⁺ ion and complex 2 could discriminately detect inorganic or aliphatic amines with high selectivity.     

Rational design of thioamide peptides as selective inhibitors of cysteine protease cathepsin L

Aberrant levels of cathepsin L (Cts L), a ubiquitously expressed endosomal cysteine protease, have been implicated in many diseases such as cancer and diabetes. Significantly, Cts L has been identified as a potential target for the treatment of COVID-19 due to its recently unveiled critical role in SARS-CoV-2 entry into the host cells. However, there are currently no clinically approved specific inhibitors of Cts L, as it is often challenging to obtain specificity against the many highly homologous cathepsin family cysteine proteases. Peptide-based agents are often promising protease inhibitors as they offer high selectivity and potency, but unfortunately are subject to degradation in vivo. Thioamide substitution, a single-atom O-to-S modification in the peptide backbone, has been shown to improve the proteolytic stability of peptides addressing this issue. Utilizing this approach, we demonstrate herein that good peptidyl substrates can be converted into sub-micromolar inhibitors of Cts L by a single thioamide substitution in the peptide backbone. We have designed and scanned several thioamide stabilized peptide scaffolds, in which one peptide, R^S_1A, was stabilized against proteolysis by all five cathepsins (Cts L, Cts V, Cts K, Cts S, and Cts B) while inhibiting Cts L with >25-fold specificity against the other cathepsins. We further showed that this stabilized R^S_1A peptide could inhibit Cts L in human liver carcinoma lysates (IC₅₀ = 19 μM). Our study demonstrates that one can rationally design a stabilized, specific peptidyl protease inhibitor by strategic placement of a thioamide and reaffirms the place of this single-atom modification in the toolbox of peptide-based rational drug design.

Information on the effects of sidechain and backbone modification on the activity of cathepsin (Cts) L, V, K, S, and B was used to design a thioamide peptide that is inert to all Cts and selectively inhibits Cts L.     

Intermolecular (1)H-(1)H Two-Dimensional Nuclear Overhauser Enhancements in the Characterization of a Rationally Designed Chiral Recognition System

Chiral stationary phases (CSPs) for liquid chromatography derived from N-(acyl)proline-3,5-dimethylanilides separate the enantiomers of N-(3,5-dinitrobenzoyl)-alpha-amino esters and amides with high levels of selectivity. These CSPs have been used to assemble a large body of chromatographic data which indirectly supports the validity of the mechanistic rationale originally used in the design of these CSPs. We herein report (1)H and (13)C chemical shift data obtained when the (S)-enantiomer of chiral solvating agent (CSA) 3, a soluble analogue of the selector used in CSP (S)-1, acts on each of the enantiomers of the dimethylamide of N-(3,5-dinitrobenzoyl)leucine, 2. The changes in chemical shift in the mixture of (S)-2 and (S)-3 support the existence of those interactions thought to be essential to chiral recognition in this system. In addition, significant intermolecular NOESY enhancements are observed in this mixture. These NOE data are consistent with the structure expected for the more stable diastereomeric adsorbate formed between (S)-2 and the (S)-proline-derived CSP 1. No intermolecular NOEs are observed for corresponding mixtures of the chiral solvating agent (S)-3 and (R)-2, the enantiomer least retained on (S)-CSP 1.     

Evaluation of several two-step scoring functions based on linear interaction energy, effective ligand size, and empirical pair potentials for prediction of protein-ligand binding geometry and free energy

The performances of several two-step scoring approaches for molecular docking were assessed for their ability to predict binding geometries and free energies. Two new scoring functions designed for "step 2 discrimination" were proposed and compared to our CHARMM implementation of the linear interaction energy (LIE) approach using the Generalized-Born with Molecular Volume (GBMV) implicit solvation model. A scoring function S1 was proposed by considering only "interacting" ligand atoms as the "effective size" of the ligand and extended to an empirical regression-based pair potential S2. The S1 and S2 scoring schemes were trained and 5-fold cross-validated on a diverse set of 259 protein-ligand complexes from the Ligand Protein Database (LPDB). The regression-based parameters for S1 and S2 also demonstrated reasonable transferability in the CSARdock 2010 benchmark using a new data set (NRC HiQ) of diverse protein-ligand complexes. The ability of the scoring functions to accurately predict ligand geometry was evaluated by calculating the discriminative power (DP) of the scoring functions to identify native poses. The parameters for the LIE scoring function with the optimal discriminative power (DP) for geometry (step 1 discrimination) were found to be very similar to the best-fit parameters for binding free energy over a large number of protein-ligand complexes (step 2 discrimination). Reasonable performance of the scoring functions in enrichment of active compounds in four different protein target classes established that the parameters for S1 and S2 provided reasonable accuracy and transferability. Additional analysis was performed to definitively separate scoring function performance from molecular weight effects. This analysis included the prediction of ligand binding efficiencies for a subset of the CSARdock NRC HiQ data set where the number of ligand heavy atoms ranged from 17 to 35. This range of ligand heavy atoms is where improved accuracy of predicted ligand efficiencies is most relevant to real-world drug design efforts.