Emergence of Dengue as one of the deadliest viral diseases prompts the need for development of effective therapeutic agents. Dengue virus (DV) exists in four different serotypes and infection caused by one serotype predisposes its host to another DV serotype heterotypic re-infection. We undertook virtual ligand screening (VLS) to filter compounds against DV that may inhibit inclusively all of its serotypes. Conserved non-structural DV protein targets such as NS1, NS3/NS2B and NS5, which play crucial role in viral replication, infection cycle and host interaction, were selected for screening of vital antiviral drug leads. A dataset of plant based natural antiviral derivatives was developed. Molecular docking was performed to estimate the spatial affinity of target compounds for the active sites of DV’s NS1, NS3/NS2B and NS5 proteins. The drug likeliness of the screened compounds was followed by ADMET analysis whereas the binding behaviors were further elucidated through molecular dynamics (MD) simulation experiments. VLS screened three potential compounds including Canthin-6-one 9-O-beta-glucopyranoside, Kushenol W and Kushenol K which exhibited optimal binding with all the three conserved DV proteins. This study brings forth novel scaffolds against DV serotypes to serve as lead molecules for further optimization and drug development against all DV serotypes with equal effect against multiple disease causing DV proteins. We therefore anticipate that the insights given in the current study could be regarded valuable towards exploration and development of a broad-spectrum natural anti-dengue therapy.

Pseudomonas aeruginosa is a Gram-negative pathogenic bacterium that is present commonly in soil and water and is responsible for causing septic shock, pneumonia, urinary tract and gastrointestinal infections, etc. The multi-drug resistance (MDR) phenomenon has increased dramatically in past years and is now considered a major threat globally, so there is an urgent need to develop new strategies to overcome drug resistance by P. aeruginosa. In P. aeruginosa, a major factor of drug resistance is associated to the formation of biofilms by the LasR enzyme, which regulates quorum sensing and has been reported as a new therapeutic target for designing novel antibacterial molecules. In this study, virtual screening and molecular docking were performed against the ligand binding domain (LBD) of LasR by employing a pharmacophore hypothesis for the screening of 2373 FDA-approved compounds to filter top-scoring hit compounds. Six inhibitors out of 2373 compounds were found to have binding affinities close to that of known LasR inhibitors. The binding modes of these compounds to the binding site in LasR-LBD were analyzed to identify the key interactions that contribute to the inhibition of LasR activity. Then, 50 ns simulations of top hit compounds were performed to elucidate the stability of their binding conformations with the LasR-LBD. This study, thus concluded that sulfamerazine showed the highest binding affinity for the LasR-LBD binding pocket exhibiting strong inhibitory binding interactions during molecular dynamics (MD) simulation.

In the cyclic guanosine monophosphate (cGMP) signaling pathway, phosphodiesterase 6 (PDE6) maintains a critical balance of the intracellular concentration of cGMP by catalyzing it to 5' guanosine monophosphate (5'-GMP). To gain insight into the mechanistic impacts of the PDE6 somatic mutations that are implicated in cancer and retinitis pigmentosa, we first defined the structure and organization of the human PDE6 heterodimer using computational comparative modelling. Each subunit of PDE6αβ possesses three domains connected through long α-helices. The heterodimer model indicates that the two chains are likely related by a pseudo two-fold axis. The N-terminal region of each subunit is comprised of two allosteric cGMP-binding domains (Gaf-A & Gaf-B), oriented in the same way and interacting with the catalytic domain present at the C-terminal in a way that would allow the allosteric cGMP-binding domains to influence catalytic activity. Subsequently, we applied an integrated knowledge-driven in silico mutation analysis approach to understand the structural and functional implications of experimentally identified mutations that cause various cancers and retinitis pigmentosa, as well as computational saturation mutagenesis of the dimer interface and cGMP-binding residues of both Gaf-A, and the catalytic domains. We studied the impact of mutations on the stability of PDE6αβ structure, subunit-interfaces and Gaf-cGMP interactions. Further, we discussed the changes in interatomic interactions of mutations that are destabilizing in Gaf-A (R93L, V141 M, F162 L), catalytic domain (D600N, F742 L, F776 L) and at the dimer interface (F426A, F248G, F424 N). This study establishes a possible link of change in PDE6αβ structural stability to the experimentally observed disease phenotypes.

