Design of new selective inhibitors of cyclooxygenase-2 by dynamic assembly of molecular building blocks

A method of dynamically assembling molecular building blocks - DycoBlock - has been proposed and tested by Liu et al. This method is based on multiple-copy stochastic dynamics simulation in the presence of a receptor molecule. In this method, a novel algorithm was used to dynamically assemble the molecular building blocks to form candidate compounds. Currently, some new improvements have been incorporated into DycoBlock to make it more efficient. In the new version of DycoBlock, the binding energy and solvent accessible surface area (SASA) can be used to screen the resulting compounds. A simple clustering algorithm based on molecular similarity was developed and used to classify the remaining compounds. The revised DycoBlock was tested by breaking SC-558 - a selective inhibitor of cyclooxygenase-2 (COX-2) - into building blocks and reassembling them in the active site of the enzyme. The accuracy of recovery grew to 58.8% while it was only 16.7% in the previous version. Then, thirty-three kinds of molecular building blocks were used in the design of novel inhibitors and the investigation of diversity. As a result, a total of 1441 compounds was generated with high diversity. After the first screening procedure, there remained 864 reasonable compounds. The results from clustering indicate that the structural motifs in the diarylheterocycle class of COX-2-selective inhibitors have been generated using the revised DycoBlock, and their binding modes were investigated.     

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

Incorporating protein flexibility in structure-based drug discovery: using HIV-1 protease as a test case

We have developed a receptor-based pharmacophore method which utilizes a collection of protein structures to account for inherent protein flexibility in structure-based drug design. Several procedures were systematically evaluated to derive the most general protocol for using multiple protein structures. Most notably, incorporating more protein flexibility improved the performance of the method. The pharmacophore models successfully discriminate known inhibitors from drug-like non-inhibitors. Furthermore, the models correctly identify the bound conformations of some ligands. We used unliganded HIV-1 protease to develop and validate this method. Drug design is always initiated with a protein-ligand structure, and such success with unbound protein structures is remarkable - particularly in the case of HIV-1 protease, which has a large conformational change upon binding. This technique holds the promise of successful computer-based drug design before bound crystal structures are even discovered, which can mean a jump-start of 1-3 years in tackling some medically relevant systems with computational methods.     

Modeling loop backbone flexibility in receptor‐ligand docking simulations

The relevance of receptor conformational change during ligand binding is well documented for many pharmaceutically relevant receptors, but is still not fully accounted for in in silico docking methods. While there has been significant progress in treatment of receptor side chain flexibility sampling of backbone flexibility remains challenging because the conformational space expands dramatically and the scoring function must balance protein–protein and protein–ligand contributions. Here, we investigate an efficient multistage backbone reconstruction algorithm for large loop regions in the receptor and demonstrate that treatment of backbone receptor flexibility significantly improves binding mode prediction starting from apo structures and in cross docking simulations. For three different kinase receptors in which large flexible loops reconstruct upon ligand binding, we demonstrate that treatment of backbone flexibility results in accurate models of the complexes in simulations starting from the apo structure. At the example of the DFG‐motif in the p38 kinase, we also show how loop reconstruction can be used to model allosteric binding. Our approach thus paves the way to treat the complex process of receptor reconstruction upon ligand binding in docking simulations and may help to design new ligands with high specificity by exploitation of allosteric mechanisms. © 2012 Wiley Periodicals, Inc.     

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors

Summary The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100–1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.     

HIV蛋白酶抑制剂——利迪链菌素的分子对接研究

用分子对接的方法, 对利迪链菌素的抗HIV蛋白酶活性进行了研究. 为了更准确地反映利迪链菌素分子与酶蛋白结合的情况, 充分考虑受体活性部位的柔性, 采用了FlexX(初步对接)和Flexidock(精确对接)分两步将配体与受体进行对接. 在初步对接中, 设计了不同的受体活性部位来考察是否有结合水分子参与抑制剂与酶的结合. 对一种作用方式已知的非肽类HIV蛋白酶抑制剂Aha006进行的对接研究显示, 分子模拟的结果与实际情况吻合得较好, 证明了本文所采用的方法的可靠性. 利迪链菌素与蛋白酶活性部位的对接结果显示, 配体分子与受体之间的结合没有结合水分子的参与, 两者通过5对氢键作用结合成为稳定的复合物. 利迪链菌素占据结合腔, 覆盖了蛋白酶的活性三联体Asp25-Thr26-Gly27, 从而起到抑制其生物活性的作用.     

