In silico studies on the interaction between bioactive ligands and ALK5, a biological target related to the cancer treatment

Studies have showed that there are many biological targets related to the cancer treatment, for example, TGF type I receptor (TGF-βRI or ALK5). The ALK5 inhibition is a strategy to treat some types of cancer, such as breast, lung, and pancreas. Here, we performed CoMFA and CoMSIA studies for 70 ligands with ALK5 inhibition. The internal validation for both models (q²_LOO = 0.887 and 0.822, respectively) showed their robustness, while the external validations showed their predictive power (CoMFA: r²_test = 0.998; CoMSIA: r²_test = 0.975). After all validations, CoMFA and CoMSIA maps indicated physicochemical evidences on the main factors involved in the interaction between bioactive ligands and ALK5. Therefore, these results suggest molecular modifications to design new ALK5 inhibitors.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Pharmacophore mapping,molecular docking and QSAR studies of structurally diverse compounds as CYP2B6 inhibitors

Pharmacophore mapping, molecular docking and quantitative structure–activity relationship (QSAR) studies were carried out for a structurally diverse set of 48 compounds as CYP2B6 inhibitors. The generated best pharmacophore hypotheses from the three methods of conformer generation (FAST, BEST and conformer algorithm based on energy screening and recursive buildup) indicate the importance of two features, namely, hydrogen bond acceptor [electron-rich centre] and ring aromaticity. The distance between the two centres of the important features for ideal inhibitors varied from 5.82 to 6.03 Å. The chemometric tools used for the QSAR analysis were genetic function approximation (GFA) and genetic partial least squares. The developed QSAR models indicate the importance of an electron-rich centre, size of molecule, impact of branching and ring system and distribution of charges in the molecular surface. The docking study confirms the importance of an electron-rich centre for binding with the iron atom of the cytochrome enzyme. A GFA model with spline option was found to be the best model based on internal validation as well as the r ² _{m (overall)} criterion (Q ² = 0.772, r ² _{m (overall)} = 0.774). According to the external prediction statistics (R ² _pred = 0.876), another GFA-derived model with spline option outperforms the remaining models.     

Deciphering the crucial molecular properties of a series of Benzothiazole Hydrazone inhibitors that targets anti-apoptotic Bcl-xL protein

The Bcl-2 family proteins are the central regulators of apoptosis. Due to its predominant role in cancer progression, the Bcl-2 family proteins act as attractive therapeutic targets. Recently, molecular series of Benzothiazole Hydrazone (BH) inhibitors that exhibits drug-likeness characteristics, which selectively targets Bcl-xL have been reported. In the present study, docking was used to explore the plausible binding mode of the highly active BH inhibitor with Bcl-xL; and Molecular Dynamics (MD) simulation was applied to investigate the stability of predicted conformation over time. Furthermore, the molecular properties of the series of BH inhibitors were extensively investigated by pharmacophore based 3D-QSAR model. The docking correctly predicted the binding mode of the inhibitor inside the Bcl-xL hydrophobic groove, whereas the MD-based free energy calculation exhibited the binding strength of the complex over the time period. Furthermore, the residue decomposition analysis revealed the major energy contributing residues – F105, L108, L130, N136, and R139 – involved in complex stability. Additionally, a six-featured pharmacophore model – AAADHR.89 – was developed using the series of BH inhibitors that exhibited high survival score. The statistically significant 3D-QSAR model exhibited high correlation co-efficient (R² = .9666) and cross validation co-efficient (Q² = .9015) values obtained from PLS regression analysis. The results obtained from the current investigation might provide valuable insights for rational drug design of Bcl-xL inhibitor synthesis.     

Molecular docking and pharmacophore model studies of Rho kinase inhibitors

Molecular docking and pharmacophore model approaches were used to characterise the binding features of four different series of Rho kinase (ROCK) inhibitors. Docking simulation of 20 inhibitors with ROCK was performed. The binding conformations and binding affinities of these inhibitors were obtained using AutoDock 4.0 software. The predicted binding affinities correlate well with the activities of these inhibitors (R ² = 0.904). 3D pharmacophore models were generated for ROCK based on highly active inhibitors implemented in Catalyst 4.11 program. The best pharmacophore model consists of one hydrogen bond acceptor feature and two hydrophobic features, and they all seemed to be essential for inhibitors in terms of their binding activities. It is anticipated that the findings reported in this paper may provide very useful information for designing new ROCK inhibitors.     

