Strong and weak hydrogen bonds in the protein-ligand interface

Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain

Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.

Molecular modelling, docking and network analysis of phytochemicals from Haritaki churna: role of protein cross-talks for their action

Phytochemicals are bioactive agents present in medicinal plants with therapeutic values. Phytochemicals isolated from plants target multiple cellular processes. In the current work, we have used fractionation techniques to identify 13 bioactive polyphenols in ayurvedic medicine Haritaki Churna. Employing the advanced spectroscopic and fractionation, structure of bioactive polyphenols was determined. Blasting the phytochemical structure allow us to identify a total of 469 protein targets from Drug bank and Binding DB. Phytochemicals with their protein targets from Drug bank was used to create a phytochemical-protein network comprising of 394 nodes and 1023 edges. It highlights the extensive cross-talk between protein target corresponding to different phytochemicals. Analysis of protein targets from Binding data bank gives a network comprised of 143 nodes and 275 edges. Taking the data together from Drug bank and binding data, seven most prominent drug targets (HSP90AA1, c-Src kinase, EGFR, Akt1, EGFR, AR, and ESR-α) were found to be target of the phytochemicals. Molecular modelling and docking experiment indicate that phytochemicals are fitting nicely into active site of the target proteins. The binding energy of the phytochemicals were better than the inhibitors of these protein targets. The strength and stability of the protein ligand complexes were further confirmed using molecular dynamic simulation studies. Further, the ADMET profiles of phytochemicals extracted from HCAE suggests that they can be potential drug targets. The phytochemical cross-talk was further proven by choosing c-Src as a model. HCAE down regulated c-Src and its downstream protein targets such as Akt1, cyclin D1 and vimentin. Hence, network analysis followed by molecular docking, molecular dynamics simulation and in-vitro studies clearly highlight the role of protein network and subsequent selection of drug candidate based on network pharmacology.Communicated by Ramaswamy H. Sarma.

Mycosporine-Like Amino Acids as a Potential Inhibitor of Tyrosinase-Related Protein 1: Computational Screening, Pharmacokinetics, and Molecular Dynamics Simulation

Melanin is the major pigment responsible for the coloring of mammalian skin, hair, and eyes to defend against ultraviolet radiation. However, excessive melanin production has resulted in numerous types of hyperpigmentation disorders. Tyrosinase-related protein 1 (TYRP1) is a transmembrane glycoprotein enzyme found in many organisms, including humans, that plays an important role in melanogenesis. Thus, controlling the enzyme activity of TYRP1 with tyrosinase inhibitors is a vital step in the treatment of hyperpigmentation problems in humans. In the present investigation, virtual screening, pharmacokinetics, drug docking, and molecular dynamics (MD) simulation were used to find the most potent drug as an inhibitor of TYRP1 to effectively treat hyperpigmentation disorder. The 3D structure of TYRP1 was retrieved from the Protein Data Bank (PDB) database (PDB ID: 5M8M) and validated by the Ramachandran plot. Pharmacokinetics and drug-likeness showed that mycosporine 2 glycine (M2G) and shinorine (SHI) were the best compounds over other ligands in the same (P-1) structural pose. However, MD simulations of the M2G showed the highest CDOCKER interaction energy (-45.182 kcal/mol) and binding affinity (-65.0529 kcal/mol) as compared to SHI and reference drugs. The molecular binding modes RMSD and RMSF plots have exhibited more relevance to the M2G ligand in comparison to other drug ligands. The bioactivity and ligand efficiency profiles revealed that M2G is the most effective compound as a TYRP1 inhibitor. Thus, M2G could be used as a most effective drug for developing valuable sunscreen products to cure hyperpigmentation-related diseases.

Computational identification of potential modulators of heme-regulated inhibitor (HRI) for pharmacological intervention against sickle cell disease

Sickle cell disease (SCD) poses a significant health challenge and therapeutic approaches often target fetal hemoglobin (HbF) to ameliorate symptoms. Hydroxyurea, a current therapeutic option for SCD, has shown efficacy in increasing HbF levels. However, concerns about myelosuppression and thrombocytopenia necessitate the exploration of alternative compounds. Heme-regulated inhibitor (HRI) presents a promising target for pharmacological intervention in SCD due to its association with HbF modulation. This study screened compounds for their potential inhibitory functions against HRI. Small-molecule compounds from 17 folkloric plants were subjected to in silico screening against HRI. Molecular docking was performed, and free binding energy calculations were determined using molecular mechanics with generalized born and surface area (MMGBSA). Lead compounds were subjected to molecular dynamics simulation at 100 ns. Computational quantum mechanical modeling of the lead compounds was subsequently performed. We further examined the pharmacodynamics, pharmacokinetic and physiological properties of the identified compounds. Five potential HRI inhibitors, including kaempferol-3-(2G-glucosyrutinoside), epigallocatechin gallate, tiliroside, myricetin-3-O-glucoside and cannabiscitrin, with respective docking scores of -16.0, -12.17, -11.37, -11.56 and 11.07 kcal/mol, were identified. The MMGBSA analysis of the complexes yielded free-binding energies of -69.76, -71.17, -60.44, -53.55 and -55 kcal/mol, respectively. The identified leads were stable within HRI binding pocket for the duration of the 100 ns simulation. The study identified five phytoligands with potential inhibitory effects on HRI. This finding holds promise for advancing SCD treatment strategies. However, additional preclinical analyses are warranted to validate the chemotherapeutic properties of the lead compounds.Communicated by Ramaswamy H. Sarma.

Elucidation of protein-ligand interactions by multiple trajectory analysis methods

The identification of interaction between protein and ligand including binding positions and strength plays a critical role in drug discovery. Molecular docking and molecular dynamics (MD) techniques have been widely applied to predict binding positions and binding affinity. However, there are few works that describe the systematic exploration of the MD trajectory evolution in this context, potentially leaving out important information. To address the problem, we build a framework, Moira (molecular dynamics trajectory analysis), which enables automating the whole process ranging from docking, MD simulations and various analyses as well as visualizations. We utilized Moira to analyze 400 MD simulations in terms of their geometric features (root mean square deviation and protein-ligand interaction profiler) and energetics (molecular mechanics Poisson-Boltzmann surface area) for these trajectories. Finally, we demonstrate the performance of different analysis techniques in distinguishing native poses among four poses.