Prediction of binding constants of protein ligands: A fast method for the prioritization of hits obtained from de novo design or 3D database search programs

A dataset of 82 protein–ligand complexes of known 3D structure and binding constant K_i was analysed to elucidate the important factors that determine the strength of protein–ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein–ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein–ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 Å in the length and up to 30 °Cs in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 × 10^-2 M to 1 × 10^-14 M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.     

Design,Synthesis and Structure—Activity Relationships of Phenylalanine-Containing Peptidomimetics as Novel HIV-1 Capsid Binders Based on Ugi Four-Component Reaction

As a key structural protein, HIV capsid (CA) protein plays multiple roles in the HIV life cycle, and is considered a promising target for anti-HIV treatment. Based on the structural information of CA modulator PF-74 bound to HIV-1 CA hexamer, 18 novel phenylalanine derivatives were synthesized via the Ugi four-component reaction. In vitro anti-HIV activity assays showed that most compounds exhibited low-micromolar-inhibitory potency against HIV. Among them, compound I-19 exhibited the best anti-HIV-1 activity (EC₅₀ = 2.53 ± 0.84 μM, CC₅₀ = 107.61 ± 27.43 μM). In addition, I-14 displayed excellent HIV-2 inhibitory activity (EC₅₀ = 2.30 ± 0.11 μM, CC₅₀ > 189.32 μM) with relatively low cytotoxicity, being more potent than that of the approved drug nevirapine (EC₅₀ > 15.02 μM, CC₅₀ > 15.2 μM). Additionally, surface plasmon resonance (SPR) binding assays demonstrated direct binding to the HIV CA protein. Moreover, molecular docking and molecular dynamics simulations provided additional information on the binding mode of I-19 to HIV-1 CA. In summary, we further explored the structure—activity relationships (SARs) and selectivity of anti-HIV-1/HIV-2 of PF-74 derivatives, which is conducive to discovering efficient anti-HIV drugs.     

Drug Design,Synthesis and Biological Evaluation of Heterocyclic Molecules as Anti-Inflammatory Agents

Non-steroidal anti-inflammatory drugs (NSAIDs) are generally utilized for numerous inflammatory ailments. The long-term utilization of NSAIDs prompts adverse reactions such as gastrointestinal ulceration, renal dysfunction and hepatotoxicity; however, selective COX-2 inhibitors prevent these adverse events. Various scientific approaches have been employed to identify safer COX-2 inhibitors, as in any case, a large portion of particular COX-2 inhibitors have been retracted from the market because of severe cardiovascular events. This study aimed to develop and synthesize a novel series of indomethacin analogues with potential anti-inflammatory properties and fewer side effects, wherein carboxylic acid moiety was substituted using DCC/DMAP coupling. This study incorporates the docking of various indomethacin analogues to detect the binding interactions with COX-2 protein (PDB ID: 3NT1). MD simulation was performed to measure the stability and flexibility of ligand–protein interactions at the atomic level, for which the top-scoring ligand–protein complex was selected. These compounds were evaluated in vitro for COX enzymes inhibition. Likewise, selected compounds were screened in vivo for anti-inflammatory potential using the carrageenan-induced rat paw oedema method and their ulcerogenic potential. The acute toxicity of compounds was also predicted using in silico tools. Most of the compounds exhibited the potent inhibition of both COX enzymes; however, 3e and 3c showed the most potent COX-2 inhibition having IC50 0.34 µM and 1.39 µM, respectively. These compounds also demonstrated potent anti-inflammatory potential without ulcerogenic liability. The biological evaluation revealed that the compound substituted with 4-nitrophenyl was most active.     

Cluster expansion models for flexible‐backbone protein energetics

Protein structure prediction and design often involve discrete modeling of side‐chain conformations on structural templates. Introducing backbone flexibility into such models has proven important in many different applications. Backbone flexibility improves model accuracy and provides access to larger sequence spaces in computational design, although at a cost in complexity and time. Here, we show that the influence of backbone flexibility on protein conformational energetics can be treated implicitly, at the level of sequence, using the technique of cluster expansion. Cluster expansion provides a way to convert structure‐based energies into functions of sequence alone. It leads to dramatic speed‐ups in energy evaluation and provides a convenient functional form for the analysis and optimization of sequence‐structure relationships. We show that it can be applied effectively to flexible‐backbone structural models using four proteins: α‐helical coiled‐coil dimers and trimers, zinc fingers, and Bcl‐x_L/peptide complexes. For each of these, low errors for the sequence‐based models when compared with structure‐based evaluations show that this new way of treating backbone flexibility has considerable promise, particularly for protein design. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Enhancing the accuracy of virtual screening: molecular dynamics with quantum-refined force fields

