Signal Propagation in the ATPase Domain of Mycobacterium tuberculosis DNA Gyrase from Dynamical-Nonequilibrium Molecular Dynamics Simulations.

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A2B2) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Ligand-Based Virtual Screening for Discovery of Indole Derivatives as Potent DNA Gyrase ATPase Inhibitors Active against Mycobacterium tuberculosis and Hit Validation by Biological Assays.

Mycobacterium tuberculosis is the single most important global infectious disease killer and a World Health Organization critical priority pathogen for development of new antimicrobials. M. tuberculosis DNA gyrase is a validated target for anti-TB agents, but those in current use target DNA breakage-reunion, rather than the ATPase activity of the GyrB subunit. Here, virtual screening, subsequently validated by whole-cell and enzyme inhibition assays, was applied to identify candidate compounds that inhibit M. tuberculosis GyrB ATPase activity from the Specs compound library. This approach yielded six compounds: four carbazole derivatives (1, 2, 3, and 8), the benzoindole derivative 11, and the indole derivative 14. Carbazole derivatives can be considered a new scaffold for M. tuberculosis DNA gyrase ATPase inhibitors. IC50 values of compounds 8, 11, and 14 (0.26, 0.56, and 0.08 µM, respectively) for inhibition of M. tuberculosis DNA gyrase ATPase activity are 5-fold, 2-fold, and 16-fold better than the known DNA gyrase ATPase inhibitor novobiocin. MIC values of these compounds against growth of M. tuberculosis H37Ra are 25.0, 3.1, and 6.2 µg/mL, respectively, superior to novobiocin (MIC > 100.0 µg/mL). Molecular dynamics simulations of models of docked GyrB:inhibitor complexes suggest that hydrogen bond interactions with GyrB Asp79 are crucial for high-affinity binding of compounds 8, 11, and 14 to M. tuberculosis GyrB for inhibition of ATPase activity. These data demonstrate that virtual screening can identify known and new scaffolds that inhibit both M. tuberculosis DNA gyrase ATPase activity in vitro and growth of M. tuberculosis bacteria.

Role of interaction mode of phenanthrene derivatives as selective PDE5 inhibitors using molecular dynamics simulations and quantum chemical calculations.

Phosphodiesterase type 5 (PDE5) inhibitors play a crucial role in blocking PDE5 to improve erectile dysfunction (ED). However, most PDE5 drugs revealed side effects including the loss of vision due to the PDE6 inhibition. Phenanthrene derivatives isolated from E. macrobulbon were previously reported as PDE5 inhibitors. Two phenanthrene derivatives (cpds 1-2) revealed better inhibition to PDE5 than PDE6 and cpd 1 is more selective to PDE5 than cpd 2. To elucidate why the phenanthrene derivatives could inhibit PDE5 and PDE6, their binding modes were investigated using molecular dynamics simulations and quantum chemical calculations, as compared to the PDE5 drugs. From the results, all four drugs and phenanthrene derivatives revealed similar π-π interactions to Phe820 in PDE5. Additional H-bond interaction to Gln817 in PDE5 resulted in better PDE5 inhibition of vardenafil and tadalafil. Moreover, cpds 1-2 were able to form the H-bond interaction with Asp764 in PDE5. In the case of the PDE6, the loss of π-π interaction to Phe776 and H-bond interaction to Gln773 indicated the important points for losing the PDE6 inhibition. In conclusion, to develop the new potent PDE5 inhibitors, not only the important interaction with PDE5 but also the interaction with PDE6 should be considered. In phenanthrene derivatives, the middle ring was significant to form π-π interactions to Phe820 in PDE5 and hydroxyl substituent was also the key part to form the H-bond interaction with Asp764 in PDE5. Principal component analysis (PCA) and free energy landscape (FEL) analysis indicated the stability of the system. The bioavailability, drug-likeness, and pharmacokinetics of phenanthrene derivatives were also predicted. These derivatives revealed good drug-likeness and GI absorption. The obtained results showed that phenanthrene derivatives could be interesting for the development of PDE5 inhibitors in the future.

