Irreversible synthesis of an ultrastrong two-dimensional polymeric material.

Polymers that extend covalently in two dimensions have attracted recent attention1,2 as a means of combining the mechanical strength and in-plane energy conduction of conventional two-dimensional (2D) materials3,4 with the low densities, synthetic processability and organic composition of their one-dimensional counterparts. Efforts so far have proven successful in forms that do not allow full realization of these properties, such as polymerization at flat interfaces5,6 or fixation of monomers in immobilized lattices7-9. Another frequently employed synthetic approach is to introduce microscopic reversibility, at the cost of bond stability, to achieve 2D crystals after extensive error correction10,11. Here we demonstrate a homogenous 2D irreversible polycondensation that results in a covalently bonded 2D polymeric material that is chemically stable and highly processable. Further processing yields highly oriented, free-standing films that have a 2D elastic modulus and yield strength of 12.7 ± 3.8 gigapascals and 488 ± 57 megapascals, respectively. This synthetic route provides opportunities for 2D materials in applications ranging from composite structures to barrier coating materials.

Identification of novel compound heterozygous variants in the SLC30A7 (ZNT7) gene in two French brothers with stunted growth, testicular hypoplasia and bone marrow failure.

Zinc is an essential trace mineral. Dietary zinc deficiency results in stunted growth, skin lesions, hypogonadism and frequent infections in humans. Mice genetically lacking Slc30a7 suffer from mild zinc deficiency and are prone to development of prostate cancer and insulin resistance. Disease-causing variants or mutations in the human SLC30A7 (ZNT7) gene have not been previously reported. Here, we describe two-boy siblings from a French family with stunted growth, testicular hypoplasia and bone marrow failure. Exome sequencing revealed compound heterozygous variants in ZNT7 consisting of NM_133496.5:c.21dup; p.Asp8ArgfsTer3 and c.842 + 15 T > C inherited from their unaffected mother and father, respectively. The c.21dup variant led to a premature stop codon generated in exon 1 of the ZNT7 coding sequence. RNA-seq analysis demonstrated that the c.842 + 15 T > C variant resulted in a leaky mRNA splicing event generating a premature stop codon right after the splicing donor site of exon 8. Moreover, the expression of ZNT7 protein was remarkably reduced by 80-96% in the affected brothers compared to the control cells. These findings strongly suggest that biallelic variants in SLC30A7 should be considered as a cause of growth retardation, testicular hypoplasia and syndromic bone marrow failure.

Properties of Configurationally Stable Atropoenantiomers in Macrocyclic Natural Products and the Chrysophaentin Family.

development of antibiotics, antineoplastics, and therapeutics for other diseases. Natural products are unique among all other small molecules in that they are produced by dedicated enzymatic assembly lines that are the protein products of biosynthetic gene clusters. As the products of chiral macromolecules, natural products have distinct three-dimensional shapes and stereochemistry is often encoded in their structures through the presence of stereocenters, or in the case of molecules that lack a stereocenter, the presence of an axis or plane of chirality. In the latter forms of chirality, if the barrier to rotation about the chiral axis or chiral plane is sufficiently high, stable conformers may exist allowing for isolation of discrete conformers, also known as atropisomers. Importantly, the diverse functions and biological activities of natural products are contingent upon their structures, stereochemistry and molecular shape. With continued innovation in methods for natural products discovery, synthetic chemistry, and analytical and computational tools, new insights into atropisomerism in natural products and related scaffolds are being made. As molecular complexity increases, more than one form of stereoisomerism may exist in a single compound (for example, point chirality, chiral axes, and chiral planes), sometimes creating atypical or noncanonical atropisomers, a term used to distinguish physically noninterconvertable atropisomers from typical atropisomers.Here we provide an account of the discovery and unusual structural and stereochemical features of the chrysophaentins, algal derived inhibitors of the bacterial cytoskeletal protein FtsZ and its associated protein partners. Eleven members of the chrysophaentin family have been discovered to date; seven of these are macrocyclic bis-bibenzyl ethers wherein the site of the ether linkage yields either a symmetrical or asymmetrical macrocyclic ring system. The asymmetrical ring system is highly strained and corresponds to the compounds having the most potent antimicrobial activity among the family. We review the structure elucidation and NMR properties that indicate restricted rotation between axes of two biaryl ethers, and the plane represented by the substituted 2-Z-butene bridge common to all of the macrocycles. Computational studies that corroborate high barriers to rotation about one representative plane, on the order of 20+ kcal/mol are presented. These barriers to rotation fix the conformation of the macrocycle into a bowl-like structure and suggest that an atropisomer should exist. Experimental evidence for atropisomerism is presented, consistent with computational predictions. These properties are discussed in the context of the total synthesis of 9-dechlorochrysophaenin A and its ring C isomers. Last, we discuss the implications for the presence of enantiomers in the biological activity and macrocyclization of the natural product.

