A Distributional Model of Bound Ligand Conformational Strain: From Small Molecules up to Large Peptidic Macrocycles.

The internal conformational strain incurred by ligands upon binding a target site has a critical impact on binding affinity, and expectations about the magnitude of ligand strain guide conformational search protocols. Estimates for bound ligand strain begin with modeled ligand atomic coordinates from X-ray co-crystal structures. By deriving low-energy conformational ensembles to fit X-ray diffraction data, calculated strain energies are substantially reduced compared with prior approaches. We show that the distribution of expected global strain energy values is dependent on molecular size in a superlinear manner. The distribution of strain energy follows a rectified normal distribution whose mean and variance are related to conformational complexity. The modeled strain distribution closely matches calculated strain values from experimental data comprising over 3000 protein-ligand complexes. The distributional model has direct implications for conformational search protocols as well as for directions in molecular design.

Discovery and Structure-Based Design of Macrocyclic Peptides Targeting STUB1.

Recent evidence suggests that deletion of STUB1─a pivotal negative regulator of interferon-γ sensing─may potentially clear malignant cells. However, current studies rely primarily on genetic approaches, as pharmacological inhibitors of STUB1 are lacking. Identifying a tool compound will be a step toward validating the target in a broader therapeutic sense. Herein, screening more than a billion macrocyclic peptides resulted in STUB1 binders, which were further optimized by a structure-enabled in silico design. The strategy to replace the macrocyclic peptides' hydrophilic and solvent-exposed region with a hydrophobic scaffold improved cellular permeability while maintaining the binding conformation. Further substitution of the permeability-limiting terminal aspartic acid with a tetrazole bioisostere retained the binding to a certain extent while improving permeability, suggesting a path forward. Although not optimal for cellular study, the current lead provides a valuable template for further development into selective tool compounds for STUB1 to enable target validation.

Denigrins and Dactylpyrroles, Arylpyrrole Alkaloids from a Dactylia sp. Marine Sponge.

Seven new arylpyrrole alkaloids (1-7), along with four known compounds, were isolated from an extract of a Dactylia sp. nov. marine sponge, and their structures were elucidated by interpretation of NMR and MS spectroscopic data. Denigrins D-G (1-4) have highly substituted pyrrole or pyrrolone rings in their core structures, while dactylpyrroles A-C (5-7) have tricyclic phenanthrene cores. Due to the proton-deficient nature of these scaffolds, key heteronuclear correlations from 1H-15N HMBC and LR-HSQMBC NMR experiments were used in the structure assignment of denigrin D (1). Dictyodendrin F (8), a previously described co-metabolite, inhibited transcription driven by the oncogenic PAX3-FOXO1 fusion gene with an IC50 value of 13 µM.

Combination of HDX-MS and in silico modeling to study enzymatic reactivity and stereo-selectivity at different solvent conditions.

The higher-order structure of a protein defines its function, and protein structural dynamics are often essential for protein binding and enzyme catalysis. Methods for protein characterization in solution are continuously being developed to understand and explore protein conformational changes with regards to function and activity. The goal of this study was to survey the use of combining HDX-MS global conformational screening with in silico modeling and continuous labeling peptide-level HDX-MS as an approach to highlight regions of interest within an enzyme required for biocatalytic processes. We surveyed in silico modeling correlated with peptide level HDX-MS experiments to characterize and localize transaminase enzyme structural dynamics at different conditions. This approach was orthogonally correlated with a global Size-Exclusion-HDX (SEC-HDX) screen for global conformational comparison and global alpha-helical content measurements by circular dichroism. Enzymatic activity and stereo-selectivity of transaminases were compared at different reaction-solution conditions that forced protein conformational changes by increasing acetonitrile concentration. The experimental peptide-level HDX-MS results demonstrated similar trends to the modeling data showing that certain regions remained folded in transaminases ATA-036 and ATA-303 with increasing acetonitrile concentration, which is also associated with shifting stereoselectivity. HDX modeling, SEC-HDX and CD experimental data showed that transaminase ATA-234 had the highest level of global unfolding with increasing acetonitrile concentration compared to the other two enzymes, which correlated with drastically reduced product conversion in transamination reaction. The combined HDX modeling/experimental workflow, based on enzymatic reactions studied at different conditions to induce changes in enzyme conformation, could be used as a tool to guide directed evolution efforts by identifying and focusing on the regions of an enzyme required for reaction product conversion and stereoselectivity.

