Fortilin interacts with TGF-ß1 and prevents TGF-ß receptor activation.

Fortilin is a 172-amino acid multifunctional protein present in both intra- and extracellular spaces. Although fortilin binds and regulates various cellular proteins, the biological role of extracellular fortilin remains unknown. Here we report that fortilin specifically interacts with TGF-ß1 and prevents it from activating the TGF-ß1 signaling pathway. In a standard immunoprecipitation-western blot assay, fortilin co-immunoprecipitates TGF-ß1 and its isoforms. The modified ELISA assay shows that TGF-ß1 remains complexed with fortilin in human serum. Both bio-layer interferometry and surface plasmon resonance (SPR) reveal that fortilin directly bind TGF-ß1. The SPR analysis also reveals that fortilin and the TGF-ß receptor II (TGFßRII) compete for TGF-ß1. Both luciferase and secreted alkaline phosphatase reporter assays show that fortilin prevents TGF-ß1 from activating Smad3 binding to Smad-binding element. Fortilin inhibits the phosphorylation of Smad3 in both quantitative western blot assays and ELISA. Finally, fortilin inhibits TGFß-1-induced differentiation of C3H10T1/2 mesenchymal progenitor cells to smooth muscle cells. A computer-assisted virtual docking reveals that fortilin occupies the pocket of TGF-ß1 that is normally occupied by TGFßRII and that TGF-ß1 can bind either fortilin or TGFßRII at any given time. These data support the role of extracellular fortilin as a negative regulator of the TGF-ß1 signaling pathway.

Using computers to ESKAPE the antibiotic resistance crisis.

Since the discovery of penicillin, the development and use of antibiotics have promoted safe and effective control of bacterial infections. However, the number of antibiotic-resistance cases has been ever increasing over time. Thus, the drug discovery process demands fast, efficient and cost-effective alternative approaches for developing lead candidates with outstanding performance. Computational approaches are appealing techniques to develop lead candidates in an in silico fashion. In this review, we provide an overview of the implementation of current in silico state-of-the-art techniques, including machine learning (ML) and deep learning (DL), in drug discovery. We also discuss the development of quantum computing and its potential benefits for antibiotics research and current bottlenecks that limit computational drug discovery advancement.

Bioactive recombinant human oncostatin M for NMR-based screening in drug discovery.

Oncostatin M (OSM) is a pleiotropic, interleukin-6 family inflammatory cytokine that plays an important role in inflammatory diseases, including inflammatory bowel disease, rheumatoid arthritis, and cancer progression and metastasis. Recently, elevated OSM levels have been found in the serum of COVID-19 patients in intensive care units. Multiple anti-OSM therapeutics have been investigated, but to date no OSM small molecule inhibitors are clinically available. To pursue a high-throughput screening and structure-based drug discovery strategy to design a small molecule inhibitor of OSM, milligram quantities of highly pure, bioactive OSM are required. Here, we developed a reliable protocol to produce highly pure unlabeled and isotope enriched OSM from E. coli for biochemical and NMR studies. High yields (ca. 10 mg/L culture) were obtained in rich and minimal defined media cultures. Purified OSM was characterized by mass spectrometry and circular dichroism. The bioactivity was confirmed by induction of OSM/OSM receptor signaling through STAT3 phosphorylation in human breast cancer cells. Optimized buffer conditions yielded 1H, 15N HSQC NMR spectra with intense, well-dispersed peaks. Titration of 15N OSM with a small molecule inhibitor showed chemical shift perturbations for several key residues with a binding affinity of 12.2 ± 3.9 µM. These results demonstrate the value of bioactive recombinant human OSM for NMR-based small molecule screening.

Ribbon α-Conotoxin KTM Exhibits Potent Inhibition of Nicotinic Acetylcholine Receptors.