In silico approaches to evaluate the molecular properties of organophosphate compounds to inhibit acetylcholinesterase activity in housefly

Organophosphate compounds (OPC) have become the primary choice as insecticides and are widely used across the world. Additionally, OPCs were also commonly used as a chemical warfare agent that triggers a great challenge to public safety. Exposure of OPCs to human causes immediate excitation of cholinergic neurotransmission through transient elevation of synaptic acetylcholine (ACh) levels and accumulations. Likewise, prolonged exposure of OPCs can affect the processes in immune response, carbohydrate metabolism, cardiovascular toxicity, and several others. Studies revealed that the toxicity of OPCs was provoked by inhibition of acetylcholinesterase (AChE). Therefore, combined in silico approaches – pharmacophore-based 3D-QSAR model; docking and Molecular Dynamics (MD) – were used to assess the precise and comprehensive effects of series of known OP-derived compounds together with its ?log LD₅₀ values. The selected five-featured pharmacophore model – AAHHR.61 – displayed the highest correlation (R² = .9166), cross-validated coefficient (Q² = .8221), F = 63.2, Pearson-R = .9615 with low RMSE = .2621 values obtained using five component PLS factors. Subsequently, the well-validated model was then used as a 3D query to search novel OPCs using a high-throughput virtual screening technique. Simultaneously, the docking studies predicted the binding pose of the most active OPC in the MdAChE binding pocket. Additionally, the stability of docking was verified using MD simulation. The results revealed that OP22 and predicted lead compounds bound tightly to S315 of MdAChE through potential hydrogen bond interaction over time. Overall, this study might provide valuable insight into binding mode of OPCs and hit compounds to inhibit AChE in housefly.     

Barnacle settlement and the adhesion of protein and diatom microfouling to xerogel films with varying surface energy and water wettability

Previous work has shown that organosilica-based xerogels have the potential to control biofouling. In this study, modifications of chemistry were investigated with respect to their resistance to marine slimes and to settlement of barnacle cyprids. Adhesion force measurements of bovine serum albumin (BSA)-coated atomic force microscopy (AFM) tips to xerogel surfaces prepared from aminopropylsilyl-, fluorocarbonsilyl-, and hydrocarbonsilyl-containing precursors, indicated that adhesion was significantly less on the xerogel surfaces in comparison to a poly(dimethylsiloxane) elastomer (PDMSE) standard. The strength of adhesion of BSA on the xerogels was highest on surfaces with the highest and the lowest critical surface tensions, γ_C and surface energies, γ_S, and duplicated the ‘Baier curve’. The attachment to and removal of cells of the diatom Navicula perminuta from a similar series of xerogel surfaces were examined. Initial attachment of cells was comparable on all of the xerogel surfaces, but the percentage removal of attached cells by hydrodynamic shear stress increased with γ_C and increased wettability as measured by the static water contact angle, θ_Ws, of the xerogel surfaces. The percentage removal of cells of Navicula was linearly correlated with both properties (R ² = 0.74 for percentage removal as a function of θ_Ws and R ² = 0.69 for percentage removal as a function of γ_C). Several of the aminopropylsilyl-containing xerogels showed significantly greater removal of Navicula compared to a PDMSE standard. Cypris larvae of the barnacle B. amphitrite showed preferred settlement on hydrophilic/higher energy surfaces. Settlement was linearly correlated with θ_Ws (R ² = 0.84) and γ_C (R ² = 0.84). Hydrophilic xerogels should prove useful as coatings for boats in regions where fouling is dominated by microfouling (protein and diatom slimes).     

Differential interactions and structural stability of chitosan oligomers with human serum albumin and α-1-glycoprotein