Bioisosteric Design Identifies Inhibitors of Mycobacterium tuberculosis DNA Gyrase ATPase Activity.

Mutations in DNA gyrase confer resistance to fluoroquinolones, second-line antibiotics for Mycobacterium tuberculosis infections. Identification of new agents that inhibit M. tuberculosis DNA gyrase ATPase activity is one strategy to overcome this. Here, bioisosteric designs using known inhibitors as templates were employed to define novel inhibitors of M. tuberculosis DNA gyrase ATPase activity. This yielded the modified compound R3-13 with improved drug-likeness compared to the template inhibitor that acted as a promising ATPase inhibitor against M. tuberculosis DNA gyrase. Utilization of compound R3-13 as a virtual screening template, supported by subsequent biological assays, identified seven further M. tuberculosis DNA gyrase ATPase inhibitors with IC50 values in the range of 0.42-3.59 µM. The most active compound 1 showed an IC50 value of 0.42 µM, 3-fold better than the comparator ATPase inhibitor novobiocin (1.27 µM). Compound 1 showed noncytotoxicity to Caco-2 cells at concentrations up to 76-fold higher than its IC50 value. Molecular dynamics simulations followed by decomposition energy calculations identified that compound 1 occupies the binding pocket utilized by the adenosine group of the ATP analogue AMPPNP in the M. tuberculosis DNA gyrase GyrB subunit. The most prominent contribution to the binding of compound 1 to M. tuberculosis GyrB subunit is made by residue Asp79, which forms two hydrogen bonds with the OH group of this compound and also participates in the binding of AMPPNP. Compound 1 represents a potential new scaffold for further exploration and optimization as a M. tuberculosis DNA gyrase ATPase inhibitor and candidate anti-tuberculosis agent.

A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase.

Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of ß-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-ß-amino acids. This enzyme has the potential to be used as a tool in the production of ß-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.

Stepwise Access of Emissive Ir(III) Complexes Bearing a Multi-Dentate Heteroaromatic Chelate: Fundamentals and Applications.

Three multi-dentate coordinated chelates LnH2 (n = 1, 2, and 3), comprising a linked 1-(pyridin-2-yl)ethylbenzene and one pyrazolyl pyridine unit and showing either tridentate or tetradentate coordination modes, are successfully designed and synthesized. Dinuclear Ir(III) complexes [Ir(κ4-Ln)(µ-Cl)]2 bearing tetradentate coordinated κ4-Ln chelate (2a, n = 1; 2b, n = 2; 2c, n = 3) were next obtained en route from the respective intermediate [Ir(κ3-LnH)Cl(µ-Cl)]2 bearing the tridentate coordinated κ3-LnH chelate (1a, n = 1; 1b, n = 2; 1c, n = 3). Next, mononuclear Ir(III) complexes Ir(κ4-Ln)(thd) (3a, n = 1; 3b, n = 2; 3c, n = 3) with the tetradentate chelate were obtained upon treatment of 2 with 2,2,6,6-tetramethyl-3,5-heptanedione (thd)H in the presence of K2CO3. Concurrently, methylation of 2c in the presence of MeI and nBu4NCl afforded tridentate Ir(κ3-L3HMe)Cl3 (4) and, next, can be converted to tetradentate Ir(κ4-L3Me)Cl2 (5) by further cyclometalation and HCl elimination in refluxing diethylene glycol monoethyl ether solution. The Ir(III) complexes 3a, 4, and 5 were unambiguously identified using spectroscopic methods, together with single-crystal X-ray structural analyses on Ir(III) derivatives 3a, 4, and 5. Their photophysical and ,electrochemical properties and device fabrication properties were also investigated and compared with results from theoretical studies.

Identification of Potent DNA Gyrase Inhibitors Active against Mycobacterium tuberculosis.