KTM is a 16 amino acid peptide with the sequence WCCSYPGCYWSSSKWC. Here, we present the nuclear magnetic resonance (NMR) structure and bioactivity of this rationally designed α-conotoxin (α-CTx) that demonstrates potent inhibition of rat α3ß2-nicotinic acetylcholine receptors (rα3ß2-nAChRs). Two bioassays were used to test the efficacy of KTM. First, a qualitative PC12 cell-based assay confirmed that KTM acts as a nAChR antagonist. Second, bioactivity evaluation by two-electrode voltage clamp electrophysiology was used to measure the inhibition of rα3ß2-nAChRs by KTM (IC50 = 0.19 ± 0.02 nM), and inhibition of the same nAChR isoform by α-CTx MII (IC50 = 0.35 ± 0.8 nM). The three-dimensional structure of KTM was determined by NMR spectroscopy, and the final set of 20 structures derived from 32 distance restraints, four dihedral angle constraints, and two disulfide bond constraints overlapped with a mean global backbone root-mean-square deviation (RMSD) of 1.7 ± 0.5 Å. The structure of KTM did not adopt the disulfide fold of α-CTx MII for which it was designed, but instead adopted a flexible ribbon backbone and disulfide connectivity of C2-C16 and C3-C8 with an estimated 12.5% α-helical content. In contrast, α-CTx MII, which has a native fold of C2-C8 and C3-C16, has an estimated 38.1% α-helical secondary structure. KTM is the first reported instance of a Framework I (CC-C-C) α-CTx with ribbon connectivity to display sub-nanomolar inhibitory potency of rα3ß2-nAChR subtypes.

Effects of Solvent Stabilization on Pharmaceutical Crystallization: Investigating Conformational Polymorphism of Probucol Using Combined Solid-State Density Functional Theory, Molecular Dynamics, and Terahertz Spectroscopy.

Solid-state density functional theory (DFT), molecular dynamics (MD), and terahertz (THz) spectroscopy were used to study the formation of enantiotropically related conformational Form I and Form II polymorphs of the pharmaceutical compound, probucol. DFT calculations were performed on the crystal systems to compare relative lattice energies and the solvent stabilization of the metastable Form II structure. The thermodynamics of solvent inclusion in the Form II·MeOH crystal system were determined from MD simulations, as was the favored conformation of molecular probucol in methanol and ethanol solutions. The findings from both solid-state DFT and MD calculations suggest that the preferred molecular orientations of the probucol molecule in solution and the probable inclusion of methanol in the crystal lattice during the crystallization process lead to the solvent selectivity of the probucol polymorph formation. The additional stabilization energy provided by the crystallization solvent facilitates the nucleation and growth of the Form II polymorph under conditions that favor this metastable crystal form over the thermodynamically stable Form I, despite the higher energy molecular and crystalline configurations of probucol Form II. We demonstrate the influence of solvent on the formation of pharmaceutical polymorphs and provide a molecular-level view of complex interactions leading to polymorphism using a combination of computational methods and THz spectral data.

Supramolecular Interactions of Group VI Metal Carbonyl Complexes: The Facilitating Role of 1,3-Bis( p-isocyanophenyl)urea.

An investigation of supramolecular phenomena involving zerovalent transition metal complexes was facilitated by the production of the ditopic isocyanide ligand 1,3-bis( p-isocyanophenyl)urea, which was synthesized via substoichiometric phosgenation of 4-isocyanophenylamine and used to coordinate group VI metal carbonyl fragments. The resulting binuclear organometallic complexes were observed to pack into ladder-like anisotropic arrays in the solid state. Crystallographic and computational evidence suggests that this packing motif can be attributed to a combination of intermolecular π-π and urea-π interactions. Similar to other N, N'-diarylureas bearing electron-withdrawing groups, 1,3-bis( p-isocyanophenyl)urea and the organometallic complexes prepared therefrom also exhibit an affinity toward anion binding in nonaqueous solution. Equilibrium constants ( K) for the formation of host-guest complexes between the organometallic derivatives of 1,3-bis( p-isocyanophenyl)urea and chloride, nitrate, and acetate anions exceed 103, 104, and 105 M-1, respectively.

SPIDR: small-molecule peptide-influenced drug repurposing.