Attenuation of Type IV pili activity by natural products

The virulence factor Type IV pili (T4P) are surface appendages used by the opportunistic pathogen Pseudomonas aeruginosa for twitching motility and adhesion in the environment and during infection. Additionally, the use of these appendages by P. aeruginosa for biofilm formation increases its virulence and drug resistance. Therefore, attenuation of the activity of T4P would be desirable to control P. aeruginosa infections. Here, a computational approach has been pursued to screen natural products that can be used for this purpose. PilB, the elongation ATPase of the T4P machinery in P. aeruginosa, has been selected as the target subunit and virtual screening of FDA-approved drugs has been conducted. Screening identified two natural compounds, ergoloid and irinotecan, as potential candidates for inhibiting this T4P-associated ATPase in P. aeruginosa. These candidate compounds underwent further rigorous evaluation through molecular dynamics (MD) simulations and then through in vitro twitching motility and biofilm inhibition assays. Notably, ergoloid emerged as a particularly promising candidate for weakening the T4P activity by inhibiting the elongation ATPases associated with T4P. This repurposing study paves the way for the timely discovery of antivirulence drugs as an alternative to classical antibiotic treatments to help combat infections caused by P. aeruginosa and related pathogens.Communicated by Ramaswamy H. Sarma.

Molecular characterization of the PhiKo endolysin from Thermus thermophilus HB27 bacteriophage phiKo and its cryptic lytic peptide RAP-29

Introduction

In the era of increasing bacterial resistance to antibiotics, new bactericidal substances are sought, and lysins derived from extremophilic organisms have the undoubted advantage of being stable under harsh environmental conditions. The PhiKo endolysin is derived from the phiKo bacteriophage infecting Gram-negative extremophilic bacterium Thermus thermophilus HB27. This enzyme shows similarity to two previously investigated thermostable type-2 amidases, the Ts2631 and Ph2119 from Thermus scotoductus bacteriophages, that revealed high lytic activity not only against thermophiles but also against Gram-negative mesophilic bacteria. Therefore, antibacterial potential of the PhiKo endolysin was investigated in the study presented here.

Methods

Enzyme activity was assessed using turbidity reduction assays (TRAs) and antibacterial tests. Differential scanning calorimetry was applied to evaluate protein stability. The Collection of Anti-Microbial Peptides (CAMP) and Antimicrobial Peptide Calculator and Predictor (APD3) were used to predict regions with antimicrobial potential in the PhiKo primary sequence. The minimum inhibitory concentration (MIC) of the RAP-29 synthetic peptide was determined against Gram-positive and Gram-negative selected strains, and mechanism of action was investigated with use of membrane potential sensitive fluorescent dye 3,3'-Dipropylthiacarbocyanine iodide (DiSC3(5)).

Results and discussion

The PhiKo endolysin is highly thermostable with melting temperature of 91.70°C. However, despite its lytic effect against such extremophiles as: T. thermophilus, Thermus flavus, Thermus parvatiensis, Thermus scotoductus, and Deinococcus radiodurans, PhiKo showed moderate antibacterial activity against mesophiles. Consequently, its protein sequence was searched for regions with potential antibacterial activity. A highly positively charged region was identified and synthetized (PhiKo105-133). The novel RAP-29 peptide lysed mesophilic strains of staphylococci and Gram-negative bacteria, reducing the number of cells by 3.7-7.1 log units and reaching the minimum inhibitory concentration values in the range of 2-31 μM. This peptide is unstructured in an aqueous solution but forms an α-helix in the presence of detergents. Moreover, it binds lipoteichoic acid and lipopolysaccharide, and causes depolarization of bacterial membranes. The RAP-29 peptide is a promising candidate for combating bacterial pathogens. The existence of this cryptic peptide testifies to a much wider panel of antimicrobial peptides than thought previously.

Computational identification of potential modulators of heme-regulated inhibitor (HRI) for pharmacological intervention against sickle cell disease

Background

Sickle cell disease (SCD) poses a significant health challenge and therapeutic approaches often target fetal hemoglobin (HbF) to ameliorate symptoms. Hydroxyurea, a current therapeutic option for SCD, has shown efficacy in increasing HbF levels. However, concerns about myelosuppression and thrombocytopenia necessitate the exploration of alternative compounds. Heme-regulated inhibitor (HRI) presents a promising target for pharmacological intervention in SCD due to its association with HbF modulation. This study systematically screened compounds for their potential inhibitory functions against HRI.

Methods

Small-molecule compounds from 17 plants commonly utilized in traditional SCD management were subjected to in silico screening against HRI. Molecular docking was performed, and free binding energy calculations were determined using molecular mechanics with generalized born and surface area (MMGBSA). The lead compounds were subjected to molecular dynamics simulation at 100 ns. Computational quantum mechanical modelling of the lead compounds was subsequently performed. We further examined the pharmacodynamics, pharmacokinetic and physiological properties of the identified compounds.

Results

Five potential HRI inhibitors, including kaempferol-3-(2G-glucosyrutinoside), epigallocatechin gallate, tiliroside, myricetin-3-O-glucoside, and cannabiscitrin, with respective docking scores of -16.0, -12.17, -11.37, -11.56 and 11.07 kcal/mol, were identified. The MMGBSA analysis of the complexes yielded free-binding energies of -69.76, -71.17, -60.44, 53.55, and - 55 kcal/mol, respectively. The identified leads were stable within HRI binding pocket for the duration of 100 ns simulation.

Conclusions

