Search

A typical drug discovery project involves identifying active compounds with significant binding potential for selected disease-specific targets. Experimental high-throughput screening (HTS) is a traditional approach to drug discovery, but is expensive and time-consuming when dealing with huge chemical libraries with billions of compounds. The search space can be narrowed down with the use of reliable computational screening approaches. In this review, we focus on various machine-learning (ML) and deep-learning (DL)-based scoring functions developed for solving classification and ranking problems in drug discovery. We highlight studies in which ML and DL models were successfully deployed to identify lead compounds for which the experimental validations are available from bioassay studies.

Machine learning (ML) research based on the quantitative structure–property relationship (QSPR) has been applied for the development of highly efficient dye-sensitized solar cells (DSSCs).

New variants of SARS-CoV-2 are being reported worldwide. The World Health Organization has reported Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2) and Omicron (B.1.1.529) as the variants of concern. There are speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19 occurs through the interaction of the spike protein's receptor-binding domain (RBD) with the peptidase domain of the human ACE-2 (hACE-2) receptor. This study aims to get a molecular-level understanding of the mechanism behind the increased infection rate in the alpha variant. We have computationally studied the spike protein interaction in both the wild-type and B.1.1.7 variant with the hACE-2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The binding free energy difference shows that the mutant variant of the spike protein has increased binding affinity for the hACE-2 receptor (

Alzheimer's disease is characterized by the accumulation of amyloid beta (Aβ) and Tau aggregates in the brain, which induces various pathological events resulting in neurodegeneration. There have been continuous efforts to develop modulators of the Aβ and Tau aggregation process to halt or modify disease progression. A few small-molecule-based inhibitors that target both Aβ and Tau pathology have been reported. Here, we report the screening of a targeted library of small molecules to modulate Aβ and Tau aggregation together with their in vitro, in silico and cellular studies. In vitro ThT fluorescence assay, dot blot assay, gel electrophoresis and transmission electron microscopy (TEM) results have shown that thiophene-based lead molecules effectively modulate Aβ aggregation and inhibit Tau aggregation. In silico studies performed by employing molecular docking, molecular dynamics and binding-free energy calculations have helped in understanding the mechanism of interaction of the lead thiophene compounds with Aβ and Tau fibril targets. In cellulo studies revealed that the lead candidate is biocompatible and effectively ameliorates neuronal cells from Aβ and Tau-mediated amyloid toxicity.

The discovery of small molecules that exhibit turn-on far-red or near-infrared (NIR) fluorescence upon DNA binding and understanding how they bind DNA are important for imaging and bioanalytical applications. Here we report the DNA-bound structure and the DNA binding mechanism of quinone cyanine dithiazole (QCy-DT), a recently reported AT-specific turn-on NIR fluorescent probe for double-stranded DNA. The nuclear magnetic resonance (NMR)-derived structure showed minor groove binding but no specific ligand-DNA interactions, consistent with an endothermic and entropy-driven binding mechanism deduced from isothermal titration calorimetry. Minor groove binding is typically fast because it minimally perturbs the DNA structure. However, QCy-DT exhibited unusually slow DNA binding. The cyanine-based probe is capable of cis-trans isomerization due to overlapping methine bridges, with 16 possible slowly interconverting cis/trans isomers. Using NMR, density functional theory, and free energy calculations, we show that the DNA-free and DNA-bound environments of QCy-DT prefer distinctly different isomers, indicating that the origin of the slow kinetics is a cis-trans isomerization and that the minor groove preferentially selects an otherwise unstable cis/trans isomer of QCy-DT. Flux analysis showed the conformational selection pathway to be the dominating DNA binding mechanism at low DNA concentrations, which switches to the induced fit pathway at high DNA concentrations. This report of cis/trans isomerization of a ligand, upon binding the DNA minor groove, expands the prevailing understanding of unique discriminatory powers of the minor groove and has an important bearing on using polymethine cyanine dyes to probe the kinetics of molecular interactions.

In this study, we investigated differential molecular interactions of crystalline and amorphous forms of telmisartan (TEL), which is a non-peptide angiotensin-II receptor antagonist commonly used in the management of hypertension. Amorphous telmisartan (AM-TEL) was prepared using quench cooling of the melt. The analysis of solid-state properties of AM-TEL using differential scanning calorimetry (DSC) and X-ray powder diffraction (XRD) confirmed formation of AM-TEL. Based on a comparative analysis of molecular interactions using spectral (FTIR and 13C solid-state NMR) and computational tools, we demonstrated that amorphous telmisartan shows altered molecular interactions. Molecular dynamics simulation of amorphous and crystalline forms demonstrate that the amorphous form retained some of the molecular interactions in its disordered molecular arrangement, with a relatively stronger (decrease in bond length) but lesser (up to only 2.6 % of the population) hydrogen bonding network as compared with the crystalline counterpart (up to 76% of the population)

With the promising advantages of the near-infrared region (NIR) emissive markers for serum albumin becoming very prominent recently, we devised CyG-NHS as the cyanine derived longest NIR-I emissive optical marker possessing albumin selective recognition ability in diverse biological milieu. Multiscale modeling involving molecular docking, molecular dynamics, and implicit solvent binding free energy calculations have been employed to gain insights into the unique binding ability of the developed probe at domain-I of albumin, in contrast to the good number of domain IIA or IIIA binding probes available in the literature reports. The binding free energy was found to be -31.8 kcal mol