Mycobacterium tuberculosis DNA gyrase manipulates the DNA topology using controlled breakage and religation of DNA driven by ATP hydrolysis. DNA gyrase has been validated as the enzyme target of fluoroquinolones (FQs), second-line antibiotics used for the treatment of multidrug-resistant tuberculosis. Mutations around the DNA gyrase DNA-binding site result in the emergence of FQ resistance in M. tuberculosis; inhibition of DNA gyrase ATPase activity is one strategy to overcome this. Here, virtual screening, subsequently validated by biological assays, was applied to select candidate inhibitors of the M. tuberculosis DNA gyrase ATPase activity from the Specs compound library (www.specs.net). Thirty compounds were identified and selected as hits for in vitro biological assays, of which two compounds, G24 and G26, inhibited the growth of M. tuberculosis H37Rv with a minimal inhibitory concentration of 12.5 µg/mL. The two compounds inhibited DNA gyrase ATPase activity with IC50 values of 2.69 and 2.46 µM, respectively, suggesting this to be the likely basis of their antitubercular activity. Models of complexes of compounds G24 and G26 bound to the M. tuberculosis DNA gyrase ATP-binding site, generated by molecular dynamics simulations followed by pharmacophore mapping analysis, showed hydrophobic interactions of inhibitor hydrophobic headgroups and electrostatic and hydrogen bond interactions of the polar tails, which are likely to be important for their inhibition. Decreasing compound lipophilicity by increasing the polarity of these tails then presents a likely route to improving the solubility and activity. Thus, compounds G24 and G26 provide attractive starting templates for the optimization of antitubercular agents that act by targeting DNA gyrase.

Virtual Screening Identifies Novel and Potent Inhibitors of Mycobacterium tuberculosis PknB with Antibacterial Activity.

Mycobacterium tuberculosis protein kinase B (PknB) is essential to mycobacterial growth and has received considerable attention as an attractive target for novel anti-tuberculosis drug development. Here, virtual screening, validated by biological assays, was applied to select candidate inhibitors of M. tuberculosis PknB from the Specs compound library (www.specs.net). Fifteen compounds were identified as hits and selected for in vitro biological assays, of which three indoles (2, AE-848/42799159; 4, AH-262/34335013; 10, AP-124/40904362) inhibited growth of M. tuberculosis H37Rv with minimal inhibitory concentrations of 6.2, 12.5, and 6.2 µg/mL, respectively. Two compounds, 2 and 10, inhibited M. tuberculosis PknB activity in vitro, with IC50 values of 14.4 and 12.1 µM, respectively, suggesting this to be the likely basis of their anti-tubercular activity. In contrast, compound 4 displayed anti-tuberculosis activity against M. tuberculosis H37Rv but showed no inhibition of PknB activity (IC50 > 128 µM). We hypothesize that hydrolysis of its ethyl ester to a carboxylate moiety generates an active species that inhibits other M. tuberculosis enzymes. Molecular dynamics simulations of modeled complexes of compounds 2, 4, and 10 bound to M. tuberculosis PknB indicated that compound 4 has a lower affinity for M. tuberculosis PknB than compounds 2 and 10, as evidenced by higher calculated binding free energies, consistent with experiment. Compounds 2 and 10 therefore represent candidate inhibitors of M. tuberculosis PknB that provide attractive starting templates for optimization as anti-tubercular agents.

Computational design, synthesis and biological evaluation of PDE5 inhibitors based on N2,N4-diaminoquinazoline and N2,N6-diaminopurine scaffolds.

We report the synthesis, and characterization of twenty-nine new inhibitors of PDE5. Structure-based design was employed to modify to our previously reported 2,4-diaminoquinazoline series. Modification include scaffold hopping to 2,6-diaminopurine core as well as incorporation of ionizable groups to improve both activity and solubility. The prospective binding mode of the compounds was determined using 3D ligand-based similarity methods to inhibitors of known binding mode, combined with a PDE5 docking and molecular dynamics based-protocol, each of which pointed to the same binding mode. Chemical modifications were then designed to both increase potency and solubility as well as validate the binding mode prediction. Compounds containing a quinazoline core displayed IC50s ranging from 0.10 to 9.39 µM while those consisting of a purine scaffold ranging from 0.29 to 43.16 µM. We identified 25 with a PDE5 IC50 of 0.15 µM, and much improved solubility (1.77 mg/mL) over the starting lead. Furthermore, it was found that the predicted binding mode was consistent with the observed SAR validating our computationally driven approach.

