14 examples of how LLMs can transform materials science and chemistry: a reflection on a large language model hackathon.

Large-language models (LLMs) such as GPT-4 caught the interest of many scientists. Recent studies suggested that these models could be useful in chemistry and materials science. To explore these possibilities, we organized a hackathon. This article chronicles the projects built as part of this hackathon. Participants employed LLMs for various applications, including predicting properties of molecules and materials, designing novel interfaces for tools, extracting knowledge from unstructured data, and developing new educational applications. The diverse topics and the fact that working prototypes could be generated in less than two days highlight that LLMs will profoundly impact the future of our fields. The rich collection of ideas and projects also indicates that the applications of LLMs are not limited to materials science and chemistry but offer potential benefits to a wide range of scientific disciplines.

Role of Steric Effects on Rates of Hydrogen Atom Transfer Reactions.

The influence of steric effects on the rates of hydrogen atom transfer (HAT) reactions between oxyradicals and alkanes is explored computationally. Quantum chemical density functional theory computations of transition states show that activation barriers and reaction enthalpies are both influenced by bulky substituents on the radical but very little by substituents on the alkane. The activation barriers remain roughly correlated with reaction enthalpies via the Evans-Polanyi relationship even when steric repulsion effects become important, although dispersion effects sometimes stabilize transition states. By making comparisons to previously developed Evans-Polanyi and modified Roberts-Steel relationships, we find that HAT reactions between bulky molecules remain well-described by these relationships.

Enzymatic Benzofuranoindoline Formation in the Biosynthesis of the Strained Bridgehead Bicyclic Dipeptide (+)-Azonazine A.

We uncovered and reconstituted a concise biosynthetic pathway of the strained dipeptide (+)-azonazine A from marine-derived Aspergillus insulicola. Formation of the hexacyclic benzofuranoindoline ring system from cyclo-(l-Trp-N-methyl-l-Tyr) is catalyzed by a P450 enzyme through an oxidative cyclization. Supplementing the producing strain with various indole-substituted tryptophan derivatives resulted in the generation of a series of azonazine A analogs.

Mitotic clustering of pulverized chromosomes from micronuclei.

Complex genome rearrangements can be generated by the catastrophic pulverization of missegregated chromosomes trapped within micronuclei through a process known as chromothripsis1-5. As each chromosome contains a single centromere, it remains unclear how acentric fragments derived from shattered chromosomes are inherited between daughter cells during mitosis6. Here we tracked micronucleated chromosomes with live-cell imaging and show that acentric fragments cluster in close spatial proximity throughout mitosis for asymmetric inheritance by a single daughter cell. Mechanistically, the CIP2A-TOPBP1 complex prematurely associates with DNA lesions within ruptured micronuclei during interphase, which poises pulverized chromosomes for clustering upon mitotic entry. Inactivation of CIP2A-TOPBP1 caused acentric fragments to disperse throughout the mitotic cytoplasm, stochastically partition into the nucleus of both daughter cells and aberrantly misaccumulate as cytoplasmic DNA. Mitotic clustering facilitates the reassembly of acentric fragments into rearranged chromosomes lacking the extensive DNA copy-number losses that are characteristic of canonical chromothripsis. Comprehensive analysis of pan-cancer genomes revealed clusters of DNA copy-number-neutral rearrangements-termed balanced chromothripsis-across diverse tumour types resulting in the acquisition of known cancer driver events. Thus, distinct patterns of chromothripsis can be explained by the spatial clustering of pulverized chromosomes from micronuclei.

Structural variation cooperates with permissive chromatin to control enhancer hijacking-mediated oncogenic transcription.

Structural variants (SVs) involving enhancer hijacking can rewire chromatin topologies to cause oncogene activation in human cancers, including hematologic malignancies; however, because of the lack of tools to assess their effects on gene regulation and chromatin organization, the molecular determinants for the functional output of enhancer hijacking remain poorly understood. Here, we developed a multimodal approach to integrate genome sequencing, chromosome conformation, chromatin state, and transcriptomic alteration for quantitative analysis of transcriptional effects and structural reorganization imposed by SVs in leukemic genomes. We identified known and new pathogenic SVs, including recurrent t(5;14) translocations that cause the hijacking of BCL11B enhancers for the allele-specific activation of TLX3 in a subtype of pediatric leukemia. Epigenetic perturbation of SV-hijacked BCL11B enhancers impairs TLX3 transcription, which are required for the growth of t(5;14) leukemia cells. By CRISPR engineering of patient-derived t(5;14) in isogenic leukemia cells, we uncovered a new mechanism whereby the transcriptional output of SV-induced BCL11B enhancer hijacking is dependent on the loss of DNA hypermethylation at the TLX3 promoter. Our results highlight the importance of the cooperation between genetic alteration and permissive chromatin as a critical determinant of SV-mediated oncogene activation, with implications for understanding aberrant gene transcription after epigenetic therapies in patients with leukemia. Hence, leveraging the interdependency of genetic alteration on chromatin variation may provide new opportunities to reprogram gene regulation as targeted interventions in human disease.

