Design of peptidomimetic inhibitors of aspartic protease of HIV-1 containing –PheΨPro– core and displaying favourable ADME-related properties

In-silico study: docking simulation and molecular dynamics of peptidomimetic fullerene-based derivatives against SARS-CoV-2 Mpro

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2, has become a global pandemic resulting in significant morbidity and mortality. This study presents 12 new peptidomimetic fullerene-based derivatives in three groups that are investigated theoretically as SARS-CoV-2 Mpro inhibitors to increase the chance of treating COVID-19. Studied compounds are designed and optimized at B88-LYP/DZVP method. Molecular descriptors results show the stability and reactivity of the compounds with Mpro, especially in the 3rd group (Ser compounds). However, Lipinski's Rule of Five values indicates that the compounds are not suitable as oral drugs. Furthermore, molecular docking simulations are carried out to investigate the binding affinity and interaction modes of the top five compounds (compounds 1, 9, 11, 2, and 10) with the Mpro protein, which have the lowest binding energy. Molecular dynamics simulations are also performed to evaluate the stability of the protein-ligand complexes with compounds 1 and 9 and compare them with natural substrate interaction. The analysis of RMSD, H-bonds, Rg, and SASA indicates that both compounds 1 (Gly-α acid) and 9 (Ser-α acid) have good stability and strong binding affinity with the Mpro protein. However, compound 9 shows slightly better stability and binding affinity compared to compound 1.

Structure-Based Design and Pharmacophore-Based Virtual Screening of Combinatorial Library of Triclosan Analogues Active against Enoyl-Acyl Carrier Protein Reductase of Plasmodium falciparum with Favourable ADME Profiles

Cost-effective therapy of neglected and tropical diseases such as malaria requires everlasting drug discovery efforts due to the rapidly emerging drug resistance of the plasmodium parasite. We have carried out computational design of new inhibitors of the enoyl-acyl carrier protein reductase (ENR) of Plasmodium falciparum (PfENR) using computer-aided combinatorial and pharmacophore-based molecular design. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) complexation QSAR model was developed for triclosan-based inhibitors (TCL) and a significant correlation was established between the calculated relative Gibbs free energies of complex formation (∆∆Gcom) between PfENR and TCL and the observed inhibitory potencies of the enzyme (IC50exp) for a training set of 20 known TCL analogues. Validation of the predictive power of the MM-PBSA QSAR model was carried out with the generation of 3D QSAR pharmacophore (PH4). We obtained a reasonable correlation between the relative Gibbs free energy of complex formation ∆∆Gcom and IC50exp values, which explained approximately 95% of the PfENR inhibition data: pIC50exp=-0.0544×∆∆Gcom+6.9336,R2=0.95. A similar agreement was established for the PH4 pharmacophore model of the PfENR inhibition (pIC50exp=0.9754×pIC50pre+0.1596, R2=0.98). Analysis of enzyme-inhibitor binding site interactions suggested suitable building blocks to be used in a virtual combinatorial library of 33,480 TCL analogues. Structural information derived from the complexation model and the PH4 pharmacophore guided us through in silico screening of the virtual combinatorial library of TCL analogues to finally identify potential new TCL inhibitors effective at low nanomolar concentrations. Virtual screening of the library by PfENR-PH4 led to a predicted IC50pre value for the best inhibitor candidate as low as 1.9 nM. Finally, the stability of PfENR-TCLx complexes and the flexibility of the active conformation of the inhibitor for selected top-ranking TCL analogues were checked with the help of molecular dynamics. This computational study resulted in a set of proposed new potent inhibitors with predicted antimalarial effects and favourable pharmacokinetic profiles that act on a novel pharmacological target, PfENR.

Computer-Aided Design of Peptidomimetic Inhibitors of Falcipain-3: QSAR and Pharmacophore Models