Chitosan is a naturally occurring deacetylated derivative of chitin with versatile biological activities. Here, we studied the interaction of chitosan oligomers with low degree of polymerization such as chitosan monomer (CM), chitosan dimer (CD), and chitosan trimer (CT) with human serum albumin (HSA) a major blood carrier protein and α-1-glycoprotein (AGP). Since, HSA and AGP are the two important plasma proteins that determine the drug disposition and affect the fate of distribution of drugs. Fluorescence emission spectra indicated that CM, CD, and CT had binding constants of K_CM = 6.2 ± .01 × 10⁵ M^?1, K_CD = 5.0 ± .01 × 10⁴ M^?1, and K_CT = 1.6 ± .01 × 10⁶ M^?1_, respectively, suggesting strong binding with HSA. However, binding of chitooligomers with AGP was insignificant. Thermodynamic and molecular docking analysis indicated that hydrogen bonds and also hydrophobic interaction played an important role in stabilizing the HSA-chitooligomer complexes with free energies of ?7.87, ?6.35, and ?8.4?Kcal/mol for CM, CD, and CT, respectively. Further, circular dichroism studies indicated a minor unfolding of HSA secondary structure, upon interaction with chitooligomers, which are supported with fluctuations of root mean square deviation (RMSD) and radius of gyration (R_g) of HSA. Docking analysis revealed that all three chitooligomers were bound to HSA within subdomain IIA (Site I). In addition, RMSD and R_g analysis depicted that HSA-chitooligomer complexes stabilized at around 4.5 ns. These results suggest that HSA might serve as a carrier in delivering chitooligomers to target tissues than AGP which has pharmacological importance.     

Structural determinants of benzodiazepinedione/peptide-based p53-HDM2 inhibitors using 3D-QSAR, docking and molecular dynamics

As a tumor suppressor, p53 protein regulates the cell cycle and is involved in preventing tumorgenesis. The protein level of p53 is under the tight control of its negative regulator human double minute 2 (HDM₂) via ubiquitination. Therefore, the design of inhibitors of HDM₂ has attracted much interest of research on developing novel anticancer drugs. Presently, two classes of molecules, i.e., the 1,4-benzodiazepine-2,5-diones (BDPs) and N-Acylpolyamine (NAPA) derivatives were studied by three-dimensional quantitative structure–activity relationship (3D-QSAR) modeling approaches including the comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) as promising p53-HDM₂ inhibitors. Based on both the ligand-based and receptor-guided (docking) alignments, two optimal 3D-QSAR models were obtained with good predictive power of q ² = 0.41, r ² _pred = 0.60 for BDPs, and q ² = 0.414, r ² _pred = 0.69 for NAPA analogs, respectively. By analysis of the model and its related contour maps, it is revealed that the electrostatic interactions contributed much larger to the compound binding affinity than the steric effects. And the contour maps intuitively suggested where to modify the molecular structures in order to improve the binding affinity. In addition, molecular dynamics simulation (MD) study was also carried out on the dataset with purpose of exploring the detailed binding modes of ligand in the HDM₂ binding pocket. Based on the CoMFA contour maps and MD-based docking analyses, some key structural aspects responsible for inhibitory activity of these two classes of compounds were concluded as follows: For BDPs, the R₁ and R₃ regions should have small electronegativity groups; substituents R₂ and R₄ should be larger, and R₃ substituent mainly involves in H-bonds forming. For NAPA derivatives, bulky and electropositive groups in ring B and ring A, small substituent at region P is favorable for the inhibitory activity. The models and related information, we hope, may provide important insight into the inhibitor-p53-HDM₂ interactions and be helpful for facilitating the design of novel potent inhibitors.     

Design,synthesis, MTT assay,DNA interaction studies of platinum(II) complexes

The square planar Pt(II) complexes of the type [Pt(Lⁿ)(Cl₂)] (where Lⁿ = L^1?3 = thiophene-2-carboxamide derivatives and L^4?6 = thiophene-2-carbothioamide derivatives) have been synthesized and characterized by physicochemical and various spectroscopic studies. MIC method was employed to inference the antibacterial potency of complexes in reference to free ligands and metal salt. Characteristic binding constant (K_b) and binding mode of complexes with calf thymus DNA (CT-DNA) were determined using absorption titration (0.76–1.61 × 10⁵ M^?1), hydrodynamic chain length assay and fluorescence quenching analysis, deducing the partial intercalative mode of binding. Molecular docking calculation displayed free energy of binding in the range of –260.06 to –219.63 kJmol^?1. The nuclease profile of complexes towards pUC19 DNA shows that the complexes cleave DNA more efficiently compared to their respective metal salt. Cytotoxicity profile of the complexes on the brine shrimp shows that all the complex exhibit noteworthy cytotoxic activity with LC₅₀ values ranging from 7.87 to 15.94 μg/mL. The complexes have been evaluated for cell proliferation potential in human colon carcinoma cells (HCT 116) and IC₅₀ value of complexes by MTT assay (IC₅₀ = 125–1000 μg/mL).     