Light-activated macromolecular phase separation modulates transcription by reconfiguring chromatin interactions.

Biomolecular condensates participate in the regulation of gene transcription, yet the relationship between nuclear condensation and transcriptional activation remains elusive. Here, we devised a biotinylated CRISPR-dCas9-based optogenetic method, light-activated macromolecular phase separation (LAMPS), to enable inducible formation, affinity purification, and multiomic dissection of nuclear condensates at the targeted genomic loci. LAMPS-induced condensation at enhancers and promoters activates endogenous gene transcription by chromatin reconfiguration, causing increased chromatin accessibility and de novo formation of long-range chromosomal loops. Proteomic profiling of light-induced condensates by dCas9-mediated affinity purification uncovers multivalent interaction-dependent remodeling of macromolecular composition, resulting in the selective enrichment of transcriptional coactivators and chromatin structure proteins. Our findings support a model whereby the formation of nuclear condensates at native genomic loci reconfigures chromatin architecture and multiprotein assemblies to modulate gene transcription. Hence, LAMPS facilitates mechanistic interrogation of the relationship between nuclear condensation, genome structure, and gene transcription in living cells.

Initiation of fatty acid biosynthesis in Pseudomonas putida KT2440.

Deciphering the mechanisms of bacterial fatty acid biosynthesis is crucial for both the engineering of bacterial hosts to produce fatty acid-derived molecules and the development of new antibiotics. However, gaps in our understanding of the initiation of fatty acid biosynthesis remain. Here, we demonstrate that the industrially relevant microbe Pseudomonas putida KT2440 contains three distinct pathways to initiate fatty acid biosynthesis. The first two routes employ conventional ß-ketoacyl-ACP synthase III enzymes, FabH1 and FabH2, that accept short- and medium-chain-length acyl-CoAs, respectively. The third route utilizes a malonyl-ACP decarboxylase enzyme, MadB. A combination of exhaustive in vivo alanine-scanning mutagenesis, in vitro biochemical characterization, X-ray crystallography, and computational modeling elucidate the presumptive mechanism of malonyl-ACP decarboxylation via MadB. Given that functional homologs of MadB are widespread throughout domain Bacteria, this ubiquitous alternative fatty acid initiation pathway provides new opportunities to target a range of biotechnology and biomedical applications.

An NmrA-like enzyme-catalysed redox-mediated Diels-Alder cycloaddition with anti-selectivity.

The Diels-Alder cycloaddition is one of the most powerful approaches in organic synthesis and is often used in the synthesis of important pharmaceuticals. Yet, strictly controlling the stereoselectivity of the Diels-Alder reactions is challenging, and great efforts are needed to construct complex molecules with desired chirality via organocatalysis or transition-metal strategies. Nature has evolved different types of enzymes to exquisitely control cyclization stereochemistry; however, most of the reported Diels-Alderases have been shown to only facilitate the energetically favourable diastereoselective cycloadditions. Here we report the discovery and characterization of CtdP, a member of a new class of bifunctional oxidoreductase/Diels-Alderase, which was previously annotated as an NmrA-like transcriptional regulator. We demonstrate that CtdP catalyses the inherently disfavoured cycloaddition to form the bicyclo[2.2.2]diazaoctane scaffold with a strict α-anti-selectivity. Guided by computational studies, we reveal a NADP+/NADPH-dependent redox mechanism for the CtdP-catalysed inverse electron demand Diels-Alder cycloaddition, which serves as the first example of a bifunctional Diels-Alderase that utilizes this mechanism.

Molecular Dynamics Simulations Guide Chimeragenesis and Engineered Control of Chemoselectivity in Diketopiperazine Dimerases.

