Templated Nucleation of Clotrimazole and Ketoprofen on Polymer Substrates.

The use of different template surfaces in crystallization experiments can directly influence the nucleation kinetics, crystal growth, and morphology of active pharmaceutical ingredients (APIs). Consequently, templated nucleation is an attractive approach to enhance crystal nucleation kinetics and preferentially nucleate desired crystal polymorphs for solid-form drug molecules, particularly large and flexible molecules that are difficult to crystallize. Herein, we investigate the effect of polymer templates on the crystal nucleation of clotrimazole and ketoprofen with both experiments and computational methods. Crystallization was carried out in toluene solvent for both APIs with a template library consisting of 12 different polymers. In complement to the experimental studies, we developed a computational workflow based on molecular dynamics (MD) and derived descriptors from the simulations to score and rank API-polymer interactions. The descriptors were used to measure the energy of interaction (EOI), hydrogen bonding, and rugosity (surface roughness) similarity between the APIs and polymer templates. We used a variety of machine learning models (14 in total) along with these descriptors to predict the crystallization outcome of the polymer templates. We found that simply rank-ordering the polymers by their API-polymer interaction energy descriptors yielded 92% accuracy in predicting the experimental outcome for clotrimazole and ketoprofen. The most accurate machine learning model for both APIs was found to be a random forest model. Using these models, we were able to predict the crystallization outcomes for all polymers. Additionally, we have performed a feature importance analysis using the trained models and found that the most predictive features are the energy descriptors. These results demonstrate that API-polymer interaction energies are correlated with heterogeneous crystallization outcomes.

Cocrystal Synthesis through Crystal Structure Prediction.

Crystal structure prediction (CSP) is an invaluable tool in the pharmaceutical industry because it allows to predict all the possible crystalline solid forms of small-molecule active pharmaceutical ingredients. We have used a CSP-based cocrystal prediction method to rank ten potential cocrystal coformers by the energy of the cocrystallization reaction with an antiviral drug candidate, MK-8876, and a triol process intermediate, 2-ethynylglyclerol. For MK-8876, the CSP-based cocrystal prediction was performed retrospectively and successfully predicted the maleic acid cocrystal as the most likely cocrystal to be observed. The triol is known to form two different cocrystals with 1,4-diazabicyclo[2.2.2]octane (DABCO), but a larger solid form landscape was desired. CSP-based cocrystal screening predicted the triol-DABCO cocrystal as rank one, while a triol-l-proline cocrystal was predicted as rank two. Computational finite-temperature corrections enabled determination of relative crystallization propensities of the triol-DABCO cocrystals with different stoichiometries and prediction of the triol-l-proline polymorphs in the free-energy landscape. The triol-l-proline cocrystal was obtained during subsequent targeted cocrystallization experiments and was found to exhibit an improved melting point and deliquescence behavior over the triol-free acid, which could be considered as an alternative solid form in the synthesis of islatravir.

Effect of Polymer Additives on the Crystal Habit of Metformin HCl.

The crystal habit can have a profound influence on the physical properties of crystalline materials, and thus controlling the crystal morphology is of great practical relevance across many industries. Herein, this work investigates the effect of polymer additives on the crystal habit of metformin HCl with both experiments and computational methods with the aim of developing a combined screening approach for crystal morphology engineering. Crystallization experiments of metformin HCl are conducted in methanol and in an isopropanol-water mixture (8:2 V/V). Polyethylene glycol, polyvinylpyrrolidone, Tween80, and hydroxypropyl methylcellulose polymer additives are used in low concentrations (1-2% w/w) in the experiments to study the effect they have on modifying the crystal habit. Additionally, this work has developed computational methods to characterize the morphology "landscape" and quantifies the overall effect of solvent and additives on the predicted crystal habits. Further analysis of the molecular dynamics simulations is used to rationalize the effect of additives on specific crystal faces. This work demonstrates that the effects of additives on the crystal habit are a result of their absorption and interactions with the slow growing {100} and {020} faces.

Absolute Configuration Determination of Chiral API Molecules by MicroED Analysis of Cocrystal Powders Formed Based on Cocrystal Propensity Prediction Calculations.