The study successfully identified five phytoligands with potential inhibitory effects on HRI, opening avenues for their use as modulators of HbF in SCD patients. This finding holds promise for advancing treatment strategies in SCD. However, additional preclinical analyses are warranted to validate the chemotherapeutic properties of the lead compounds.

Norsesquiterpenes from the Latex of Euphorbia dentata and Their Chemical Defense Mechanisms against Helicoverpa armigera

Euphorbia dentata (Euphorbiaceae), an invasive weed, is rarely eaten by herbivorous insects and could secrete a large amount of white latex, causing a serious threat to local natural vegetation, agricultural production and human health. In order to prevent this plant from causing more negative effects on humans, it is necessary to understand and utilize the chemical relationships between the latex of E. dentata and herbivorous insects. In this study, three new norsesquiterpenes (1-3), together with seven known analogues (4-10), were isolated and identified from the latex of E. dentata. All norsesquiterpenes (1-10) showed antifeedant and growth-inhibitory effects on H. armigera with varying levels, especially compounds 1 and 2. In addition, the action mechanisms of active compounds (1-3) were revealed by detoxifying enzyme (AchE, CarE, GST and MFO) activities and corresponding molecular docking analyses. Our findings provide a new idea for the development and utilization of the latex of E. dentata, as well as a potential application of norsesquiterpenes in botanical insecticides.

Al distribution and structural stability of H-BEA zeolites at different Si/Al ratios and temperatures: a first-principles study

Beta zeolites have been widely used in acid-catalyzed reactions because of their excellent properties. An in-depth study of the position, quantity, and distribution of beta zeolites substituted by Al is significant to understand the catalytic performance of the active site of zeolite catalysts. The distribution of Al in H-BEA and the structure of silanol nests in dealuminated BEA at different Si/Al ratios and synthesis temperatures were studied by the DFT method. T1, T2, T7, and T9 sites were chosen to be simulated. The synthesis temperature can change the distribution of Al and the proportion of T sites at different Si/Al ratios. The proportion of T7 and T9 is more than 70% at different Si/Al ratios of H-BEA and decreases with the synthesis temperature. T1 and T2 sites begin to appear when Si/Al < 20 and the proportion of T1 and T2 sites is less than 20%. When Si/Al < 8, the substitution energy of the AlSiAl structure, which has Si(2Al, 2Si) species, is obviously lower than that of the normal structure, which indicates that the Al-O-Si-O-Al species will appear in H-BEA. The Al(T7)Si(T5)Al(T9)Si(T5)Al(T7) and Al(T1)Si(T1)Al(T9) groups can not only stabilize H-BEA but also play an essential role in the formation of Si(2Al, 2Si) species. For dealuminated BEA zeolites, the silanol nest forms four hydrogen bonds through four silanols. The orientation of silanol groups in the silanol nest formed after dealumination at different T sites is different. The T7 and T9 sites in H-BEA are more likely to undergo dealumination. By contrast, the dealumination of the T1 and T2 sites is a challenge.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Piperidine-based natural products targeting Type IV pili antivirulence: A computational approach

Type IV (T4) pilus is among the virulence factors with a key role in serious bacterial diseases. Specifically, in Neisseria meningitidis and Pseudomonas aeruginosa, it determines pathogenicity and causes infection. Here, a computational approach has been pursued to find piperidine-based inhibitor molecules against the elongation ATPase of T4 pili in these two selected pathogens. Using the modeled structures of the PilF and PilB ATPases of N. meningitidis and P. aeruginosa, virtual library screening via molecular docking has returned inhibitor molecule candidates. The dynamics of the best three binders have further been investigated in detail via molecular dynamic simulations. Among these, ligands with COCONUT IDs CNP0030078 and CNP0051517 were found to have higher potential in the inhibition of ATPases based on molecular dynamic simulation analysis and biological activity information. The obtained results will guide future efforts in antivirulence drug development against T4 pili of N. meningitidis and P. aeruginosa.

Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses

The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.

Phenylpropanoid amides from Solanum rostratum and their phytotoxic activities against Arabidopsis thaliana

Introduction

Solanum rostratum, an annual malignant weed, has seriously damaged the ecological environment and biodiversity of invasion area. This alien plant gains a competitive advantage by producing some new phytotoxic substances to inhibit the growth of native plants, thus achieving successful invasion. However, the chemical structures, inhibitory functions and action mechanisms of phytotoxic substances of S. rostratum remain unclear.

Methods

In this study, to clarify the chemical structures of phytotoxic substances from S. rostratum, we isolated phenylpropanoid amides from the plant. Their structures were identified by comprehensive HR-ESIMS, NMR and ECD data. And the inhibitory functions of isolated phenylpropanoid amides on one model plant (Arabidopsis thaliana) were also investigated. In addition, the action mechanisms of active phenylpropanoid amides were revealed by antioxidant-related enzymes [Catalase (CAT), Peroxidase (POD), Superoxide dismutase (SOD)] activities and corresponding molecular docking analyses.

Results and Discussion

Phytochemical research on the whole plant of S. rostratum led to the isolation and identification of four new phenylpropanoid amides (1-4), together with two known analogues (5-6). All the compounds showed phytotoxic effects with varying levels on the seed germination and root elongation of one model plant (Arabidopsis thaliana), especially compound 2 and 4. Likewise, compounds 2 and 4 displayed potent inhibitory effects on antioxidant-related enzyme (POD). In addition, compounds 2 and 4 formed common conventional hydrogen bonds with residues Ala34 and Ser35 in POD revealed by molecular docking analyses. These findings not only helped to reveal the invasion mechanism of S. rostratum from the perspective of "novel weapons hypothesis", but also opened up new ways for the exploitation and utilization of S. rostratum.

fingeRNAt—A novel tool for high-throughput analysis of nucleic acid-ligand interactions