Crystal Structure and Conformation of 8-(2-Hydroxyethylamino) and 8-(Pyrrolidin-1-yl) Adenosines

In the course of investigation of 8-alkylamino substituted adenosines, the title compounds were synthesized as potential partial agonists for adenosine receptors. The structure determination of these compounds was carried out with the X-ray crystallography study. Crystals of 8-(2-hydroxyethylamino)adenosine are monoclinic, space group P 2₁; a = 7.0422(2), b = 11.2635(3), c = 8.9215(2) Å, β = 92.261(1)°, V = 707.10(3) ų, Z = 2; R-factor is 0.0339. The nucleoside is characterized by the anti conformation; the ribose ring has the C(2′)-endo conformation and gauche–gauche form across C(4′)–C(5′) bond. The molecular structure is stabilized by intramolecular hydrogen bond of N–H·O type. Crystals of 8-(pyrrolidin-1-yl)adenosine are monoclinic, space group C 2; a = 19.271(1), b = 7.3572(4), c = 11.0465(7) Å, β = 103.254(2)°, V = 1524.4(2) ų, Z = 4; R-factor is 0.0498. In this compound, there is syn conformation of the nucleoside; the ribose has the C(2′)-endo conformation and gauche–gauche form across C(4′)–C(5′) bond. The molecular structure is stabilized by intramolecular hydrogen bond of O–H·N type. For both compounds, the branching net of intermolecular hydrogen bonds occur in the crystal structures.     

The molecular basis for the inhibition of phosphodiesterase-4D by three natural resveratrol analogs. Isolation,molecular docking,molecular dynamics simulations,binding free energy,and bioassay

The phosphodiesterase-4 (PDE4) enzyme is a promising therapeutic target for several diseases. Our previous studies found resveratrol and moracin M to be natural PDE4 inhibitors. In the present study, three natural resveratrol analogs [pterostilbene, (E)-2′,3,5′,5-tetrahydroxystilbene (THSB), and oxyresveratrol] are structurally related to resveratrol and moracin M, but their inhibition and mechanism against PDE4 are still unclear. A combined method consisting of molecular docking, molecular dynamics (MD) simulations, binding free energy, and bioassay was performed to better understand their inhibitory mechanism. The binding pattern of pterostilbene demonstrates that it involves hydrophobic/aromatic interactions with Phe340 and Phe372, and forms hydrogen bond(s) with His160 and Gln369 in the active site pocket. The present work also reveals that oxyresveratrol and THSB can bind to PDE4D and exhibits less negative predicted binding free energies than pterostilbene, which was qualitatively validated by bioassay (IC₅₀ = 96.6, 36.1, and 27.0 μM, respectively). Additionally, a linear correlation (R² = 0.953) is achieved for five PDE4D/ligand complexes between the predicted binding free energies and the experimental counterparts approximately estimated from their IC₅₀ values (≈RT ln IC₅₀). Our results imply that hydrophobic/aromatic forces are the primary factors in explaining the mechanism of inhibition by the three products. Results of the study help to understand the inhibitory mechanism of the three natural products, and thus help the discovery of novel PDE4 inhibitors from resveratrol, moracin M, and other natural products.     

In silico studies on p21-activated kinase 4 inhibitors: comprehensive application of 3D-QSAR analysis,molecular docking,molecular dynamics simulations,and MM-GBSA calculation

Abstract

P21-activated kinase 4 (PAK4) is a serine/threonine protein kinase, which is associated with many cancer diseases, and thus being considered as a potential drug target. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations were performed to explore the structure-activity relationship of a series of pyrropyrazole PAK4 inhibitors. The statistical parameters of comparative molecular field analysis (CoMFA, Q ² = 0.837, R ² = 0.990, and R ² _pred = 0.967) and comparative molecular similarity indices analysis (CoMSIA, Q ² = 0.720, R ² = 0.972, and R ² _pred = 0.946) were obtained from 3D-QSAR model, which exhibited good predictive ability and significant statistical reliability. The binding mode of PAK4 with its inhibitors was obtained through molecular docking study, which indicated that the residues of GLU396, LEU398, LYS350, and ASP458 were important for activity. Molecular mechanics generalized born surface area (MM-GBSA) method was performed to calculate the binding free energy, which indicated that the coulomb, lipophilic and van der Waals (vdW) interactions made major contributions to the binding affinity. Furthermore, through 100?ns MD simulations, we obtained the key amino acid residues and the types of interactions they participated in. Based on the constructed 3D-QSAR model, some novel pyrropyrazole derivatives targeting PAK4 were designed with improved predicted activities. Pharmacokinetic and toxicity predictions of the designed PAK4 inhibitors were obtained by the pkCSM, indicating these compounds had better absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. Above research provided a valuable insight for developing novel and effective pyrropyrazole compounds targeting PAK4.     