Complex macrocycle exploration: parallel, heuristic, and constraint-based conformer generation using ForceGen.

ForceGen is a template-free, non-stochastic approach for 2D to 3D structure generation and conformational elaboration for small molecules, including both non-macrocycles and macrocycles. For conformational search of non-macrocycles, ForceGen is both faster and more accurate than the best of all tested methods on a very large, independently curated benchmark of 2859 PDB ligands. In this study, the primary results are on macrocycles, including results for 431 unique examples from four separate benchmarks. These include complex peptide and peptide-like cases that can form networks of internal hydrogen bonds. By making use of new physical movements ("flips" of near-linear sub-cycles and explicit formation of hydrogen bonds), ForceGen exhibited statistically significantly better performance for overall RMS deviation from experimental coordinates than all other approaches. The algorithmic approach offers natural parallelization across multiple computing-cores. On a modest multi-core workstation, for all but the most complex macrocycles, median wall-clock times were generally under a minute in fast search mode and under 2 min using thorough search. On the most complex cases (roughly cyclic decapeptides and larger) explicit exploration of likely hydrogen bonding networks yielded marked improvements, but with calculation times increasing to several minutes and in some cases to roughly an hour for fast search. In complex cases, utilization of NMR data to constrain conformational search produces accurate conformational ensembles representative of solution state macrocycle behavior. On macrocycles of typical complexity (up to 21 rotatable macrocyclic and exocyclic bonds), design-focused macrocycle optimization can be practically supported by computational chemistry at interactive time-scales, with conformational ensemble accuracy equaling what is seen with non-macrocyclic ligands. For more complex macrocycles, inclusion of sparse biophysical data is a helpful adjunct to computation.

Optical Chirality Sensing with a Stereodynamic Aluminum Biphenolate Probe.

The determination of the enantiopurity and the concentration of chiral compounds by chiroptical sensing with molecular probes is increasingly attractive for high-throughput screening applications including streamlined asymmetric reaction development. In this study, we use stereodynamic aluminum biphenolate complexes for quantitative ee and concentration analysis of amino alcohols and α-hydroxy acids. An important feature of the tropos biphenolate ligand used is the presence of a phenylacetylene antenna for optimal chirality recognition and CD/UV responses at high wavelengths. The complexation-driven chirality amplification yields strong CD signals, which allows quantitative chiroptical sensing with good accuracy. We show that aluminate biphenolate sensors can exhibit linear and nonlinear correlations between the induced CD signals and the enantiomeric composition or concentration of the chiral substrate.

A multifunctional catalyst that stereoselectively assembles prodrugs.

The catalytic stereoselective synthesis of compounds with chiral phosphorus centers remains an unsolved problem. State-of-the-art methods rely on resolution or stoichiometric chiral auxiliaries. Phosphoramidate prodrugs are a critical component of pronucleotide (ProTide) therapies used in the treatment of viral disease and cancer. Here we describe the development of a catalytic stereoselective method for the installation of phosphorus-stereogenic phosphoramidates to nucleosides through a dynamic stereoselective process. Detailed mechanistic studies and computational modeling led to the rational design of a multifunctional catalyst that enables stereoselectivity as high as 99:1.

Mechanistic Study of the Gas-Phase In-Source Hofmann Elimination of Doubly Quaternized Cinchona-Alkaloid Based Phase-Transfer Catalysts by (+)-Electrospray Ionization/Tandem Mass Spectrometry.