Computational methods play a pivotal role in drug discovery and are widely applied in virtual screening, structure optimization, and compound activity profiling. Over the last decades, almost all the attention in medicinal chemistry has been directed to protein-ligand binding, and computational tools have been created with this target in mind. With novel discoveries of functional RNAs and their possible applications, RNAs have gained considerable attention as potential drug targets. However, the availability of bioinformatics tools for nucleic acids is limited. Here, we introduce fingeRNAt—a software tool for detecting non-covalent interactions formed in complexes of nucleic acids with ligands. The program detects nine types of interactions: (i) hydrogen and (ii) halogen bonds, (iii) cation-anion, (iv) pi-cation, (v) pi-anion, (vi) pi-stacking, (vii) inorganic ion-mediated, (viii) water-mediated, and (ix) lipophilic interactions. However, the scope of detected interactions can be easily expanded using a simple plugin system. In addition, detected interactions can be visualized using the associated PyMOL plugin, which facilitates the analysis of medium-throughput molecular complexes. Interactions are also encoded and stored as a bioinformatics-friendly Structural Interaction Fingerprint (SIFt)—a binary string where the respective bit in the fingerprint is set to 1 if a particular interaction is present and to 0 otherwise. This output format, in turn, enables high-throughput analysis of interaction data using data analysis techniques. We present applications of fingeRNAt-generated interaction fingerprints for visual and computational analysis of RNA-ligand complexes, including analysis of interactions formed in experimentally determined RNA-small molecule ligand complexes deposited in the Protein Data Bank. We propose interaction fingerprint-based similarity as an alternative measure to RMSD to recapitulate complexes with similar interactions but different folding. We present an application of interaction fingerprints for the clustering of molecular complexes. This approach can be used to group ligands that form similar binding networks and thus have similar biological properties. The fingeRNAt software is freely available at https://github.com/n-szulc/fingeRNAt.

Investigation of antidiabetic properties of shikonin by targeting aldose reductase enzyme: In silico and in vitro studies

Diabetes is a complicated multifactorial disorder in which the patient generally observes polyphagia, polydipsia, and polyuria due to uncontrolled growth in blood sugar levels. For its management, the pharmaceutical industry is working day and night to find a better drug with no or least toxicity. That's why nowadays a more focused branch is to use herbal phytoconstituents for its prevention. Shikonin is a naphthoquinone natural dye that is isolated from the plants of the Boraginaceae family and has proven its role as an anti-cancer, anti-inflammatory, and anti-gonadotrophic agent. In our previous study, we have published its anti-diabetic action by inhibiting the enzyme protein tyrosine phosphatase 1B. In this study, we were more focused on finding out the role of Shikonin and its pharmacophores by inhibiting the action of aldose reductase (AR) enzyme. The study was conducted using pharmacophore modeling, molecular docking, and molecular dynamics simulation studies. The absorption, distribution, metabolism, excretion (ADME), and toxicity profile were also evaluated in this study. Along with all the computational biology parameters we also focused on the in vitro activity and kinetic study of inhibitory activity of Shikonin against aldose reductase.

Forchlorfenuron and Novel Analogs Cause Cytotoxic Effects in Untreated and Cisplatin-Resistant Malignant Mesothelioma-Derived Cells

Malignant mesothelioma (MM) is a currently incurable, aggressive cancer derived from mesothelial cells, most often resulting from asbestos exposure. The current first-line treatment in unresectable MM is cisplatin/pemetrexed, which shows very little long-term effectiveness, necessitating research for novel therapeutic interventions. The existing chemotherapies often act on the cytoskeleton, including actin filaments and microtubules, but recent advances indicate the ‘fourth’ form consisting of the family of septins, representing a novel target. The septin inhibitor forchlorfenuron (FCF) and FCF analogs inhibit MM cell growth in vitro, but at concentrations which are too high for clinical applications. Based on the reported requirement of the chloride group in the 2-position of the pyridine ring of FCF for MM cell growth inhibition and cytotoxicity, we systematically investigated the importance (cell growth-inhibiting capacity) of the halogen atoms fluorine, chlorine, bromine and iodine in the 2- or 3-position of the pyridine ring. The MM cell lines ZL55, MSTO-211H, and SPC212, and—as a control—immortalized Met-5A mesothelial cells were used. The potency of the various halogen substitutions in FCF was mostly correlated with the atom size (covalent radius); the small fluoride analogs showed the least effect, while the largest one (iodide) most strongly decreased the MTT signals, in particular in MM cells derived from epithelioid MM. In the latter, the strongest effects in vitro were exerted by the 2-iodo and, unexpectedly, the 2-trifluoromethyl (2-CF3) FCF analogs, which were further tested in vivo in mice. However, FCF-2-I and, more strongly, FCF-2-CF3 caused rapidly occurring strong symptoms of systemic toxicity at doses lower than those previously obtained with FCF. Thus, we investigated the effectiveness of FCF (and selected analogs) in vitro in MM cells which were first exposed to cisplatin. The slowly appearing population of cisplatin-resistant cells was still susceptible to the growth-inhibiting/cytotoxic effect of FCF and its analogs, indicating that cisplatin and FCF target non-converging pathways in MM cells. Thus, a combination therapy of cisplatin and FCF (analogs) might represent a new avenue for the treatment of repopulating chemo-resistant MM cells in this currently untreatable cancer.

Computational investigation of the alkaloids of Pilocarpus microphyllus species as phytopharmaceuticals for the inhibition of sterol 14α-demethylase protease of Trypanosoma cruzi

Trypanosoma cruzi is a protozoan transmitted by the insect Triatoma infestans, popularly known as kissing bug. This protozoan causes the Chagas disease, a Neglected Tropical Disease. This study aimed to investigate, through DFT method and B3LYP hybrid functional, the physicochemical, pharmacokinetic, and pharmacodynamic properties of the alkaloids present in the leaves of the species Pilocarpus microphyllus (jaborandi) as a potential inhibitory activity on the protease sterol 14α-demethylase of T. cruzi associated with the techniques of molecular docking, molecular dynamics, MM-PBSA and ADMET predictions. The molecules of isopilosine, epiisopiloturine, epiisopilosine, and pilosine showed up the lowest binding energies by molecular docking, good human intestinal absorption, low penetration in the blood-brain barrier, antiprotozoal and anticarcinogenic activities in ADMET studies. It has been observed a better binding affinity of the sterol 14α-demethylase protease with isopilosine in molecular dynamics and MM-PBSA studies, which indicates it as a potential drug candidate for Chagas disease.Communicated by Ramaswamy H. Sarma.