We present fluorogenic cationic organo chalcogens that are highly selective to RNA. We have demonstrated that the conformational dynamics and subsequently the optical properties of these dyes can be controlled to facilitate efficient bioimaging. We report the application of organoselenium and organosulfur-based cell-permeable red-emissive probes bearing a favorable cyclic sidearm for selective and high contrast imaging of cell nucleoli. The probes exhibit high quantum yield upon interacting with RNA in an aqueous solution. An in-depth multiscale simulation study reveals that the prominent rotational freezing of the electron-donating sidearm of the probes in the microenvironment of RNA helps in attaining more planar conformation when compared to DNA. It exerts a greater extent of intramolecular charge transfer and hence leads to enhanced fluorescence emission. A systematic structure-interaction relationship study highlighted the impact of heavy-chalcogens toward the improved emissive properties of the probes.

The selective recognition and imaging of oncogene specific G-quadruplex (GQ) structures holds great promise in the development of diagnostic therapy (theranostics) for cancer and has been challenging due to their structural dynamics and diversity. We report selective recognition of GQ by a small molecule through unique hybrid loop stacking and groove binding mode with turn on far-red fluorescence response and anticancer activity demonstrating the potential implications for GQ-targeted cancer theranostics. Methods: Biophysical investigation reveal the turn on far-red emission property of TGP18 for selective recognition of GQ. In cellulo studies including DNA damage and oxidative stress evaluation guided us to perform in vitro (3D spheroid) and in vivo (xenograft mice model) anti-cancer activity, and tumor tissue imaging to assess the theranostic potential of TGP18. Results: Neocuproine-based far-red turn on fluorescence probe TGP18 shows GQ-to-duplex selectivity and specifically recognizes BCL-2 GQ with high affinity through a unique hybrid binding mode involving loop-stacking and groove interactions. Our study reveals that the selective recognition originating from the distinct loop structure of GQ that alters the overall probe interaction and binding affinity. TGP18 binding to anti-apoptotic BCL-2 GQ ablates the pro-survival function and elicit anti-cancer activity by inducing apoptosis in cancer cells. We deciphered that inhibition of BCL-2 transcription synergized with signaling cascade of nucleolar stress, DNA damage and oxidative stress in triggering apoptosis signaling pathway. Conclusion: Intervention of GQ mediated lethality by TGP18 has translated into anti-cancer activity in both in vitro 3D spheroid culture and in vivo xenograft models of lung and breast cancer with superior efficacy for the former. In vivo therapeutic efficacy supplemented with tumor 3D spheroid and tissue imaging potential define the role of TGP18 in GQ-targeted cancer theranostics.

Abnormal deposition of hyperphosphorylated tau as neurofibrillary tangles (NFTs) is an important pathological hallmark of Alzheimer's disease (AD) and of other neurodegenerative disorders. A noninvasive positron emission tomography (PET) tracer that quantifies neurofibrillary tangles in vivo can enhance the clinical diagnosis of AD and can also be used to evaluate the efficacy of therapeutics aimed at reducing the abnormal aggregation of the tau fibril in the brain. In this paper, we study the binding profile of fibrillar tau aggregates with a PET tracer PI2620, which is a new second generation tau PET tracer that is presently experimentally and clinically studied. The target structure for the tau fibril is based on cryo-electron microscopy (cryo-EM) structure. A multiscale simulation workflow including molecular docking, molecular dynamics simulation, metadynamics simulation, and free energy calculations was implemented. We find that PI2620 can bind to eight surface binding sites, three core binding sites, and one entry site. The binding at the core sites and entry site is found to be much more favorable than that on the surface sites due to stronger hydrophobic interactions and less solvent exposure. Furthermore, the entry site which is formed by the terminal β-sheets of the fibril is found to have the highest binding affinity to PI2620. Importantly, the binding capacity at the entry site can be much higher than that at other core sites, due to its easy accessibility. Therefore, the entry site is believed to be the major binding site for PI2620. A previous computational study on tracers with tau fibrils reports a maximum of four binding sites. Through use of methods that allow us to locate "cryptic binding sites", we report here additional core sites available for binding and we address the limitation of using the cryo-EM structure alone for structure-based tracer design. Our results could be helpful for elucidating the binding mechanism of imaging tracers with the fibrillar form of tau, a knowledge that in turn can be used to guide the development of compounds with higher affinity and selectivity for tau using structure-based design strategies.

The new variant of SARS-CoV-2, Omicron, has been quickly spreading in many countries worldwide. Compared to the original virus, Omicron is characterized by several mutations in its genomic region, including spike protein’s receptor-binding domain (RBD). We have computationally investigated the interaction between RBD of both wild-type and omicron variants with hACE2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The mode of the interaction between Omicron’s RBD to the human ACE2 (hACE2) receptor is similar to the original SARS-CoV-2 RBD except for a few key differences. The binding free energy difference shows that the spike protein of Omicron has increased binding affinity for the hACE-2 receptor. The mutated residues in the RBD showed strong interactions with a few amino acid residues of the hACE2. More specifically, strong electrostatic interactions (salt bridges) and hydrogen bonding were observed between R493 and R498 residues of the Omicron RBD with D30/E35 and D38 residues of the hACE2, respectively. Other mutated amino acids in the Omicron RBD, e.g. S496 and H505, also exhibited hydrogen bonding with the hACE2 receptor. The pi-stacking interaction was also observed between tyrosine residues (RBD-Tyr501: hACE2-Tyr41) in the complex, which contributes majorly to binding free energies suggesting this as one of the key interactions stabilizing the complex formation. The structural insights of RBD:hACE2 complex, their binding mode information and residue wise contributions to binding free energy provide insight on the increased transmissibility of Omicron and pave the way to design and optimize novel antiviral agents.