Five Hypotheses on the Origins of Temperature Dependence of 77Se NMR Shifts in Diselenides.

Notable thermal shifts in diselenides have been documented in 77Se NMR for more than 50 years, but no satisfactory explanation has been found. Here, five hypotheses are considered as possible explanations for the large temperature dependence of the 77Se chemical shifts of diaryl and dialkyl diselenides compared to monoselenides and selenols. Density functional theory calculations are provided to bolster hypotheses and better understand the effects of barrier height and dipole energies. It is proposed that the temperature dependence of diselenide 77Se NMR chemical shifts is due to rotation around the Se-Se bond and sampling of twisted conformers at higher temperatures. The molecular twisting is solvent dependent; here, DMSO-d6 and toluene-d8 were evaluated. No correlation was established between para-substituents on diaryl diselenides and the magnitude of the change in the 77Se NMR shift (Δδ) with temperature.

EnzyHTP Computational Directed Evolution with Adaptive Resource Allocation.

Directed evolution facilitates enzyme engineering via iterative rounds of mutagenesis. Despite the wide applications of high-throughput screening, building "smart libraries" to effectively identify beneficial variants remains a major challenge in the community. Here, we developed a new computational directed evolution protocol based on EnzyHTP, a software that we have previously reported to automate enzyme modeling. To enhance the throughput efficiency, we implemented an adaptive resource allocation strategy that dynamically allocates different types of computing resources (e.g., GPU/CPU) based on the specific need of an enzyme modeling subtask in the workflow. We implemented the strategy as a Python library and tested the library using fluoroacetate dehalogenase as a model enzyme. The results show that compared to fixed resource allocation where both CPU and GPU are on-call for use during the entire workflow, applying adaptive resource allocation can save 87% CPU hours and 14% GPU hours. Furthermore, we constructed a computational directed evolution protocol under the framework of adaptive resource allocation. The workflow was tested against two rounds of mutational screening in the directed evolution experiments of Kemp eliminase (KE07) with a total of 184 mutants. Using folding stability and electrostatic stabilization energy as computational readout, we identified all four experimentally observed target variants. Enabled by the workflow, the entire computation task (i.e., 18.4 µs MD and 18,400 QM single-point calculations) completes in 3 days of wall-clock time using ∼30 GPUs and ∼1000 CPUs.

LassoHTP: A High-Throughput Computational Tool for Lasso Peptide Structure Construction and Modeling.

Lasso peptides are a subclass of ribosomally synthesized and post-translationally modified peptides with a slipknot conformation. With superior thermal stability, protease resistance, and antimicrobial activity, lasso peptides are promising candidates for bioengineering and pharmaceutical applications. To enable high-throughput computational prediction and design of lasso peptides, we developed a software, LassoHTP, for automatic lasso peptide structure construction and modeling. LassoHTP consists of three modules, including the scaffold constructor, mutant generator, and molecular dynamics (MD) simulator. With a user-provided sequence and conformational annotation, LassoHTP can either generate the structure and conformational ensemble as is or conduct random mutagenesis. We used LassoHTP to construct eight known lasso peptide structures de novo and to simulate their conformational ensembles for 100 ns MD simulations. For benchmarking, we calculated the root mean square deviation (RMSD) of these ensembles with reference to their experimental crystal or NMR PDB structures; we also compared these RMSD values against those of the MD ensembles that are initiated from the PDB structures. Dihedral principal component analysis was also conducted. The results show that the LassoHTP-initiated ensembles are similar to those of the PDB-initiated ensembles. LassoHTP offers a computational platform to develop strategies for lasso peptide prediction and design.

