Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor

Structure-based identification of novel scaffolds as potential HIV-1 entry inhibitors involving CCR5

C-C chemokine receptor 5 (CCR5), which is part of the chemokine receptor family, is a member of the G protein-coupled receptor superfamily. The interactions of CCR5 with HIV-1 during viral entry position it as an effective therapeutic target for designing potent antiviral therapies. The small-molecule Maraviroc was approved by the FDA as a CCR5 drug in 2007, while clinical trials failure has characterised many of the other CCR5 inhibitors. Thus, the continual identification of potential CCR5 inhibitors is, therefore, warranted. In this study, a structure-based discovery approach has been utilised to screen and retrieved novel potential CCR5 inhibitors from the Asinex antiviral compound (∼ 8,722) database. Explicit lipid-bilayer molecular dynamics simulation, in silico physicochemical and pharmacokinetic analyses, were further performed for the top compounds. A total of 23 structurally diverse compounds with binding scores higher than Maraviroc were selected. Subsequent molecular dynamics (MD) simulations analysis of the top four compounds LAS 51495192, BDB 26405401, BDB 26419079, and LAS 34154543, maintained stability at the CCR5 binding site. Furthermore, these compounds made pertinent interactions with CCR5 residues critical for the HIV-1 gp120-V3 loop binding such as Trp86, Tyr89, Phe109, Tyr108, Glu283 and Tyr251. Additionally, the predicted in silico physicochemical and pharmacokinetic descriptors of the selected compounds were within the acceptable range for drug-likeness. The results suggest positive indications that the identified molecules may represent promising CCR5 entry inhibitors. Further structural optimisations and biochemical testing of the proposed compounds may assist in the discovery of effective HIV-1 therapy.Communicated by Ramaswamy H. Sarma.

Molecular modeling studies of dihydro-alkyloxy-benzyl-oxopyrimidines (DABOs) as non-nucleoside inhibitors of HIV-1 reverse transcriptase using 3D-QSAR, Topomer CoMFA and molecular docking simulations

The human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) is generally regarded as a target for the treatment of acquired immune deficiency syndrome (AIDS).

Isolation, purification and molecular mechanism of a peanut protein-derived ACE-inhibitory peptide

Although a number of bioactive peptides are capable of angiotensin I-converting enzyme (ACE) inhibitory effects, little is known regarding the mechanism of peanut peptides using molecular simulation. The aim of this study was to obtain ACE inhibiting peptide from peanut protein and provide insight on the molecular mechanism of its ACE inhibiting action. Peanut peptides having ACE inhibitory activity were isolated through enzymatic hydrolysis and ultrafiltration. Further chromatographic fractionation was conducted to isolate a more potent peanut peptide and its antihypertensive activity was analyzed through in vitro ACE inhibitory tests and in vivo animal experiments. MALDI-TOF/TOF-MS was used to identify its amino acid sequence. Mechanism of ACE inhibition of P8 was analyzed using molecular docking and molecular dynamics simulation. A peanut peptide (P8) having Lys-Leu-Tyr-Met-Arg-Pro amino acid sequence was obtained which had the highest ACE inhibiting activity of 85.77% (half maximal inhibitory concentration (IC50): 0.0052 mg/ml). This peanut peptide is a competitive inhibitor and show significant short term (12 h) and long term (28 days) antihypertensive activity. Dynamic tests illustrated that P8 can be successfully docked into the active pocket of ACE and can be combined with several amino acid residues. Hydrogen bond, electrostatic bond and Pi-bond were found to be the three main interaction contributing to the structural stability of ACE-peptide complex. In addition, zinc atom could form metal-carboxylic coordination bond with Tyr, Met residues of P8, resulting into its high ACE inhibiting activity. Our finding indicated that the peanut peptide (P8) having a Lys-Leu-Tyr-Met-Arg-Pro amino acid sequence can be a promising candidate for functional foods and prescription drug aimed at control of hypertension.

Integrated computational tools for identification of CCR5 antagonists as potential HIV-1 entry inhibitors: homology modeling, virtual screening, molecular dynamics simulations and 3D QSAR analysis