Our previous efforts have proved that modifications targeting the 150-cavity of influenza neuraminidase can achieve more potent and more selective inhibitors. In this work, four subseries of C5-NH

The current global health emergency in the form of the Coronavirus 2019 (COVID-19) pandemic has highlighted the need for fast, accurate, and efficient drug discovery pipelines. Traditional drug discovery projects relying on in vitro high-throughput screening (HTS) involve large investments and sophisticated experimental set-ups, affordable only to big biopharmaceutical companies. In this scenario, application of efficient state-of-the-art computational methods and modern artificial intelligence (AI)-based algorithms for rapid screening of repurposable chemical space [approved drugs and natural products (NPs) with proven pharmacokinetic profiles] to identify the initial leads is a powerful option to save resources and time. Structure-based drug repurposing is a popular in silico repurposing approach. In this review, we discuss traditional and modern AI-based computational methods and tools applied at various stages for structure-based drug discovery (SBDD) pipelines. Additionally, we highlight the role of generative models in generating molecules with scaffolds from repurposable chemical space.

Three types of chemical entities, namely, small organic molecules (organics), peptides, and biologics, are mainly used as drug candidates for the treatment of various diseases. Even though the peptide drugs are known since 1920 in association with the clinical use of insulin, only a limited number of peptides are currently used for therapeutics due to various disadvantages associated with them such as limited serum and blood stability, oral bioavailability, and permeability. Since, through chemical modifications and structure tuning, many of these limitations can be overcome, peptide-based drugs are gaining attention in pharmaceutical research. As of today, there are more than 60 peptide-based drugs approved by FDA, and over 150 peptides are in the advanced clinical studies. In this book chapter, the peptide-based lead compounds and drugs available for treating various viral diseases and their advantages and disadvantages when compared to small molecules drugs are discussed.

The discovery of small molecules that exhibit turn-on far-red or near-infrared (NIR) fluorescence upon DNA binding and understanding how they bind DNA are important for imaging and bioanalytical applications. Here we report the DNA-bound structure and the DNA binding mechanism of quinone cyanine dithiazole (QCy-DT), a recently reported AT-specific turn-on NIR fluorescent probe for double-stranded DNA. The nuclear magnetic resonance (NMR)-derived structure showed minor groove binding but no specific ligand-DNA interactions, consistent with an endothermic and entropy-driven binding mechanism deduced from isothermal titration calorimetry. Minor groove binding is typically fast because it minimally perturbs the DNA structure. However, QCy-DT exhibited unusually slow DNA binding. The cyanine-based probe is capable of

To address drug resistance to HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs), a series of novel diarylpyrimidine (DAPY) derivatives targeting “tolerant region I” and “tolerant region II” of the NNRTIs binding pocket (NNIBP) were designed utilizing a structure-guided scaffold-hopping strategy. The dihydrofuro[3,4-d]pyrimidine derivatives 13c2 and 13c4 proved to be exceptionally potent against a wide range of HIV-1 strains carrying single NNRTI-resistant mutations (EC(50) = 0.9–8.4 nM), which were remarkably superior to that of etravirine (ETV). Meanwhile, both compounds exhibited comparable activities with ETV toward the virus with double mutations F227L+V106A and K103N+Y181C. Furthermore, the most active compound 13c2 showed favorable pharmacokinetic properties with an oral bioavailability of 30.96% and a half-life of 11.1 h, which suggested that 13c2 is worth further investigation as a novel NNRTI to circumvent drug resistance.

Understanding the initial events of aggregation of amyloid β monomers to form β-sheet rich fibrils is useful for the development of therapeutics for Alzheimer's disease. In this context, the changes in energetics involved in the aggregation of helical amyloid β monomers into β-sheet rich dimers have been investigated using umbrella sampling simulations and density functional theory calculations. The results from umbrella sampling simulations for the free energy profile for the interconversion closely agree with the results of density functional theory calculations. The results reveal that helical peptides converted to β-sheet structures through coil-like conformations as intermediates that are mostly stabilized by intramolecular hydrogen bonds. The stabilization of intermediate structures could be a possible way to inhibit fibril formation. Mutations substantially decrease the height of the energy barrier for interconversion from α-helix to β-sheet structure when compared to that of the wild type, something that is attributed to an increase in the number of intramolecular hydrogen bonds between backbone atoms in the coil structures that correspond to a maximum value on the free energy surface. The reduction of the energy barrier leads to an enhancement of the rate of aggregation of amyloid β monomers upon introduction of various familial mutations, which is consistent with previous experimental reports.

Our previous efforts have proved that modifications targeting the 150-cavity of influenza neuraminidase can achieve more potent and more selective inhibitors. In this work, four subseries of C5-NH2 modified oseltamivir derivatives were designed and synthesized to explore every region inside the 150-cavity. Among them, compound 23d was exceptionally potent against the whole panel of Group-1 NAs with IC50 values ranging from 0.26 to 0.73 nM, being 15-53 times better than oseltamivir carboxylate (OSC) and 7-11 times better than zanamivir. In cellular assays, 23d showed more potent or equipotent antiviral activities against corresponding virus strains compared to OSC with no cytotoxicity. Furthermore, 23d exhibited high metabolic stability in human liver microsomes (HLM) and low inhibitory effect on main cytochrome P450 enzymes. Notably, 23d displayed favorable druggability in vivo and potent antiviral efficacy in the embryonated egg model and mice model. Overall, 23d appears to be a promising candidate for the treatment of influenza virus infection.