In this work, antiparasitic peptidomimetics inhibitors (PEP) of falcipain-3 (FP3) of Plasmodium falciparum (Pf) are proposed using structure-based and computer-aided molecular design. Beginning with the crystal structure of PfFP3-K11017 complex (PDB ID: 3BWK), three-dimensional (3D) models of FP3-PEPx complexes with known activities ( IC50exp) were prepared by in situ modification, based on molecular mechanics and implicit solvation to compute Gibbs free energies (GFE) of inhibitor-FP3 complex formation. This resulted in a quantitative structure–activity relationships (QSAR) model based on a linear correlation between computed GFE (ΔΔGcom) and the experimentally measured  IC50exp. Apart from the structure-based relationship, a ligand-based quantitative pharmacophore model (PH4) of novel PEP analogues where substitutions were directed by comparative analysis of the active site interactions was derived using the proposed bound conformations of the PEPx. This provided structural information useful for the design of virtual combinatorial libraries (VL), which was virtually screened based on the 3D-QSAR PH4. The end results were predictive inhibitory activities falling within the low nanomolar concentration range.

Adhesion of fibroblast cells on thin films representing surfaces of polymeric scaffolds of human urethra rationalized by molecular models of integrin binding: cell adhesion on polymeric scaffolds for regenerative medicine

This work combines experimental and computational study of Balb/3T3 clone A31 mouse embryo fibroblasts cell line adhesion and proliferation on fourteen different polymeric surfaces prepared from poly(dioxanone) (PDO), poly(glycolic acid) (PGA), poly(hydroxybutyrate) (PHB), and poly(L-lactic acid) (PLA), and their 1:1 mixtures. The study was done with the aim to explore the attractive interactions between various synthetic biomaterials and simple model of the cell attachment mechanism involving the trans-membrane protein integrin. The considered polymeric biodegradable biomaterials can be used as scaffolds for tissue engineering and regenerative urology. During the growth of new tissue, the polymer scaffold is replaced by the extracellular matrix (ECM) synthetized by the proliferating cells. The adhesion and proliferation experiments were done on thin polymer films produced by solvent casting. The computational approach used 3D molecular models of two layers of ordered parallel polymeric fibres, which formed quasi-planar nanosized models of the scaffold surface. Experimental data showed that PGA based polymer films promote the cell adhesion. Cell proliferation testing, performed by incubating the fibroblast cells with the studied polymer films, disclosed that PLA, PHB/PLA and PHB/PGA systems are able to support proliferation of Balb/3T3 clone A31 cells equal to the plain glass. Relative interaction energies between 3D models of polymeric films and the α₂ I domain of the cell adhesion receptor integrin α₂β₁ computed by molecular mechanics suggest that plain polymers PGA, PDO and mixtures PDO/PGA, PHB/PGA, and especially PGA/PLA display elevated affinity to the cell-attachment protein, which confirms the experimental observations. The combination of experimental and modelling approach can assist rational design of synthetic polymeric biomaterial for scaffolds of artificial human urethra that can be efficiently colonized by cells.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

Antiviral agents against COVID-19: structure-based design of specific peptidomimetic inhibitors of SARS-CoV-2 main protease

Despite the intense development of vaccines and antiviral therapeutics, no specific treatment of coronavirus disease 2019 (COVID-19), caused by the new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently available. Recently, X-ray crystallographic structures of a validated pharmacological target of SARS-CoV-2, the main protease (M^pro also called 3CL^pro) in complex with peptide-like irreversible inhibitors have been published. We have carried out computer-aided structure-based design and optimization of peptidomimetic irreversible α-ketoamide M^pro inhibitors and their analogues using MM, MD and QM/MM methodology, with the goal to propose lead compounds with improved binding affinity to SARS-CoV-2 M^pro, enhanced specificity for pathogenic coronaviruses, decreased peptidic character, and favourable drug-like properties. The best inhibitor candidates designed in this work show largely improved interaction energies towards the M^pro and enhanced specificity due to 6 additional hydrogen bonds to the active site residues. The presented results on new SARS-CoV-2 M^pro inhibitors are expected to stimulate further research towards the development of specific anti-COVID-19 drugs.

Structure-based design of SARS-CoV-2 main protease inhibitors identified hydantoin, benzothiazine and cresol moieties as promising residues of new peptidomimetic inhibitors.

Pathogenicity of new BEST1 variants identified in Italian patients with best vitelliform macular dystrophy assessed by computational structural biology

Background

Best vitelliform macular dystrophy (BVMD) is an autosomal dominant macular degeneration. The typical central yellowish yolk-like lesion usually appears in childhood and gradually worsens. Most cases are caused by variants in the BEST1 gene which encodes bestrophin-1, an integral membrane protein found primarily in the retinal pigment epithelium.