Biophysical Characterization of Novel DNA Aptamers against K103N/Y181C Double Mutant HIV-1 Reverse Transcriptase.

The human immunodeficiency virus type-1 Reverse Transcriptase (HIV-1 RT) plays a pivotal role in essential viral replication and is the main target for antiviral therapy. The anti-HIV-1 RT drugs address resistance-associated mutations. This research focused on isolating the potential specific DNA aptamers against K103N/Y181C double mutant HIV-1 RT. Five DNA aptamers showed low IC50 values against both the KY-mutant HIV-1 RT and wildtype (WT) HIV-1 RT. The kinetic binding affinity forms surface plasmon resonance of both KY-mutant and WT HIV-1 RTs in the range of 0.06-2 µM and 0.15-2 µM, respectively. Among these aptamers, the KY44 aptamer was chosen to study the interaction of HIV-1 RTs-DNA aptamer complex by NMR experiments. The NMR results indicate that the aptamer could interact with both WT and KY-mutant HIV-1 RT at the NNRTI drug binding pocket by inducing a chemical shift at methionine residues. Furthermore, KY44 could inhibit pseudo-HIV particle infection in HEK293 cells with nearly 80% inhibition and showed low cytotoxicity on HEK293 cells. These together indicated that the KY44 aptamer could be a potential inhibitor of both WT and KY-mutant HIV-RT.

Characterization of New DNA Aptamers for Anti-HIV-1 Reverse Transcriptase.

HIV-1 RT is a necessary enzyme for retroviral replication, which is the main target for antiviral therapy against AIDS. Effective anti-HIV-1 RT drugs are divided into two groups; nucleoside inhibitors (NRTI) and non-nucleoside inhibitors (NNRTI), which inhibit DNA polymerase. In this study, new DNA aptamers were isolated as anti-HIV-1 RT inhibitors. The selected DNA aptamer (WT62) presented with high affinity and inhibition against wild-type (WT) HIV-1 RT and gave a KD value of 75.10±0.29 nM and an IC50 value of 84.81±8.54 nM. Moreover, WT62 decreased the DNA polymerase function of K103 N/Y181 C double mutant (KY) HIV-1 RT by around 80 %. Furthermore, the ITC results showed that this aptamer has small binding enthalpies with both WT and KY HIV-1 RTs through which the complex might form a hydrophobic interaction or noncovalent bonding. The NMR result also suggested that the WT62 aptamer could bind with both WT and KY mutant HIV-1 RTs at the connection domain.

Preparation, biological & cheminformatics-based assessment of N2,N4-diphenylpyrimidine-2,4-diamine as potential Kinase-targeted antimalarials.

Twenty eight new N2,N4-diphenylpyrimidine-2,4-diamines have been prepared in order to expand our understanding of the anti-malarial SAR of the scaffold. The aim of the study was to make structural modifications to improve the overall potency, selectivity and solubility of the series by varying the anilino groups attached to the 2- and 4-position. We evaluated the activity of the compounds against Plasmodium falciparum (Pf) 3D7, cytotoxicity against HepG2, % inhibition at a panel of 10 human kinases, solubility, permeability and lipophilicity, and human and rat in vitro clearance. 11 was identified as a potent anti-malarial with an IC50 of 0.66 µM at the 3D7 strain and a selectivity (SI) of ~ 40 in terms of cytotoxicity against the HepG2 cell line. It also displayed low experimental logD7.4 (2.27), reasonable solubility (124 µg/ml), good metabolic stability, but low permeability. A proteo-chemometric workflow was employed to identify putative Pf targets of the most promising compounds. Ligand-based similarity searching of the ChEMBL database led to the identification of most probable human targets. These were then used as input for sequence-based searching of the Pf proteome. Homology modelling and molecular docking were used to evaluate whether compounds could indeed bind to these targets with valid binding modes. In vitro biological testing against close human analogs of these targets was subsequently undertaken. This allowed us to identify potential Pf targets and human anti-targets that could be exploited in future development.