Tauopathies are a subclass of neurodegenerative diseases characterized by an accumulation of microtubule binding tau fibrils in brain regions. Diseases such as Alzheimer's (AD), chronic traumatic encephalopathy (CTE), Pick's disease (PiD), and corticobasal degeneration (CBD) belong to this subclass. Development of tracers which can visualize and discriminate between different tauopathies is of clinical importance in the diagnosis of various tauopathies. Currently, several tau tracers are available for in vivo imaging using a positron emission tomography (PET) technique. Among these tracers, PBB3 is reported to bind to various types of tau fibrils with comparable binding affinities. In contrast, tau tracer AV-1451 is reported to bind to specific types of tau fibrils (in particular to AD-associated and CTE) with higher binding affinity and only show nonspecific or weaker binding toward tau fibrils dominant with 3R isoforms (associated with PiD). The tau fibrils associated with different tauopathies can adopt different microstructures with different binding site microenvironments. By using detailed studies of the binding profiles of tau tracers for different types of tau fibrils, it may be possible to design tracers with high selectivity toward a specific tauopathy. The microstructures for the tau fibrils from patients with AD, PiD, and CTE have recently been demonstrated by cryogenic electron microscopy (cryo-EM) measurements allowing structure-based in silico simulations. In the present study, we have performed a multiscale computational study involving molecular docking, molecular dynamics, free energy calculations, and QM fragmentation calculations to understand the binding profiles of tau tracer AV-1451 and its potential use for diagnosis of AD, CTE, and PiD tauopathies. Our computational study reveals that different affinity binding sites exist for AV-1451 in the tau fibrils associated with different tauopathies. The binding affinity of this tracer toward different tau fibrils goes in this order: PiD > AD > CTE. The interaction energies for different tau fibril-tracer complexes using the QM fragmentation scheme also showed the same trend. However, by carrying out molecular dynamics simulations for the AD-derived tau fibrils in organic solvents, we found additional high affinity binding sites for AV-1451. The AV-1451 binding profile in these cryptic sites correctly explains the preferential binding of this tracer toward the AD fibrils when compared with the PiD fibrils. This study clearly demonstrates having a cryo-EM structure is still not sufficient for the structure-based tracer discovery for certain targets, as they may have "potential but hidden" high affinity binding sites, and we need additional strategies to identify them.

Three types of chemical entities, namely, small organic molecules (organics), peptides, and biologics, are mainly used as drug candidates for the treatment of various diseases. Even though the peptide drugs are known since 1920 in association with the clinical use of insulin, only a limited number of peptides are currently used for therapeutics due to various disadvantages associated with them such as limited serum and blood stability, oral bioavailability, and permeability. Since, through chemical modifications and structure tuning, many of these limitations can be overcome, peptide-based drugs are gaining attention in pharmaceutical research. As of today, there are more than 60 peptide-based drugs approved by FDA, and over 150 peptides are in the advanced clinical studies. In this book chapter, the peptide-based lead compounds and drugs available for treating various viral diseases and their advantages and disadvantages when compared to small molecules drugs are discussed.

We present here a new insight into the chemical reactivity and bioactivity of dopamine (DA) and gallic acid (GA) and their hyaluronic acid (HA) conjugates. Our data suggest that HA-GA scaffolds are superior to HA-DA, with higher oxidation kinetics, improved tissue adhesive properties, and radical scavenging ability with a lower pro-inflammatory response. This journal is publishedVersion

To assess a tumor biomarker like the cyclooxygenase-2 enzyme (COX-2), non-invasive imaging techniques are powerful tools. The (non-) linear optical properties of activatable fluorescent probes which are selectively bound to the biomarker can therefore be exploited. The here presented molecular modelling results based on multi-scale modelling techniques highlight the importance of the conformational versatility and of changes in the electronic interactions of such probes when they are embedded in water or in the COX-2 homodimer enzyme. The ANQ-IMC-6 probe, which combines the binding domain/scaffold of indomethacin (IMC) on COX-2 with the optical properties of acenaphtho[1,2–b]quinoxaline (ANQ), is found to be folded in the solvent and unfolded in the enzyme. A concerted movement of the probe and the protein is seen, while the rotational autocorrelation function exhibits also the intrinsic properties of the probe. Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) calculations are used to simulate the one-photon and two-photon absorption spectra along with the first hyperpolarizability. The transition has a local character in vacuum, but changes to a charge transfer one in the presence of the microenvironment of the enzyme. This is also visible through a change of the shape of the absorption spectrum, while at the same time the simulated signals of second harmonic generation experiments are strongly enhanced. The results of this work prove that an environment sensitive probe with an anchoring group and an optical active part can be constructed for use in absorption spectroscopy, without the need to revert to fluorescence experiments.