Methods

Here we describe the spectrum of BEST1 variants identified in a cohort of 57 Italian patients analyzed by Sanger sequencing. In 13 cases, the study also included segregation analysis in affected and unaffected relatives. We used molecular mechanics to calculate two quantitative parameters related to calcium-activated chloride channel (CaCC composed of 5 BEST1 subunits) stability and calcium-dependent activation and related them to the potential pathogenicity of individual missense variants detected in the probands.

Results

Thirty-six out of 57 probands (63% positivity) and 16 out of 18 relatives proved positive to genetic testing. Family study confirmed the variable penetrance and expressivity of the disease. Six of the 27 genetic variants discovered were novel: p.(Val9Gly), p.(Ser108Arg), p.(Asn179Asp), p.(Trp182Arg), p.(Glu292Gln) and p.(Asn296Lys). All BEST1 variants were assessed in silico for potential pathogenicity. Our computational structural biology approach based on 3D model structure of the CaCC showed that individual amino acid replacements may affect channel shape, stability, activation, gating, selectivity and throughput, and possibly also other features, depending on where the individual mutated amino acid residues are located in the tertiary structure of BEST1. Statistically significant correlations between mean logMAR best-corrected visual acuity (BCVA), age and modulus of computed BEST1 dimerization energies, which reflect variations in the in CaCC stability due to amino acid changes, permitted us to assess the pathogenicity of individual BEST1 variants.

Conclusions

Using this computational approach, we designed a method for estimating BCVA progression in patients with BEST1 variants.

Structure-Based Design and in Silico Screening of Virtual Combinatorial Library of Benzamides Inhibiting 2-trans Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Predicted Pharmacokinetic Profiles

Background: During the previous decade a new class of benzamide-based inhibitors of 2-trans enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (Mt) with unusual binding mode have emerged. Here we report in silico design and evaluation of novel benzamide InhA-Mt inhibitors with favorable predicted pharmacokinetic profiles. Methods: By using in situ modifications of the crystal structure of N-benzyl-4-((heteroaryl)methyl) benzamide (BHMB)-InhA complex (PDB entry 4QXM), 3D models of InhA-BHMBx complexes were prepared for a training set of 19 BHMBs with experimentally determined inhibitory potencies (half-maximal inhibitory concentrations IC50exp). In the search for active conformation of the BHMB1-19, linear QSAR model was prepared, which correlated computed gas phase enthalpies of formation (∆∆HMM) of InhA-BHMBx complexes with the IC50exp. Further, taking into account the solvent effect and entropy changes upon ligand, binding resulted in a superior QSAR model correlating computed complexation Gibbs free energies (∆∆Gcom). The successive pharmacophore model (PH4) generated from the active conformations of BHMBs served as a virtual screening tool of novel analogs included in a virtual combinatorial library (VCL) of compounds containing benzamide scaffolds. The VCL filtered by Lipinski’s rule-of-five was screened by the PH4 model to identify new BHMB analogs. Results: Gas phase QSAR model: −log₁₀(IC₅₀^exp) = pIC₅₀^exp = −0.2465 × ∆∆H_MM + 7.95503, R² = 0.94; superior aqueous phase QSAR model: pIC₅₀^exp = −0.2370 × ∆∆G_com + 7.8783, R² = 0.97 and PH4 pharmacophore model: pIC50exp = 1.0013 × pIC50exp − 0.0085, R² = 0.95. The VCL of more than 114 thousand BHMBs was filtered down to 73,565 analogs Lipinski’s rule. The five-point PH4 screening retained 90 new and potent BHMBs with predicted inhibitory potencies IC₅₀^pre up to 65 times lower than that of BHMB1 (IC₅₀^exp = 20 nM). Predicted pharmacokinetic profile of the new analogs showed enhanced cell membrane permeability and high human oral absorption compared to current anti-tuberculotics. Conclusions: Combined use of QSAR models that considered binding of the BHMBs to InhA, pharmacophore model, and ADME properties helped to recognize bound active conformation of the benzamide inhibitors, permitted in silico screening of VCL of compounds sharing benzamide scaffold and identification of new analogs with predicted high inhibitory potencies and favorable pharmacokinetic profiles.