Identification of MRTX1133, a Noncovalent, Potent, and Selective KRASG12D Inhibitor

KRAS^G12D, the most common oncogenic KRAS mutation, is a promising target for the treatment of solid tumors. However, when compared to KRAS^G12C, selective inhibition of KRAS^G12D presents a significant challenge due to the requirement of inhibitors to bind KRAS^G12D with high enough affinity to obviate the need for covalent interactions with the mutant KRAS protein. Here, we report the discovery and characterization of the first noncovalent, potent, and selective KRAS^G12D inhibitor, MRTX1133, which was discovered through an extensive structure-based activity improvement and shown to be efficacious in a KRAS^G12D mutant xenograft mouse tumor model.

Structural basis of reactivation of oncogenic p53 mutants by a small molecule: methylene quinuclidinone (MQ)

In response to genotoxic stress, the tumor suppressor p53 acts as a transcription factor by regulating the expression of genes critical for cancer prevention. Mutations in the gene encoding p53 are associated with cancer development. PRIMA-1 and eprenetapopt (APR-246/PRIMA-1^MET) are small molecules that are converted into the biologically active compound, methylene quinuclidinone (MQ), shown to reactivate mutant p53 by binding covalently to cysteine residues. Here, we investigate the structural basis of mutant p53 reactivation by MQ based on a series of high-resolution crystal structures of cancer-related and wild-type p53 core domains bound to MQ in their free state and in complexes with their DNA response elements. Our data demonstrate that MQ binds to several cysteine residues located at the surface of the core domain. The structures reveal a large diversity in MQ interaction modes that stabilize p53 and its complexes with DNA, leading to a common global effect that is pertinent to the restoration of non-functional p53 proteins.

The tumor suppressor p53 is mutated in more than half of human cancers and the compound methylene quinuclidinone (MQ) was shown to reactivate p53 mutants by binding covalently to cysteine residues. Here, the authors present crystal structures of wild-type and cancer related p53 mutant core domains bound to MQ alone and in complex with their DNA response elements and observe that MQ is bound to several cysteine residues located at the surface of the core domain.

Computational details of molecular structure, spectroscopic properties, topological studies and SARS-Cov-2 enzyme molecular docking simulation of substituted triazolo pyrimidine thione heterocycles

Investigation the molecular structure of the system requires a detailed experience in dealing with theoretical computational guides to highlight its important role. Molecular structure of three heterocyclic compounds 8,10-diphenylpyrido[3,2-e][1,2,4]triazolo[4,3-c]pyrimidine-3(2H)-thione (HL), 8-phenyl-10-(p-tolyl)pyrido[3,2-e][1,2,4]triazolo[4,3-c]pyrimidine-3(2H)-thione (CH3L) and10-(4-nitrophenyl)-8-phenylpyrido[3,2-e][1,2,4]triazolo[4,3-c]pyrimidine-3(2H)-thione (NO2L) was studied at DFT/B3LYP/6-31G (d,p) level in ethanol solvent. Spectroscopic properties such Infrared (IR, 1H NMR and 13C NMR) and ultraviolet-visible (UV-VIS) analyses were computed. Some quantum and reactivity parameters (HOMO energy, LUMO energy, energy gap, ionization potential, electron affinity, chemical potential, global softness, lipophelicity) were studied, also molecular electrostatic potential (MEP) was performed to indicate the reactive nucleophilic and electrophilic sites. The effects of H-, CH3- and NO2- substituents on heterocyclic ligands were studied and it was found that the electron donation sites concerned with hydrogen and methyl substituents over nitro substituent. Topological analysis using reduced density gradient (RDG) was discussed in details. To predict the relevant antiviral activity of the reported heterocyclic compounds, molecular docking simulation was applied to the crystal structure of SARS-Cov-2 viral Mpro enzyme with 6WTT code and PLpro with 7JRN code. The enzymatic viral protein gives an image about the binding affinity between the target protein receptor and the heterocyclic ligands entitled. The hydrogen bonding interactions were evaluated from molecular docking with different strength for each ligand compound to discuss the efficiency of heterocyclic ligands toward viral inhibition.

Anomalous amide proton chemical shifts as signatures of hydrogen bonding to aromatic sidechains

Hydrogen bonding between an amide group and the p-π cloud of an aromatic ring was first identified in a protein in the 1980s. Subsequent surveys of high-resolution X-ray crystal structures found multiple instances, but their preponderance was determined to be infrequent. Hydrogen atoms participating in a hydrogen bond to the p-π cloud of an aromatic ring are expected to experience an upfield chemical shift arising from a shielding ring current shift. We surveyed the Biological Magnetic Resonance Data Bank for amide hydrogens exhibiting unusual shifts as well as corroborating nuclear Overhauser effects between the amide protons and ring protons. We found evidence that Trp residues are more likely to be involved in p-π hydrogen bonds than other aromatic amino acids, whereas His residues are more likely to be involved in in-plane hydrogen bonds, with a ring nitrogen acting as the hydrogen acceptor. The p-π hydrogen bonds may be more abundant than previously believed. The inclusion in NMR structure refinement protocols of shift effects in amide protons from aromatic sidechains, or explicit hydrogen bond restraints between amides and aromatic rings, could improve the local accuracy of sidechain orientations in solution NMR protein structures, but their impact on global accuracy is likely be limited.

Metabolic profiling, ADME pharmacokinetics, molecular docking studies and antibacterial potential of Phyllantus muellerianus leaves