Previous efforts in our lab have led to the development of human immunodeficiency virus type 1 (HIV-1) non-nucleoside reverse transcriptase inhibitor (NNRTI) thiophene[3,2-d]pyrimidine compound 1 (K-5a2) with promising activity against wild-type and mutant HIV-1 strains. In this work, a series of novel diarylpyrimidines derivatives carrying a structurally diverse motif at the right wing of the lead K-5a2 was designed and synthesized as potential anti-HIV-1 agents. The results demonstrated that 8a yielded exceptionally potent activity against HIV-1 wild-type (50% effective concentration (EC50) = 3.30 nM) and mutant strain RES056 (EC50 = 22.6 nM) in MT-4 cells; in the reverse transcriptase inhibitory assay, 8a (half maximal inhibitory concentration (IC50) = 0.028 μM) was remarkably superior to that of K-5a2 (IC50 = 0.300 μM) and comparable to that of etravirine (ETR; IC50 = 0.011 μM). Notably, 8a exhibited better druggability than that of K-5a2, including significantly reduced CYP enzymatic inhibitory activity (IC50 > 50 μM), lower human ether-à-go-go related gene (hERG) inhibition (IC50 > 30 μM), and improved metabolic stability (short half-life, T1/2 = 77.5 min) in vitro. doi: 10.1021/acsinfecdis.0c00327 ispartof: Acs Infectious Diseases vol:6 issue:8 pages:2225-2234 ispartof: location:United States status: published

Summary Multiple lines of evidence indicate that amyloid beta (Aβ) peptide is responsible for the pathological devastation caused in Alzheimer's disease (AD). Aβ aggregation species predominantly contribute to multifaceted toxicity observed in neuronal cells including generation of reactive oxygen species (ROS), mitochondrial dysfunction, interfering with synaptic signaling, and activation of premature apoptosis. Herein, we report a natural product berberine-derived (Ber-D) multifunctional inhibitor to ameliorate in cellulo multifaceted toxicity of AD. The structural attributes of polyphenolic Ber-D have contributed to its efficient Cu chelation and arresting the redox cycle to prevent the generation of ROS and rescue biomacromolecules from oxidative damage. Ber-D inhibits metal-dependent and -independent Aβ aggregation, which is supported by in silico studies. Ber-D treatment averts mitochondrial dysfunction and corresponding neuronal toxicity contributing to premature apoptosis. These key multifunctional attributes make Ber-D a potential therapeutic candidate to ameliorate multifaceted Aβ toxicity in AD., Graphical Abstract, Highlights • Multifaceted Aβ toxicity activates premature apoptosis and neuronal death in AD • Berberine is modified to nontoxic and polyphenolic multifunctional derivative (Ber-D) • Ber-D modulates Aβ aggregation, metal toxicity, ROS, oxidative and biomolecules damage • Study revealed mitochondrial protection and downregulation of apoptosis markers, Drugs; Applied Chemistry; Molecular Biology

Correction for ‘Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7’ by N. Arul Murugan et al., Phys. Chem. Chem. Phys., 2022, 24, 20371–20380, https://doi.org/10.1039/D2CP00469K.

Parkinson's disease (PD) is a serious neurodegenerative disease and is characterized by abnormal α-synuclein (α-syn) accumulation in Lewy bodies (LB) and Lewy neurites (LN), which makes α-syn an important imaging target for PD. An imaging probe that quantifies fibrillar α-syn can enhance the clinical diagnosis of PD and can also be used to evaluate the efficacy of therapeutics aimed at reducing the abnormal aggregation of the α-syn fibril in the brain. In this paper, we study the binding profile of fibrillar α-syn with a fluorescent probe 4-(dicyanovinyl)julolidine (DCVJ), which is being explored for identifying α-syn imaging agents. A multiscale simulation workflow including molecular docking, molecular dynamics, metadynamics, and QM/MM calculations was implemented. We find that DCVJ can bind to multiple sites of α-syn which are located either at the surface or in the core. Free energy calculations using implicit solvent models reveal that the most favorable binding mode for DCVJ is associated with the core binding site and is further confirmed by metadyamics simulation. Besides, a dynamic binding pathway is discovered, which reveals that DCVJ binds gradually into the core of the fibril passing through several intermediate states. The conformational arrest of the dicyano vinyl group in the fibrillar environment could explain the reason behind the fibril-specific fluorescence of DCVJ. Furthermore, based on hybrid QM/MM calculations, the molecular geometry of the dicyano vinyl group is found to be environment specific which explains why DCVJ serves as a staining agent for such fibrillar-like environments. Our results could be helpful for elucidating the binding mechanism of imaging tracers with the fibrillar form of α-syn and explain their fibrillar-specific optical properties, a knowledge that in turn can be used to guide the design and development of compounds with higher affinity and selectivity for α-syn using structure-based strategies.

International audience; A series of novel BF2 complexes, bearing a five-membered heterocyclic ring (with X = NMe, O, and S), were synthesized and characterized with a focus on the influence of atom exchange on the photophysical properties of both unsubstituted and dimethylamino derivatives. The experimental results show that the optical spectra substantially differ in both sets of dyes. In particular, the dimethylamino series are more strongly affected by heteroatom substitution, i.e., the insertion of X = O or X = S in lieu of X = NMe causes substantial bathochromic shifts of the absorption and emission bands, as well as marked changes in their topologies. In contrast, the optical spectra of the unsubstituted compounds undergo only relatively small red-shifts, and no variation of band shapes is observed. Moreover, the measured absorption spectra of the unsubstituted compounds bearing X = NMe and X = O are almost identical. Interestingly, the fluorescence yields of the dimethylamino derivatives are much larger (up to one order of magnitude) than those of the corresponding unsubstituted compounds. The experimental analyses are supported by state-of-the-art quantum chemistry calculations, which satisfactorily reproduced the experimental trends and provided further insights into the observed optical signatures