BACKGROUND: Conventional de novo drug design is costly and time consuming, making it accessible to only the best resourced research organizations. An emergent approach to new drug development is drug repurposing, in which compounds that have already gone through some level of clinical testing are examined for efficacy against diseases divergent than their original application. Repurposing of existing drugs circumvents the time and considerable cost of early stages of drug development, and can be accelerated by using software to screen existing chemical databases to identify suitable drug candidates. RESULTS: Small-molecule Peptide-Influenced Drug Repurposing (SPIDR) was developed to identify small molecule drugs that target a specific receptor by exploring the conformational binding space of peptide ligands. SPIDR was tested using the potent and selective 16-amino acid peptide α-conotoxin MII ligand and the α3ß2-nicotinic acetylcholine receptor (nAChR) isoform. SPIDR incorporates a genetic algorithm-based, heuristic search procedure, which was used to explore the ligand binding domain of the α3ß2-nAChR isoform using a library consisting of 640,000 α-conotoxin MII peptide analogs. The peptides that exhibited the highest affinity for α3ß2-nAChR were used as models for a small-molecule structure similarity search of the PubChem Compound database. SPIDR incorporates the SimSearcher utility, which generates shape distribution signatures of molecules and employs multi-level K-means clustering to insure fast database queries. SPIDR identified non-peptide drugs with estimated binding affinities nearly double that of the native α-conotoxin MII peptide. CONCLUSIONS: SPIDR has been generalized and integrated into DockoMatic v 2.1. This software contains an intuitive graphical interface for peptide mutant screening workflow and facilitates mapping, clustering, and searching of local molecular databases, making DockoMatic a valuable tool for researchers in drug design and repurposing.

Fortilin binds IRE1α and prevents ER stress from signaling apoptotic cell death.

The endoplasmic reticulum, the cytoplasmic organelle that matures a massive amount of nascent secretory polypeptides, is particularly sensitive to stress. Endoplasmic reticulum stress causes unfolded proteins to populate the organelle, eliciting the unfolded protein response. During the unfolded protein response, GRP78-an endoplasmic reticulum master stress regulator-detaches from three endoplasmic reticulum stress sensors (IRE1α, PERK, and ATF6) and allows them to activate the apoptotic signaling pathway. Fortilin, a pro-survival molecule, is known to inhibit apoptosis by binding and inhibiting p53, but its role in endoplasmic reticulum stress-induced apoptosis remains unknown. Here, we report that fortilin directly interacts with the cytoplasmic domain of IRE1α, inhibits both kinase and endoribonuclease (RNase) activities of the stress sensor, and protects cells against apoptotic cell death at both cellular and whole animal levels. Our data support a role of fortilin in the unfolded protein response and its potential participation in human diseases caused by unfolded protein response.IRE1α is an ER stress sensor, whose activity induces apoptosis. Here, the authors report that fortilin, a pro-survival factor, with yet unknown roles in ER stress, interacts with active IRE1α, inhibits both its kinase end RNase activities, and protects cells from apoptosis both in vitro and in vivo.

Genetic Algorithm Managed Peptide Mutant Screening: Optimizing Peptide Ligands for Targeted Receptor Binding.

This study demonstrates the utility of genetic algorithms to search exceptionally large and otherwise intractable mutant libraries for sequences with optimal binding affinities for target receptors. The Genetic Algorithm Managed Peptide Mutant Screening (GAMPMS) program was used to search an α-conotoxin (α-CTx) MII mutant library of approximately 41 billion possible peptide sequences for those exhibiting the greatest binding affinity for the α3ß2-nicotinic acetylcholine receptor (nAChR) isoform. A series of top resulting peptide ligands with high sequence homology was obtained, with each mutant having an estimated ΔGbind approximately double that of the potent native α-CTx MII ligand. A consensus sequence from the top GAMPMS results was subjected to more rigorous binding free energy calculations by molecular dynamics and compared to α-CTx MII and other related variants for binding with α3ß2-nAChR. In this study, the efficiency of GAMPMS to substantially reduce the sample population size through evolutionary selection criteria to produce ligands with higher predicted binding affinity is demonstrated.

Quantitative Determination of Vinpocetine in Dietary Supplements.

Current United States regulatory policies allow for the addition of pharmacologically active substances in dietary supplements if derived from a botanical source. The inclusion of certain nootropic drugs, such as vinpocetine, in dietary supplements has recently come under scrutiny due to the lack of defined dosage parameters and yet unproven short- and long-term benefits and risks to human health. This study quantified the concentration of vinpocetine in several commercially available dietary supplements and found that a highly variable range of 0.6-5.1 mg/serving was present across the tested products, with most products providing no specification of vinpocetine concentrations.

Cyclopamine bioactivity by extraction method from Veratrum californicum.