Global increase in the level of antimicrobial resistance among bacterial pathogens has prompted the search for alternative treatment from medicinal plants. Phyllantus muellerianus leaves has been used traditionally against microorganisms of medical importance, hence the need to evaluate the pharmacological pathways and mode of actions using in vitro and in silico approaches. Clinical isolates of eight (8) microorganisms associated with urinary tract infections were obtained and identified using morphological and biochemical methods. Phyllantus muellerianus leaves were extracted and purified by solvent partitioning. Ethyl acetate fraction of PM had the highest yield and zone diameter range from 13.5 ± 1.00 to 28 ± 1.53 mm. The rate of protein leakage per time interval of Staphylococcus aureus increased from 9.29 μg/ml at 0 min to 17.43 μg/ml at 120 min while leakage in Candida albicans also increased from 8.57 μg/ml at 0 min to 70.43 μg/ml at 120 min. GCMS fingerprints, pharmacodynamics and pharmacokinetic studies revealed the active agent as quindoline, an azaindole and isotere of indoles having a binding energy of −9.1 kcal/mol. Analyses of the structural and atomic orientations of quindoline, and superimposition on ciprofloxacin, a common antibiotic revealed an interesting comparison, effecting a stronger binding affinity of Quindoline-HMG-CoA complex.

HIV-1-RT inhibition activity of Satureja spicigera (C.KOCH) BOISS. Aqueous extract and docking studies of phenolic compounds identified by RP-HPLC-DAD

AIDS is a global disease caused by HIV, affecting millions of people and causing death. The current limitations of antiretroviral therapy used in the therapy of HIV/AIDS have led to the need to search for new and effective drugs from natural products, especially plants. Herewith, using the present study, the detection of HIV-1-RT inhibition of aqueous extract of Satureja spicigera (C.KOCH) BOISS. was performed for the first time. Besides, total phenolic content (TPC), analysis of phenolic constituents by RP-HPLC-DAD and antioxidant capacity by DPPH and Ferric reducing antioxidant power (FRAP) methods were determined for the first time. In addition, molecular docking studies were carried out between HIV-1-RT and phenolic substances, the presence of which was determined in the aqueous extract, for the determination of the phenolics that may be responsible for HIV-1-RT activity. HIV-1-RT inhibition was defined as IC₅₀ : 22.83 μg/ml. Benzoic acid, vanillin, rutin, and chlorogenic acid were present as main phenolics in quantities of 621.96, 505.87, 349.33, and 323.23 µg phenolic/g extract, respectively. Further, TPC, DPPH, and FRAP were calculated as in the order of 151.69 mg GAE/g extract, 23.77 µg/ml, and 445.7 µmol TE/g extract. Chlorogenic acid (-8.48 kcal/mol) was found to be the most effective ligand in docking studies, with a value close to positive standard nevirapine (-9.35 kcal/mol). Hereby, although the aqueous extract of S. spicigera can be used as a natural antioxidant, the crude extract or its phenolics have the potential to be used in the treatment of AIDS due to its high HIV-1-RT activity. PRACTICAL APPLICATIONS: In this study, anti-HIV-1-RT and antioxidant activity and total phenolic content of Satureja spicigera aqueous extract were determined. In addition, HPLC analysis of some phytochemicals and the activities of these phytochemicals against HIV-1-RT enzyme was determined by molecular docking studies. The results showed that the aqueous extract of S. spicigera and some of the phytochemicals it contains have the potential to be used as a natural product against HIV infection or in the treatment of AIDS.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Supramolecular cancer nanotheranostics

Among the many challenges in medicine, the treatment and cure of cancer remains an outstanding goal given the complexity and diversity of the disease. Nanotheranostics, the integration of therapy and diagnosis in nanoformulations, is the next generation of personalized medicine to meet the challenges in precise cancer diagnosis, rational management and effective therapy, aiming to significantly increase the survival rate and improve the life quality of cancer patients. Different from most conventional platforms with unsatisfactory theranostic capabilities, supramolecular cancer nanotheranostics have unparalleled advantages in early-stage diagnosis and personal therapy, showing promising potential in clinical translations and applications. In this review, we summarize the progress of supramolecular cancer nanotheranostics and provide guidance for designing new targeted supramolecular theranostic agents. Based on extensive state-of-the-art research, our review will provide the existing and new researchers a foundation from which to advance supramolecular cancer nanotheranostics and promote translationally clinical applications.

Carveol a Naturally-Derived Potent and Emerging Nrf2 Activator Protects Against Acetaminophen-Induced Hepatotoxicity

Acetaminophen (N-acetyl p-aminophenol or APAP) is used worldwide for its antipyretic and anti-inflammatory potential. However, APAP overdose sometimes causes severe liver damage. In this study, we elucidated the protective effects of carveol in liver injury, using molecular and in silico approaches. Male BALB/c mice were divided into two experimental cohorts, to identify the best dose and to further assess the role of carveol in the nuclear factor E2-related factor; nuclear factor erythroid 2; p45-related factor 2 (Nrf2) pathway. The results demonstrated that carveol significantly modulated the detrimental effects of APAP by boosting endogenous antioxidant mechanisms, such as nuclear translocation of Nrf2 gene, a master regulator of the downstream antioxidant machinery. Furthermore, an inhibitor of Nrf2, called all-trans retinoic acid (ATRA), was used, which exaggerated APAP toxicity, in addition to abrogating the protective effects of carveol; this effect was accompanied by overexpression of inflammatory mediators and liver = 2ltoxicity biomarkers. To further support our notion, we performed virtual docking of carveol with Nrf2-keap1 target, and the resultant drug-protein interactions validated the in vivo findings. Together, our findings suggest that carveol could activate the endogenous master antioxidant Nrf2, which further regulates the expression of downstream antioxidants, eventually ameliorating the APAP-induced inflammation and oxidative stress.

A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach

The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.

Inhibitory activity of quercetin, its metabolite, and standard antiviral drugs towards enzymes essential for SARS-CoV-2: the role of acid-base equilibria