Establishing the absolute configuration of chiral active pharmaceutical ingredients (APIs) is of great importance. Single crystal X-ray diffraction (scXRD) has traditionally been the method of choice for such analysis, but scXRD requires the growth of large crystals, which can be challenging. Here, we present a method for determining absolute configuration that does not rely on the growth of large crystals. By examining microcrystals formed with chiral probes (small chiral compounds such as amino acids), absolute configuration can be unambiguously determined by microcrystal electron diffraction (MicroED). Our streamlined method employs three steps: (1) virtual screening to identify promising chiral probes, (2) experimental cocrystal screening and (3) structure determination by MicroED and absolute configuration assignment. We successfully applied this method to analyze two chiral API molecules currently on the market for which scXRD was not used to determine absolute configuration.

Selecting a stable solid form of remdesivir using microcrystal electron diffraction and crystal structure prediction.

Therapeutic options in response to the coronavirus disease 2019 (COVID-19) outbreak are urgently needed. In this communication, we demonstrate how to support selection of a stable solid form of an antiviral drug remdesivir in quick time using the microcrystal electron diffraction (MicroED) technique and a cloud-based and artificial intelligence implemented crystal structure prediction platform. We present the MicroED structures of remdesivir forms II and IV and conclude that form II is more stable than form IV at ambient temperature in agreement with experimental observations. The combined experimental and theoretical study can serve as a template for formulation scientists in the pharmaceutical industry.

Theoretical Insights into the Mechanism of Wavelength Regulation in Blue-Absorbing Proteorhodopsin.

Proteorhodopsin (PR) is a light-driven proton pump that is most notable for ushering in the discovery of an ever-increasing number of microbial retinal proteins that are at the forefront of fields such as optogenetics. Two variants, blue (BPR) and green (GPR) proteorhodopsin, have evolved to harvest light at different depths of the ocean. The color-tuning mechanism in PR is controlled by a single residue at position 105: in BPR it is a glutamine, whereas in GPR it is a leucine. Although the majority of studies on the spectral tuning mechanism in PR have focused on GPR, detailed understanding of the electronic environment responsible for spectral tuning in BPR is lacking. In this work, several BPR models were investigated using quantum mechanics/molecular mechanics (QM/MM) calculations to obtain fundamental insights into the color tuning mechanism of BPR. We find that the molecular mechanism of spectral tuning in BPR depends on two geometric parameters, the bond length alternation and the torsion angle deviation of the all-trans-retinyl chromophore. Both parameters are influenced by the strength of the hydrogen-bonded networks in the chromophore-binding pocket, which shows how BPR is different from other microbial rhodopsins.

Structure-Based Macrocycle Design in Small-Molecule Drug Discovery and Simple Metrics To Identify Opportunities for Macrocyclization of Small-Molecule Ligands.

Interest is growing in the use of macrocycles in pharmaceutical discovery. Macrocylization may provide a gateway to an expanded chemical space for small-molecule drug discovery, and this could be beneficial in prosecuting difficult targets, e.g., protein-protein interactions. Most, but not all, macrocycle drugs are derived from natural products. Studies on synthetic drug-like small-molecule macrocycles are limited, and our current understanding of macrocycle drugs is similarly limited. Following some background discussion, we review several examples of the structure-based design of synthetic macrocycles. Our opinion is that in conformationally suitable systems macrocycles are an analog class worthy of consideration. We then summarize an approach for the initial evaluation of molecules as candidates for macrocyclization.

Activation of the ß-common receptor by erythropoietin impairs acetylcholine-mediated vasodilation in mouse mesenteric arterioles.