Using integrated in-silico computational techniques, including homology modeling, structure-based and pharmacophore-based virtual screening, molecular dynamic simulations, per-residue energy decomposition analysis and atom-based 3D-QSAR analysis, we proposed ten novel compounds as potential CCR5-dependent HIV-1 entry inhibitors. Via validated docking calculations, binding free energies revealed that novel leads demonstrated better binding affinities with CCR5 compared to maraviroc, an FDA-approved HIV-1 entry inhibitor and in clinical use. Per-residue interaction energy decomposition analysis on the averaged MD structure showed that hydrophobic active residues Trp86, Tyr89 and Tyr108 contributed the most to inhibitor binding. The validated 3D-QSAR model showed a high cross-validated r_cv² value of 0.84 using three principal components and non-cross-validated r² value of 0.941. It was also revealed that almost all compounds in the test set and training set yielded a good predicted value. Information gained from this study could shed light on the activity of a new series of lead compounds as potential HIV entry inhibitors and serve as a powerful tool in the drug design and development machinery.

Structural insights from binding poses of CCR2 and CCR5 with clinically important antagonists: a combined in silico study

Chemokine receptors are G protein-coupled receptors that contain seven transmembrane domains. In particular, CCR2 and CCR5 and their ligands have been implicated in the pathophysiology of a number of diseases, including rheumatoid arthritis and multiple sclerosis. Based on their roles in disease, they have been attractive targets for the pharmaceutical industry, and furthermore, targeting both CCR2 and CCR5 can be a useful strategy. Owing to the importance of these receptors, information regarding the binding site is of prime importance. Structural studies have been hampered due to the lack of X-ray crystal structures, and templates with close homologs for comparative modeling. Most of the previous models were based on the bovine rhodopsin and β2-adrenergic receptor. In this study, based on a closer homolog with higher resolution (CXCR4, PDB code: 3ODU 2.5 Å), we constructed three-dimensional models. The main aim of this study was to provide relevant information on binding sites of these receptors. Molecular dynamics simulation was done to refine the homology models and PROCHECK results indicated that the models were reasonable. Here, binding poses were checked with some established inhibitors of high pharmaceutical importance against the modeled receptors. Analysis of interaction modes gave an integrated interpretation with detailed structural information. The binding poses confirmed that the acidic residues Glu291 (CCR2) and Glu283 (CCR5) are important, and we also found some additional residues. Comparisons of binding sites of CCR2/CCR5 were done sequentially and also by docking a potent dual antagonist. Our results can be a starting point for further structure-based drug design.

Theoretical studies on the interactions and interferences of HIV-1 glycoprotein gp120 and its coreceptor CCR5

The interaction between the HIV gp120 protein and coreceptor CCR5 or CXCR4 of the host cell is critical in mediating the HIV entry process. A model for the CCR5-gp120 complex has been developed. In the model, the N-terminus of CCR5 binds to three discontinuous domains of gp120, including the fourth conserved (C4) region, β19/β20 connecting loop, and V3 loop. The second extra-cellular loop (ECL2) of CCR5 also interacts with the crown part of the gp120 V3 loop. The bindings of the three CCR5 antagonists, maraviroc, aplaviroc, and vicriviroc, to the trans-membrane domain of CCR5 have been modeled. The bindings are found to affect the conformation of the ECL2 domain, which in turn drives the N-terminus of CCR5 to an altered state. Aplaviroc is more hydrophilic than maraviroc and vicriviroc, and its binding is more interfered by solvent, resulting in a quite different effect to the structure of CCR5 compared with those of the other two molecules. The above results are in accord with experimental observations and provide a structural basis for further design of CCR5 antagonists.

Efficient synthesis and identification of novel propane-1,3-diamino bridged CCR5 antagonists with variation on the basic center carrier

By employing pharmacophore-based design and the privileged fragments reassembly, a series of piperidine-/tropane-/piperazine-bridged CCR5 antagonists were designed and synthesized via an efficient convergent synthesis strategy, with focus on the optimal choice of the basic center carrier structure. Significantly, the 4-amino-4-methylpiperidine bridged 1-acyl-1,3-propanediamine compounds were identified as a new class of nanomolar CCR5 antagonists, providing an efficient approach and novel scaffolds for further development of potent CCR5 inhibitors.