The recently declared global pandemic of a new human coronavirus called SARS-CoV-2, which causes respiratory tract disease COVID-19, has reached worldwide resonance and global efforts are being made to look for possible cures. Sophisticated molecular docking software, as well as available protein sequence and structure information, offer the ability to test the inhibition of two important targets of SARS-CoV-2, furin (FUR) enzyme, and spike glycoprotein, or spike protein (SP), that are key to host cell adhesion and hijacking. The potential inhibitory effect and mechanism of action of acid-base forms of different antiviral drugs, dominant at physiological pH, chloroquine (CQ), hydroxychloroquine (HCQ), and cinanserin (CIN), which have been shown to be effective in the treatment of SARS-CoV-2 virus, is reported with the special emphasis on their relative abundances. On the other hand, the potential inhibitory effect of the dominant acid-base forms of quercetin (Q) and its oxidative metabolite 2-(3,4-dihydroxybenzoyl)-2,4,6-trihydroxy-3(2H) benzofuranone (BZF), which are constituents of traditional food products believed to exhibit antiviral effects, was also examined. The undertaken study includes the determination of the major energy contributions to the binding energy as well as in-depth analysis of amino acid residues at the active pocket and possible interactions. The approach that we propose here may be an additional strategy for combating the deadly virus by preventing the first step of the virus replication cycle. Preliminary research has shown that the investigated compounds exert an inhibitory effect against the SARS-CoV-2 furin enzyme and spiked glycoprotein through different acid-base forms. These investigations may be helpful in creating potential therapeutic agents in the fight against the SARS-CoV-2 virus. On the other hand, the results we predicted in this computational study may be the basis for new experimental in vitro and in vivo studies.

A Geometric Definition of Short to Medium Range Hydrogen-Mediated Interactions in Proteins

We present a method to rapidly identify hydrogen-mediated interactions in proteins (e.g., hydrogen bonds, hydrogen bonds, water-mediated hydrogen bonds, salt bridges, and aromatic π-hydrogen interactions) through heavy atom geometry alone, that is, without needing to explicitly determine hydrogen atom positions using either experimental or theoretical methods. By including specific real (or virtual) partner atoms as defined by the atom type of both the donor and acceptor heavy atoms, a set of unique angles can be rapidly calculated. By comparing the distance between the donor and the acceptor and these unique angles to the statistical preferences observed in the Protein Data Bank (PDB), we were able to identify a set of conserved geometries (15 for donor atoms and 7 for acceptor atoms) for hydrogen-mediated interactions in proteins. This set of identified interactions includes every polar atom type present in the Protein Data Bank except OE1 (glutamate/glutamine sidechain) and a clear geometric preference for the methionine sulfur atom (SD) to act as a hydrogen bond acceptor. This method could be readily applied to protein design efforts.

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Small molecule therapeutics for COVID-19: repurposing of inhaled furosemide

The novel coronavirus SARS-CoV-2 has become a global health concern. The morbidity and mortality of the potentially lethal infection caused by this virus arise from the initial viral infection and the subsequent host inflammatory response. The latter may lead to excessive release of pro-inflammatory cytokines, IL-6 and IL-8, as well as TNF-α ultimately culminating in hypercytokinemia (“cytokine storm”). To address this immuno-inflammatory pathogenesis, multiple clinical trials have been proposed to evaluate anti-inflammatory biologic therapies targeting specific cytokines. However, despite the obvious clinical utility of such biologics, their specific applicability to COVID-19 has multiple drawbacks, including they target only one of the multiple cytokines involved in COVID-19’s immunopathy. Therefore, we set out to identify a small molecule with broad-spectrum anti-inflammatory mechanism of action targeting multiple cytokines of innate immunity. In this study, a library of small molecules endogenous to the human body was assembled, subjected to in silico molecular docking simulations and a focused in vitro screen to identify anti-pro-inflammatory activity via interleukin inhibition. This has enabled us to identify the loop diuretic furosemide as a candidate molecule. To pre-clinically evaluate furosemide as a putative COVID-19 therapeutic, we studied its anti-inflammatory activity on RAW264.7, THP-1 and SIM-A9 cell lines stimulated by lipopolysaccharide (LPS). Upon treatment with furosemide, LPS-induced production of pro-inflammatory cytokines was reduced, indicating that furosemide suppresses the M1 polarization, including IL-6 and TNF-α release. In addition, we found that furosemide promotes the production of anti-inflammatory cytokine products (IL-1RA, arginase), indicating M2 polarization. Accordingly, we conclude that furosemide is a reasonably potent inhibitor of IL-6 and TNF-α that is also safe, inexpensive and well-studied. Our pre-clinical data suggest that it may be a candidate for repurposing as an inhaled therapy against COVID-19.

Opportunities and challenges in neutron crystallography

Neutron and X-ray crystallography are complementary to each other. While X-ray scattering is directly proportional to the number of electrons of an atom, neutrons interact with the atomic nuclei themselves. Neutron crystallography therefore provides an excellent alternative in determining the positions of hydrogens in a biological molecule. In particular, since highly polarized hydrogen atoms (H+) do not have electrons, they cannot be observed by X-rays. Neutron crystallography has its own limitations, mainly due to inherent low flux of neutrons sources, and as a consequence, the need for much larger crystals and for different data collection and analysis strategies. These technical challenges can however be overcome to yield crucial structural insights about protonation states in enzyme catalysis, ligand recognition, as well as the presence of unusual hydrogen bonds in proteins.

On the evolution of protein-adenine binding

How proteins evolved to recognize and bind their ligands is a key mystery in protein function evolution. To explore this mystery, we study how proteins bind adenine, an ancient fragment. We characterize physicochemical patterns of protein–adenine interactions and link these to proteins’ evolutionary origins. In conflict with previous findings, we see that all of adenine’s hydrogen donors and acceptors have been used to bind proteins, and that adenine binding is likely to have emerged multiple times in evolution. To identify adenine-binding sites of shared origin, we use “themes”: short amino acid segments suggested to constitute evolutionary building blocks. We detect specific themes that are engaged in adenine binding; the detection of these in a protein’s sequence might reveal its function.

