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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Effect of waters of crystallization on terahertz spectra: anhydrous oxalic acid and its dihydrate.

Oxalic acid and oxalic acid dihydrate were studied using terahertz spectroscopy and solid-state density functional theory (DFT) in the spectral range 10-100 cm(-1). The size of the oxalic acid molecule and its limited internal degrees of freedom make it ideal for evaluating the performance of computational methods for the structural and dynamical simulation of strongly hydrogen-bonded solids. Calculations of the solid-state structures and terahertz spectra of oxalic acid and oxalic acid dihydrate were performed using the hybrid B3LYP and B3PW91 and the nonhybrid BLYP and PW91 density functionals employing the 6-311G(2d,2p) basis set. When these simulations were compared to the experimental spectra of the oxalic acid solids, a constant overprediction of the dihydrate frequencies was observed in contrast to the results of the anhydrous system. This change in behavior is connected to the nature of the vibrational motions being accessed. The primary molecular motion contributions to the terahertz vibrations of oxalic acid dihydrate were found to originate in the external motions of the cocrystallized H(2)O molecules. The observed overestimation of the vibrational energies in the simulated terahertz spectra is attributed to increased anharmonicity of the vibrational motions in the dihydrate system versus the anhydrous, resulting from weaker hydrogen bonding through the networked water molecules.

Investigating the anharmonicity of lattice vibrations in water-containing molecular crystals through the terahertz spectroscopy of L-serine monohydrate.

The influence of cocrystallized H(2)O molecules on the terahertz (THz) spectra and corresponding computational treatment of hydrated molecular crystals was investigated in the study of protonated and deuterium-substituted l-serine.H(2)O. The THz spectra of both solids have been measured in the range of 10 to 90 cm(-1), with simulations of the crystalline structure and THz vibrational modes performed using solid-state density functional theory. Significant and systematic overestimations of the predicted vibrational frequencies were observed in all calculations. Evidence provided by the comparison of the experimental and calculated vibrational frequencies for both the protonated and deuterated l-serine.H(2)O solids indicates the presence of significant anharmonicity in the observed lattice vibrations. The results suggest that vibrational anharmonicity may play a much larger role in the interpretation of the THz spectra of hydrates in contrast to their corresponding anhydrous forms.

Discrimination of chiral solids: a terahertz spectroscopic investigation of L- and DL-serine.

The terahertz (THz, far-infrared) spectra of enantiomerically pure and racemic crystalline serine were investigated using time-domain THz spectroscopy and solid-state density functional theory (DFT) in the spectral range of 10-90 cm(-1). The experimental THz spectra of L- and DL-serine at 78 K appear quite similar despite the significant differences in arrangement of the molecules in their crystal structures. Structural analyses of the two systems and calculation of the vibrational modes and intensities were performed using DFT with periodic boundary conditions employing the B3LYP and PW91 density functionals with the 6-31G(d,p) and 6-311G(d,p) basis sets. The applied computational methods produced simulations of the THz spectra in good agreement with experiment, with accurate predictions of the subtle differences in the THz spectra of the two chiral solid-state mixtures of serine. The observed spectral features are assigned as primarily external lattice translations and rotations with lesser contributions due to intramolecular torsions of the -NH(3)(+) and -COO(-) groups modified by intermolecular hydrogen bonding.