An unusual in-source fragmentation pattern observed for 14 doubly quaternized cinchona alkaloid-based phase-transfer catalysts (PTC) was studied using (+)-ESI high resolution mass spectrometry. Loss of the substituted benzyl cation (R1 or R2) was found to be the major product ion [M2+ - R1+ or R2+]+ in MS spectra of all PTC compounds. A Hofmann elimination product ion [M - H]+ was also observed. Only a small amount of the doubly charged M2+ ions were observed in the MS spectra, likely due to strong Columbic repulsion between the two quaternary ammonium cations in the gas phase. The positive voltage in the MS inlet but not the ESI probe was found to induce this extensive fragmentation for all PTC diboromo-salts. Compound 1 was used as an example to illustrate the proposed in-source fragmentation mechanism. The mechanism of formation of the Hofmann elimination product ion [M - H]+ was further investigated using HRMS/MS, H/D exchange, and DFT calculations. The proposed formation of 2b as the major Hofmann elimination product ion was supported both by HRMS/MS and DFT calculations. Formation of product ion 2b through a concerted unimolecular Ei elimination pathway is proposed rather than a bimolecular E2 elimination pathway for common solution Hofmann eliminations. Graphical Abstract ᅟ.

1,3,8-Triazaspiro[4.5]decane-2,4-diones as efficacious pan-inhibitors of hypoxia-inducible factor prolyl hydroxylase 1-3 (HIF PHD1-3) for the treatment of anemia.

The discovery of 1,3,8-triazaspiro[4.5]decane-2,4-diones (spirohydantoins) as a structural class of pan-inhibitors of the prolyl hydroxylase (PHD) family of enzymes for the treatment of anemia is described. The initial hit class, spirooxindoles, was identified through affinity selection mass spectrometry (AS-MS) and optimized for PHD2 inhibition and optimal PK/PD profile (short-acting PHDi inhibitors). 1,3,8-Triazaspiro[4.5]decane-2,4-diones (spirohydantoins) were optimized as an advanced lead class derived from the original spiroindole hit. A new set of general conditions for C-N coupling, developed using a high-throughput experimentation (HTE) technique, enabled a full SAR analysis of the spirohydantoins. This rapid and directed SAR exploration has resulted in the first reported examples of hydantoin derivatives with good PK in preclinical species. Potassium channel off-target activity (hERG) was successfully eliminated through the systematic introduction of acidic functionality to the molecular structure. Undesired upregulation of alanine aminotransferese (ALT) liver enzymes was mitigated and a robust on-/off-target margin was achieved. Spirohydantoins represent a class of highly efficacious, short-acting PHD1-3 inhibitors causing a robust erythropoietin (EPO) upregulation in vivo in multiple preclinical species. This profile deems spirohydantoins as attractive short-acting PHDi inhibitors with the potential for treatment of anemia.

Efficient and accurate characterization of the Bergman cyclization for several enediynes including an expanded substructure of esperamicin A1.

Incorporation of enediynes into anticancer drugs remains an intriguing yet elusive strategy for the design of therapeutically active agents. Density functional theory was used to locate reactants, products, and transition states along the Bergman cyclization pathways connecting enediynes to reactive para-biradicals. Sum method correction to low-level calculations confirmed B3LYP/6-31G(d,p) as the method of choice in investigating enediynes. Herein described as MI:Sum, calculated reaction enthalpies differed from experiment by an average of 2.1 kcal x mol(-1) (mean unsigned error). A combination of strain energy released across the reaction coordinate and the critical intramolecular distance between reacting diynes explains reactivity differences. Where experimental and calculated barrier heights are in disagreement, higher level multireference treatment of the enediynes confirms lower level estimates. Previous work concerning the chemically reactive fragment of esperamcin, MTC, is expanded to our model system MTC2.