Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify.

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.

Active Learning Exploration of Transition-Metal Complexes to Discover Method-Insensitive and Synthetically Accessible Chromophores.

Transition-metal chromophores with earth-abundant transition metals are an important design target for their applications in lighting and nontoxic bioimaging, but their design is challenged by the scarcity of complexes that simultaneously have well-defined ground states and optimal target absorption energies in the visible region. Machine learning (ML) accelerated discovery could overcome such challenges by enabling the screening of a larger space but is limited by the fidelity of the data used in ML model training, which is typically from a single approximate density functional. To address this limitation, we search for consensus in predictions among 23 density functional approximations across multiple rungs of "Jacob's ladder". To accelerate the discovery of complexes with absorption energies in the visible region while minimizing the effect of low-lying excited states, we use two-dimensional (2D)efficient global optimization to sample candidate low-spin chromophores from multimillion complex spaces. Despite the scarcity (i.e., ∼0.01%) of potential chromophores in this large chemical space, we identify candidates with high likelihood (i.e., >10%) of computational validation as the ML models improve during active learning, representing a 1000-fold acceleration in discovery. Absorption spectra of promising chromophores from time-dependent density functional theory verify that 2/3 of candidates have the desired excited-state properties. The observation that constituent ligands from our leads have demonstrated interesting optical properties in the literature exemplifies the effectiveness of our construction of a realistic design space and active learning approach.

Emergence of a proton exchange-based isomerization and lactonization mechanism in the plant coumarin synthase COSY.

Plants contain rapidly evolving specialized enzymes that support the biosynthesis of functionally diverse natural products. In coumarin biosynthesis, a BAHD acyltransferase-family enzyme COSY was recently discovered to accelerate coumarin formation as the only known BAHD enzyme to catalyze an intramolecular acyl transfer reaction. Here we investigate the structural and mechanistic basis for COSY's coumarin synthase activity. Our structural analyses reveal an unconventional active-site configuration adapted to COSY's specialized activity. Through mutagenesis studies and deuterium exchange experiments, we identify a unique proton exchange mechanism at the α-carbon of the o-hydroxylated trans-hydroxycinnamoyl-CoA substrates during the catalytic cycle of COSY. Quantum mechanical cluster modeling and molecular dynamics further support this key mechanism for lowering the activation energy of the rate-limiting trans-to-cis isomerization step in coumarin production. This study unveils an unconventional catalytic mechanism mediated by a BAHD-family enzyme, and sheds light on COSY's evolutionary origin and its recruitment to coumarin biosynthesis in eudicots.

Influence of the Greater Protein Environment on the Electrostatic Potential in Metalloenzyme Active Sites: The Case of Formate Dehydrogenase.

The Mo/W-containing metalloenzyme formate dehydrogenase (FDH) is an efficient and selective natural catalyst that reversibly converts CO2 to formate under ambient conditions. In this study, we investigate the impact of the greater protein environment on the electrostatic potential (ESP) of the active site. To model the enzyme environment, we used a combination of classical molecular dynamics and multiscale quantum-mechanical (QM)/molecular-mechanical (MM) simulations. We leverage charge shift analysis to systematically construct QM regions and analyze the electronic environment of the active site by evaluating the degree of charge transfer between the core active site and the protein environment. The contribution of the terminal chalcogen ligand to the ESP of the metal center is substantial and dependent on the chalcogen identity, with similar, less negative ESPs for Se and S terminal chalcogens in comparison to O regardless of whether the metal is Mo or W. The orientation of the side chains and conformations of the cofactor also affect the ESP, highlighting the importance of sampling dynamic fluctuations in the protein. Overall, our observations suggest that the terminal chalcogen ligand identity plays an important role in the enzymatic activity of FDH, suggesting opportunities for a rational bioinspired catalyst design.

Synthesis and evaluation of potent yaku'amide A analogs.

Two full-length analogs of the anticancer peptide yaku'amide A (1a) and four partial structures have been synthesized. These analogs were identified by computational studies in which the three E- and Z-ΔIle residues of the natural product were replaced by the more accessible dehydroamino acids ΔVal and ΔEnv. Of the eight possible analogs, modeling showed that the targeted structures 2a and 2b most closely resembled the three-dimensional structure of 1a. Synthesis of 2a and 2b followed a convergent route that was streamlined by the absence of ΔIle in the targets. Screening of the compounds against various cancer cell lines revealed that 2a and 2b mimic the potent anticancer activity of 1a, thereby validating the computational studies.

MOFSimplify, machine learning models with extracted stability data of three thousand metal-organic frameworks.