Therapeutic strategies underpinning the development of novel techniques for the treatment of HIV infection

The HIV replication cycle offers multiple targets for chemotherapeutic intervention, including the viral exterior envelope glycoprotein, gp120; viral co-receptors CXCR4 and CCR5; transmembrane glycoprotein, gp41; integrase; reverse transcriptase; protease and so on. Most currently used anti-HIV drugs are reverse transcriptase inhibitors or protease inhibitors. The expanding application of simulation to drug design combined with experimental techniques have developed a large amount of novel inhibitors that interact specifically with targets besides transcriptase and protease. This review presents details of the anti-HIV inhibitors discovered with computer-aided approaches and provides an overview of the recent five-year achievements in the treatment of HIV infection and the application of computational methods to current drug design.

Ligand and structure-based methodologies for the prediction of the activity of G protein-coupled receptor ligands

Accurate in silico models for the quantitative prediction of the activity of G protein-coupled receptor (GPCR) ligands would greatly facilitate the process of drug discovery and development. Several methodologies have been developed based on the properties of the ligands, the direct study of the receptor-ligand interactions, or a combination of both approaches. Ligand-based three-dimensional quantitative structure-activity relationships (3D-QSAR) techniques, not requiring knowledge of the receptor structure, have been historically the first to be applied to the prediction of the activity of GPCR ligands. They are generally endowed with robustness and good ranking ability; however they are highly dependent on training sets. Structure-based techniques generally do not provide the level of accuracy necessary to yield meaningful rankings when applied to GPCR homology models. However, they are essentially independent from training sets and have a sufficient level of accuracy to allow an effective discrimination between binders and nonbinders, thus qualifying as viable lead discovery tools. The combination of ligand and structure-based methodologies in the form of receptor-based 3D-QSAR and ligand and structure-based consensus models results in robust and accurate quantitative predictions. The contribution of the structure-based component to these combined approaches is expected to become more substantial and effective in the future, as more sophisticated scoring functions are developed and more detailed structural information on GPCRs is gathered.

Discovery of a novel CCR5 antagonist lead compound through fragment assembly

CCR5, as the major co-receptor for HIV-1 entry, is an attractive novel target for the pharmaceutical industry in the HIV-1 therapeutic area. In this study, based on the structures of maraviroc and 1,4-bis(4-(7-chloroquinolin-4-yl)piperazin-1-yl)butane-1,4- dione (1), which was identified using structure-based virtual screening in conjunction with a calcium mobilization assay, a series of novel small molecule CCR5 antagonists have been designed and synthesized through fragment assembly. Preliminary SARs were obtained, which are in good agreement with the molecular binding model and should prove helpful for future antagonist design. The novel scaffold presented here might also be useful in the development of maraviroc-derived second generation CCR5 antagonists.

Synthesis and biological evaluation of 1,3,3,4-tetrasubstituted pyrrolidine CCR5 receptor antagonists. Discovery of a potent and orally bioavailable anti-HIV agent

A series of 1,3,3,4-tetrasubstituted pyrrolidine containing CCR5 receptor antagonists were designed, which were elaborated either by condensation of a lithium salt of 3-(N,N-dibenzyl)aminopropionic acid methyl ester with ethyl benzoformate or by Baylis-Hillman reaction of ethyl acrylate with ethyl benzoformate and subsequent 1,4-addition of benzylamine, in the key steps. These compounds bearing 4-(N,N-disubstituted)amino piperidine units showed low nanomolar potency against the CCR5 receptor, whereas molecules with a 4-phenylpiperidine moiety displayed poor activity. Asymmetric synthesis of the most potent compound 23 a gave rise to the (3R,4S)-enantiomer 30 and the (3S,4R)-enantiomer 31, which showed IC(50) values of 2.9 and 385.9 nM, respectively. These results indicated that (3R,4S)-configuration in the series of compounds is favored for their interaction with the CCR5 receptor. The possible binding mode of these antagonists with the CCR5 receptor was discussed using a computer-modeling method. Compound 30 displayed excellent replication inhibition of seven genetically diverse R5 HIV-1 strains in the PBMC model, in a concentration-dependent manner with EC(50) values ranging from 0.3 nM to 30 nM. This molecule showed oral bioavailabilities of 41.2 % and 21.6 % in rats and dogs, respectively. Thus, compound 30 is a promising candidate for the treatment of HIV-1 infection.

