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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Color vision: "OH-site" rule for seeing red and green.

Eyes gather information, and color forms an extremely important component of the information, more so in the case of animals to forage and navigate within their immediate environment. By using the ONIOM (QM/MM) (ONIOM = our own N-layer integrated molecular orbital plus molecular mechanics) method, we report a comprehensive theoretical analysis of the structure and molecular mechanism of spectral tuning of monkey red- and green-sensitive visual pigments. We show that interaction of retinal with three hydroxyl-bearing amino acids near the ß-ionone ring part of the retinal in opsin, A164S, F261Y, and A269T, increases the electron delocalization, decreases the bond length alternation, and leads to variation in the wavelength of maximal absorbance of the retinal in the red- and green-sensitive visual pigments. On the basis of the analysis, we propose the "OH-site" rule for seeing red and green. This rule is also shown to account for the spectral shifts obtained from hydroxyl-bearing amino acids near the Schiff base in different visual pigments: at site 292 (A292S, A292Y, and A292T) in bovine and at site 111 (Y111) in squid opsins. Therefore, the OH-site rule is shown to be site-specific and not pigment-specific and thus can be used for tracking spectral shifts in any visual pigment.

QM/MM study of the structure, energy storage, and origin of the bathochromic shift in vertebrate and invertebrate bathorhodopsins.