Discovery of New and Potent InhA Inhibitors as Antituberculosis Agents: Structure-Based Virtual Screening Validated by Biological Assays and X-ray Crystallography.

The enoyl-acyl carrier protein reductase InhA of Mycobacterium tuberculosis is an attractive, validated target for antituberculosis drug development. Moreover, direct inhibitors of InhA remain effective against InhA variants with mutations associated with isoniazid resistance, offering the potential for activity against MDR isolates. Here, structure-based virtual screening supported by biological assays was applied to identify novel InhA inhibitors as potential antituberculosis agents. High-speed Glide SP docking was initially performed against two conformations of InhA differing in the orientation of the active site Tyr158. The resulting hits were filtered for drug-likeness based on Lipinski's rule and avoidance of PAINS-like properties and finally subjected to Glide XP docking to improve accuracy. Sixteen compounds were identified and selected for in vitro biological assays, of which two (compounds 1 and 7) showed MIC of 12.5 and 25 µg/mL against M. tuberculosis H37Rv, respectively. Inhibition assays against purified recombinant InhA determined IC50 values for these compounds of 0.38 and 0.22 µM, respectively. A crystal structure of the most potent compound, compound 7, bound to InhA revealed the inhibitor to occupy a hydrophobic pocket implicated in binding the aliphatic portions of InhA substrates but distant from the NADH cofactor, i.e., in a site distinct from those occupied by the great majority of known InhA inhibitors. This compound provides an attractive starting template for ligand optimization aimed at discovery of new and effective compounds against M. tuberculosis that act by targeting InhA.

Design, synthesis and evaluation of N2,N4-diaminoquinazoline based inhibitors of phosphodiesterase type 5.

We describe the design, synthesis and evaluation of a series of N2,N4-diaminoquinazoline analogs as PDE5 inhibitors. Twenty compounds were prepared and these were assessed in terms of their PDE5 and PDE6 activity, ex-vivo vasodilation response, mammalian cytotoxicity and aqueous solubility. Molecular docking was used to determine the binding mode of the series and this was demonstrated to be consistent with the observed SAR. Compound 15 was the most active PDE5 inhibitor (IC50 = 0.072 ±â€¯0.008 µM) and exhibited 4.6-fold selectivity over PDE6. Ex-vivo assessment of 15 and 22 in a rat pulmonary artery vasodilation model demonstrated EC50s of 1.63 ±â€¯0.72 µM and 2.28 ±â€¯0.74 µM respectively.

Simulations of Shikimate Dehydrogenase from Mycobacterium tuberculosis in Complex with 3-Dehydroshikimate and NADPH Suggest Strategies for MtbSDH Inhibition.