Binding modes of CCR5-targetting HIV entry inhibitors: partial and full antagonists

Since the CC-chemokine receptor 5 (CCR5) was identified as a major co-receptor for human immunodeficiency virus type 1 (HIV-1) entry into a host cell, CCR5-targetting HIV entry inhibitors have been developed and some of them are currently in clinical trials. Most of these inhibitors also inhibit the physiological chemokine reaction function of CCR5, which is so far considered to be safe to patients based on the observation that individuals that naturally lack CCR5 do not show apparent health problems. Nevertheless, to minimize the toxicity and side effects, it would be ideal to preserve the chemokine receptor activity. In this work, we simulated the flexible docking of two small molecule inhibitors to CCR5 in a solvated phospholipid bilayer environment. One of the inhibitors, aplaviroc has a unique feature of preserving two of the natural chemokine ligands binding to CCR5 and subsequent activation whereas the other one, SCH-C fully blocks chemokine-CCR5 interactions. Our results revealed significantly different binding modes of these two inhibitors although both established extensive interaction networks with CCR5. Comparison of the different binding modes suggests that avoiding the deep insertion of inhibitors into the transmembrane helix bundle may be able to preserve chemokine-CCR5 interactions. These results could help design HIV co-receptor activity-specific inhibitors.

3D-QSAR studies of substituted 1-(3, 3-diphenylpropyl)-piperidinyl amides and ureas as CCR5 receptor antagonists

3D-QSAR studies on the derivatives of 1-(3,3-diphenylpropyl)-piperidinyl amide and urea as CCR5 receptor antagonists were performed by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA) methods to rationalize the structural requirements responsible for the inhibitory activity of these compounds. The global minimum energy conformer of the template molecule, the most active and pharmacokinetically stable molecule of the series, was obtained by systematic search and used to build structures of the molecules in the dataset. The best predictions for the CCR5-receptor were obtained with the CoMFA standard model (q (2) = 0.787, r (2) = 0.962) and CoMSIA model combined steric, electrostatic and hydrophobic fields (q (2) = 0.809, r (2) = 0.951). The predictive ability of CoMFA and CoMSIA were determined using a test set of 12 compounds giving predictive correlation coefficients of 0.855 and 0.83, respectively, indicating good predictive power. Further, the robustness of the model was verified by bootstrapping analysis. The contour maps produced by the CoMFA and CoMSIA models were used to identify the structural features relevant to the biological activity in this series. Based on the CoMFA and CoMSIA analysis, we have identified some key features in the series that are responsible for CCR5 antagonistic activity which may be used to design more potent 1-(3,3-diphenylpropyl)-piperidinyl derivatives and predict their activity prior to synthesis.

Understanding the regulation mechanisms of PAF receptor by agonists and antagonists: Molecular modeling and molecular dynamics simulation studies

Platelet-activating factor receptor (PAFR) is a member of G-protein coupled receptor (GPCR) superfamily. Understanding the regulation mechanisms of PAFR by its agonists and antagonists at the atomic level is essential for designing PAFR antagonists as drug candidates for treating PAF-mediated diseases. In this study, a 3D model of PAFR was constructed by a hierarchical approach integrating homology modeling, molecular docking and molecular dynamics (MD) simulations. Based on the 3D model, regulation mechanisms of PAFR by agonists and antagonists were investigated via three 8-ns MD simulations on the systems of apo-PAFR, PAFR-PAF and PAFR-GB. The simulations revealed that binding of PAF to PAFR triggers the straightening process of the kinked helix VI, leading to its activated state. In contrast, binding of GB to PAFR locks PAFR in its inactive state.