In silico studies on potential MCF-7 inhibitors: a combination of pharmacophore and 3D-QSAR modeling,virtual screening,molecular docking,and pharmacokinetic analysis

Gallic acid and its derivatives exhibit a diverse range of biological applications, including anti-cancer activity. In this work, a data-set of forty-six molecules containing the galloyl moiety, and known to show anticarcinogenic activity against the MCF-7 human cancer cell line, have been chosen for pharmacophore modeling and 3D-Quantitative Structure Activity Relationship (3D-QSAR) studies. A tree-based partitioning algorithm has been used to find common pharmacophore hypotheses. The QSAR model was generated for three, four, and five featured hypotheses with increasing PLS factors and analyzed. Results for five featured hypotheses with three acceptors and two aromatic rings were the best out of all the possible combinations. On analyzing the results, the most robust (R^2?=?.8990) hypothesis with a good predictive power (Q^2?=?.7049) was found to be AAARR.35. A good external validation (R² = .6109) was also obtained. In order to design new MCF-7 inhibitors, the QSAR model was further utilized in pharmacophore-based virtual screening of a large database. The predicted IC₅₀ values of the identified potential MCF-7 inhibitors were found to lie in the micromolar range. Molecular docking into the colchicine domain of tubulin was performed in order to examine one of the probable mechanisms. This revealed various interactions between the ligand and the active site protein residues. The present study is expected to provide an effective guide for methodical development of potent MCF-7 inhibitors.     

Studies of H4R antagonists using 3D-QSAR,molecular docking and molecular dynamics

Three-dimensional quantitative structure–activity relationship studies were performed on a series of 88 histamine receptor 4 (H4R) antagonists in an attempt to elucidate the 3D structural features required for activity. Several in silico modeling approaches, including comparative molecular field analysis (CoMFA), comparative similarity indices analysis (CoMSIA), molecular docking, and molecular dynamics (MD), were carried out. The results show that both the ligand-based CoMFA model (Q ² = 0.548, R _ncv² = 0.870, R _pre² = 0.879, SEE = 0.410, SEP = 0.386) and the CoMSIA model (Q ² = 0.526, R _ncv² =0.866, R _pre² = 0.848, SEE = 0.416, SEP = 0.413) are acceptable, as they show good predictive capabilities. Furthermore, a combined analysis incorporating CoMFA, CoMSIA contour maps and MD results shows that (1) compounds with bulky or hydrophobic substituents at positions 4–6 in ring A (R2 substituent), positively charged or hydrogen-bonding (HB) donor groups in the R1 substituent, and hydrophilic or HB acceptor groups in ring C show enhanced biological activities, and (2) the key amino acids in the binding pocket are TRP67, LEU71, ASP94, TYR95, PHE263 and GLN266. To our best knowledge, this work is the first to report the 3D-QSAR modeling of these H4R antagonists. The conclusions of this work may lead to a better understanding of the mechanism of antagonism and aid in the design of new, more potent H4R antagonists.     

3D-QSAR (CoMFA,CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q² = 0.802, r²_ncv = 0.979, and the best CoMSIA model has q² = 0.799, r²_ncv = 0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. All the results can provide us more useful information for our further drug design.     

Substituted α-mercaptoketones,new types of specific neprilysin inhibitors

New neprilysin inhibitors containing an α-mercaptoketone HSC(R¹R²)CO group, as zinc ligand were designed. Two parameters were explored for potency optimization: the size of the inhibitor which could interact with the S₁, S₁′ or S₂′ domain of the enzyme and the nature of the substituents R¹, R² of the mercaptoketone group. Introduction of a cyclohexyl chain in R¹, R² position and a (3-thiophen)benzyl group in position R³ (compound 12n) yielded to the most potent inhibitor of this series with a Ki value of 2 ± 0.3 nM. This result suggests that this new inhibitor interacts within the S₁, S₁′ domain of NEP allowing a pentacoordination of the catalytic Zn²⁺ ion by the mercaptoketone moiety.