In the biosynthesis of the tryptophan-linked dimeric diketopiperazines (DKPs), cytochromes P450 selectively couple DKP monomers to generate a variety of intricate and isomeric frameworks. To determine the molecular basis for selectivity of these biocatalysts we obtained a high-resolution crystal structure of selective Csp2 -N bond forming dimerase, AspB. Overlay of the AspB structure onto C-C and C-N bond forming homolog NzeB revealed no significant structural variance to explain their divergent chemoselectivities. Molecular dynamics (MD) simulations identified a region of NzeB with increased conformational flexibility relative to AspB, and interchange of this region along with a single active site mutation led to a variant that catalyzes exclusive C-N bond formation. MD simulations also suggest that intermolecular C-C or C-N bond formation results from a change in mechanism, supported experimentally through use of a substrate mimic.

Transcriptional profiling of Hutchinson-Gilford Progeria syndrome fibroblasts reveals deficits in mesenchymal stem cell commitment to differentiation related to early events in endochondral ossification.

The expression of a mutant Lamin A, progerin, in Hutchinson-Gilford Progeria Syndrome leads to alterations in genome architecture, nuclear morphology, epigenetic states, and altered phenotypes in all cells of the mesenchymal lineage. Here, we report a comprehensive analysis of the transcriptional status of patient derived HGPS fibroblasts, including nine cell lines not previously reported, in comparison with age-matched controls, adults, and old adults. We find that Progeria fibroblasts carry abnormal transcriptional signatures, centering around several functional hubs: DNA maintenance and epigenetics, bone development and homeostasis, blood vessel maturation and development, fat deposition and lipid management, and processes related to muscle growth. Stratification of patients by age revealed misregulated expression of genes related to endochondral ossification and chondrogenic commitment in children aged 4-7 years old, where this differentiation program starts in earnest. Hi-C measurements on patient fibroblasts show weakening of genome compartmentalization strength but increases in TAD strength. While the majority of gene misregulation occurs in regions which do not change spatial chromosome organization, some expression changes in key mesenchymal lineage genes coincide with lamin associated domain misregulation and shifts in genome compartmentalization.

Loops, topologically associating domains, compartments, and territories are elastic and robust to dramatic nuclear volume swelling.

Layers of genome organization are becoming increasingly better characterized, but less is known about how these structures respond to perturbation or shape changes. Low-salt swelling of isolated chromatin fibers or nuclei has been used for decades to investigate the structural properties of chromatin. But, visible changes in chromatin appearance have not been linked to known building blocks of genome structure or features along the genome sequence. We combine low-salt swelling of isolated nuclei with genome-wide chromosome conformation capture (Hi-C) and imaging approaches to probe the effects of chromatin extension genome-wide. Photoconverted patterns on nuclei during expansion and contraction indicate that global genome structure is preserved after dramatic nuclear volume swelling, suggesting a highly elastic chromosome topology. Hi-C experiments before, during, and after nuclear swelling show changes in average contact probabilities at short length scales, reflecting the extension of the local chromatin fiber. But, surprisingly, during this large increase in nuclear volume, there is a striking maintenance of loops, TADs, active and inactive compartments, and chromosome territories. Subtle differences after expansion are observed, suggesting that the local chromatin state, protein interactions, and location in the nucleus can affect how strongly a given structure is maintained under stress. From these observations, we propose that genome topology is robust to extension of the chromatin fiber and isotropic shape change, and that this elasticity may be beneficial in physiological circumstances of changes in nuclear size and volume.

Chromosome compartmentalization alterations in prostate cancer cell lines model disease progression.

Prostate cancer aggressiveness and metastatic potential are influenced by gene expression and genomic aberrations, features that can be influenced by the 3D structure of chromosomes inside the nucleus. Using chromosome conformation capture (Hi-C), we conducted a systematic genome architecture comparison on a cohort of cell lines that model prostate cancer progression, from normal epithelium to bone metastasis. We describe spatial compartment identity (A-open versus B-closed) changes with progression in these cell lines and their relation to gene expression changes in both cell lines and patient samples. In particular, 48 gene clusters switch from the B to the A compartment, including androgen receptor, WNT5A, and CDK14. These switches are accompanied by changes in the structure, size, and boundaries of topologically associating domains (TADs). Further, compartment changes in chromosome 21 are exacerbated with progression and may explain, in part, the genesis of the TMPRSS2-ERG translocation. These results suggest that discrete 3D genome structure changes play a deleterious role in prostate cancer progression. .