I want a new drug: G-protein-coupled receptors in drug development

Huey Lewis and the News summed it up nicely in their 1980s hit record: 'I want a new drug, one that won't make me sick, one that won't make me crash my car, or make me feel three feet thick'. The song could be an anthem for drug discovery in the pharmaceutical industry. We all want new and better drugs with fewer side effects, which are effective for combating the major diseases of our time: cancer, heart disease, obesity and autoimmune diseases. How do we get these new drugs? There are currently some new ideas in drug discovery, centered on that staple diet of the pharmaceutical industry, the G-protein-coupled receptor (GPCR) superfamily. In silico methods, employing receptor-based modeling, offer a more rational approach in the design of drugs targeting GPCRs. These approaches can be used to understand receptor selectivity and species specificity of drugs that interact with GPCRs. In addition, there are various novel approaches, such as the design and potential utility of drugs that target more than one GPCR ('dual specificity' drugs).

Investigation of substituent effect of 1-(3,3-diphenylpropyl)-piperidinyl phenylacetamides on CCR5 binding affinity using QSAR and virtual screening techniques

A linear quantitative-structure activity relationship model is developed in this work using Multiple Linear Regression Analysis as applied to a series of 51 1-(3,3-diphenylpropyl)-piperidinyl phenylacetamides derivatives with CCR5 binding affinity. For the selection of the best variables the Elimination Selection-Stepwise Regression Method (ES-SWR) is utilized. The predictive ability of the model is evaluated against a set of 13 compounds. Based on the produced QSAR model and an analysis on the domain of its applicability, the effects of various structural modifications on biological activity are investigated. The study leads to a number of guanidine derivatives with significantly improved predicted activities.

In silico modelling of the interaction of flavonoids with human P-glycoprotein nucleotide-binding domain

A three-dimensional model of human ABCB1 nucleotide-binding domain (NBD) was developed by homology modelling using the high-resolution human TAP1 transporter structure as template. Interactions between NBD and flavonoids were investigated using in silico docking studies. Ring-A of unmodified flavonoid was located within the NBD P-loop with the 5-hydroxyl group involved in hydrogen bonding with Lys1076. Ring-B was stabilised by hydrophobic stacking interactions with Tyr1044. The 3-hydroxyl group and carbonyl oxygen were extensively involved in hydrogen bonding interactions with amino acids within the NBD. Addition of prenyl, benzyl or geranyl moieties to ring-A (position-6) and hydrocarbon substituents (O-n-butyl to O-n-decyl) to ring-B (position-4) resulted in a size-dependent decrease in predicted docking energy which reflected the increased binding affinities reported in vitro.

QSAR analyses of 3-(4-benzylpiperidin-1-yl)-N-phenylpropylamine derivatives as potent CCR5 antagonists