Veratrum californicum, commonly referred to as corn lily or Californian false hellebore, grows in high mountain meadows and produces the steroidal alkaloid cyclopamine, a potent inhibitor of the Hedgehog (Hh) signaling pathway. The Hh pathway is a crucial regulator of many fundamental processes during vertebrate embryonic development. However, constitutive activation of the Hh pathway contributes to the progression of various cancers. In the present study, a direct correlation was made between the extraction efficiency for cyclopamine from root and rhizome by eight methods, and the associated biological activity in Shh-Light II cells using the Dual-Glo® Luciferase Assay System. Alkaloid recovery ranged from 0.39 to 8.03mg/g, with ethanol soak being determined to be the superior method for obtaining biologically active cyclopamine. Acidic ethanol and supercritical extractions yielded degraded or contaminated cyclopamine with lower antagonistic activity towards Hh signaling.

Computational exploration of a protein receptor binding space with student proposed peptide ligands.

Computational molecular docking is a fast and effective in silico method for the analysis of binding between a protein receptor model and a ligand. The visualization and manipulation of protein to ligand binding in three-dimensional space represents a powerful tool in the biochemistry curriculum to enhance student learning. The DockoMatic tutorial described herein provides a framework by which instructors can guide students through a drug screening exercise. Using receptor models derived from readily available protein crystal structures, docking programs have the ability to predict ligand binding properties, such as preferential binding orientations and binding affinities. The use of computational studies can significantly enhance complimentary wet chemical experimentation by providing insight into the important molecular interactions within the system of interest, as well as guide the design of new candidate ligands based on observed binding motifs and energetics. In this laboratory tutorial, the graphical user interface, DockoMatic, facilitates docking job submissions to the docking engine, AutoDock 4.2. The purpose of this exercise is to successfully dock a 17-amino acid peptide, α-conotoxin TxIA, to the acetylcholine binding protein from Aplysia californica-AChBP to determine the most stable binding configuration. Each student will then propose two specific amino acid substitutions of α-conotoxin TxIA to enhance peptide binding affinity, create the mutant in DockoMatic, and perform docking calculations to compare their results with the class. Students will also compare intermolecular forces, binding energy, and geometric orientation of their prepared analog to their initial α-conotoxin TxIA docking results.

Broadband nonlinear vibrational spectroscopy by shaping a coherent fiber supercontinuum.

Vibrational spectroscopy has been widely applied in different fields due to its label-free chemical-sensing capability. Coherent anti-Stokes Raman scattering (CARS) provides stronger signal and faster acquisition than spontaneous Raman scattering, making it especially suitable for molecular imaging. Coherently-controlled single-beam CARS simplifies the conventional multi-beam setup, but the vibrational bandwidth and non-trivial spectrum retrieval have been limiting factors. In this work, a coherent supercontinuum generated in an all-normal-dispersion nonlinear fiber is phase-shaped within a narrow bandwidth for broadband vibrational spectroscopy. The Raman spectra can be directly retrieved from the CARS measurements, covering the fingerprint regime up to 1750 cm(-1). The retrieved spectra of several chemical species agree with their spontaneous Raman data. The compact fiber supercontinuum source offers broad vibrational bandwidth with high stability and sufficient power, showing the potential for spectroscopic imaging in a wide range of applications.

Modified corrections for London forces in solid-state density functional theory calculations of structure and lattice dynamics of molecular crystals.

Dispersion forces are critical for defining the crystal structures and vibrational potentials of molecular crystals. It is, therefore, important to include corrections for these forces in periodic density functional theory (DFT) calculations of lattice vibrational frequencies. In this study, DFT was augmented with a correction term for London-type dispersion forces in the simulations of the structures and terahertz (THz) vibrational spectra of the dispersion-bound solids naphthalene and durene. The parameters of the correction term were modified to best reproduce the experimental crystal structures and THz spectra. It was found that the accurate reproduction of the lattice dimensions by adjusting the magnitude of the applied dispersion forces resulted in the highest-quality fit of the calculated vibrational modes with the observed THz absorptions. The method presented for the modification of the dispersion corrections provides a practical approach to accurately simulating the THz spectra of molecular crystals, accounting for inherent systematic errors imposed by computational and experimental factors.

Revealing the true crystal structure of L-phenylalanine using solid-state density functional theory.

Solid-state density functional theory can be used for crystal structure determination from powder X-ray diffraction data of molecular crystals that are too large and complex for conventional refinement methods.

Noncovalent interactions between modified cytosine and guanine DNA base pair mimics investigated by terahertz spectroscopy and solid-state density functional theory.