The need of circumventing life-threatening cardiovascular disorders (CVDs) and pulmonary hypertension (PHT) worldwide prompts researchers to develop effective therapeutic agents. Crucial role of cyclic nucleotide phosphodiesterase-5 (PDE5A) and cyclic nucleotide phosphodiesterase-3 (PDE3A) in cardiovascular signaling makes them potential drug targets for the treatment of CVDs and PHT. In this study, one-drug-multiple-target strategy has been employed to screen inhibitors exhibiting dual specificity through Phase-generated and statistically validated e-pharmacophore models of PDE5A and PDE3A. An extensive CoCoCo database of 7 million compounds with similar to 150,000,000 conformations was virtually screened by sequential e-pharmacophore mapping followed by Lipinski Rule of Five (RO5) evaluation and hierarchical docking simulations. Finally, docked hits were subjected to rigorous MMGBSA analysis to estimate the relative spatial affinity of the drug-like compounds. The hits (354 and 366 ligands against PDE5A and PDE3A, respectively) were further optimized through 2D clustering followed by a comprehensive 2D and 3D interaction analysis. Five structurally diverse hits mapped equally well with the e-pharmacophore models and showed promising inhibitory interactions with conserved four catalytic features of PDE5A and PDE3A, thus exhibiting dual specificity. Proposed lead compounds exhibited the lowest MMGBSA binding energies and were found to be in agreement with Lipinski Rule of Five (RO5) and ADME/Tox criteria as compared to sildenafil. The proposed dual inhibitors could thus provide promising outcomes for the discovery of dual as well as multipotent drug like compounds after lead optimization and primary therapeutic interventions.

Cardiovascular disorders are the world’s major cause of death nowadays. To treat cardiovascular diseases especially coronary artery diseases and hypertension, researchers found potential ROCK2 (Rho-associated coiled-coil-containing protein kinase 2) target due to its substantial role in NO-cGMP and RhoA/ROCK pathway. Available drugs for ROCK2 are less effective and some of them depict side effects. Therefore, a set of novel compounds were screened that can potentially inhibit the activity of ROCK2 and help to treat cardiovascular diseases by employing In-silico techniques. In this study, we undertook ligand based virtual screening of 50 million compound’s library, to that purpose shape and features (contain functional groups) based pharmacophore query was modelled and validated by Area Under Curve graph (AUC). 2000 best hits were screened for Lipinski’s rule of 5 compliance. Subsequently, these selected compounds were docked into the binding site of ROCK2 to gain insights into the interactions between hit compounds and the target protein. Based on binding affinity and RMSD scores, a final cohort of 15 compounds were chosen which were further refined by pharmacokinetics, ADMET and bioactivity scores. 2 potential hits were screened using density functional theory, revealing remarkable biological and chemical activity. Potential inhibitors (F847-0007 and 9543495) underwent rigorous examination through MD Simulations and MMGBSA analysis, elucidating their stability and strong binding affinities. Results of current study unveil the potential of identified novel hits as promising lead compounds for ROCK2 associated with cardiovascular diseases. These findings will further investigate via In-vitro and In-vivo studies to develop novel druglike molecules against ROCK2.

Soluble guanylate cyclase (sGC) is a key enzyme implicated in various physiological processes such as vasodilation, thrombosis and platelet aggregation. The enzyme’s Heme-Nitric oxide/Oxygen (H-NOX) binding domain is the only sensor of nitric oxide (NO) in humans, which on binding with NO activates sGC to produce the second messenger cGMP. H-NOX is thus a hot target for drug design programs. BAY60-2770 and BAY58-2667 are two widely studied activators of sGC. Here we present comparative molecular dynamics studies to understand the molecular details characterizing the binding of BAY60-2770 and BAY58-2667 with the human H-NOX (hH-NOX) and bacterial H-NOX (bH-NOX) domains. HartreeFock method was used for parametrization of both the activators. A 50 ns molecular dynamics (MD) simulation was run to identify the functionally critical regions of the H-NOX domains. The CPPTRAJ module was used for analysis. BAY60-2770 on binding with bH-NOX, triggered rotational movement in signaling helix F and significant dynamicity in loops α and β, but in hH-NOX domain the compound showed relatively lesser aforementioned structural fluctuations. Conversely, hH-NOX ligated BAY58-2667 experienced highest transitions in its helix F due to electrostatic interactions with D84, T85 and R88 residues which are not conserved in bH-NOX. These conformational transformations might be essential to communicate with downstream PAS, CC and cyclase domains of sGC. Comparative MD studies revealed that BAY bound bHNOX dynamics varied from that of hH-NOX, plausibly due to some key residues such as R40, F74 and Y112 which are not conserved in bacteria. These findings will help to the design of novel drug leads to cure diseases associated to human sGC.