Using the recently reported cryo-EM structure for the tau fibril [ Fitzpatrick et al. (2017) Nature 547, 185-190 ], which is a potential target concerning Alzheimer's disease, we present the first molecular modeling studies on its interaction with various positron emission tomography (PET) tracers. Experimentally, based on the binding assay studies, at least three different high-affinity binding sites have been reported for tracers in the tau fibril. Herein, through integrated modeling using molecular docking, molecular dynamics, and binding free energy calculations, we provide insight into the binding patterns of various tracers to the tau fibril. We suggest that there are four different high-affinity binding sites available for many of the studied tracers showing varying binding affinity to different binding sites. Thus, PBB3 binds most strongly to site 4, and interestingly, this site is not a preferable site for any other tracers. For THK5351, our data show that it strongly binds to sites 3 and 1, the former one being more preferable. We also find that MK6240 and T807 bind to site 1 specifically. The modeling data also give some insight into whether a tracer bound to a specific site can be replaced by others or not. For example, the displacement of T807 by PBB3 as reported experimentally can also be explained and attributed to the larger binding affinity of the latter compound in all binding sites. The binding free energy results explain very well the small binding affinity of THK523 compared to all the aryl quinoline moieties containing THK tracers. The ability of certain tau tracers, like FDDNP and THK523, to bind to amyloid fibrils has also been investigated. Furthermore, such off-target interaction of tau tracers with amyloid beta fibrils has been validated using a quantum mechanical fragmentation approach.

The thrombin-binding aptamer (TBA), which shows anticoagulant properties, is one of the most studied G-quadruplex-forming aptamers. In this study, we investigated the impact of different chemical modifications such as a three-carbon spacer (spacer-C3), unlocked nucleic acid (UNA) and 3′-amino-modified UNA (amino-UNA) on the structural dynamics and stability of TBA. All three modifications were incorporated at three different loop positions (T3, T7, T12) of the TBA G-quadruplex structure to result in a series of TBA variants and their stability was studied by thermal denaturation; folding was studied by circular dichroism spectroscopy and thrombin clotting time. The results showed that spacer-C3 introduction at the T7 loop position (TBA-SP7) significantly improved stability and thrombin clotting time while maintaining a similar binding affinity as TBA to thrombin. Detailed molecular modelling experiments provided novel insights into the experimental observations, further supporting the efficacy of TBA-SP7. The results of this study could provide valuable information for future designs of TBA analogues with superior thrombin inhibition properties.

Sirtuins are a family of nicotinamide adenine dinucleotide (NAD+)-dependent enzymes, which undergo robust deacetylase activity, resulting in the production of nicotinamide. It is well known that nicotinamide, which is one of the products, can also act as an inhibitor for further deacetylation process by forming NAD+ again. Hence, the removal of nicotinamide from sirtuins is a demanding process, and the mechanistic understanding of the process remains elusive. In this investigation, we have made an attempt to unravel the unbinding pathways of nicotinamide from SIRT1, SIRT2, and SIRT3 (SIRT1-3) using Random Acceleration Molecular Dynamics (RAMD) Simulations, and we have successfully identified various unbinding channels. The selectivity of the egression channel is determined by using a thorough analysis of the frequency of egression trajectories. Similarly, various inhibitors have been docked with the active sites of SIRT1-3, and their egression pathways have been investigated to understand whether they follow the same egression pathway as that of nicotinamide. The residues that are responsible for the unbinding pathways have been determined from the analysis of RAMD trajectories. From these results, it is clear that phenylalanine and histidine residues play major roles in the egression of inhibitors. Additionally, the key residues Leu, Pro, Met, Phe, Tyr, and Ile are found to control the release by acting as gateway residues. The role of these residues from different egression channels has been studied by carrying out mutations with alanine residue. This is the first report on sirtuins, which demonstrates the novel unbinding pathways for nicotinamide/inhibitors. This work provides new insights for developing more promising SIRT1-3 inhibitors.

In recent times, our healthcare system is being challenged by many drug-resistant microorganisms and ageing-associated diseases for which we do not have any drugs or drugs with poor therapeutic profile. With pharmaceutical technological advancements, increasing computational power and growth of related biomedical fields, there have been dramatic increase in the number of drugs approved in general, but still way behind in drug discovery for certain class of diseases. Now, we have access to bigger genomics database, better biophysical methods, and knowledge about chemical space with which we should be able to easily explore and predict synthetically feasible compounds for the lead optimization process. In this chapter, we discuss the limitations and highlights of currently available computational methods used for protein–ligand binding affinities estimation and this includes force-field, ab initio electronic structure theory and machine learning approaches. Since the electronic structure-based approach cannot be applied to systems of larger length scale, the free energy methods based on this employ certain approximations, and these have been discussed in detail in this chapter. Recently, the methods based on electronic structure theory and machine learning approaches also are successfully being used to compute protein–ligand binding affinities and other pharmacokinetic and pharmacodynamic properties and so have greater potential to take forward computer-aided drug discovery to newer heights.