Structural diversification of hapalindole and fischerindole natural products via cascade biocatalysis.

Hapalindoles and related compounds (ambiguines, fischerindoles, welwitindolinones) are a diverse class of indole alkaloid natural products. They are typically isolated from the Stigonemataceae order of cyanobacteria and possess a broad scope of biological activities. Recently the biosynthetic pathway for assembly of these metabolites has been elucidated. In order to generate the core ring system, L-tryptophan is converted into the cis-indole isonitrile subunit before being prenylated with geranyl pyrophosphate at the C-3 position. A class of cyclases (Stig) catalyzes a three-step process including a Cope rearrangement, 6-exo-trig cyclization and electrophilic aromatic substitution to create a polycyclic core. Formation of the initial alkaloid is followed by diverse late-stage tailoring reactions mediated by additional biosynthetic enzymes to give rise to the wide array of structural variations observed in this compound class. Herein, we demonstrate the versatility and utility of the Fam prenyltransferase and Stig cyclases toward core structural diversification of this family of indole alkaloids. Through synthesis of cis-indole isonitrile subunit derivatives, and aided by protein engineering and computational analysis, we have employed cascade biocatalysis to generate a range of derivatives, and gained insights into the basis for substrate flexibility in this system.

Overcoming Kinetic and Thermodynamic Challenges of Classic Cope Rearrangements.

Systematic evaluation of 1,5-dienes bearing 3,3-electron-withdrawing groups and 4-methylation results in the discovery of a Cope rearrangement for Meldrum's acid-containing substrates that have unexpectedly favorable kinetic and thermodynamic profiles. The protocol is quite general due to a concise and convergent synthesis from abundant starting materials. Furthermore, products with an embedded Meldrum's acid moiety are prepared, which, in turn, can yield complex amides under neutral conditions. We have now expanded the scope of the reductive Cope rearrangement, which, via chemoselective reduction, can promote thermodynamically unfavorable [3,3] sigmatropic rearrangements of 3,3-dicyano-1,5-dienes to form reduced Cope rearrangement products. The Cope rearrangement is found to be stereospecific and can yield enantioenriched building blocks when chiral, nonracemic 1,3-disubstituted allylic electrophiles are utilized. We expand further the use of Cope rearrangements for the synthesis of highly valuable building blocks for complex- and drug-like molecular synthesis.

A thermodynamic atlas of carbon redox chemical space.

Redox biochemistry plays a key role in the transduction of chemical energy in living systems. However, the compounds observed in metabolic redox reactions are a minuscule fraction of chemical space. It is not clear whether compounds that ended up being selected as metabolites display specific properties that distinguish them from nonbiological compounds. Here, we introduce a systematic approach for comparing the chemical space of all possible redox states of linear-chain carbon molecules to the corresponding metabolites that appear in biology. Using cheminformatics and quantum chemistry, we analyze the physicochemical and thermodynamic properties of the biological and nonbiological compounds. We find that, among all compounds, aldose sugars have the highest possible number of redox connections to other molecules. Metabolites are enriched in carboxylic acid functional groups and depleted of ketones and aldehydes and have higher solubility than nonbiological compounds. Upon constructing the energy landscape for the full chemical space as a function of pH and electron-donor potential, we find that metabolites tend to have lower Gibbs energies than nonbiological molecules. Finally, we generate Pourbaix phase diagrams that serve as a thermodynamic atlas to indicate which compounds are energy minima in redox chemical space across a set of pH values and electron-donor potentials. While escape from thermodynamic equilibrium toward kinetically driven states is a hallmark of life and its origin, we envision that a deeper quantitative understanding of the environment-dependent thermodynamic landscape of putative prebiotic molecules will provide a crucial reference for future origins-of-life models.

Radiation-induced DNA damage and repair effects on 3D genome organization.