Shikimate dehydrogenase (SDH) from Mycobacterium tuberculosis ( MtbSDH), encoded by the aroE gene, is essential for viability of M. tuberculosis but absent from humans. Therefore, it is a potentially promising target for antituberculosis agent development. Molecular-level understanding of the interactions of MtbSDH with its 3-dehydroshikimate (DHS) substrate and NADPH cofactor will help in the design of novel and effective MtbSDH inhibitors. However, this is limited by the lack of relevant crystal structures for MtbSDH complexes. Here, molecular dynamics (MD) simulations were performed to generate these MtbSDH complexes and investigate interactions of MtbSDH with substrate and cofactor and the role of MtbSDH dynamics within these. The results indicate that, while structural rearrangements are not necessary for DHS binding, reorientation of individual side chains in the NADPH binding pocket is involved in ternary complex formation. The mechanistic roles for Lys69, Asp105, and Ala213 were investigated by generating Lys69Ala, Asp105Asn, and Ala213Leu mutants in silico and investigating their complexes with DHS and NADPH. Our results show that Lys69 plays a dual role, in positioning NADPH and in catalysis. Asp105 plays a crucial role in positioning both the ε-amino group of Lys69 and nicotinamide ring of NADPH for MtbSDH catalysis but makes no direct contribution to DHS binding. Ala213 is the selection key for NADPH binding with the nicotinamide ring in the proS, rather than proR, conformation in the MtbSDH complex. Our results identify three strategies for MtbSDH inhibition: prevention of MtbSDH binary and ternary complex formation by blocking DHS and NADPH binding (first and second strategies, respectively) and the prevention of MtbSDH complex formation with either DHS or NADPH by blocking both DHS and NADPH binding (third strategy). Further, based on this third strategy, we propose guidelines for the rational design of "hybrid" MtbSDH inhibitors able to bind in both the substrate (DHS) and cofactor (NADPH) pockets, providing a new avenue of exploration in the search for anti-TB therapeutics.

Elucidation of Vasodilation Response and Structure Activity Relationships of N²,N4-Disubstituted Quinazoline 2,4-Diamines in a Rat Pulmonary Artery Model.

Pulmonary arterial hypertension (PAH) is a rare and progressive disease arising from various etiologies and pathogenesis. PAH decreases life expectancy due to pulmonary vascular remodeling, elevation of mean pulmonary arterial pressure, and ultimately progresses to heart failure. While clinical treatments are available to reduce the associated symptoms, a complete cure has yet to be found. Phosphodiesterase-5 (PDE-5) inhibition has been identified as a possible intervention point in PAH treatment. The functional vasodilation response to N²,N4-diamino quinazoline analogues with differing PDE-5 inhibitory activities and varying physicochemical properties were assessed in both endothelium-intact and denuded rat pulmonary arteries to gain greater insight into their mode of action. All analogues produced vasorelaxant effects with EC50s ranging from 0.58 ± 0.22 µM to ˃30 µM. It was observed that vasodilation response in intact vessels was highly correlated with that of denuded vessels. The ~10% drop in activity is consistent with a loss of the nitric oxide mediated cyclic guanosine monophosphate (NO/cGMP) pathway in the latter case. A moderate correlation between the vasodilation response and PDE-5 inhibitory activity in the intact vessels was observed. Experimental protocol using the alpha-adrenergic (α1) receptor agonist, phenylephrine (PE), was undertaken to assess whether quinazoline derivatives showed competitive behavior similar to the α1 receptor blocker, prazosin, itself a quinazoline derivative, or to the PDE-5 inhibitor, sildenafil. Competitive experiments with the α1-adrenergic receptor agonist point to quinazoline derivatives under investigation here act via PDE-5 inhibition and not the former. The pre-incubation of pulmonary arterial rings with quinazoline test compounds (10 µM) reduced the contractile response to PE around 40⁻60%. The most promising compound (9) possessed ~32 folds higher selectivity in terms of vasodilation to its mammalian A549 cell cytotoxicity. This study provides experi0 0mental basis for PDE-5 inhibition as the mode of action for vasodilation by N²,N4-diamino quinazoline analogues along with their safety studies that may be beneficial in the treatment of various cardiovascular pathologies.

Synthesis and evaluation of the NSCLC anti-cancer activity and physical properties of 4-aryl-N-phenylpyrimidin-2-amines.

Reported herein are efforts to profile 4-aryl-N-phenylpyrimidin-2-amines in terms of their anti-cancer activity towards non small-cell lung carcinoma (NSCLC) cells. We have synthesized new 4-aryl-N-phenylpyrimidin-2-amines and assessed them in terms of their cytotoxicity (A549, NCI-H187, MCF7, Vero & KB) and physicochemical properties (logD7.4 and solubility). 13f and 13c demonstrated potent anti-cancer activity in A549 cells (0.2µM), compared to 0.4µM for the NSCLC drug Doxorubicin. 13f also displayed low experimental logD7.4 (2.9) and the best solubility (∼40µM). Compounds 13b and 13d showed the best balance of A549 anti-cancer activity and selectivity. 13g showed good activity and selectivity comparable with the anti-cancer drug Doxorubicin.