Epithelial cell adhesion on films mimicking surface of polymeric scaffolds of artificial urethra compared to molecular modeling of integrin binding

In this study, a combined experimental and computational study of long-term human bladder epithelial cell line HBLAK adhesion and proliferation on five different polymeric surfaces, namely hyaluronic acid, amylose, collagen, polyhydroxybutyrate, and polylactic acid, was performed with the goal to understand the nature of the attraction between various surface materials and a simplified model of the cell surface (transmembrane protein integrin). These biodegradable polymers are frequently used as scaffolds for tissue engineering. During formation of the new tissue, the scaffold polymers are gradually replaced by the natural extracellular matrix of the proliferating cells. Cell adhesion and proliferation experiments were carried out employing thin polymer films prepared by solvent casting while for the computational approach three-dimensional molecular models of layers of ordered polymeric fibers were used as quasi-planar nano-sized models of polymeric surface patches. Experimental results indicated a good capability of amylose, polyhydroxybutyrate, and hyaluronic acid polymer films to foster cell adhesion. Proliferation experiment, carried out by incubating cells with the investigated polymer films for 72 h, showed that all the investigated polymers are able to sustain a good proliferation of HBLAK cells almost comparable to plain glass. Computational estimate of molecular mechanic interaction energies between three-dimensional models of polymeric films and the collagen-binding α2 I domain of the cell adhesion receptor integrin α2β1 confirmed elevated affinity to amylose and polyhydroxybutyrate that is related to higher polarity of function groups on the film surface as documented by the maps of molecular electrostatic potential. This combined experimental and modeling approach can contribute to rational design and surface modifications of polymeric material suitable for forming the scaffolds of human urethra that can be effectively colonized by stem cells.

Virtual design of novel Plasmodium falciparum cysteine protease falcipain-2 hybrid lactone-chalcone and isatin-chalcone inhibitors probing the S2 active site pocket