Purpose Several tracers have been designed for tracking the abnormal accumulation of tau pathology in vivo. Recently, concerns have been raised about the sources of off-target binding for these tracers; inconclusive data propose binding for some tracers to monoamine oxidase B (MAO-B). Methods Molecular docking and dynamics simulations were used to estimate the affinity and free energy for the binding of several tau tracers (FDDNP, THK523, THK5105, THK5317, THK5351, T807 [aka AV-1451, flortaucipir], T808, PBB3, RO-948, MK-6240, JNJ-311 and PI-2620) to MAO-B. These values were then compared with those for safinamide (MAO-B inhibitor). PET imaging was used with the tau tracer [18F]THK5317 and the MAO-B tracer [11C]DED in five patients with Alzheimer’s disease to investigate the MAO-B binding component of this first generation tau tracer in vivo. Results The computational modelling studies identified a binding site for all the tau tracers on MAO-B; this was the same site as that for safinamide. The binding affinity and free energy of binding for the tau tracers to MAO-B was substantial and in a similar range to those for safinamide. The most recently developed tau tracers MK-6240, JNJ-311 and PI-2620 appeared, in silico, to have the lowest relative affinity for MAO-B. The in vivo investigations found that the regional distribution of binding for [18F]THK5317 was different from that for [11C]DED, although areas of suspected off-target [18F]THK5317 binding were detected. The binding relationship between [18F]THK5317 and [11C]DED depended on the availability of the MAO-B enzyme. Conclusions The developed tau tracers show in silico and in vivo evidence of cross-interaction with MAO-B; the MAO-B component of the tracer binding was dependent on the regional concentration of the enzyme. Electronic supplementary material The online version of this article (10.1007/s00259-019-04305-8) contains supplementary material, which is available to authorized users.

We study the linear and nonlinear optical properties of a well-known acid-base indicator, bromophenol blue (BPB), in aqueous solution by employing static and integrated approaches. In the static approach, optical properties have been calculated using time-dependent density functional theory (TD-DFT) on the fully relaxed geometries of the neutral and different unprotonated forms of BPB. Moreover, both closed and open forms of BPB were considered. In the integrated approach, the optical properties have been computed over many snapshots extracted from molecular dynamics simulation using a hybrid time-dependent density functional theory/molecular mechanics approach. The static approach suggests closed neutral ⇒ anionic interconversion as the dominant mechanism for the red shift in the absorption spectra of BPB due to a change from acidic to basic pH. It is found by employing an integrated approach that the two interconversions, namely open neutral ⇒ anionic and open neutral ⇒ dianionic, can contribute to the pH-dependent shift in the absorption spectra of BPB. Even though both static and integrated approaches reproduce the pH-dependent red shift in the absorption spectra of BPB, the latter one is suitable to determine both the spectra and spectral broadening. Finally, the computed static first hyperpolarizability for various protonated and deprotonated forms of BPB reveals that this molecule can be used as a nonlinear optical probe for pH sensing in addition to its highly exploited use as an optical probe.

The design of molecular probes and tracer molecules with specificity toward amyloid beta (Aβ) fibrils is of paramount importance for the selective diagnosis of Alzheimer's disease. This requires a detailed understanding of the binding sites in amyloid targets, their number, and their binding mechanism for various tracer molecules. We adopt an integrated approach including molecular docking, molecular dynamics, and generalized Born-based free energy calculations to investigate site-specific interactions of different amyloid binding molecules. Our study reproduces the experimental results on the relative binding affinity of the tracers and amyloid binders and explains the feature of "multiple binding sites" in amyloid targets as probed by competition binding experiments. A major outcome of this study is that it is the core sites of the Aβ fibrils that are responsible for the experimentally reported binding profiles of tracers in amyloid targets rather than the surface sites that received much focus in earlier investigations.

A new series of piperidine derivatives namely 3,5-dimethyl-2,6-diaryl-1-propionyl-piperidin-4-ones were synthesized. The structure of the compound, 3, 5-dimethyl-2,6-phenyl-1-propionyl-piperidin-4-one was determined unambiguously using NMR, IR and single crystal X-ray diffraction techniques. DFT and its time dependent version based calculations have been carried out to analyze its ground state electronic structure and to interpret the experimental spectroscopic data. The current study also explores into the spectroscopic property of these piperidine-4-one compounds and proposes these molecules as potential sunscreens and UV filters.

Understanding the mechanism of fibrillization of amyloid forming peptides could be useful for the development of therapeutics for Alzheimer's disease (AD). Taking this standpoint, we have explored in this work the free energy profile for the interconversion of monomeric and dimeric forms of amyloid forming peptides into different secondary structures namely beta-sheet, helix, and random coil in aqueous solution using umbrella sampling simulations and density functional theory calculations. We show that the helical structures of amyloid peptides can form β sheet rich aggregates through random coil conformations in aqueous condition. Recent experiments ( Chem. Eur. J. 2018, 24, 3397-3402 and ACS Appl. Mater. Interfaces 2017, 9, 21116-21123) show that molybdenum disulfide nanosurface and nanoparticles can reduce the fibrillization process of amyloid beta peptides. We have unravelled the free energy profile for the interconversion of helical forms of amyloid forming peptides into beta-sheet and random coil in the presence of a two-dimensional nanosurface of MoS2. Results indicate that the monomer and dimeric forms of the peptides adopt the random coil conformation in the presence of MoS2 while the helical form is preferable for the monomeric form and that the beta-sheet and helix forms are the preferable forms for dimers in aqueous solution. This is due to strong interaction with MoS2 and intramolecular hydrogen bonds of random coil conformation. The stabilization of random coil conformation does not lead to a β sheet like secondary structure for the aggregate. Thus, the confinement of MoS2 promotes deaggregation of amyloid beta peptides rather than aggregation, something that could be useful for the development of therapeutics for AD.