Clinically, erythropoietin (EPO) is known to increase systemic vascular resistance and arterial blood pressure. However, EPO stimulates the production of the potent vasodilator, nitric oxide (NO), in culture endothelial cells. The mechanism by which EPO causes vasoconstriction despite stimulating NO production may be dependent on its ability to activate two receptor complexes, the homodimeric EPO (EPOR2 ) and the heterodimeric EPOR/ß-common receptor (ßCR). The purpose of this study was to investigate the contribution of each receptor to the vasoactive properties of EPO. First-order, mesenteric arteries were isolated from 16-week-old male C57BL/6 mice, and arterial function was studied in pressure arteriographs. To determine the contribution of each receptor complex, EPO-stimulating peptide (ESP), which binds and activates the heterodimeric EPOR/ßCR complex, and EPO, which activates both receptors, were added to the arteriograph chamber 20 min prior to evaluation of endothelium-dependent (acetylcholine, bradykinin, A23187) and endothelium-independent (sodium nitroprusside) vasodilator responses. Only ACh-induced vasodilation was impaired in arteries pretreated with EPO or ESP. EPO and ESP pretreatment abolished ACh-induced vasodilation by 100% and 60%, respectively. EPO and ESP did not affect endothelium-independent vasodilation by SNP. Additionally, a novel ßCR inhibitory peptide (ßIP), which was computationally developed, prevented the impairment of acetylcholine-induced vasodilation by EPO and ESP, further implicating the EPOR/ßCR complex. Last, pretreatment with either EPO or ESP did not affect vasoconstriction by phenylephrine and KCl. Taken together, these findings suggest that acute activation of the heterodimeric EPOR/ßCR in endothelial cells leads to a selective impairment of ACh-mediated vasodilator response in mouse mesenteric resistance arteries.

Molecular mechanism of activation of human musk receptors OR5AN1 and OR1A1 by (R)-muscone and diverse other musk-smelling compounds.

Understanding olfaction at the molecular level is challenging due to the lack of crystallographic models of odorant receptors (ORs). To better understand the molecular mechanism of OR activation, we focused on chiral (R)-muscone and other musk-smelling odorants due to their great importance and widespread use in perfumery and traditional medicine, as well as environmental concerns associated with bioaccumulation of musks with estrogenic/antiestrogenic properties. We experimentally and computationally examined the activation of human receptors OR5AN1 and OR1A1, recently identified as specifically responding to musk compounds. OR5AN1 responds at nanomolar concentrations to musk ketone and robustly to macrocyclic sulfoxides and fluorine-substituted macrocyclic ketones; OR1A1 responds only to nitromusks. Structural models of OR5AN1 and OR1A1 based on quantum mechanics/molecular mechanics (QM/MM) hybrid methods were validated through direct comparisons with activation profiles from site-directed mutagenesis experiments and analysis of binding energies for 35 musk-related odorants. The experimentally found chiral selectivity of OR5AN1 to (R)- over (S)-muscone was also computationally confirmed for muscone and fluorinated (R)-muscone analogs. Structural models show that OR5AN1, highly responsive to nitromusks over macrocyclic musks, stabilizes odorants by hydrogen bonding to Tyr260 of transmembrane α-helix 6 and hydrophobic interactions with surrounding aromatic residues Phe105, Phe194, and Phe207. The binding of OR1A1 to nitromusks is stabilized by hydrogen bonding to Tyr258 along with hydrophobic interactions with surrounding aromatic residues Tyr251 and Phe206. Hydrophobic/nonpolar and hydrogen bonding interactions contribute, respectively, 77% and 13% to the odorant binding affinities, as shown by an atom-based quantitative structure-activity relationship model.

Computational design and experimental characterization of a novel ß-common receptor inhibitory peptide.

In short-term animal models of ischemia, erythropoietin (EPO) signaling through the heterodimeric EPO receptor (EPOR)/ß-common receptor (ßCR) is believed to elicit tissue protective effects. However, large, randomized, controlled trials demonstrate that targeting a higher hemoglobin level by administering higher doses of EPO, which are more likely to activate the heterodimeric EPOR/ßCR, is associated with an increase in adverse cardiovascular events. Thus, inhibition of long-term activation of the ßCR may have therapeutic implications. This study aimed to design and evaluate the efficacy of novel computationally designed ßCR inhibitory peptides (ßIP). These novel ßIPs were designed based on a truncated portion of Helix-A from EPO, specifically residues 11-26 (VLERYLLEAKEAEKIT). Seven novel peptides (P1 to P7) were designed. Peptide 7 (P7), VLERYLHEAKHAEKIT, demonstrated the most robust inhibitory activity. We also report here the ability of P7 to inhibit ßCR-induced nitric oxide (NO) production and angiogenesis in human umbilical vein endothelial cells (HUVECs). Specifically, we found that P7 ßIP completely abolished EPO-induced NO production. The inhibitory effect could be overcome with super physiological doses of EPO, suggesting a competitive inhibition. ßCR-induced angiogenesis in HUVEC's was also abolished with treatment of P7 ßIP, but P7 ßIP did not inhibit vascular endothelial growth factor (VEGF)-induced angiogenesis. In addition, we demonstrate that the novel P7 ßIP does not inhibit EPO-induced erythropoiesis with use of peripheral blood mononuclear cells (PBMCs). These results, for the first time, describe a novel, potent ßCR peptide inhibitor that inhibit the actions of the ßCR without affecting erythropoiesis.