Summary A methodology aimed at improving the accuracy of current docking–scoring procedures is proposed, and validated through detailed tests of its performance in predicting the activity of HIV-1 protease inhibitors. This methodology is based on molecular dynamics simulations using a force field whose effective charges are refined by means of a novel procedure that relies on quantum-mechanical calculations and preserves the internal consistency of the parameterization scheme.     

Tuning of hydrogen bond strength using substituents on phenol and aniline: A possible ligand design strategy

Using Density Functional Theory, the hydrogen bonding energy is calculated for the interaction of phenol and aniline with four model compounds representing the protein backbone and various amino acid site chain residues. The models are methanol, protonated methylamine, formaldehyde and acetate anion. The H-bond energies for the uncharged species are 2.5kcalmol^–1, whereas the charged model compounds bind with much higher energies of 20kcalmol^–1. The effect of para-substitution on the hydrogen bond energies is determined. Substitution has little effect on the H-bond energy of the neutral complexes (<2kcalmol^–1), but for the positively and negatively charged systems substitution drastically alters the binding energies, e.g., 14.3kcalmol^–1 for para-NO₂. In the context of protein–ligand binding, relatively small changes in binding energy can cause large changes in affinity due to their exponential relationship. This means that for –NO₂ an enormous change of 10 orders of magnitude for the affinity constant is predicted. These calculations allow prediction of H-bonds, using different substituents, in order to fine-tune and optimize ligand–protein interactions in the search for drug candidates.     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.     

The Glu143 Residue Might Play a Significant Role in T20 Peptide Binding to HIV-1 Receptor gp41: An In Silico Study

Despite the enormous efforts made to develop other fusion inhibitors for HIV, the enfuvirtide (known as T20) peptide is the only approved HIV-1 inhibitory drug so far. Investigating the role of potential residues of the T20 peptide’s conformational dynamics could help us to understand the role of potential residues of the T20 peptide. We investigated T20 peptide conformation and binding interactions with the HIV-1 receptor (i.e., gp41) using MD simulations and docking techniques, respectively. Although the mutation of E143 into alanine decreased the flexibility of the E143A mutant, the conformational compactness of the mutant was increased. This suggests a potential role of E143 in the T20 peptide’s conformation. Interestingly, the free energy landscape showed a significant change in the wild-type T20 minimum, as the E143A mutant produced two observed minima. Finally, the docking results of T20 to the gp41 receptor showed a different binding interaction in comparison to the E143A mutant. This suggests that E143 residue can influence the binding interaction with the gp41 receptor. Overall, the E143 residue showed a significant role in conformation and binding to the HIV-1 receptor. These findings can be helpful in optimizing and developing HIV-1 inhibitor peptides.     

Drug-resistant molecular mechanism of CRF01_AE HIV-1 protease due to V82F mutation

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery. The circulating recombinant form 01 A/E (CRF01_AE, abbreviated AE) subtype is one of the most common HIV-1 subtypes, which is infecting more humans and is expanding rapidly throughout the world. It is, therefore, necessary to develop inhibitors against subtype AE HIV-1 PR. In this work, we have performed computer simulation of subtype AE HIV-1 PR with the drugs lopinavir (LPV) and nelfinavir (NFV), and examined the mechanism of resistance of the V82F mutation of this protease against LPV both structurally and energetically. The V82F mutation at the active site results in a conformational change of 79′s loop region and displacement of LPV from its proper binding site, and these changes lead to rotation of the side-chains of residues D25 and I50′. Consequently, the conformation of the binding cavity is deformed asymmetrically and some interactions between PR and LPV are destroyed. Additionally, by comparing the interactive mechanisms of LPV and NFV with HIV-1 PR we discovered that the presence of a dodecahydroisoquinoline ring at the P1′ subsite, a [2-(2,6-dimethylphenoxy)acetyl]amino group at the P2′ subsite, and an N2 atom at the P2 subsite could improve the binding affinity of the drug with AE HIV-1 PR. These findings are helpful for promising drug design. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

EGFR and COX-2 Dual Inhibitor: The Design,Synthesis, and Biological Evaluation of Novel Chalcones