The three-dimensional structure of chromosomes plays an important role in gene expression regulation and also influences the repair of radiation-induced DNA damage. Genomic aberrations that disrupt chromosome spatial domains can lead to diseases including cancer, but how the 3D genome structure responds to DNA damage is poorly understood. Here, we investigate the impact of DNA damage response and repair on 3D genome folding using Hi-C experiments on wild type cells and ataxia telangiectasia mutated (ATM) patient cells. We irradiate fibroblasts, lymphoblasts, and ATM-deficient fibroblasts with 5 Gy X-rays and perform Hi-C at 30 minutes, 24 hours, or 5 days after irradiation. We observe that 3D genome changes after irradiation are cell type-specific, with lymphoblastoid cells generally showing more contact changes than irradiated fibroblasts. However, all tested repair-proficient cell types exhibit an increased segregation of topologically associating domains (TADs). This TAD boundary strengthening after irradiation is not observed in ATM deficient fibroblasts and may indicate the presence of a mechanism to protect 3D genome structure integrity during DNA damage repair.

Computational Investigation of the Mechanism of Diels-Alderase PyrI4.

We studied the mechanisms of activation and stereoselectivity of a monofunctional Diels-Alderase (PyrI4)-catalyzed intramolecular Diels-Alder reaction that leads to formation of the key spiro-tetramate moiety in the biosynthesis of the pyrroindomycin family of natural products. Key activation effects of PyrI4 include acid catalysis and an induced-fit mechanism that cooperate with the unique "lid" feature of PyrI4 to stabilize the Diels-Alder transition state. PyrI4 enhances the intrinsic Diels-Alder stereoselectivity of the substrate and leads to stereospecific formation of the product.

A Qualitative Assessment of Patient Satisfaction with Radical Cystectomy for Bladder Cancer at a Single Institution: How Can We Improve?

PURPOSE: To evaluate patient satisfaction (with emphasis on preoperative education) with radical cystectomy for bladder cancer at our institution, the University of Missouri Hospital, qualitatively in order to identify specific areas where improvements can be made. MATERIALS AND METHODS: We developed a patient survey that used open-ended questions to identify positive and negative experiences that contributed to patient satisfaction. We administered the survey to radical cystectomy patients who met inclusion criteria and agreed to participate. We recorded, transcribed and qualitatively coded the responses. We identified four themes under which both positive and negative responses were placed, and constructed two diagrams to better illustrate contributors to patient experience and satisfaction. RESULTS: We identified 25 patients who met inclusion criteria. Of those, 13 participated in the survey. Regarding overall experience, 92.3% of patients rated their care as excellent or good. Regarding preoperative education, 76.9% of patients reported they definitely or somewhat received enough information on what to expect after surgery, and 76.9% definitely received enough guidance on how to care for themselves after surgery. From qualitative coding of patient responses to open-ended questions, we identified preoperative preparation, delivery of care, caregiver availability, and patient-centered care as themes that contributed positively and negatively to patient experience. CONCLUSION: Although the overall patient satisfaction could be perceived as high (92.3%), qualitative analysis revealed several areas where improvements can be made to improve patient experience with radical cystectomy at our institution. As previously expected, preoperative preparation was a contributor.

Mechanism of an Organocatalytic Cope Rearrangement Involving Iminium Intermediates: Elucidating the Role of Catalyst Ring Size.

The mechanism of the organocatalytic Cope rearrangement is elucidated through a combined computational and experimental approach. As reported previously, hydrazides catalyze the Cope rearrangement of 1,5-hexadiene-2-carboxaldehydes via iminium ion formation, and seven- and eight-membered ring catalysts are more active than smaller ring sizes. In the present work, quantum mechanical computations and kinetic isotope effect experiments demonstrate that the Cope rearrangement step, rather than iminium formation, is rate-limiting. The computations further explain how the hydrazide catalyst lowers the free-energy barrier of the Cope rearrangement via an associative transition state that is stabilized by enehydrazine character. The computations also explain the catalyst ring size effect, as larger hydrazide rings are able to accommodate optimal transition-state geometries that minimize the unfavorable lone-pair repulsion between neighboring nitrogen atoms and maximize the favorable hyperconjugative donation from each nitrogen atom into neighboring electron-poor sigma bonds, with the seven-membered catalyst achieving a nearly ideal transition-state geometry that is comparable to that of an unconstrained acyclic catalyst. Experimental kinetics studies support the computations, showing that the seven-membered and acyclic hydrazide catalysts react 10 times faster than the six-membered catalyst. Unraveling the mechanism of this reaction is an important step in understanding other reactions catalyzed by hydrazides, and explaining the ring size effect is critical because cyclic catalysts provide a constrained scaffold, enabling the development of asymmetric variants of these reactions.