In Silico Prediction of Ligand Binding Energies in Multiple Therapeutic Targets and Diverse Ligand Sets-A Case Study on BACE1, TYK2, HSP90, and PERK Proteins.

We present here the use of QM/MM LIE (linear interaction energy) based binding free energy calculations that greatly improve the precision and accuracy of predicting experimental binding affinities, in comparison to most current binding free energy methodologies, while maintaining reasonable computational times. Calculations are done for four sets of ligand-protein complexes, chosen on the basis of diversity of protein types and availability of experimental data, totaling 140 ligands binding to therapeutic protein targets BACE1, TYK2, HSP90, and PERK. This method allows calculations for a diverse set of ligands and multiple protein targets without the need for parametrization or extra calculations. The accuracy achieved with this method can be used to guide small molecule computational drug design programs.

Probing the remarkable thermal kinetics of visual rhodopsin with E181Q and S186A mutants.

We recently reported a very unusual temperature dependence of the rate of thermal reaction of wild type bovine rhodopsin: the Arrhenius plot exhibits a sharp "elbow" at 47 °C and, in the upper temperature range, an unexpectedly large activation energy (114 ± 8 kcal/mol) and an enormous prefactor (1072±5 s-1). In this report, we present new measurements and a theoretical model that establish convincingly that this behavior results from a collective, entropy-driven breakup of the rigid hydrogen bonding networks (HBNs) that hinder the reaction at lower temperatures. For E181Q and S186A, two rhodopsin mutants that disrupt the HBNs near the binding pocket of the 11-cis retinyl chromophore, we observe significant decreases in the activation energy (∼90 kcal/mol) and prefactor (∼1060 s-1), consistent with the conclusion that the reaction rate is enhanced by breakup of the HBN. The results provide insights into the molecular mechanism of dim-light vision and eye diseases caused by inherited mutations in the rhodopsin gene that perturb the HBNs.

Implausibility of the vibrational theory of olfaction.

The vibrational theory of olfaction assumes that electron transfer occurs across odorants at the active sites of odorant receptors (ORs), serving as a sensitive measure of odorant vibrational frequencies, ultimately leading to olfactory perception. A previous study reported that human subjects differentiated hydrogen/deuterium isotopomers (isomers with isotopic atoms) of the musk compound cyclopentadecanone as evidence supporting the theory. Here, we find no evidence for such differentiation at the molecular level. In fact, we find that the human musk-recognizing receptor, OR5AN1, identified using a heterologous OR expression system and robustly responding to cyclopentadecanone and muscone, fails to distinguish isotopomers of these compounds in vitro. Furthermore, the mouse (methylthio)methanethiol-recognizing receptor, MOR244-3, as well as other selected human and mouse ORs, responded similarly to normal, deuterated, and (13)C isotopomers of their respective ligands, paralleling our results with the musk receptor OR5AN1. These findings suggest that the proposed vibration theory does not apply to the human musk receptor OR5AN1, mouse thiol receptor MOR244-3, or other ORs examined. Also, contrary to the vibration theory predictions, muscone-d30 lacks the 1,380- to 1,550-cm(-1) IR bands claimed to be essential for musk odor. Furthermore, our theoretical analysis shows that the proposed electron transfer mechanism of the vibrational frequencies of odorants could be easily suppressed by quantum effects of nonodorant molecular vibrational modes. These and other concerns about electron transfer at ORs, together with our extensive experimental data, argue against the plausibility of the vibration theory.