We report a workflow and the output of a natural language processing (NLP)-based procedure to mine the extant metal-organic framework (MOF) literature describing structurally characterized MOFs and their solvent removal and thermal stabilities. We obtain over 2,000 solvent removal stability measures from text mining and 3,000 thermal decomposition temperatures from thermogravimetric analysis data. We assess the validity of our NLP methods and the accuracy of our extracted data by comparing to a hand-labeled subset. Machine learning (ML, i.e. artificial neural network) models trained on this data using graph- and pore-geometry-based representations enable prediction of stability on new MOFs with quantified uncertainty. Our web interface, MOFSimplify, provides users access to our curated data and enables them to harness that data for predictions on new MOFs. MOFSimplify also encourages community feedback on existing data and on ML model predictions for community-based active learning for improved MOF stability models.

Probing the Mechanism of Isonitrile Formation by a Non-Heme Iron(II)-Dependent Oxidase/Decarboxylase.

The isonitrile moiety is an electron-rich functionality that decorates various bioactive natural products isolated from diverse kingdoms of life. Isonitrile biosynthesis was restricted for over a decade to isonitrile synthases, a family of enzymes catalyzing a condensation reaction between l-Trp/l-Tyr and ribulose-5-phosphate. The discovery of ScoE, a non-heme iron(II) and α-ketoglutarate-dependent dioxygenase, demonstrated an alternative pathway employed by nature for isonitrile installation. Biochemical, crystallographic, and computational investigations of ScoE have previously been reported, yet the isonitrile formation mechanism remains obscure. In the present work, we employed in vitro biochemistry, chemical synthesis, spectroscopy techniques, and computational simulations that enabled us to propose a plausible molecular mechanism for isonitrile formation. Our findings demonstrate that the ScoE reaction initiates with C5 hydroxylation of (R)-3-((carboxymethyl)amino)butanoic acid to generate 1, which undergoes dehydration, presumably mediated by Tyr96 to synthesize 2 in a trans configuration. (R)-3-isocyanobutanoic acid is finally generated through radical-based decarboxylation of 2, instead of the common hydroxylation pathway employed by this enzyme superfamily.

ASCL1 represses a SOX9+ neural crest stem-like state in small cell lung cancer.

ASCL1 is a neuroendocrine lineage-specific oncogenic driver of small cell lung cancer (SCLC), highly expressed in a significant fraction of tumors. However, ∼25% of human SCLC are ASCL1-low and associated with low neuroendocrine fate and high MYC expression. Using genetically engineered mouse models (GEMMs), we show that alterations in Rb1/Trp53/Myc in the mouse lung induce an ASCL1+ state of SCLC in multiple cells of origin. Genetic depletion of ASCL1 in MYC-driven SCLC dramatically inhibits tumor initiation and progression to the NEUROD1+ subtype of SCLC. Surprisingly, ASCL1 loss promotes a SOX9+ mesenchymal/neural crest stem-like state and the emergence of osteosarcoma and chondroid tumors, whose propensity is impacted by cell of origin. ASCL1 is critical for expression of key lineage-related transcription factors NKX2-1, FOXA2, and INSM1 and represses genes involved in the Hippo/Wnt/Notch developmental pathways in vivo. Importantly, ASCL1 represses a SOX9/RUNX1/RUNX2 program in vivo and SOX9 expression in human SCLC cells, suggesting a conserved function for ASCL1. Together, in a MYC-driven SCLC model, ASCL1 promotes neuroendocrine fate and represses the emergence of a SOX9+ nonendodermal stem-like fate that resembles neural crest.

MYC Drives Temporal Evolution of Small Cell Lung Cancer Subtypes by Reprogramming Neuroendocrine Fate.

Small cell lung cancer (SCLC) is a neuroendocrine tumor treated clinically as a single disease with poor outcomes. Distinct SCLC molecular subtypes have been defined based on expression of ASCL1, NEUROD1, POU2F3, or YAP1. Here, we use mouse and human models with a time-series single-cell transcriptome analysis to reveal that MYC drives dynamic evolution of SCLC subtypes. In neuroendocrine cells, MYC activates Notch to dedifferentiate tumor cells, promoting a temporal shift in SCLC from ASCL1+ to NEUROD1+ to YAP1+ states. MYC alternatively promotes POU2F3+ tumors from a distinct cell type. Human SCLC exhibits intratumoral subtype heterogeneity, suggesting that this dynamic evolution occurs in patient tumors. These findings suggest that genetics, cell of origin, and tumor cell plasticity determine SCLC subtype.

Impact of Dehydroamino Acids on the Structure and Stability of Incipient 310-Helical Peptides.