We report computer-aided design of new lactone–chalcone and isatin–chalcone (HLCIC) inhibitors of the falcipain-2 (PfFP-2). 3D models of 15 FP-2:HLCIC1-15 complexes with known observed activity (IC₅₀^exp) were prepared to establish a quantitative structure–activity (QSAR) model and linear correlation between relative Gibbs free energy of enzyme:inhibitor complex formation (ΔΔG_com) and IC₅₀^exp: pIC₅₀^exp = −0.0236 × ΔΔG_com+5.082(#); R² = 0.93. A 3D pharmacophore model (PH4) derived from the QSAR directed our effort to design novel HLCIC analogues. During the design, an initial virtual library of 2621440 HLCIC was focused down to 18288 drug-like compounds and finally, PH4 screened to identify 81 promising compounds. Thirty-three others were added from an intuitive substitution approach intended to fill better the enzyme S2 pocket. One hundred and fourteen theoretical IC₅₀ (IC₅₀^pre) values were predicted by means of (#) and their pharmacokinetics (ADME) profiles. More than 30 putative HLCICs display IC₅₀^pre 100 times superior to that of the published most active training set inhibitor HLCIC1.

Design, Synthesis, and Neuroprotective Effects of a Series of Pyrazolines against 6-Hydroxydopamine-Induced Oxidative Stress

Parkinson's disease (PD) is a chronic, progressive, and age-related neurodegenerative disorder characterized by the loss of midbrain dopaminergic neurons caused by the accumulation of free radicals and oxidative stress. Based on the neuroprotective properties of 2-pyrazoline derivatives, in the current work, 1-(phenyl/4-substituted phenyl)-3-(2-furanyl/thienyl)-5-aryl-2-pyrazolines (3a⁻i, 4a⁻i) were synthesized via the cyclization of the chalcones (1, 2) with suitable phenylhydrazine hydrochloride derivatives. All these compounds were investigated for their neuroprotective effects using an in vitro 6-hydroxydopamine (6-OHDA)-induced neurotoxicity model of PD in the rat pheochromocytoma (PC-12) Adh cell line. In addition, some different pharmacokinetic parameters of all compounds were in silico predicted by the QikProp module of Schrödinger's Maestro molecular modeling package. 4-Methylsulfonylphenyl substituted compounds 3h (20%) and 4h (23%) were determined as the most promising neuroprotective agents related to their inductive roles in cell viability when compared with the 6-OHDA-positive control group (43% and 42%, respectively). Moreover, in silico pharmacokinetic results indicated that all compounds were within the acceptable range intended for human use. According to both in vitro and in silico studies, compounds 3h and 4h draw attention as potential orally bioavailable therapeutic drug candidates against neurodegeneration in PD.

Diarylcyclopropane hydroxamic acid inhibitors of histone deacetylase 4 designed by combinatorial approach and QM/MM calculations

Inhibitors of histone deacetylase superfamily (HDAC), which induce cell cycle arrest, trigger cell death and reduce angiogenesis appear as promising anti-cancer drugs targeting the epigenetic regulation of gene expression. Approved HDAC inhibitors were found effective against haematological and solid malignancies, other HDACIs are currently in clinical trials for the treatment of neurological diseases or immune disorders. Among those, diarylcyclopropane hydroxamic acids (DCHA) were found to be potent and selective inhibitors of the class IIa HDACs, specifically HDAC4, a pharmacological target for the treatment of Huntington's disease and muscular atrophy. Crystallographic analysis revealed that one of the aryl groups of the DCHA fills the lower specificity pocket of the HDAC4 catalytic site that is specific for the class IIa HDACs. We have used computer-assisted combinatorial chemistry, hybrid quantum mechanics/molecular mechanics (QM/MM) with implicit solvation and QSAR models to optimize DCHA inhibitors and propose more potent DCHA analogues. The QM/MM approach has been selected since the process of inhibitor binding to the catalytic zinc and polar amino acid residues of the deacetylase active site induces considerable rearrangement of electron density of the inhibitor. Virtual combinatorial library consisting of 12180 DCHA analogues was focused by means of structure-based evaluation to form a small combinatorial subset enriched in potentially interesting inhibitor candidates. Two validated QSAR models making use of computed relative binding affinities of the DCHA inhibitors to the HDAC4 (ΔΔG_com^QM/MM) were utilized to estimate the inhibitory potencies of the new analogues. The predicted half-maximal inhibitory concentrations (IC₅₀^pre) of the designed analogues fall into the low nanomolar concentration range and their predicted ADME properties are also favourable. The best designed DCHA analogues contain indazole, phenylpiperidine, phenyloxazole or hydroxypyridine moieties and stabilize bound inhibitors by hydrogen bonds to the catalytic water molecule and backbone carbonyl groups of the deacetylase active site residues. This makes them more potent and more specific inhibitors towards the HDAC4 isoform than the known diarylcyclopropane hydroxamic acids. The analogues are recommended for synthesis and experimental verification of inhibitory potencies in medicinal chemistry laboratories.

In silico Molecular Docking and ADME Studies of 1,3,4-Thiadiazole Derivatives in Relation to in vitro PON1 Activity

BACKGROUND

Paraoxonase 1 (PON1) is a paraoxonase, arylesterase and lactonase associated with protection of lipoproteins and cell membranes against oxidative modification.

OBJECTIVE

Based on antioxidative properties of PON1 and significance of 1,3,4-thiadiazoles in pharmaceutical chemistry, herein we aimed to evaluate the potentials of 1,3,4-thiadiazole derivatives as PON1 activators.

METHODS

2-[[5-(2,4-Difluoro/dichlorophenylamino)-1,3,4-thiadiazol-2-yl]thio]acetophenone derivatives (1-18) were in vitro evaluated for their activator effects on PON1 which was purified using ammonium sulfate precipitation (60-80%) and DEAE-Sephadex anion exchange chromatography. Molecular docking studies were performed for the detection of affinities of all compounds to the active site of PON1. Moreover, Absorption, Distribution, Metabolism and Excretion (ADME) properties of all compounds were also in silico predicted. In silico molecular docking and ADME studies were carried out according to modules of Schrodinger's Maestro molecular modeling package.

RESULTS

All compounds, particularly compounds 10, 13 and 17, were determined as promising PON1 activators and apart from compound 1, all of them were detected in the active site of PON1. Besides, ADME results indicated that all compounds were potential orally bioavailable drug-like molecules.

CONCLUSION

PON1 activators, compounds 10, 13 and 17 stand out as potential drug candidates for further antioxidant studies and these compounds can be investigated for their therapeutic effects in many disorders such as atherosclerosis, diabetes mellitus, obesity, chronic liver inflammation and many more.

Stereoisomers of oseltamivir – synthesis, in silico prediction and biological evaluation

Three diastereoisomers of oseltamivir were synthesized, their properties predicted by quantum-chemical calculations and their antiviral activities evaluated.

Computer-Aided Design of Orally Bioavailable Pyrrolidine Carboxamide Inhibitors of Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Pharmacokinetic Profiles

We have carried out a computational structure-based design of new potent pyrrolidine carboxamide (PCAMs) inhibitors of enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (MTb). Three-dimensional (3D) models of InhA-PCAMx complexes were prepared by in situ modification of the crystal structure of InhA-PCAM1 (Protein Data Bank (PDB) entry code: 4U0J), the reference compound of a training set of 20 PCAMs with known experimental inhibitory potencies (IC50(exp)). First, we built a gas phase quantitative structure-activity relationships (QSAR) model, linearly correlating the computed enthalpy of the InhA-PCAM complex formation and the IC50(exp). Further, taking into account the solvent effect and loss of inhibitor entropy upon enzyme binding led to a QSAR model with a superior linear correlation between computed Gibbs free energies (ΔΔGcom) of InhA-PCAM complex formation and IC50(exp) (pIC50(exp) = -0.1552·ΔΔGcom + 5.0448, R² = 0.94), which was further validated with a 3D-QSAR pharmacophore model generation (PH4). Structural information from the models guided us in designing of a virtual combinatorial library (VL) of more than 17 million PCAMs. The VL was adsorption, distribution, metabolism and excretion (ADME) focused and reduced down to 1.6 million drug like orally bioavailable analogues and PH4 in silico screened to identify new potent PCAMs with predicted IC50(pre) reaching up to 5 nM. Combining molecular modeling and PH4 in silico screening of the VL resulted in the proposed novel potent antituberculotic agent candidates with favorable pharmacokinetic profiles.

Advances in computationally modeling human oral bioavailability

Although significant progress has been made in experimental high throughput screening (HTS) of ADME (absorption, distribution, metabolism, excretion) and pharmacokinetic properties, the ADME and Toxicity (ADME-Tox) in silico modeling is still indispensable in drug discovery as it can guide us to wisely select drug candidates prior to expensive ADME screenings and clinical trials. Compared to other ADME-Tox properties, human oral bioavailability (HOBA) is particularly important but extremely difficult to predict. In this paper, the advances in human oral bioavailability modeling will be reviewed. Moreover, our deep insight on how to construct more accurate and reliable HOBA QSAR and classification models will also discussed.

Virtually Designed Triclosan-Based Inhibitors of Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis and of Plasmodium falciparum

We report here new chemical structures of predicted nanomolar triclosan-based inhibitors (TCLs) of Mycobacterium tuberculosis enoyl-acyl carrier protein reductase (InhA) virtually proposed by computer-assisted molecular design. 3D models of InhA-TCL complexes were prepared by in situ modifications of the reference crystal structure (PDB entry 1P45) for a training set of 15 TCLs with known InhA inhibitory activities. A QSAR model was built leading to linear correlation between the calculated free energies of complexation (ΔΔGcom ) and experimental values IC50 (exp) : pIC50 =-0.0657×ΔΔGcom +3.0502, R(2) =0.96. In addition, ligand-based quantitative pharmacophore model (PH4) was built from bound conformations of the training set compounds and confirmed the correlation between molecular models and observed activities: pIC50 (exp=) 0.8929×pIC50 (pre) -0.441, R(2) =0.95. Structural information from both models helped us to propose new TCL analogues. A virtual library of TCLs with known predicted activities against enoyl-acyl carrier protein reductase of Plasmodium falciparum (PfENR) was evaluated, revealing dual target TCLs. Moreover, analysis of binding site interactions suggested enriching substitutions, which led to more potent TCLs with predicted pIC50 (pre) as low as 7 nM. The computational approach, which used both free energy estimated from molecular modeling and 3D-QSAR pharmacophore model, was helpful in virtually proposing the dual-targeted drugs and provided valuable information for the design of novel potential antituberculotic agents.

Quantitative structure–activity relationships and design of thymine-like inhibitors of thymidine monophosphate kinase of Mycobacterium tuberculosis with favourable pharmacokinetic profiles

We have designed new potent inhibitors of thymidine monophosphate kinase of Mycobacterium tuberculosis (TMPK_mt) using structure-based molecular design.

Design of Thymidine Analogues Targeting Thymidilate Kinase of Mycobacterium tuberculosis

We design here new nanomolar antituberculotics, inhibitors of Mycobacterium tuberculosis thymidine monophosphate kinase (TMPKmt), by means of structure-based molecular design. 3D models of TMPKmt-inhibitor complexes have been prepared from the crystal structure of TMPKmt cocrystallized with the natural substrate deoxythymidine monophosphate (dTMP) (1GSI) for a training set of 15 thymidine analogues (TMDs) with known activity to prepare a QSAR model of interaction establishing a correlation between the free energy of complexation and the biological activity. Subsequent validation of the predictability of the model has been performed with a 3D QSAR pharmacophore generation. The structural information derived from the model served to design new subnanomolar thymidine analogues. From molecular modeling investigations, the agreement between free energy of complexation (ΔΔG_com) and K_i values explains 94% of the TMPKmt inhibition (pK_i = −0.2924ΔΔG_com + 3.234; R² = 0.94) by variation of the computed ΔΔG_com and 92% for the pharmacophore (PH4) model (pK_i = 1.0206 × pK_i^pred − 0.0832, R² = 0.92). The analysis of contributions from active site residues suggested substitution at the 5-position of pyrimidine ring and various groups at the 5′-position of the ribose. The best inhibitor reached a predicted K_i of 0.155 nM. The computational approach through the combined use of molecular modeling and PH4 pharmacophore is helpful in targeted drug design, providing valuable information for the synthesis and prediction of activity of novel antituberculotic agents.

Synthesis and biological activity of potent HIV-1 protease inhibitors based on Phe-Pro dihydroxyethylene isosteres

Peptidomimetic inhibitors of HIV-1 PR are still a key resource in the fight against AIDS. Here we describe the synthesis and biological activity of HIV-1 PR inhibitors based on four novel dihydroxyethylene isosteres of the Phe-Pro and Pro-Pro dipeptides. The isosteres, containing four stereogenic centers, were synthesized in high yield and excellent stereoselectivity via the cyclization of epoxy amines derived from α-amino acids. The inhibitors were assembled by coupling the isosteres with suitable flanking groups and were screened against recombinant HIV PR showing activities in the subnanomolar to micromolar range. Two Phe-Pro-based inhibitors active at the nanomolar level were further investigated: both inhibitors combine the ability to suppress HIV-1 replication in infected MT-2 cells with low cytotoxicity against the same cells, thereby displaying a high therapeutic index. These results demonstrate the potential of the new Phe-Pro dihydroxyethylene isostere as a core unit of powerful HIV-1 PR inhibitors.

Computer-assisted combinatorial design of bicyclic thymidine analogs as inhibitors of Mycobacterium tuberculosis thymidine monophosphate kinase

Thymidine monophosphate kinase (TMPK(mt)) is an essential enzyme for nucleotide metabolism in Mycobacterium tuberculosis, and thus an attractive target for novel antituberculosis agents. In this work, we have explored the chemical space around the 2',3'-bicyclic thymidine nucleus by designing and in silico screening of a virtual focused library selected via structure based methods to identify more potent analogs endowed with favorable ADME-related properties. In all the library members we have exchanged the ribose ring of the template with a cyclopentane moiety that is less prone to enzymatic degradation. In addition, we have replaced the six-membered 2',3'-ring by a number of five-membered and six-membered heterocyclic rings containing alternative proton donor and acceptor groups, to exploit the interaction with the carboxylate groups of Asp9 and Asp163 as well as with several cationic residues present in the vicinity of the TMPK(mt) binding site. The three-dimensional structure of the TMPK(mt) complexed with 5-hydroxymethyl-dUMP, an analog of dTMP, was employed to develop a QSAR model, to parameterize a scoring function specific for the TMPK(mt) target and to select analogues which display the highest predicted binding to the target. As a result, we identified a small highly focused combinatorial subset of bicyclic thymidine analogues as virtual hits that are predicted to inhibit the mycobacterial TMPK in the submicromolar concentration range and to display favorable ADME-related properties.

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.