N. S. Hari Narayana Moorthy‡
*,
Natércia F. Brás‡,
Maria J. Ramos and
Pedro A. Fernandes*
REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, 687, Rua do Campo Alegre, 4169-007, Porto, Portugal. E-mail: hari.moorthy@fc.up.pt; hari.nmoorthy@gmail.com; pafernan@fc.up.pt; Fax: +351-220 402 506; Tel: +351-220 402 506
First published on 15th April 2014
Virtual screening studies were performed with docking and pharmacophore methods on a natural products (NPs) data set to investigate its inhibitory effect on cardiovascular targets such as renin (REN) and angiotensin-converting enzymes (ACE). Conformers obtained from the ligand–protein complex (from pdb) and the flexialigned structure of the reference compounds were used to generate pharmacophore query models. The results derived from the analysis on natural product data set showed that the compounds such as Nat-4, Nat-59, Nat-99 and Nat-141 provided good results against REN and Nat-7, Nat-6, Nat-31, Nat-59 and Nat-61 are considered as good HITs against ACE. The present studies revealed that the compound Nat-59 has remarkable interaction with both the studied targets (ACE and REN). Also we have observed better effects for those compounds on the hERG target through the pharmacophore based virtual screening method, revealing that these can be used as antiarrythmic agents. A new lead compound (NLC-1) was designed on the basis of the pharmacophoric features of the selected HITs, provided good HIT properties by pharmacophore and docking studies. Molecular dynamics simulations were performed on the compound NLC-1 in order to find out the binding features of this molecule on the targets (ACE and REN). These results concluded that the natural compound identified from these studies and the designed compound (NLC-1) can have multiple activities on cardiovascular targets.
The RAAS pathway has primary role in controlling blood pressure through the production of angiotensin-II (AT-II) from angiotensin-I (AT-I) by the action of angiotensin-I converting enzyme (ACE).5 Initially, the RAAS proteolytic cascade starts with the release of the enzyme renin (REN) (EC 3.4.23.15) by the juxtaglomerular cells located in the kidney. REN activates a plasma protein named angiotensinogen (452 residues) by hydrolyzing the peptidic bond between 10th and 11th positions, releasing a decapeptide known as AT-I. Finally, the inactive AT-I is further cleaved into the octapeptide AT-II by the metalloprotease ACE.2,6 ACE hydrolyzes both the inactive AT-I into vasoconstrictor AT-II and the vasodilator bradykinin into an inactive metabolite leading to blood pressure up-regulation.5,7,8
The REN and ACE play important role for the production of the vasoconstrictor AT-II and the inhibitors discovered for these targets are used for the treatment of hypertension. However, the inhibition of ACE causes REN overexpression, which in part compensates the effect of inhibiting ACE and AT-I functions.4 The inhibition of REN is much more attractive, it catalyzes the rate-limiting step of this enzymatic cascade and angiotensinogen is the only known substrate. Therefore, drugs that inhibit REN are expected to have several advantages and they should be an attractive antihypertensive strategy with fewer side effects.5,9–11 The peptidomimetic REN inhibitor, for example, remikiren (RO 42-5892) was discovered during the 1980s. More recent research has lead to the discovery of several classes of non-peptidic REN inhibitors.12,13 In 2007, the Food and Drug Administration (FDA) approved the first and unique direct REN inhibitor, aliskiren (commercialized by Novartis as Tekturna). The first breakthrough in the area of ACE inhibitor was captopril, which was launched in 1981 for hypertension management.14–16
The development of novel class of antihypertensive drugs include the development of single compound capable of inhibiting more than one enzymatic activity simultaneously involved in hypertension pathophysiology is important nowadays.5,17 Hence, in the present investigation, we have used some computational methods to discover novel antihypertensive agents from natural molecules (natural products (NPs)). NPs from different structural classes have been investigated for its action on REN and ACE enzymes. NPs traditionally have played an important role in drug discovery and were the basis of most early medicines.18–20
However, in the drug discovery research, still there is a shortage of lead compounds progressing into clinical trials. This is especially in the therapeutic areas such as oncology, immunosuppression and metabolic diseases where NPs have played a central role in lead discovery. In the world market, NPs and NP-derived drugs were well represented in the top 35 selling ethical drugs in between 2000 and 2002.18,21 Hence, the virtual screening analysis of NPs for the therapeutic targets provides significant lead for the design of novel molecules. In the present investigation, we have performed virtual screening of NPs on REN and ACE targets using computational techniques (pharmacophore, docking, molecular dynamic (MD) simulations and correlation analysis) (Table S1†). On the basis of the study results, we have designed a novel lead compound for further development.
The details of the pharmacophore queries used for the pharmacophore based virtual screening studies are given below. The pharmacophore contours and its radius used in the pharmacophore query models development are provided in Table 1.
Pharmacophores | Code | Contour | Radius |
---|---|---|---|
Pharmacophore-1 | F1 | Acc | 2 |
F2 | Acc | 1.5 | |
F3 | ML|Acc2 | 1.5 | |
F4 | Aro|Hyd | 2 | |
F5 | Don | 1 | |
Pharmacophore-2 | F1 | Acc | 1 |
F2 | Acc | 1 | |
F3 | ML|Acc2 | 1.5 | |
F4 | Acc | 1.3 | |
F5 | Aro|Hyd | 1.4 | |
Pharmacophore-3 | F1 | Ani|Acc|Acc2|Don2 | 2.5 |
F2 | Cat|Acc2|Don | 1.5 | |
F3 | Acc|Don2 | 2.5 | |
F4 | Acc2|Don | 1.5 | |
F5 | Hyd|Acc2|Don2 | 2.5 | |
+V1 h | Excluded volume | 1.5 | |
Pharmacophore-4 | F1 | Ani|Acc|Acc2|Don2 | 2 |
F2 | Hyd|Acc | 1 | |
F3 | Hyd | 1.5 | |
F4 | Aro|Hyd|Acc|Don | 1.5 | |
F5 | Hyd|Acc2|Don2 | 1.5 | |
+V1 h | Exterior volume | 1.35 | |
Pharmacophore-5 | F1 | Don | 1.5 |
F2 | Don | 1.5 | |
F3 | Hyd|Aro | 1.3 | |
F4 | Hyd | 1.5 |
(1) In complex based pharmacophore analysis, the protein–ligand complexes obtained from protein data bank were used. The PDB structures such as 1O86 and 2V0Z were utilized for ACE (pharmacophore-1) and REN (pharmacophore-3) respectively.
(2) The flexialigned structure of the reference compounds such as captopril, enalapril, lisinopril and rimipril for ACE (pharmacophore-2) and compounds such as aliskiren, enalkiren, rimikiren and zankiren for REN (pharmacophore-4) were used to create pharmacophore query models.
(3) In order to analyze the effect of the molecules on hERG target, a pharmacophore query was generated with the inhibitor present in pdb 3O0U (pharmacophore-5).29
Subsequently, the docking parameters used in the initial studies executed correctly the binding mode and affinity of these known compounds. Hence, we have applied the same docking protocol to evaluate the binding interactions of all NPs data set compounds studied. The binding energy and inhibitory constant values obtained for the first twenty ligands are shown in Table 2.
REN | ACE | ||||||
---|---|---|---|---|---|---|---|
Inhibitor | Ki (μM) | ΔGbind (kcal mol−1) | Inhibitor | Ki (μM) | ΔGbind (kcal mol−1) | ||
Docking | Correlation analysis | Docking | Correlation analysis | ||||
Nat-20 | 1.46 × 10−2 | −10.69 | −9.1704 | Nat-34 | 7.10 × 10−4 | −12.48 | −10.3527 |
Nat-38 | 1.94 × 10−2 | −10.52 | −9.1704 | Nat-7 | 7.72 × 10−4 | −12.43 | −10.3527 |
Nat-39 | 4.15 × 10−2 | −10.07 | −9.1703 | Nat-6 | 1.20 × 10−3 | −12.17 | −10.3527 |
Nat-59 | 3.06 × 10−2 | −10.25 | −9.1703 | Nat-31 | 6.82 × 10−3 | −11.14 | −10.3523 |
NLC-1 | 6.33 × 10−2 | −9.82 | −9.1702 | Nat-59 | 7.54 × 10−3 | −11.08 | −10.3523 |
Nat-165 | 6.43 × 10−2 | −9.81 | −9.1702 | NLC-1 | 1.38 × 10−2 | −10.72 | −10.352 |
Nat-7 | 7.36 × 10−2 | −9.73 | −9.1702 | Nat-61 | 2.11 × 10−2 | −10.47 | −10.3516 |
Nat-19 | 8.43 × 10−2 | −9.65 | −9.1701 | Nat-69 | 2.26 × 10−2 | −10.43 | −10.3515 |
Nat-132 | 8.72 × 10−2 | −9.63 | −9.1701 | Nat-72 | 2.34 × 10−2 | −10.41 | −10.3515 |
Nat-14 | 1.63 × 10−1 | −9.26 | −9.1697 | Nat-10 | 2.46 × 10−2 | −10.38 | −10.3515 |
Nat-12 | 1.77 × 10−1 | −9.21 | −9.1697 | Nat-47 | 2.81 × 10−2 | −10.3 | −10.3512 |
Nat-77 | 2.21 × 10−1 | −9.08 | −9.1695 | Nat-165 | 2.86 × 10−2 | −10.29 | −10.3512 |
Nat-15 | 2.28 × 10−1 | −9.06 | −9.1694 | Nat-38 | 2.96 × 10−2 | −10.27 | −10.3512 |
Nat-99 | 2.28 × 10−1 | −9.06 | −9.1694 | Nat-39 | 3.45 × 10−2 | −10.18 | −10.3509 |
Nat-141 | 2.32 × 10−1 | −9.05 | −9.1694 | Nat-40 | 3.50 × 10−2 | −10.17 | −10.3509 |
Nat-23 | 2.53 × 10−1 | −9 | −9.1693 | Nat-67 | 3.69 × 10−2 | −10.14 | −10.3508 |
Nat-24 | 2.66 × 10−1 | −8.97 | −9.1693 | Nat-19 | 4.08 × 10−2 | −10.08 | −10.3506 |
Nat-81 | 2.89 × 10−1 | −8.92 | −9.1692 | Nat-24 | 4.08 × 10−2 | −10.08 | −10.3506 |
Nat-22 | 2.89 × 10−1 | −8.92 | −9.1692 | Nat-68 | 5.25 × 10−2 | −9.93 | −10.35 |
Nat-163 | 2.99 × 10−1 | −8.9 | −9.1691 | Nat-99 | 5.25 × 10−2 | −9.93 | −10.35 |
Nat-82 | 3.25 × 10−1 | −8.85 | −9.1690 | Nat-15 | 5.43 × 10−2 | −9.91 | −10.3499 |
![]() | ||
Fig. 2 Representation of REN enzyme, with both catalytic Asp residues in ball and sticks, the flap residues in green and the main non-polar residues around the active pocket in dark blue. |
The database used in the present study has 165 natural compounds, which are described in ESI in Table S1.† We have considered only the five top-ranked inhibitor molecules for further discussion to simplify our results. The natural compounds Nat-20, Nat-38, Nat-39, Nat-59 and Nat-165 are the best five inhibitors bound to REN active site (structures provided in Table 3). Their representations and the most important interactions established within the active site of REN are shown in Fig. 3. It has observed that the major interactions within the pocket site of REN occurred by hydrophobic and van der Waals (vdW) interactions, as well as strong hydrogen bridges established with both catalytic aspartate residues. The residues that interact with the natural compounds are the following: Val36, Asp38, Tyr83, Thr85, Pro118, Phe119, Phe124, Val127, Asp226, Ala229, Tyr231 and Met303 (residue number provided as per the human REN sequence number). The compound Nat-20 establishes a very short H-bridge (1.88 Å) with the negative charged Asp226 by its hydroxyl group, whilst the compound Nat-38 establishes two important H-bonds with the carboxylic group of catalytic Asp38 (2.37 Å) and the carbonyl side chain group of Thr85 (2.28 Å). Other strong hydrogen bridges are also established between a carbonyl and hydroxyl groups of compound Nat-39 to the carboxyl group of Asp38 (1.84 Å) and the OH group of Tyr231. Compounds Nat-59 and NLC-1 possess several hydroxyl groups, which establish several relevant H-bonds with the carboxylic groups of Asp38 and Asp226, the carbonyl side chain of Thr85 and the carbonyl backbone of Pro118. The oxygen atom in Nat-165 interacts with the carboxylic side chain of Asp38 by an H-bonding (3.17 Å). This oxygen atom is less reactive than the hydroxyl and carbonyl groups present in other molecules, which is probably the reason for its weakest binding energy when compared with the previously described compounds. In addition, the aromatic side chains of residues Tyr83, Phe119, Phe124 and Tyr231 make several vdW interactions and T-shaped π–π stacking interactions with different aromatic moieties present in these five NPs. Similar hydrophobic contacts also occurred between the Val36, Val127, Pro118 and Met303 residues to some non-polar groups of all five-ranked inhibitors.
![]() | ||
Fig. 3 Representation of the interactions established by the residues of REN enzyme to the best five-ranked inhibitors: compounds Nat-20, Nat-38, Nat-39, Nat-59, Nat-165 and NLC-1. |
In fact, the binding free energy values and binding pose of these five-ranked inhibitors are very similar (a difference of 0.88 kcal mol−1 between the ΔGbind of the first and fifth compound), which indicates the crucial importance of REN residues such as Val36, Asp38, Tyr83, Thr85, Pro118, Phe119, Phe124, Val127, Asp226, Ala229, Tyr231 and Met303 interact with these natural compounds.
The binding interaction derived for the NPs identical with the reported compounds in the literatures. In the reported literatures, Asp38 and Asp226 residues in REN make hydrogen bonding interactions with the hydrogen bond acceptor or donor groups present in the molecules. Depending upon the nature of the inhibitors present in the active site, the residues such as Tyr14, Arg74, Tyr75, Ser76, Tyr77, Phe112, Gln128, Ala218, etc. are predominantly interacted with the molecules by hydrogen bonds or π–CH interactions or CH–CH or π–π stacking interactions.41–44
Fig. 5 shows the interactions established between these compounds and the residues present in the binding pocket of this enzyme. The compound Nat-34 is the one that establish more interactions with ACE protein. As seen, the NH2 group makes a hydrogen bond with the hydroxyl group of Tyr360, a carbonyl group establish hydrogen bonds with the Tyr394, Glu403, His410 and a NH group interacts with the backbone carbonyl of Ala356. Its carboxylate group strongly interacts with the Zn2+ ion as well as with the OH group of Tyr523. In addition, the benzene ring possesses π–π stacking contacts with the aromatic side chain of Phe391 and Phe512 residues. The second-ranked inhibitor establishes H-bridges with CO and NH backbone groups of Ser355 and Ala356, respectively, as well as hydrophobic contacts with the Phe391. This compound also has a carboxylate group that coordinates in a bidentated form to Zn2+ ion. Compound Nat-6 interacts with Asp358, Ala356 and Tyr523 by H-bonds, as well as its COO− group is bidentated to Zn2+. The binding mode of compound Nat-31 is very similar to the compound Nat-34. Thus, it interacts with Tyr360, Tyr394, His410, Glu403 and Tyr523 by hydrophilic contacts, while the aromatic residues Phe391, Phe512, Tyr523 and His410 make hydrophobic contacts. Similarly to the REN enzyme, the compounds Nat-59 and NLC-1 also interact well with this metalloprotease (ACE). It makes T-shaped π–π stacking interactions with Phe512 and Tyr523; hydrophilic contacts with Asn70, Ser355, Lys368, Glu384 and Tyr520; and bidentated coordination with Zn2+ ion by one hydroxyl group and an oxygen-ring atom. All these five compounds show similar interactions with the ACE active site, which agrees with the proximity of their binding mode conformations and binding energies values (a difference of 1.76 kcal mol−1 between the binding energies of the first and fifth compounds). Therefore, the main hydrophilic amino acids those contribute to an efficient inhibitor binding are the following: Asn70, Ser355, Asp358, Lys368, His383, Glu384, His387, Glu403, His410, Glu411 and Glu485; whilst the crucial aromatic residues are the Tyr360, Phe391, Tyr394, Phe512, Tyr520 and Tyr523.
![]() | ||
Fig. 5 Representation of the interactions established by the residues of ACE enzyme to the best five-ranked inhibitors: compounds Nat-34, Nat-7, Nat-6, Nat-31, Nat-59 and NLC-1. |
The binding interaction analysis of different ACE inhibitors on the catalytic site of N and C-domain of the ACE reported on the literature showed that the residues such as Gln281, His353, His513, Tyr520 and Tyr523 are predominantly interacted with the inhibitors in C-domain ACE. However, in the N-domain ACE, the residues such as Gln259, His331, His491, Tyr498 and Tyr501 are present in the active site for the interaction. This shows that mainly the same amino acid residues present in both the domain of the enzyme for the interaction.45 The study performed by Tzakos showed that the residues such as Gln281, His353, His513, Tyr520 and Tyr523 importantly present for the ligand interaction. They reported that the interactions governing ACE_C domain and the ligand are: a salt bridge between Asp377 and Glu162, and the NH2 group, a hydrogen bond of the inhibitor with Gln281, and a π–π stacking interaction of Phe391 with aromatic group. In ACE_N, inhibitors interact with Glu431, Tyr369 and Arg381 by hydrogen bonds and a salt bridge between the carboxy group in the inhibitors and Arg500. Some specific residues also present in different compounds for the interaction.45–47 The results derived from the present study also coincided with the reported results in the literature. In our studies His383, His387, Glu384, Phe391, Tyr394, His410, Glu411, Phe512, Tyr520 and Tyr523 are preferably interacted with the inhibitors.
The pharmacophore-1 (receptor ACE) and the pharmacophore-2 (reference ACE) developed with 5 pharmacophore contours features. The pharmacophore-3 (receptor REN) and the pharmacophore-4 (reference REN) also possessed 5 contour features and a volume parameter such as excluded volume in pharmacophore-3 and exterior volume in pharmacophore-4 models. The generated pharmacophore queries (pharmacophores 1–4) have been used to perform virtual screening, in order to derive HIT molecules from the NPs data set. The detail of the significant HIT compounds obtained from the pharmacophore based virtual screening study is given in Tables 4 and 5. Among the 165 NP compounds present in the data set, 6 compounds are identified as significant compounds through pharmacophore-1 and 16 compounds identified through pharmacophore-2. It is interesting that the HITs identified through the pharmacophore-1 are also present in the pharmacophore-2. The compounds Nat-61, Nat-59, Nat-33, Nat-82, Nat-3 and Nat-31 are considered as significant HITs from pharmacophore pharmacophore-1 and 2. These compounds exhibited the RMSD values ≤1. The HITs identified through pharmacophore-3 and 4 (REN) are provide in Table 5. Those pharmacophore query models provided four compounds each as significant HITs. It is interesting that the compound Nat-59 is selected as one of the significant compounds against the REN target. This compound possessed the RMSD value <2 for both the queries (pharmacophore-3 and 4). This compound (Nat-59) has been selected as significant HIT through all the pharmacophore queries (1–4) for REN and ACE.
Comp. code | RMSD | |
---|---|---|
Pharmacophore-1 | Pharmacophore-2 | |
Nat-61 | 0.8804 | 0.8541 |
Nat-59 | 0.8867 | 0.8470 |
Nat-33 | 0.9612 | 0.7533 |
Nat-82 | 0.9814 | 0.7771 |
Nat-3 | 1.0569 | 0.7465 |
Nat-31 | 1.0586 | 0.6173 |
Nat-4 | — | 0.4865 |
Nat-6 | — | 0.5959 |
Nat-34 | — | 0.6091 |
Nat-60 | — | 0.7440 |
Nat-27 | — | 0.7932 |
Nat-133 | — | 0.8294 |
Nat-24 | — | 0.9210 |
Nat-45 | — | 0.9335 |
Nat-7 | — | 0.9600 |
Nat-44 | — | 0.9788 |
NLC-1 | 0.8743 | 0.6832 |
Comp. code | RMSD | |
---|---|---|
Pharmacophore-3 | Pharmacophore-4 | |
Nat-10 | 1.6178 | — |
Nat-4 | 1.6857 | — |
Nat-59 | 1.7001 | 1.2119 |
Nat-99 | 1.4345 | — |
Nat-141 | — | 1.0065 |
Nat-46 | — | 1.1549 |
Nat-71 | — | 1.0943 |
NLC-1 | 1.7600 | 1.3456 |
The virtual screening results (docking and pharmacophore studies) on the REN and ACE showed that some compounds possessed significant docking score and those compounds also identified as good HITs through the pharmacophore analysis. Compounds such as Nat-4, Nat-59, Nat-99 and Nat-141 are considered as good compounds through both methods against REN target. The compounds Nat-7, Nat-6, Nat-31, Nat-59 and Nat-61 are considered as good HITs against ACE. These derived virtual screening results on REN and ACE by docking and pharmacophore analyses describe that the compound Nat-59 is a significant compound against both targets. Structure of the significant compounds selected from the studies is provided in Table 3.
The compound especially Nat-59 has significant effect on the active site of REN and ACE. This has been confirmed by the docking and pharmacophore analyses. Also the compounds such as Nat-6, Nat-7, Nat-38, Nat-39 and Nat-99 also possessed interaction with REN and ACE targets. The pharmacophore analysis also showed that Nat-31 has been identified as good HIT for ACE. On the basis of the structural description and pharmacophore features of Nat-59, Nat-31, Nat-6, Nat-38 and Nat-34, we have designed a new lead compound against ACE and REN targets. According to docking results, the developed new compound NLC-1 has significant binding free energy (−10.25 kcal mol−1) comparing to the reference compound aliskiren (−11.49 kcal mol−1) on REN target. In ACE target, NLC-1 possessed the binding free energy of −11.08 kcal mol−1, which is comparable to the reference ACE inhibitors. The pharmacophore analysis (pharmacophore 1–4) also selected this compound as significant HIT, the RMSD value of the compound is provided in Tables 4 and 5. This confirms that NLC-1 has significant interaction with REN and ACE targets. The interaction of NLC-1 revealed that mode of interactions is same as Nat-59 against both targets (Fig. 2 and 4).
We have also performed pharmacophore analysis of the NPs against one of the cardiovascular target human ether-a-go-go related gene (hERG). The pharmacophore analysis performed on the compounds showed that 72 compounds were identified as HITs for the target. It is interesting that those compounds identified as significant compounds from docking and pharmacophore analysis for ACE and REN inhibitory effects also considered as significant hit. The structural features of these molecules are comparatively similar to the features needed for the hERG blocking activities. Earlier reports on hERG blockers showed that the flexibility of the molecules orients the substituent towards the active site. Low polarizable groups and aromatic rings (hydrophobicity) provide optimum vdW surface properties required for the hERG blocking activity. The identified HITs also possessed flexible bonds, aromatic rings and polarizable groups in their structure.48–50 Hence, these compounds have significant interactions with the hERG protein and can used as antiarrythmic agents. The compound NLC-1 has been considered as significant hit for all cardiovascular targets. The conformers selected through the pharmacophore based virtual screening analysis are provided in Table 6.
![]() | ||
Fig. 7 RMSD values obtained for the backbone of REN (black) and ACE (red) proteins, as well as for the compound NLC-1 bound to REN (green) and ACE (blue). |
Fig. 8 shows the RMSF values obtained for the REN and ACE bound enzymes by residue. Concerning to REN, the residues such as Leu1-Thr6, Gln150-Lys154, Glu167-Ser171, Glu187-Phe190, Val210-Thr214, Gly262-Pro266 and Tyr290-Lys294 show the highest RMSF values and possess high mobility, which probably may promote an easier inhibitor binding. It is expected that all these key residues will assist the flap-opening and closure to allow access for the substrate to the REN active pocket, enabling its catalytic reaction. The compound NLC-1 occupies the active site and the flap is stably closed over part of it, the RMSF values for the flap (residues 75–85) are indeed low (approximately 0.6 Å). This fact agrees with the higher stabilization provided by the establishment of several and strong interactions between the NLC-1 and the enzymatic active site.
Similarly analyzing the RMSF values for ACE, it was observed higher values for residues Val1-Lys9, Asn68-Asn72, His116-Ser120, Lys137-Trp143, Ala212-Pro234, Val254-Pro260, Ser398-Asn408, Thr440-Asn443, Cys459-Asp468, Gln502-Pro509, Ser518-Glu520 and Lys560-Pro586. Since, these groups are the ones that mostly modify their position during the simulation studied, that probably help for the binding of substrate/inhibitor to the enzymes.
In order to examine the parts (atom or groups) of compound (NLC-1) have significant flexibility, the RMSF values were also calculated for this molecule when it is bound to REN and ACE. Fig. 9 shows these RMSF values and the numeration of all heavy atoms of NLC-1. It is possible to observe several peaks that have RMSF values above 0.5 Å. The regions of NLC-1 with the highest RMSF values correspond to several hydroxyl groups, which have higher flexibility during the MD simulation, due to their rotatable bonds that have high conformational flexibility to establish H-bonds with the active pocket residues of REN and ACE.
![]() | ||
Fig. 9 RMSF values obtained for NLC-1 compound by atom bound to REN (blue) and ACE (red), as well as the representation of its 2D structure with the respective numeration of the heavy atoms. |
Structural information on the REN and ACE binding sites, in particular the H-bonds and vdW contacts established between the NLC-1 and the neighbour residues during the MD simulations were also analyzed. It is noteworthy that the NLC-1 compound has mainly hydrophobic character due to its aromatic rings, but it also possesses several hydroxyl groups. Fig. 10 shows the structures that are closest to the average structures of the REN:NCL-1 and ACE:NLC-1 complexes. In relation to REN, it is verified that several hydrophobic residues (Val36, Pro118, Ala122, Val127, Ala229, Leu224, Ile305 and Ala314) establish dispersive contacts with the NLC-1 and may be responsible for its efficient binding. It is also observed a T-shaped π–π stacking contacts between the NLC-1 rings and the aromatic side chains of Tyr83, Phe119 and Phe124 residues, as well as H-bonds between some hydroxyl groups of this molecule and the backbone NH groups of Ser84 and Gly86 (5% and 9% of occupancies), the backbone NH and OH groups of Thr85 (46% and 61% of occupancies), the OH groups of Ser230 and Thr309 (28% and 11% of occupancies) and the unprotonated carboxylic group of Asp226 with an average distance of 1.83 ± 0.45 Å and 62% of occupancy.
![]() | ||
Fig. 10 Representation of the structures that is closest to the average structures of the REN:NLC-1 (left) and ACE:NLC-1 (right) complexes. |
The binding pocket of ACE is a narrow tunnel highly populated by aromatic residues such as Trp318, Tyr484, His371, Phe418, Phe473 and His474, which establish strong π–π stacking contacts with the NLC-1 aromatic rings. These non-polar groups of NLC-1 also made dispersive contacts with the non-polar Val341 and Val479 residues of ACE. However, the most important contact of this compound with the ACE active site has its direct and strong interaction with the metallic Zn2+ cation with an average distance of 2.00 ± 0.38 Å. Furthermore, several H-bonds are established between the OH groups of NLC-1 compound and the amino groups of Gln242 (14% of occupancy), Arg483 (28% of occupancy), His371 (10% of occupancy) and Lys472 (5% of occupancy), the negatively charged carboxylic group of Glu337 (89% of occupancy), the OH groups of Tyr480 and Tyr484 (8% and 7% of occupancies) and the backbone carbonyl group of Ala315 (86% of occupancy).
Our results highlight the importance, for the inhibitor efficiency, with the short and strong H-bonds, dispersive and π–π stacking contacts, demonstrating that all these interactions contribute to an efficient binding of NLC-1 to the REN and ACE active sites, and allow the stabilization of the active site scaffold of these two complexes.
Statistical Parameters | REN | ACE | ||||
---|---|---|---|---|---|---|
Model 1 | Model 2 | Model 3 | Model A | Model B | Model C | |
a R = correlation coefficient, Q2 = crossvalidated correlation coefficient, LOO = leave one out. | ||||||
Number of compounds | 86 | 55 | 63 | 63 | 46 | 43 |
R values | 0.6405 | 0.7027 | 0.6825 | 0.6003 | 0.7038 | 0.6710 |
Ftest | 58.4207 | 51.6971 | 53.1902 | 34.3710 | 43.1951 | 33.5855 |
Ttest | −35.38 | −32.24 | −30.16 | −38.5000 | −39.85 | −35.02 |
Q2LOO | 0.5116 | 0.4981 | 0.5139 | 0.6309 | 0.4872 | 0.4658 |
Q2test | — | 0.5024 | 0.4625 | — | 0.5455 | 0.4165 |
![]() | ||
Fig. 14 Scatterplot of ΔGbind (kcal mol−1) (experimental vs. predicted) of REN and ACE inhibitors from the correlation studies. |
In summary, it has verified that the active sites of both REN and ACE enzymes possess predominantly hydrophilic charged amino acids (mainly negatively) such as Glu and Asp and aromatic residues (Phe, Tyr, etc.) with hydrophobic characteristics. Therefore, the main interactions between all inhibitors and both enzymatic classes occurred by hydrogen bonds, hydrophobic π–π stacking and very close contacts with the Zn2+ ion (for ACE). The pharmacophore queries developed for the pharmacophore based virtual screening also showed predominantly hydrophilic contours such as Acc, Acc2, Don, Don2, etc. and aromatic contour (Aro/Hyd). These query contours also coincided with the binding mode studies performed with the docking analysis. Herein, comprehensive molecular docking and pharmacophore studies about the binding mode of several NPs into the REN and ACE active sites were analysed.
The results obtained indicate that the scoring function of AutoDock is adequate in predicting the correct binding modes of well-known inhibitors in relation to their available crystallographic structures. Subsequently, similar protocol was applied to correctly describe the binding mode of several NPs into the REN and ACE active sites. Hence these analysis results such as pharmacophore and docking analysis provide active site characters of REN and ACE. Our results suggest that the best compound to inhibit both enzymatic targets is the compound Nat-59. The designed compound NLC-1 also has the same binding interaction with both targets as Nat-59. The MD simulations performed on the NLC-1 with both the enzymes also confirm that they bind with the same residues in the enzymes. The new compound shares the pharmacophore of Nat-59, Nat-6, Nat-7, Nat-31 and Nat-34. Therefore, these findings are particularly relevant for the design of novel compounds and these ligand and structure-based approaches may lead to improved inhibitors.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4ra00856a |
‡ Both authors are equally contributed in this work. |
This journal is © The Royal Society of Chemistry 2014 |