Activation of OR1A1 suppresses PPAR-γ expression by inducing HES-1 in cultured hepatocytes.

Olfactory receptors (ORs) comprise the largest G protein-coupled receptor gene superfamily. Recent studies indicate that ORs are also expressed in non-olfactory organs, including metabolically active tissues, although their biological functions in these tissues are largely unknown. In this study, OR1A1 expression was detected in HepG2 liver cells. OR1A1 activation by (-)-carvone, a known OR1A1 ligand, increased the cyclic adenosine monophosphate (cAMP), but not intracellular Ca(2+) concentration, thereby inducing protein kinase A (PKA) activity with subsequent phosphorylation of cAMP response element-binding protein (CREB) and upregulation of the CREB-responsive gene hairy and enhancer of split (HES)-1, a corepressor of peroxisome proliferator-activated receptor-γ (PPAR-γ) in hepatocytes. In (-)-carvone-stimulated cells, the repression of PPAR-γ reduced the expression of the target gene, mitochondrial glycerol-3-phosphate acyltransferase, which encodes a key enzyme involved in triglyceride synthesis. Intracellular triglyceride level and lipid accumulation were reduced in cells stimulated with (-)-carvone, effects that were diminished following the loss of OR1A1 function. These results indicate that OR1A1 may function as a non-redundant receptor in hepatocytes that regulates the PKA-CREB-HES-1 signaling axis and thereby modulates hepatic triglyceride metabolism.

Kinetics of thermal activation of an ultraviolet cone pigment.

Visual pigments can be thermally activated via isomerization of the retinyl chromophore and hydrolysis of the Schiff base (SB) through which the retinyl chromophore is bound to the opsin protein. Here, we present the first combined experimental and theoretical study of the thermal activation of a Siberian hamster ultraviolet (SHUV) pigment. We measured the rates of thermal isomerization and hydrolysis in the SHUV pigment and bovine rhodopsin. We found that these rates were significantly faster in the UV pigment than in rhodopsin due to the difference in the structural and electrostatic effects surrounding the unprotonated Schiff base (USB) retinyl chromophore in the UV pigment. Theoretical (DFT-QM/MM) calculations of the cis-trans thermal isomerization revealed a barrier of ∼23 kcal/mol for the USB retinyl chromophore in SHUV compared to ∼40 kcal/mol for protonated Schiff base (PSB) chromophore in rhodopsin. The lower barrier for thermal isomerization in the SHUV pigment is attributed to the (i) lessening of the steric restraints near the ß-ionone ring and SB ends of the chromophore, (ii) displacement of the transmembrane helix 6 (TM6) away from the binding pocket toward TM5 due to absence of the salt bridge between the USB and the protonated E113 residue, and (iii) change in orientation of the hydrogen-bonding networks (HBNs) in the extracellular loop 2 (EII). The results in comparing thermal stability of UV cone pigment and rhodopsin provide insight into molecular evolution of vertebrate visual pigments in achieving low discrete dark noise and high photosensitivity in rod pigments for dim-light vision.

QM/MM model of the mouse olfactory receptor MOR244-3 validated by site-directed mutagenesis experiments.

Understanding structure/function relationships of olfactory receptors is challenging due to the lack of x-ray structural models. Here, we introduce a QM/MM model of the mouse olfactory receptor MOR244-3, responsive to organosulfur odorants such as (methylthio)methanethiol. The binding site consists of a copper ion bound to the heteroatoms of amino-acid residues H105, C109, and N202. The model is consistent with site-directed mutagenesis experiments and biochemical measurements of the receptor activation, and thus provides a valuable framework for further studies of the sense of smell at the molecular level.

Unusual kinetics of thermal decay of dim-light photoreceptors in vertebrate vision.