Modified cytosine and guanine nucleobases cocrystallize in a hydrogen bonding configuration similar to that observed in native DNA. The noncovalent interactions binding these base pairs in the crystalline solid were investigated using terahertz (THz) spectroscopy and solid-state density functional theory (DFT). While stronger hydrogen bonding interactions are responsible for the general molecular orientations in the crystalline state, it is the weaker dipole-dipole and dispersion forces that determine the overall packing arrangement. The inclusion of dispersion interactions in the DFT calculations was found to be necessary to accurately simulate the unit cell structure and THz vibrational spectrum. Using properly modeled intermolecular potentials, the lattice vibrational motions of the cytosine and guanine derivatives were calculated. The vibrational characters of the modes exhibited by the DNA base pair mimic in the THz region were primarily rotational motions and are indicative of the energies and the nature of vibrations that would likely be observed between similar base pairs in DNA molecules.

Importance of accurate spectral simulations for the analysis of terahertz spectra: citric acid anhydrate and monohydrate.

The terahertz (THz) spectra of crystalline solids are typically uniquely sensitive to the molecular packing configurations, allowing for the detection of polymorphs and hydrates by THz spectroscopic techniques. It is possible, however, that coincident absorptions may be observed between related crystal forms, in which case careful assessment of the lattice vibrations of each system must be performed. Presented here is a THz spectroscopic investigation of citric acid in its anhydrous and monohydrate phases. Remarkably similar features were observed in the THz spectra of both systems, requiring the accurate calculation of the low-frequency vibrational modes by solid-state density functional theory to determine the origins of these spectral features. The results of the simulations demonstrate the necessity of reliable and rigorous methods for THz vibrational modes to ensure the proper evaluation of the THz spectra of molecular solids.

Identification and quantification of polymorphism in the pharmaceutical compound diclofenac acid by terahertz spectroscopy and solid-state density functional theory.

Polymorph detection and quantification in crystalline materials is a principle interest of the pharmaceutical industry. Terahertz (THz) spectroscopy can be used for such analytical applications since this technique is sensitive to the intermolecular interactions of molecules in the solid state. Understanding the fundamental nature of the lattice vibrational motions leading to absorptions in THz spectra is challenging, but may be achieved through computational approaches. In this study, the THz spectra of two diclofenac acid polymorphs were obtained by THz spectroscopy, and the vibrational characters of the observed absorptions were analyzed using solid-state density functional theory (DFT). The results demonstrate the quantitative capacity of THz spectroscopy and the reliability and utility of solid-state DFT in the calculation of low-frequency vibrational motions.

Noncovalent interactions in paired DNA nucleobases investigated by terahertz spectroscopy and solid-state density functional theory.

Cocrystallized adenine and thymine derivatives, along with the pure monomeric crystals, were investigated by terahertz spectroscopy and solid-state density functional theory (DFT). The methylated nucleobase derivatives crystallize in planar hydrogen-bonded adenine-thymine pairs similar to the manner found in DNA. The spectra obtained for 1-methylthymine, 9-methyladenine, and the 1:1 cocrystal in the range of 10-100 cm(-1) clearly demonstrate that absorptions in this spectral range originate from the uniquely ordered assembly and the intermolecular interactions found in each individual crystal system. The quality of spectral reproduction for the DFT simulations of each system was clearly improved by the inclusion of an empirical correction term for London-type dispersion forces to the calculations. Notably, it was found that these weak dispersion forces in the adenine-thymine cocrystal were necessary to produce a properly converged crystal structure and meaningful simulation of the terahertz vibrational spectrum.

Application of London-type dispersion corrections to the solid-state density functional theory simulation of the terahertz spectra of crystalline pharmaceuticals.

The effects of applying an empirical dispersion correction to solid-state density functional theory methods were evaluated in the simulation of the crystal structure and low-frequency (10 to 90 cm(-1)) terahertz spectrum of the non-steroidal anti-inflammatory drug, naproxen. The naproxen molecular crystal is bound largely by weak London force interactions, as well as by more prominent interactions such as hydrogen bonding, and thus serves as a good model for the assessment of the pair-wise dispersion correction term in systems influenced by intermolecular interactions of various strengths. Modifications to the dispersion parameters were tested in both fully optimized unit cell dimensions and those determined by X-ray crystallography, with subsequent simulations of the THz spectrum being performed. Use of the unmodified PBE density functional leads to an unrealistic expansion of the unit cell volume and the poor representation of the THz spectrum. Inclusion of a modified dispersion correction enabled a high-quality simulation of the THz spectrum and crystal structure of naproxen to be achieved without the need for artificially constraining the unit cell dimensions.