SAM-dependent enzyme-catalysed pericyclic reactions in natural product biosynthesis.

Pericyclic reactions-which proceed in a concerted fashion through a cyclic transition state-are among the most powerful synthetic transformations used to make multiple regioselective and stereoselective carbon-carbon bonds. They have been widely applied to the synthesis of biologically active complex natural products containing contiguous stereogenic carbon centres. Despite the prominence of pericyclic reactions in total synthesis, only three naturally existing enzymatic examples (the intramolecular Diels-Alder reaction, and the Cope and the Claisen rearrangements) have been characterized. Here we report a versatile S-adenosyl-l-methionine (SAM)-dependent enzyme, LepI, that can catalyse stereoselective dehydration followed by three pericyclic transformations: intramolecular Diels-Alder and hetero-Diels-Alder reactions via a single ambimodal transition state, and a retro-Claisen rearrangement. Together, these transformations lead to the formation of the dihydropyran core of the fungal natural product, leporin. Combined in vitro enzymatic characterization and computational studies provide insight into how LepI regulates these bifurcating biosynthetic reaction pathways by using SAM as the cofactor. These pathways converge to the desired biosynthetic end product via the (SAM-dependent) retro-Claisen rearrangement catalysed by LepI. We expect that more pericyclic biosynthetic enzymatic transformations remain to be discovered in naturally occurring enzyme 'toolboxes'. The new role of the versatile cofactor SAM is likely to be found in other examples of enzyme catalysis.

Ionization behavior of nanoporous polyamide membranes.

Escalating global water scarcity necessitates high-performance desalination membranes, for which fundamental understanding of structure-property-performance relationships is required. In this study, we comprehensively assess the ionization behavior of nanoporous polyamide selective layers in state-of-the-art nanofiltration (NF) membranes. In these films, residual carboxylic acids and amines influence permeability and selectivity by imparting hydrophilicity and ionizable moieties that can exclude coions. We utilize layered interfacial polymerization to prepare physically and chemically similar selective layers of controlled thickness. We then demonstrate location-dependent ionization of carboxyl groups in NF polyamide films. Specifically, only surface carboxyl groups ionize under neutral pH, whereas interior carboxyl ionization requires pH >9. Conversely, amine ionization behaves invariably across the film. First-principles simulations reveal that the low permittivity of nanoconfined water drives the anomalous carboxyl ionization behavior. Furthermore, we report that interior carboxyl ionization could improve the water-salt permselectivity of NF membranes over fourfold, suggesting that interior charge density could be an important tool to enhance the selectivity of polyamide membranes. Our findings highlight the influence of nanoconfinement on membrane transport properties and provide enhanced fundamental understanding of ionization that could enable novel membrane design.

Mechanistic Studies of a Skatole-Forming Glycyl Radical Enzyme Suggest Reaction Initiation via Hydrogen Atom Transfer.

Gut microbial decarboxylation of amino acid-derived arylacetates is a chemically challenging enzymatic transformation which generates small molecules that impact host physiology. The glycyl radical enzyme (GRE) indoleacetate decarboxylase from Olsenella uli (Ou IAD) performs the non-oxidative radical decarboxylation of indole-3-acetate (I3A) to yield skatole, a disease-associated metabolite produced in the guts of swine and ruminants. Despite the importance of IAD, our understanding of its mechanism is limited. Here, we characterize the mechanism of Ou IAD, evaluating previously proposed hypotheses of: (1) a Kolbe-type decarboxylation reaction involving an initial 1-e- oxidation of the carboxylate of I3A or (2) a hydrogen atom abstraction from the α-carbon of I3A to generate an initial carbon-centered radical. Site-directed mutagenesis, kinetic isotope effect experiments, analysis of reactions performed in D2O, and computational modeling are consistent with a mechanism involving initial hydrogen atom transfer. This finding expands the types of radical mechanisms employed by GRE decarboxylases and non-oxidative decarboxylases, more broadly. Elucidating the mechanism of IAD decarboxylation enhances our understanding of radical enzymes and may inform downstream efforts to modulate this disease-associated metabolism.

Self-Immolative Hydroxybenzylamine Linkers for Traceless Protein Modification.