Proteins’ interactions with ancient ligands may reveal how molecular recognition emerged and evolved. We explore how proteins recognize adenine: a planar rigid fragment found in the most common and ancient ligands. We have developed a computational pipeline that extracts protein–adenine complexes from the Protein Data Bank, structurally superimposes their adenine fragments, and detects the hydrogen bonds mediating the interaction. Our analysis extends the known motifs of protein–adenine interactions in the Watson–Crick edge of adenine and shows that all of adenine’s edges may contribute to molecular recognition. We further show that, on the proteins' side, binding is often mediated by specific amino acid segments (“themes”) that recur across different proteins, such that different proteins use the same themes when binding the same adenine-containing ligands. We identify numerous proteins that feature these themes and are thus likely to bind adenine-containing ligands. Our analysis suggests that adenine binding has emerged multiple times in evolution.

The Brain Exposure Efficiency (BEE) Score

The blood-brain barrier (BBB), composed of microvascular tight junctions and glial cell sheathing, selectively controls drug permeation into the central nervous system (CNS) by either passive diffusion or active transport. Computational techniques capable of predicting molecular brain penetration are important to neurological drug design. A novel prediction algorithm, termed the Brain Exposure Efficiency Score (BEE), is presented. BEE addresses the need to incorporate the role of trans-BBB influx and efflux active transporters by considering key brain penetrance parameters, namely, steady state unbound brain to plasma ratio of drug (K_p,uu) and dose normalized unbound concentration of drug in brain (C_u,b). BEE was devised using quantitative structure-activity relationships (QSARs) and molecular modeling studies on known transporter proteins and their ligands. The developed algorithms are provided as a user-friendly open source calculator to assist in optimizing a brain penetrance strategy during the early phases of small molecule molecular therapeutic design.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

The phytohormone forchlorfenuron decreases viability and proliferation of malignant mesothelioma cells in vitro and in vivo

Malignant mesothelioma (MM) is one of the most aggressive cancer types with a patient’s life expectancy of typically less than one year upon diagnosis. The urgency of finding novel therapeutic approaches to treat mesothelioma is evident. Here we tested the effect of the plant-growth regulator forchlorfenuron (FCF), an inhibitor of septin function(s) in mammalian cells, on the viability and proliferation of MM cell lines, as well as other tumor cell lines derived from lung, prostate, colon, ovary, cervix and breast. Exposure to FCF strongly inhibited proliferation of human and mouse (most efficiently epithelioid) MM cells and all other tumor cells in a concentration-dependent manner and led to cell cycle arrest and cell death. The role of septin 7 (SEPT7), a presumably essential target of FCF in MM cells was confirmed by an shRNA strategy. FCF was robustly inhibiting tumor cell growth in vitro at low micromolar (IC₅₀: ≈20-60µM) concentrations and more promisingly also in vivo. Initial experiments with FCF analogous revealed the importance of FCF’s chloride group for efficient cell growth inhibition. FCF’s rather low systemic toxicity might warrant for an extended search for other related and possibly more potent FCF analogues to target MM and putatively other septin-dependent tumors.

Crystallographic structure versus homology model: a case study of molecular dynamics simulation of human and zebrafish histone deacetylase 10

Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100 ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design.Communicated by Ramaswamy H. Sarma.

Cooperative nature of the sulfur centered hydrogen bond: investigation of (H2S)n (n = 2-4) clusters using an affordable yet accurate level of theory

Existing studies have shown that appreciably high level quantum chemical calculations are required to reproduce experimental energetic and geometric features of a H₂S dimer. This condition severely restricts any practical possibility of obtaining reliable results for higher order H₂S clusters. We have shown here that the binding energies calculated at the CCSD(T)/CBS level with counterpoise corrected geometries calculated at the MP2/aug-cc-pV(Q+d)Z level of theory excellently match with the experimental results for the H₂S dimer. Subsequently, the above mentioned levels of theory were used for trimers and tetramers. (H₂S)_n (n = 2-4) clusters were found to show cooperative strengthening of S-HS hydrogen bonds, which is clearly evident from the evolution of binding energies and hydrogen bond lengths, with increasing cluster size. Localized molecular orbital energy decomposition analyses have been carried out to understand how the contributions of various energy components modulate with the size of the clusters and what are their relative contributions towards the overall stabilization of the clusters. Natural bond orbital and atoms in molecules analyses were also carried out in order to look into the evolution of the electronic charge transfer and electron density topology with cluster size.

Computational Analysis of Non-covalent Interactions in Phycocyanin Subunit Interfaces

Protein-protein interactions are an important phenomenon in biological processes and functions. We used the manually curated non-redundant dataset of 118 phycocyanin interfaces to gain additional insight into this phenomenon using a robust inter-atomic non-covalent interaction analyzing tool PPCheck. Our observations indicate that there is a relatively high composition of hydrophobic residues at the interfaces. Most of the interface residues are clustered at the middle of the range which we call "standard-size" interfaces. Furthermore, the multiple interaction patterns founded in the present study indicate that more than half of the residues involved in these interactions participate in multiple and water-bridged hydrogen bonds. Thus, hydrogen bonds contribute maximally towards the stability of protein-protein complexes. The analysis shows that hydrogen bond energies contribute to about 88 % to the total energy and it also increases with interface size. Van der Waals (vdW) energy contributes to 9.3 %±1.7 % on average in these complexes. Moreover, there is about 1.9 %±1.5 % contribution by electrostatic energy. Nevertheless, the role by vdW and electrostatic energy could not be ignored in interface binding. Results show that the total binding energy is more for large phycocyanin interfaces. The normalized energy per residue was less than -16 kJ mol^-1 , while most of them have energy in the range from -6 to -14 kJ mol^-1 . The non-covalent interacting residues in these proteins were found to be highly conserved. Obtained results might contribute to the understanding of structural stability of this class of evolutionary essential proteins with increased practical application and future designs of novel protein-bioactive compound interactions.