CCR5 receptor binding affinity of a series of 3-(4-benzylpiperidin-1-yl)propylamine congeners was subjected to QSAR study using the linear free energy related (LFER) model of Hansch. Appropriate indicator variables encoding different group contributions and different physicochemical variables such as hydrophobicity (pi), electronic (Hammett sigma), and steric (molar refractivity, STERIMOL values) parameters of phenyl ring substituents of the compounds were used as predictor variables. The Hansch analysis explores the importance of the lipophilicity and electron-donating substituents for the binding affinity. However, this method could not give more insight into the structure-activity relationships because of the diverse molecular features in the data set. 3D-QSAR analyses of the same data set using Molecular Shape Analysis (MSA), Receptor Surface Analysis (RSA), and Molecular Field Analysis (MFA) techniques were also performed. The best model with acceptable statistical quality was derived from the MSA, which showed the importance of the relative negative charge (RNCG): substituents with a high RNCG value have more binding affinity than the unsubstituted piperidine and phenyl (R1 position) congeners. The relative negative charge surface area (RNCS) is detrimental (e.g. R2 = 3,4-Cl2) for the activity. An increase in the length of the molecule in the Z dimension (Lz) is conducive (e.g. R3 = sulfonylmorpholino), while an increase in the area of the molecular shadow in the XZ plane (Sxz) is detrimental (e.g. R1 = N-c-hexylmethyl-5-oxopyrrolidin-3-yl) for the binding affinity. The presence of a chiral center makes the molecule less active (e.g. R1 = N-methyl-5-oxopyrrolidin-3-yl). An increase in the van der Waals area, the molecular volume, and the difference between the volume of the individual molecule and the shape reference compound are conducive (e.g. R3 = (CH3)2NSO2-) for the binding affinity. Substituents with higher JursFPSA_2 values (fractional charged partial surface area) like the N-methylsulfonylpiperidin-4-yl (R1 position) group have better binding affinity than the substituents such as 4-chlorophenylamino (R1 position). Unsubstituted piperidines (R1 position) with less JursFNSA_1 values have lower binding affinity than the 4-chlorophenyl substituted compounds. The MFA derived equation shows interaction energies at different grid points, while the RSA model shows the importance of hydrophobicity and charge at different regions of the molecules. The models were validated through the leave-one-out, leave-15%-out, and leave-25%-out cross-validation techniques. The developed models were also subjected to a randomization test (99% confidence level). Although the MSA derived models had excellent statistical qualities both for the training as well as test sets, RSA and MFA results for the test sets are not comparable statistically with the MSA derived models.

Quantitative structure-activity relationships for small non-peptide antagonists of CXCR2: indirect 3D approach using the frontal polygon method

The chemokine receptor, CXCR2, plays an important role in recruiting granulocytes to sites of inflammation and has been proposed as an important therapeutic target. A number of CXCR2 antagonists have been synthesized and evaluated; however, quantitative structure-activity relationship (QSAR) models have not been developed for these molecules. Most CXCR2 antagonists can be grouped into four related categories: N,N'-diphenylureas, nicotinamide N-oxides, quinoxalines, and triazolethiols. Based on these categories, we developed a QSAR model for 59 nonpeptide antagonists of CXCR2 using a partial 3D comparison of the antagonists with local fingerprints obtained from rigid and flexible fragments of the molecules. Each compound was represented by calculated structural descriptors that encoded atomic charge, molar refraction, hydrophobicity, and geometric features. We obtained good conventional R(2) coefficients, high leave-one-out cross-validated values for the whole dataset (R(cv)(2)=0.785), as well as for the dataset divided into subsets of triazolethiol derivatives (R(cv)(2)=0.821) and joint subset of N'-diphenylureas, nicotinamide N-oxides, N,N'-diphenylureas, and quinoxaline derivatives and quinoxalines derivatives (R(cv)(2)=0.766), indicating a good predictive ability and robustness of the model. Additionally, charge distribution was found to be a significant contributor in modeling whole dataset. Using our model, structural fragments (submolecules) responsible for the antagonist activity were also identified. These data suggest the QSAR models developed here may be useful in guiding the design of CXCR2 antagonists from molecular fragments.

Molecular docking and 3D-QSAR studies on the binding mechanism of statine-based peptidomimetics with β-secretase

beta-Secretase is an important protease in the pathogenesis of Alzheimer's disease. Some statine-based peptidomimetics show inhibitory activities to the beta-secretase. To explore the inhibitory mechanism, molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies on these analogues were performed. The Lamarckian Genetic Algorithm (LGA) was applied to locate the binding orientations and conformations of the peptidomimetics with the beta-secretase. A good correlation between the calculated binding free energies and the experimental inhibitory activities suggests that the identified binding conformations of these potential inhibitors are reliable. Based on the binding conformations, highly predictive 3D-QSAR models were developed with q(2) values of 0.582 and 0.622 for CoMFA and CoMSIA, respectively. The predictive abilities of these models were validated by some compounds that were not included in the training set. Furthermore, the 3D-QSAR models were mapped back to the binding site of the beta-secretase, to get a better understanding of vital interactions between the statine-based peptidomimetics and the protease. Both the CoMFA and the CoMSIA field distributions are in well agreement with the structural characteristics of the binding groove of the beta-secretase. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction would be useful in developing new drug leads against Alzheimer's disease.