Traceless self-immolative linkers are widely used for the reversible modification of proteins and peptides. This article describes a new class of traceless linkers based on ortho- or para-hydroxybenzylamines. The introduction of electron-donating substituents on the aromatic core stabilizes the quinone methide intermediate, thus providing a platform for payload release that can be modulated. To determine the extent to which the electronics affect the rate of release, we prepared a small library of hydroxybenzylamine linkers with varied electronics in the aromatic core, resulting in half-lives ranging from 20 to 144 h. Optimization of the linker design was carried out with mechanistic insights from density functional theory (DFT) and the in silico design of an intramolecular trapping agent through the use of DFT and intramolecular distortion energy calculations. This resulted in the development of a faster self-immolative linker with a half-life of 4.6 h. To demonstrate their effectiveness as traceless linkers for bioconjugation, reversible protein-polyethylene glycol conjugates with a model protein lysozyme were prepared, which had reduced protein activity but recovered ≥94% activity upon traceless release of the polymer. This new class of linkers with tunable release rates expands the traceless linkers toolbox for a variety of bioconjugation applications.

EnzyHTP: A High-Throughput Computational Platform for Enzyme Modeling.

Molecular simulations, including quantum mechanics (QM), molecular mechanics (MM), and multiscale QM/MM modeling, have been extensively applied to understand the mechanism of enzyme catalysis and to design new enzymes. However, molecular simulations typically require specialized, manual operation ranging from model construction to data analysis to complete the entire life cycle of enzyme modeling. The dependence on manual operation makes it challenging to simulate enzymes and enzyme variants in a high-throughput fashion. In this work, we developed a Python software, EnzyHTP, to automate molecular model construction, QM, MM, and QM/MM computation, and analyses of modeling data for enzyme simulations. To test the EnzyHTP, we used fluoroacetate dehalogenase (FAcD) as a model system and simulated the enzyme interior electrostatics for 100 FAcD mutants with a random single amino acid substitution. For each enzyme mutant, the workflow involves structural model construction, 1 ns molecular dynamics (MD) simulations, and quantum mechanical calculations in 100 MD-sampled snapshots. The entire simulation workflow for 100 mutants was completed in 7 h with 10 GPUs and 160 CPUs. EnzyHTP improves the efficiency of computational enzyme modeling, setting a basis for high-throughput identification of function-enhancing enzymes and enzyme variants. The software is expected to facilitate the fundamental understanding of catalytic origins across enzyme families and to accelerate the optimization of biocatalysts for non-native substrates.

IntEnzyDB: an Integrated Structure-Kinetics Enzymology Database.

Data-driven modeling has emerged as a new paradigm for biocatalyst design and discovery. Biocatalytic databases that integrate enzyme structure and function data are in urgent need. Here we describe IntEnzyDB as an integrated structure-kinetics database for facile statistical modeling and machine learning. IntEnzyDB employs a relational database architecture with a flattened data structure, which allows rapid data operation. This architecture also makes it easy for IntEnzyDB to incorporate more types of enzyme function data. IntEnzyDB contains enzyme kinetics and structure data from six enzyme commission classes. Using 1050 enzyme structure-kinetics pairs, we investigated the efficiency-perturbing propensities of mutations that are close or distal to the active site. The statistical results show that efficiency-enhancing mutations are globally encoded and that deleterious mutations are much more likely to occur in close mutations than in distal mutations. Finally, we describe a web interface that allows public users to access enzymology data stored in IntEnzyDB. IntEnzyDB will provide a computational facility for data-driven modeling in biocatalysis and molecular evolution.

Influence of Water and Enzyme on the Post-Transition State Bifurcation of NgnD-Catalyzed Ambimodal [6+4]/[4+2] Cycloaddition.

The enzyme NgnD catalyzes an ambimodal cycloaddition that bifurcates to [6+4]- and [4+2]-adducts. Both products have been isolated in experiments, but it remains unknown how enzyme and water influence the bifurcation selectivity at the femtosecond time scale. Here, we study the impact of water and enzyme on the post-transition state bifurcation of NgnD-catalyzed [6+4]/[4+2] cycloaddition by integrating quantum mechanics/molecular mechanics quasiclassical dynamics simulations and biochemical assays. The ratio of [6+4]/[4+2] products significantly differs in the gas phase, water, and enzyme. Biochemical assays were employed to validate computational predictions. The study informs how water and enzyme affect the bifurcation selectivity through perturbation of the reaction dynamics in the femtosecond time scale, revealing the fundamental roles of condensed media in dynamically controlling the chemical selectivity for biosynthetic reactions.