A comparative study of the impact of small, medium-sized, and bulky α,ß-dehydroamino acids (ΔAAs) on the structure and stability of Balaram's incipient 310-helical peptide (1) is reported. Replacement of the N-terminal Aib residue of 1 with a ΔAA afforded peptides 2a-c that maintained the 310-helical shape of 1. In contrast, installation of a ΔAA in place of Aib-3 yielded peptides 3a-c that preferred a ß-sheet-like conformation. The impact of the ΔAA on peptide structure was independent of size, with small (ΔAla), medium-sized (Z-ΔAbu), and bulky (ΔVal) ΔAAs exerting similar effects. The proteolytic stabilities of 1 and its analogs were determined by incubation with Pronase. Z-ΔAbu and ΔVal increased the resistance of peptides to proteolysis when incorporated at the 3-position and had negligible impact on stability when placed at the 1-position, whereas ΔAla-containing peptides degraded rapidly regardless of position. Exposure of peptides 2a-c and 3a-c to the reactive thiol cysteamine revealed that ΔAla-containing peptides underwent conjugate addition at room temperature, while Z-ΔAbu- and ΔVal-containing peptides were inert even at elevated temperatures. These results suggest that both bulky and more accessible medium-sized ΔAAs should be valuable tools for bestowing rigidity and proteolytic stability on bioactive peptides.

Humanin induces conformational changes in the apoptosis regulator BAX and sequesters it into fibers, preventing mitochondrial outer-membrane permeabilization.

The mitochondrial, or intrinsic, apoptosis pathway is regulated mainly by members of the B-cell lymphoma 2 (BCL-2) protein family. BCL-2-associated X apoptosis regulator (BAX) plays a pivotal role in the initiation of mitochondria-mediated apoptosis as one of the factors causing mitochondrial outer-membrane permeabilization (MOMP). Of current interest are endogenous BAX ligands that inhibit its MOMP activity. Mitochondrial-derived peptides (MDPs) are a recently identified class of mitochondrial retrograde signaling molecules and are reported to be potent apoptosis inhibitors. Among them, humanin (HN) has been shown to suppress apoptosis by inhibiting BAX translocation to the mitochondrial outer membrane, but the molecular mechanism of this interaction is unknown. Here, using recombinant protein expression, along with light-scattering, CD, and fluorescence spectroscopy, we report that HN and BAX can form fibers together in vitro Results from negative stain EM experiments suggest that BAX undergoes secondary and tertiary structural rearrangements and incorporates into the fibers, and that its membrane-associating C-terminal helix is important for the fibrillation process. Additionally, HN mutations known to alter its anti-apoptotic activity affect fiber morphology. Our findings reveal for the first time a potential mechanism by which BAX can be sequestered by fibril formation, which can prevent it from initiating MOMP and committing the cell to apoptosis.

Bulky Dehydroamino Acids Enhance Proteolytic Stability and Folding in ß-Hairpin Peptides.

The bulky dehydroamino acids dehydrovaline (ΔVal) and dehydroethylnorvaline (ΔEnv) can be inserted into the turn regions of ß-hairpin peptides without altering their secondary structures. These residues increase proteolytic stability, with ΔVal at the (i + 1) position having the most substantial impact. Additionally, a bulky dehydroamino acid can be paired with a d-amino acid (i.e., d-Pro) to synergistically enhance resistance to proteolysis. A link between proteolytic stability and peptide structure is established by the finding that a stabilized ΔVal-containing ß-hairpin is more highly folded than its Asn-containing congener.

Potential involvement of mi-RNA 574-3p in progression of prostate cancer: A bioinformatic study.

Aberrant gene expression is a hallmark of prostate cancer (PCa), the second deadliest disease affecting males worldwide. Dysregulation of miRNA has been associated with the progression of PCa and in recent studies, miRNA 574-3p was found to be upregulated in cancerous prostate tissue. In this study, we characterize the effects of upregulated miRNA 574-3p on gene expression in the tumor microenvironment through different bioinformatic tools such as Diana-Tools, the KEGG Pathway Database, and the Reactome Database. We have identified nine regulatory genes that are targeted by miRNA 574-3p and downregulated in prostate cells. Pathway analysis of these genes shows that they are involved in the regulation of the Notch signaling pathway, Wnt signaling pathway, apoptosis, DNA damage response, G1 to S cell-cycle control, inflammatory response pathway, angiogenesis, translation factors, and the expression of oncogenes. Our results show the oncogenic potential of miRNA 574-3p in PCa progression and metastasis. Moreover, this study highlights the complex molecular mechanisms and pathways affected by the upregulation of miRNA 574-3p in prostate cells. In future studies, the presented data may aid in designing new therapies for PCa with improved efficacy.