For most researchers, discovering new anticancer drugs to avoid the adverse effects of current ones, to improve therapeutic benefits and to reduce resistance is essential. Because the COX-2 enzyme plays an important role in various types of cancer leading to malignancy enhancement, inhibition of apoptosis, and tumor-cell metastasis, an indispensable objective is to design new scaffolds or drugs that possess combined action or dual effect, such as kinase and COX-2 inhibition. The start compounds A1 to A6 were prepared through the diazo coupling of 3-aminoacetophenone with a corresponding phenol and then condensed with two new chalcone series, C7–18. The newly synthesized compounds were assessed against both COX-2 and epidermal growth factor receptor (EGFR) for their inhibitory effect. All novel compounds were screened for cytotoxicity against five cancer cell lines. Compounds C9 and G10 exhibited potent EGFR inhibition with IC₅₀ values of 0.8 and 1.1 µM, respectively. Additionally, they also displayed great COX-2 inhibition with IC₅₀ values of 1.27 and 1.88 µM, respectively. Furthermore, the target compounds were assessed for their cytotoxicity against pancreatic ductal cancer (Panc-1), lung cancer (H-460), human colon cancer (HT-29), human malignant melanoma (A375) and pancreatic cancer (PaCa-2) cell lines. Interestingly, compounds C10 and G12 exhibited the strongest cytotoxic effect against PaCa-2 with average IC₅₀ values of 0.9 and 0.8 µM, respectively. To understand the possible binding modes of the compounds under investigation with the receptor cites of EGFR and COX-2, a virtual docking study was conducted.     

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein,Mpro and PLpro of SARS-CoV-2: Protein–Protein Interactions,Dynamics Simulations and Free Energy Calculations

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (−16.8 ± 0.02 kcal/mol, −12.3 ± 0.03 kcal/mol and −13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (M^pro) and the papainlike protease (PL^pro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.     

SPROUT: A program for structure generation

Summary SPROUT is a new computer program for constrained structure generation that is designed to generate molecules for a range of applications in molecular recognition. It uses artificial intelligence techniques to moderate the combinatorial explosion that is inherent in structure generation. The program is presented here for the design of enzyme inhibitors. Structure generation is divided into two phases: (i) primary structure generation to produce molecular graphs to fit the steric constraints; and (ii) secondary structure generation which is the process of introducing appropriate functionality to the graphs to produce molecules that satisfy the secondary constraints, e.g., electrostatics and hydrophobicity. Primary structure generation has been tested on two enzyme receptor sites; the p-amidino-phenyl-pyruvate binding site of trypsin and the acetyl pepstatin binding site of HIV-1 protease. The program successfully generates structures that resemble known substrates and, more importantly, the predictive power of the program has been demonstrated by its ability to suggest novel structures.     

A hybrid method of molecular dynamics and harmonic dynamics for docking of flexible ligand to flexible receptor

We have developed a new docking method to consider receptor flexibility, a hybrid method of molecular dynamics and harmonic dynamics. The global motions of the whole receptor were approximately introduced into those of the receptor in the docking simulation as harmonic dynamics. On the other hand, the local flexibility of the side chains was also considered by conventional molecular dynamics. We confirmed that this new method can reproduce the fluctuations of the whole receptor by making a comparison of the directions and amplitudes of the global fluctuations. Then this method was applied to the docking of HIV-1 protease and its ligand. As a result, we observed a docking process where the ligand enters into the binding pocket well, which implies that this method is effective enough to reproduce a molecular complex formation.     

Antiinflammation Derived Suzuki-Coupled Fenbufens as COX-2 Inhibitors: Minilibrary Construction and Bioassay

A small fenbufen library comprising 18 compounds was prepared via Suzuki Miyara coupling. The five-step preparations deliver 9–17% biphenyl compounds in total yield. These fenbufen analogs exert insignificant activity against the IL-1 release as well as inhibiting cyclooxygenase 2 considerably. Both the para-amino and para-hydroxy mono substituents display the most substantial COX-2 inhibition, particularly the latter one showing a comparable activity as celecoxib. The most COX-2 selective and bioactive disubstituted compound encompasses one electron-withdrawing methyl and one electron-donating fluoro groups in one arene. COX-2 is selective but not COX-2 to bioactive compounds that contain both two electron-withdrawing groups; disubstituted analogs with both resonance-formable electron-donating dihydroxy groups display high COX-2 activity but inferior COX-2 selectivity. In silico simulation and modeling for three COX-2 active—p-fluoro, p-hydroxy and p-amino—fenbufens show a preferable docking to COX-2 than COX-1. The most stabilization by the p-hydroxy fenbufen with COX-2 predicted by theoretical simulation is consistent with its prominent COX-2 inhibition resulting from experiments.