Rapid activity prediction of HIV-1 integrase inhibitors: harnessing docking energetic components for empirical scoring by chemometric and artificial neural network approaches.

Improving performance of scoring functions for drug docking simulations is a challenging task in the modern discovery pipeline. Among various ways to enhance the efficiency of scoring function, tuning of energetic component approach is an attractive option that provides better predictions. Herein we present the first development of rapid and simple tuning models for predicting and scoring inhibitory activity of investigated ligands docked into catalytic core domain structures of HIV-1 integrase (IN) enzyme. We developed the models using all energetic terms obtained from flexible ligand-rigid receptor dockings by AutoDock4, followed by a data analysis using either partial least squares (PLS) or self-organizing maps (SOMs). The models were established using 66 and 64 ligands of mercaptobenzenesulfonamides for the PLS-based and the SOMs-based inhibitory activity predictions, respectively. The models were then evaluated for their predictability quality using closely related test compounds, as well as five different unrelated inhibitor test sets. Weighting constants for each energy term were also optimized, thus customizing the scoring function for this specific target protein. Root-mean-square error (RMSE) values between the predicted and the experimental inhibitory activities were determined to be <1 (i.e. within a magnitude of a single log scale of actual IC50 values). Hence, we propose that, as a pre-functional assay screening step, AutoDock4 docking in combination with these subsequent rapid weighted energy tuning methods via PLS and SOMs analyses is a viable approach to predict the potential inhibitory activity and to discriminate among small drug-like molecules to target a specific protein of interest.

Mass Spectrometric Characterization of HIV-1 Reverse Transcriptase Interactions with Non-nucleoside Reverse Transcriptase Inhibitors.

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) of human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) have been developed for the treatment of acquired immunodeficiency syndrome. HIV-1 RT binding to NNRTIs has been characterized by various biophysical techniques. However, these techniques are often hampered by the low water solubility of the inhibitors, such as the current promising diarylpyrimidine-based inhibitors rilpivirine and etravirine. Hence, a conventional and rapid method that requires small sample amounts is desirable for studying NNRTIs with low water solubility. Here we successfully applied a recently developed mass spectrometric technique under non-denaturing conditions to characterize the interactions between the heterodimeric HIV-1 RT enzyme and NNRTIs with different inhibitory activities. Our data demonstrate that mass spectrometry serves as a semi-quantitative indicator of NNRTI binding affinity for HIV-1 RT using low and small amounts of samples, offering a new high-throughput screening tool for identifying novel RT inhibitors as anti-HIV drugs.

Skin sensitization prediction using quantum chemical calculations: a theoretical model for the S(N)Ar domain.

It is widely accepted that skin sensitization begins with the sensitizer in question forming a covalent adduct with a protein electrophile or nucleophile. We investigate the use of quantum chemical methods in an attempt to rationalize the sensitization potential of chemicals of the S(N)Ar reaction domain. We calculate the full reaction profile for 23 chemicals with experimental sensitization data. For all quantitative measurements, we find that there is a good correlation between the reported pEC3 and the calculated barrier to formation of the low energy product or intermediate (r(2) = 0.64, N = 12) and a stronger one when broken down by specific subtype (r(2) > 0.9). Using a barrier cutoff of ∼10 kcal/mol allows us to categorize 100% (N = 12) of the sensitizers from the nonsensitizers (N = 11), with just 1 nonsensitizer being mispredicted as a weak sensitizer (9%). This model has an accuracy of ∼96%, with a sensitivity of 100% and a specificity of ∼91%. We find that the kinetic and thermodynamic information provided by the complete profile can help in the rationalization process, giving additional insight into a chemical's potential for skin sensitization.