Influence of water and enzyme SpnF on the dynamics and energetics of the ambimodal [6+4]/[4+2] cycloaddition.

SpnF is the first monofunctional Diels-Alder/[6+4]-ase that catalyzes a reaction leading to both Diels-Alder and [6+4] adducts through a single transition state. The environment-perturbed transition-state sampling method has been developed to calculate free energies, kinetic isotope effects, and quasi-classical reaction trajectories of enzyme-catalyzed reactions and the uncatalyzed reaction in water. Energetics calculated in this way reproduce the experiment and show that the normal Diels-Alder transition state is stabilized by H bonds with water molecules, while the ambimodal transition state is favored in the enzyme SpnF by both intramolecular hydrogen bonding and hydrophobic binding. Molecular dynamics simulations show that trajectories passing through the ambimodal transition state bifurcate to the [6+4] adduct and the Diels-Alder adduct with a ratio of 1:1 in the gas phase, 1:1.6 in water, and 1:11 in the enzyme. This example shows how an enzyme acts on a vibrational time scale to steer post-transition state trajectories toward the Diels-Alder adduct.

Evolution of the Diels-Alder Reaction Mechanism since the 1930s: Woodward, Houk with Woodward, and the Influence of Computational Chemistry on Understanding Cycloadditions.

This review article describes the evolution of Woodward's mechanistic thinking, beginning in the late 1930s and early 1940s with his proposal of a charge-transfer mechanism for the Diels-Alder reaction, eventually leading to the Woodward-Katz two-stage concerted mechanism in 1959, and then to its mechanistic solution in terms of orbital symmetry control. Houk's research in the Woodward labs, testing the predictions of this theory, is described. Subsequent modern calculations with quantum mechanics and molecular dynamics simulations have shown that Woodward indeed had perfectly described not only the cyclopentadiene dimerization mechanism, but a new class of transition states now known as ambimodal or bis-pericyclic transition states. In recent years, the Houk group has found that ambimodal reactions are operative in the [6+4] cycloaddition. Molecular dynamics simulations of many Diels-Alder and ambimodal cycloadditions provide a time-resolved picture of how these reactions occur. Lastly, Roald Hoffmann provides a Coda in which he describes his joy in "being taken along the journey" of the cycloaddition story from Woodward's youth to today's trajectory simulations.

Ambimodal Trispericyclic Transition State and Dynamic Control of Periselectivity.

We report an ambimodal trispericyclic transition state leading to [6+4]-, [4+6]-, and [8+2]-cycloadducts in the reactions of 8,8-disubstituted heptafulvenes with 6,6-dimethylfulvene. The potential energy surfaces for these reactions were explored with ωB97X-D density functional theory. Quasi-classical direct molecular dynamics simulations gave information on the ratios of products expected in these reactions.

Carcinogenic Metabolic Activation Process of Naphthalene by the Cytochrome P450 Enzyme 1B1: A Computational Study.

The metabolic activation and transformation of naphthalene by the cytochrome P450 enzyme (CYP 1B1) plays an important role in its potential carcinogenicity. The process has been explored by a quantum mechanics/molecular mechanics (QM/MM) computational method. Molecular dynamic simulations were performed to explore the interaction between naphthalene and CYP 1B1. Naphthalene involves α- and ß-carbon, the electrophilic addition of which would result in different reaction pathways. Our computational results show that both additions on α- and ß-carbon can generate naphthalene 1,2-oxide. The activation barrier for the addition on ß-carbon is higher than that for the α-carbon by 2.6 kcal·mol-1, which is possibly caused by the proximity between ß-carbon and the iron-oxo group of Cpd I in the system. We also found that naphthalene 1,2-oxide is unstable and the O-C bond cleavage easily occurs via cellular hydronium ion, hydroxyl radical/anion; then it will convert to the potential ultimate carcinogen 1,2-naphthoquinone. The results demonstrate and inform a detailed process of generating naphthalene 1,2-oxide and new predictions for its conversion.