The widespread interest in phase recognition of lipid membranes has led to the use of different optical techniques to enable differentiation of healthy and not fully functional cells. In this work, we show how the combination of different (non)linear optical methods such as One Photon Absorption (OPA), Two Photon Absorption (TPA) and Second Harmonic Generation (SHG) as well as the study of the fluorescence decay time leads to an enhanced screening of membrane phases using a fluorescent dioctadecyltetramethylindocarbocyanine (DiD) probe. In the current study we consider the pure liquid disordered phases of DOPC (room temperature) and DPPC (323 K), the solid gel phase of DPPC (298 K), and the liquid ordered phase of a 2:1 binary mixture of Sphingomyelin and Cholesterol. By means of extensive hybrid quantum mechanics molecular mechanics calculations and based upon the (non) linear absorption of the embedded probes, it is found that DiD can be used to identify the lipid bilayer phase. The joint TPA and SHG as well as fluorescence analyses qualifies DiD as versatile probe for phase recognition. In particular, the SHG data obtained by means of HyperRayleigh Scattering and by Electric Field Induced Second Harmonic Generation reveal differences in polarization of the probe in the different environments. The TPA results finally confirm the particular location of the probe in between the polar head group region of the 2:1 SM:Chol mixture in the liquid ordered phase., 33 pages, 4 figures

Optical imaging of amyloid fibrils serves as a cost-effective route for the diagnosis of Alzheimer-like conformational diseases. However{,} the challenge here is to optimize the binding affinity and photophysical properties of the optical imaging agents in a way specific to certain types of amyloids. In a few occasions it is shown that novel optical imaging agents can be designed to bind to a particular type of amyloid fibril with larger binding affinity and specificity. There is also a recent report on photoluminescent polythiophenes which display photophysical properties that can be used to distinguish the variants or subtypes of amyloids (J. Rasmussen et al.{,} Proc. Natl. Acad. Sci. U. S. A.{,} 2017{,} 114(49){,} 13018–13023). Based on a multiscale modeling approach{,} here{,} we report on the complementary aspect that the photophysical properties of a benzothiazole based optical probe (referred to as BTA-3) can be specific to the binding sites in the same amyloid fibrils and we attribute this to its varying electronic structure in different sites. As reported experimentally from competitive binding assay studies for many amyloid staining molecules and tracers{,} we also show multiple binding sites in amyloid fibrils for this probe. In particular{,} BTA-3 displayed a red-shift in its low-frequency absorption band only in site-4{,} a surface site of amyloid fibrils when compared to the spectra in water solvent. In the remaining sites{,} it exhibited a less significant blue shift for the same absorption band.

The extraordinary ligand binding properties of albumin makes it a key player in the pharmacokinetics and pharmacodynamics of many vital drugs. Albumin is highly susceptible for nonenzymatic glycation mediated structural modifications, and there is a need to determine structural and functional impact of specific AGEs modifications. The present study was aimed toward determining the AGE mediated structure and function changes, primarily looking into the effect on binding affinity of drugs in the two major drug binding sites of albumin. The impact of the two most predominant AGEs modifications, i.e., carboxyethyllysine (CEL) and argpyrimidine (Arg-P), was studied on the basis of the combination of in vitro and in silico experiments. In vitro studies were carried out by AGEs modification of bovine serum albumin (BSA) for the formation of Arg-P and CEL followed by drug interaction studies. In silico studies involved molecular dynamics (MD) simulations and docking studies for native and AGEs modified BSAs. In particular the side chain modification was specifically carried out for the residues in the drug binding sites, i.e., Arg-194, Arg-196, Arg-198, and Arg-217, and Lys-204 (site I) and Arg-409 and Lys-413 (site II). The equilibrated structures of native BSA (n-BSA) and glycated BSA (G-BSA) as obtained from MD were used for drug binding studies using molecular docking approach. It was evident from the results of both in vitro and in silico drug interaction studies that AGEs modification results in the reduced drug binding affinity for tolbutamide (TLB) and ibuprofen (IBP) in sites I and II. Moreover, the AGEs modification mediated conformational changes resulted in the shallow binding pockets with reduced accessibility for drugs.

Understanding the interaction of gold nanoclusters with proteins has important ramifications in various fields. We present a study of the interaction between gold nanoclusters and lysozyme investigated using classical molecular dynamics and center-of-mass pulling simulations. The results reveal that the gold nanoclusters induce significant structural changes in lysozyme. Because the internal motions of lysozyme are related to its function, the changes in these internal motions have been quantified using principal component analysis of the molecular dynamics trajectories. The internal motions of lysozyme that are important for its function have been altered because of the interaction with the gold nanocluster. We have also explored how these induced changes in the lysozyme structure affect specific optical properties of the gold nanocluster using the complex polarization propagator method within the time-dependent density functional theory framework, which is of relevance for studies of the optical imaging of lysozyme using gold nanoclusters as molecular probes.

This study demonstrates that a hybrid density functional theory/molecular mechanics approach can be successfully combined with time-dependent wavepacket approach to predict the shape of optical bands for molecules in solutions, including vibrational fine structure. A key step in this treatment is the estimation of the inhomogeneous broadening based on the hybrid approach, where the polarization between solute and atomically decomposed solvent is taken into account in a self-consistent manner. The potential of this approach is shown by predicting optical absorption bands for three heterocyclic ketoimine difluoroborates in solution.