By comparing the results from a hybrid quantum mechanics/molecular mechanics method (SORCI+Q//B3LYP/6-31G*:Amber) between vertebrate (bovine) and invertebrate (squid) visual pigments, the mechanism of molecular rearrangements, energy storage, and origin of the bathochromic shift accompanying the transformation of rhodopsin to bathorhodopsin have been evaluated. The analysis reveals that, in the presence of an unrelaxed binding site, bathorhodopsin was found to carry almost 27 kcal/mol energy in both visual pigments and absorb (λ(max)) at 528 nm in bovine and 554 nm in squid. However, when the residues within 4.0 Å radius of the retinal are relaxed during the isomerization event, almost ∼16 kcal/mol energy is lost in squid compared to only ∼8 kcal/mol in bovine. Loss of a larger amount of energy in squid is attributed to the presence of a flexible binding site compared to a rigid binding site in bovine. Structure of the squid bathorhodopsin is characterized by formation of a direct H-bond between the Schiff base and Asn87.

Why 11-cis-retinal? Why not 7-cis-, 9-cis-, or 13-cis-retinal in the eye?

One of the basic and unresolved puzzles in the chemistry of vision concerns the natural selection of 11-cis-retinal as the light-sensing chromophore in visual pigments. A detailed computational examination of the structure, stability, energetics, and spectroscopy of 7-cis-, 9-cis-, 11-cis-, and 13-cis-retinal isomers in vertebrate (bovine, monkey) and invertebrate (squid) visual pigments was carried out using a hybrid quantum mechanics/molecular mechanics (QM/MM) method. The results show that the electrostatic interaction between retinal and opsin dominates the natural selection of 11-cis-retinal over other cis isomers in the dark state. In all of the pigments, 9-cis-retinal was found to be only slightly higher in energy than 11-cis-retinal, which provides strong evidence for the presence of 9-cis-rhodopsin in nature. 7-cis-Retinal is suggested to be an "upside-down" version of the all-trans isomer because the structural rearrangements observed for 7-cis-rhodopsin from squid were found to be very similar to those for squid bathorhodopsin. The progressive red shift in the calculated absorption wavelength (λ(max)) (431, 456, 490, and 508 nm for the 7-cis-, 9-cis-, 11-cis-, and 13-cis-retinal isomers) is due to the decrease in bond length alternation of the retinal.

Self-assembly of dendronized perylene bisimides into complex helical columns.

The synthesis of perylene 3,4:9,10-tetracarboxylic acid bisimides (PBIs) dendronized with first-generation dendrons containing 0 to 4 methylenic units (m) between the imide group and the dendron, (3,4,5)12G1-m-PBI, is reported. Structural analysis of their self-organized arrays by DSC, X-ray diffraction, molecular modeling, and solid-state (1)H NMR was carried out on oriented samples with heating and cooling rates of 20 to 0.2 °C/min. At high temperature, (3,4,5)12G1-m-PBI self-assemble into 2D-hexagonal columnar phases with intracolumnar order. At low temperature, they form orthorhombic (m = 0, 2, 3, 4) and monoclinic (m = 1) columnar arrays with 3D periodicity. The orthorhombic phase has symmetry close to hexagonal. For m = 0, 2, 3, 4 ,they consist of tetramers as basic units. The tetramers contain a pair of two molecules arranged side by side and another pair in the next stratum of the column, turned upside-down and rotated around the column axis at different angles for different m. In contrast, for m = 1, there is only one molecule in each stratum, with a four-strata 2(1) helical repeat. All molecules face up in one column, and down in the second column, of the monoclinic cell. This allows close and extended π-stacking, unlike in the disruptive up-down alteration from the case of m = 0, 2, 3, 4. Most of the 3D structures were observed only by cooling at rates of 1 °C/min or less. This complex helical self-assembly is representative for other classes of dendronized PBIs investigated for organic electronics and solar cells.

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.

QM/MM study of dehydro and dihydro ß-ionone retinal analogues in squid and bovine rhodopsins: implications for vision in salamander rhodopsin.

Visual pigment rhodopsin provides a decisive crossing point for interaction between organisms and environment. Naturally occurring visual pigments contain only PSB11 and 3,4-dehydro-PSB11 as chromophores. Therefore, the ability of visual opsin to discriminate between the retinal geometries is investigated by means of QM/MM incorporation of PSB11, 6-s-cis and 6-s-trans forms of 3,4-dehydro-PSB11, and 3,4-dehydro-5,6-dihydro-PSB11 and 5,6-dihydro-PSB11 analogues into squid and bovine rhodopsin environments. The analogue-protein interaction reveals the binding site of squid rhodopsin to be malleable and ductile, while that of bovine rhodopsin is rigid and stiff. On the basis of these studies, a tentative model of the salamander rhodopsin binding site is also proposed.

Columnar packing motifs of functionalized perylene derivatives: local molecular order despite long-range disorder.

We elucidate local packing motifs and dynamical order parameters in a perylene tetracarboxydiimide derivative (C(8,7)-PDI), one of the most promising candidates for rationally designed, self-assembling, and self-healing molecular wires. Spectroscopic fingerprints obtained from solid-state NMR spectroscopy are interpreted by means of first-principles calculations and molecular dynamics simulations. The interplay of steric repulsion, H bonding, and pi-pi packing effects leads to a specific relative molecular pitch angle of approximately 35 +/- 10 degrees between successive molecules in the stack. Dynamical order parameters, determined from NMR sideband patterns as a measure of molecular motion, yield values of S approximately = 1.0 in the core of the columnar stack, corresponding to an almost frozen molecular dynamics at ambient temperature. This rigidity is compatible with characteristic intermolecular distances obtained from dipolar couplings between specific hydrogens via double-quantum NMR experiments and further supported by ab initio calculations.

Water-mediated spectral shifts in rhodopsin and bathorhodopsin.

The role of water molecules in spectral tuning of proteins has been left largely unexplored. This topic is important because changing hydrogen bond patterns during the activation process may lead to spectral shifts which can be of diagnostic value for the underlying structures. Arguments put forward in this article are based on spectral shift calculations of the rhodopsin and bathorhodopsin chromophore due to wat2a and 2b in the presence and absence of the counterion and of the amino acids lining the rhodopsin binding pocket. They show, among others, that a single water molecule can shift the absorbance by up to 0.1 eV or 34 nm depending on the environment of the chromophore.