Relationships between Product Ratios in Ambimodal Pericyclic Reactions and Bond Lengths in Transition Structures.

Ambimodal reactions involve a single transition state leading to multiple products. In such reactions, transition state theory gives no information about the ratio of products that are formed, and molecular dynamics must be performed to predict this ratio. Understanding the relationship between the transition structure and the product ratio is a long-standing problem in molecular dynamics. We have studied 15 ambimodal pericyclic reactions and investigated the relationship between the TS bond lengths in the saddle points and the product ratios from trajectory simulations. A linear correlation, ln(B:A) = -9.4(Bond 3 - Bond 2), is found with R2 = 0.92, where A and B refer to the products formed upon formation of bonds 2 and 3, respectively. The correlation shows that the ratio of products formed after the bifurcation is related to the partial bond lengths, and corresponding bond orders, in the transition state.

The Dynamics of Chemical Reactions: Atomistic Visualizations of Organic Reactions, and Homage to van†'t Hoff.

Jacobus Henricus van 't Hoff was the first Nobel Laureate in Chemistry. He pioneered in the study of chemical dynamics, which referred at that time to chemical kinetics and thermodynamics. The term has evolved in modern times to refer to the exploration of chemical transformations in a time-resolved fashion. Chemical dynamics has been driven by the development of molecular dynamics trajectory simulations, which provide atomic visualization of chemical processes and illuminate how dynamic effects influence chemical reactivity and selectivity. In homage to the legend of van 't Hoff, we review the development of the chemical dynamics of organic reactions, our area of research. We then discuss our trajectory simulations of pericyclic reactions, and our development of dynamic criteria for concerted and stepwise reaction mechanisms. We also describe a method that we call environment-perturbed transition state sampling, which enables trajectory simulations in condensed-media using quantum mechanics and molecular mechanics (QM/MM). We apply the method to reactions in solvent and in enzyme. Jacobus Henricus van 't Hoff (1852, Rotterdam-1911, Berlin) received the Nobel Prize for Chemistry in 1901 "in recognition of the extraordinary services he has rendered by the discovery of the laws of chemical dynamics and osmotic pressure in solutions". van 't Hoff was born the Netherlands, and earned his doctorate in Utrecht in 1874. In 1896 he moved to Berlin, where he was offered a position with more research and less teaching. van 't Hoff is considered one of the founders of physical chemistry. A key step in establishing this new field was the start of Zeitschrift für Physikalische Chemie in 1887.

Theoretical study of radiative and nonradiative decay rates for Cu(i) complexes with double heteroleptic ligands.

In this paper, we conducted DFT and TDDFT calculations on three double heteroleptic Cu(i) complexes to understand how different substituents on N^N ligands influence the phosphorescence quantum yield (PLQY). Both radiative and nonradiative decay processes were thoroughly investigated. Factors that determine the rate of radiative process (kr) were considered, including the lowest triplet excited state E(T1), the transition dipole moment MSm,j of the Sm → S0 transition, the spin-coupled matrix element SOC, and the singlet-triplet splitting energies ΔE(Sm-T1). The results indicate that E(T1), MSm,j and SOC increase and ΔE(Sm-T1) decreases upon introducing -Ph and -CH2- groups on the N^N ligands. The net results lead to a gradual increase of kr in the three Cu(i) complexes, from 1 (0.48 × 104 s-1) to 2 (0.64 × 104 s-1) and then to 3 (1.61 × 104 s-1). The rate of nonradiative decay process (knr) was computed by a convolution method. We explored how knr is determined by SOC between T1 and S0 states (T1|SOC|S02), effective energy gap ΔE' and the Huang-Rhys factor (S). We found that T1|SOC|S02 and ΔE' contribute significantly to knr, but S does not determine the order of knr. knr gradually decreases from complex 1 (2.51 × 106 s-1) to 2 (0.32 × 106 s-1) and then to 3 (0.14 × 106 s-1) after introducing -Ph and -CH2- groups on the N^N ligands. The computed PLQYs for the three complexes are 1: 0.0019, 2: 0.0198, and 3: 0.1011. These are quantitatively consistent with the experimental observation (1: 0.0028, 2: 0.0061, and 3: 0.1000).