We present measurements of rate constants for thermal-induced reactions of the 11-cis retinyl chromophore in vertebrate visual pigment rhodopsin, a process that produces noise and limits the sensitivity of vision in dim light. At temperatures of 52.0-64.6 °C, the rate constants fit well to an Arrhenius straight line with, however, an unexpectedly large activation energy of 114 ± 8 kcal/mol, which is much larger than the 60-kcal/mol photoactivation energy at 500 nm. Moreover, we obtain an unprecedentedly large prefactor of 10(72±5) s(-1), which is roughly 60 orders of magnitude larger than typical frequencies of molecular motions! At lower temperatures, the measured Arrhenius parameters become more normal: Ea = 22 ± 2 kcal/mol and Apref = 10(9±1) s(-1) in the range of 37.0-44.5 °C. We present a theoretical framework and supporting calculations that attribute this unusual temperature-dependent kinetics of rhodopsin to a lowering of the reaction barrier at higher temperatures due to entropy-driven partial breakup of the rigid hydrogen-bonding network that hinders the reaction at lower temperatures.

Spectral tuning of ultraviolet cone pigments: an interhelical lock mechanism.

Ultraviolet (UV) cone pigments can provide insights into the molecular evolution of vertebrate vision since they are nearer to ancestral pigments than the dim-light rod photoreceptor rhodopsin. While visible-absorbing pigments contain an 11-cis retinyl chromophore with a protonated Schiff-base (PSB11), UV pigments uniquely contain an unprotonated Schiff-base (USB11). Upon F86Y mutation in model UV pigments, both the USB11 and PSB11 forms of the chromophore are found to coexist at physiological pH. The origin of this intriguing equilibrium remains to be understood at the molecular level. Here, we address this phenomenon and the role of the USB11 environment in spectral tuning by combining mutagenesis studies with spectroscopic (UV-vis) and theoretical [DFT-QM/MM (SORCI+Q//B3LYP/6-31G(d): Amber96)] analysis. We compare structural models of the wild-type (WT), F86Y, S90A and S90C mutants of Siberian hamster ultraviolet (SHUV) cone pigment to explore structural rearrangements that stabilize USB11 over PSB11. We find that the PSB11 forms upon F86Y mutation and is stabilized by an "inter-helical lock" (IHL) established by hydrogen-bonding networks between transmembrane (TM) helices TM6, TM2, and TM3 (including water w2c and amino acid residues Y265, F86Y, G117, S118, A114, and E113). The findings implicate the involvement of the IHL in constraining the displacement of TM6, an essential component of the activation of rhodopsin, in the spectral tuning of UV pigments.

Spectral tuning in halorhodopsin: the chloride pump photoreceptor.

The spectral tuning of halorhodopsin from Halobacterium salinarum (shR) during anion transport was analyzed at the molecular level using DFT-QM/MM [SORCI+Q//B3LYP/6-31G(d):Amber96] hybrid methods. Insights into the influence of Cl(-) depletion, Cl(-) substitution by N3(-) or NO3(-), and mutation of key amino acid residues along the ion translocation pathway (H95A, H95R, Q105E, R108H, R108I, R108K, R108Q, T111V, R200A, R200H, R200K, R200Q, and T203V) were analyzed for the first time in a fully atomistic model of the shR photoreceptor. We found evidence that structural rearrangements mediated by specific hydrogen bonds of internal water molecules and counterions (D238 and Cl(-)) in the active site induce changes in the bond-length alternation of the all-trans retinyl chromophore and affect the wavelength of maximal absorption in shR.

The active site of melanopsin: the biological clock photoreceptor.

The nonvisual ocular photoreceptor melanopsin, found in the neurons of vertebrate inner retina, absorbs blue light and triggers the "biological clock" of mammals by activating the suprachiasmatic nuclei (a small region of the brain that regulates the circadian rhythms of neuronal and hormonal activities over 24 h cycles). The structure of melanopsin, however, has yet to be established. Here, we propose for the first time a structural model of the active site of mouse melanopsin. The homology model is based on the crystal structure of squid rhodopsin (λ(max) = 490 nm) and shows a maximal absorbance (λ(max) = 447 nm) consistent with the observed absorption of the photoreceptor. The 43 nm spectral shift is due to an increased bond-length alternation of the protonated Schiff base of 11-cis-retinal chromophore, induced by N87Q mutation and water-mediated H-bonding interactions with the Schiff base linkage. These findings, analogous to spectral changes observed in the G89Q bovine rhodopsin mutant, suggest that single site mutations can convert photopigments into visual light sensors or nonvisual sensory photoreceptors.