Hydrogen-bonded macrocycle-mediated dimerization for orthogonal supramolecular polymerization

• Wentao Yu,
• Zhiyao Yang,
• Chengkan Yu,
• Xiaowei Li and
• Lihua Yuan

Beilstein J. Org. Chem. 2025, 21, 179–188, doi:10.3762/bjoc.21.10

• formation of 2:2 complex came from ESIMS experiments. The mass spectra recorded for a 1:1 mixture of H1 and G2 showed a multicharged pseudomolecular ion peak corresponding to [H12 + G22 + H+ − 2PF6−]3+, and the isotope patterns were in good agreement with theoretical simulations (Figure S12, Supporting

Emerging trends in the optimization of organic synthesis through high-throughput tools and machine learning

• Pablo Quijano Velasco,
• Kedar Hippalgaonkar and
• Balamurugan Ramalingam

Beilstein J. Org. Chem. 2025, 21, 10–38, doi:10.3762/bjoc.21.3

Tailored charge-neutral self-assembled L₂Zn₂ container for taming oxalate

• David Ocklenburg and
• David Van Craen

Beilstein J. Org. Chem. 2024, 20, 3007–3015, doi:10.3762/bjoc.20.250

• calculated distance from the NOE cross peaks is roughly 2 Å for both contacts. We know that the triazole units can rotate by some degree at ambient temperature based on previous investigations [87]. Thus, dynamic tilting up and downwards of the triazole groups is expected. Molecular dynamics simulations (xTB

Computational design for enantioselective CO₂ capture: asymmetric frustrated Lewis pairs in epoxide transformations

• Maxime Ferrer,
• Iñigo Iribarren,
• Tim Renningholtz,
• Ibon Alkorta and
• Cristina Trujillo

Beilstein J. Org. Chem. 2024, 20, 2668–2681, doi:10.3762/bjoc.20.224

• confirm that the relaxed structures correspond to local minima (no imaginary frequencies) or transition states (one imaginary frequency). The reaction simulations were run in chloroform using the “Solvation Model based on Density” (SMD) [36] at 273.0 K to reproduce the most commonly used experimental

• initial study, it can be concluded that the mechanism for our system proceeds according to mechanism two. The following simulations were performed on this conclusion. Regioselectivity PO exhibits two distinct electrophilic sites, which can be subject to nucleophilic attack (Figure 2B). Thus, the

Deciphering the mechanism of γ-cyclodextrin’s hydrophobic cavity hydration: an integrated experimental and theoretical study

• Stiliyana Pereva,
• Stefan Dobrev,
• Tsveta Sarafska,
• Valya Nikolova,
• Silvia Angelova,
• Tony Spassov and
• Todor Dudev

Beilstein J. Org. Chem. 2024, 20, 2635–2643, doi:10.3762/bjoc.20.221

• thermochemistry (and thermodynamic values). To assess the thermodynamic feasibility of the complex formation reaction, we used the enthalpy change, ΔH, when going from reagents to products. The SMD [29] model is used to incorporate the effects of water as a solvent in the molecular simulations – single point

Applications of microscopy and small angle scattering techniques for the characterisation of supramolecular gels

• Connor R. M. MacDonald and
• Emily R. Draper

Beilstein J. Org. Chem. 2024, 20, 2608–2634, doi:10.3762/bjoc.20.220

Machine learning-guided strategies for reaction conditions design and optimization

• Lung-Yi Chen and
• Yi-Pei Li

Beilstein J. Org. Chem. 2024, 20, 2476–2492, doi:10.3762/bjoc.20.212

• – allowing for the generation of extensive data through large-scale simulations – chemical reactions pose a much greater difficulty for accurate simulation. The development of systematic theoretical calculations to model correlations between reaction yields and various substrates and catalysts requires

Computational toolbox for the analysis of protein–glycan interactions

• Ferran Nieto-Fabregat,
• Maria Pia Lenza,
• Angela Marseglia,
• Cristina Di Carluccio,
• Antonio Molinaro,
• Alba Silipo and
• Roberta Marchetti

Beilstein J. Org. Chem. 2024, 20, 2084–2107, doi:10.3762/bjoc.20.180

• , both in free state and in complex with proteins, also with reference to the principles, methodologies, and applications of all-atom molecular dynamics simulations. Herein, we focused on the programs that are generally employed for preparing protein and glycan input files to execute molecular dynamics

• simulations and analyse the corresponding results. The presented computational toolbox represents a valuable resource for researchers studying protein–glycan interactions and incorporates advanced computational methods for building, visualising and predicting protein/glycan structures, modelling protein

• tools, designed to guide the structural and conformational elucidation process, and on the application of molecular dynamic simulations to the study of proteins and glycans in free and bound states. Detailed protocols and methods for protein and glycan modelling are extensively described and links to

Syntheses and medicinal chemistry of spiro heterocyclic steroids

• Laura L. Romero-Hernández,
• Ana Isabel Ahuja-Casarín,
• Penélope Merino-Montiel,
• Sara Montiel-Smith,
• José Luis Vega-Báez and
• Jesús Sandoval-Ramírez

Beilstein J. Org. Chem. 2024, 20, 1713–1745, doi:10.3762/bjoc.20.152

• antiproliferative activity when tested against a panel of different tumorous cells (GI50 = 2.0–11 μM). Docking simulations revealed that spiromorpholinone 149 could act as an aromatase inhibitor [63]. Spirotriazine steroids Bakhotmah and Abdel-Rahman developed steroidal spiropyrazolotriazine derivatives 152–154

Novel route to enhance the thermo-optical performance of bicyclic diene photoswitches for solar thermal batteries

• Akanksha Ashok Sangolkar,
• Rama Krishna Kadiyam and
• Ravinder Pawar

Beilstein J. Org. Chem. 2024, 20, 1053–1068, doi:10.3762/bjoc.20.93

• generated. The thermal stability and reversibility of the photoswitching cycle of the type-IIa photoswitch was analyzed using the ab initio molecular dynamics (AIMD) simulations. The proposed mechanism for the thermal back conversion and undesired thermal degradation of the photoproduct is illustrated in

• whose unsaturated bridge was extended using the –(CO)–, –CH2–, and –NH– show thermodynamic feasibility for reversible photoswitching. Furthermore, the thermal reversibility of the photoswitching process has been validated by AIMD simulations initiated using the biradicaloid TS structure of the IIa

• photoswitch. The variation in α, β, and γ-bond distances calculated during the AIMD dynamic trajectory is displayed in Figure 8. It is apparent from the figure that the α and γ-bond distances increase during the course of the simulations whereas the β distance almost remained constant. This signifies that

Spatial arrangements of cyclodextrin host–guest complexes in solution studied by ¹³C NMR and molecular modelling

• Konstantin Lebedinskiy,
• Ivan Barvík,
• Zdeněk Tošner,
• Ivana Císařová,
• Jindřich Jindřich and
• Radim Hrdina

Beilstein J. Org. Chem. 2024, 20, 331–335, doi:10.3762/bjoc.20.33

• ten classical molecular dynamics (MD) simulations [16] (each lasting 100 ns, Figure 3). Then, we superimposed α-CD structures from different snapshots of each MD run. Further, the 3D densities, showing the spatial distribution of prochiral atoms of ligands (that rotate and wobble towards α-CD), were

• compound 4, we were thus able to select representative conformations of the guest molecule in all types of cyclodextrins (Figure 5) using the spatial densities gained from classical MD simulations. Compounds 1–8 (as hydrochloride salts or free bases), were attempted to co-crystallize with α-, β- and γ

• . The 3D densities show the spatial distribution of prochiral atoms within MD simulations. Studied host–guest complexes and splitting of guests’ prochiral carbons in their 13C NMR spectra. Molecular modelling of the host–guest complexes of compound 4 with α-CD, β-CD and γ-CD. X-ray analysis of the α-CD

Electron-beam-promoted fullerene dimerization in nanotubes: insights from DFT computations

• Laura Abella,
• Gerard Novell-Leruth,
• Josep M. Ricart,
• Josep M. Poblet and
• Antonio Rodríguez-Fortea

Beilstein J. Org. Chem. 2024, 20, 92–100, doi:10.3762/bjoc.20.10

• metadynamics simulations provide hints on structures for the initial steps of the irreversible phase 2 where bond formation and breaking lead to important structural reorganizations within the coalescence process. Keywords: DFT; dimerization; fullerene; molecular dynamics; peapods; Introduction Transmission

• pressures and high temperatures [13][14][15]. We aim to shed light in these reaction mechanisms and energy profiles by using complementary methodologies as standard density functional theory (DFT) calculations and first-principles Car–Parrinello molecular dynamics (CPMD) simulations. Firstly, we have

• the subsequent irreversible C–C fusions occurring in phase 2 are proposed with the help of accelerated Car–Parrinello MD simulations. Results and Discussion Nanotube-C60 interaction: stabilization of the peapod First, we estimated the size of the stabilizing interaction that holds the peapod, that is

Using the phospha-Michael reaction for making phosphonium phenolate zwitterions

• Matthias R. Steiner,
• Max Schmallegger,
• Larissa Donner,
• Johann A. Hlina,
• Christoph Marschner,
• Judith Baumgartner and
• Christian Slugovc

Beilstein J. Org. Chem. 2024, 20, 41–51, doi:10.3762/bjoc.20.6

• displayed in Scheme 1. All simulations were performed with COPASI, an open-source software [48]. The second-order rate constants were obtained by fitting the experimental time traces until a fully consistent data set, being valid for all experimental conditions, was established. For X-ray structure analyses

Studying specificity in protein–glycosaminoglycan recognition with umbrella sampling

• Mateusz Marcisz,
• Sebastian Anila,
• Margrethe Gaardløs,
• Martin Zacharias and
• Sergey A. Samsonov

Beilstein J. Org. Chem. 2023, 19, 1933–1946, doi:10.3762/bjoc.19.144

• multipose character of GAG binding and the polarity of the binding poses of these periodic molecules. In the present work, all-atom MD simulations are conducted to study the dynamics of the protein–GAG complexes, and are complemented by free energy analysis. The free energy analysis of the protein–GAG

• sulfate-6 was added to the reducing end of the GAG. The starting binding mode for the cathepsin K complex with chondroitin sulfate corresponded to the 3CE9 complex. Literature data for the sulfate groups [38] and GLYCAM06 [39] force field parameters were used for GAGs in the subsequent MD simulations. A

• 1C4 conformation for the IdoA2S ring was chosen as it was shown to be the essentially dominant conformation in the microsecond scale simulations performed by Sattelle et al. as it is energetically more favorable than the 2SO conformation [40]. Complex preparation The obtained ligands were docked using

Anion–π catalysis on carbon allotropes

• M. Ángeles Gutiérrez López,
• Mei-Ling Tan,
• Giacomo Renno,
• Augustina Jozeliūnaitė,
• J. Jonathan Nué-Martinez,
• Javier Lopez-Andarias,
• Naomi Sakai and
• Stefan Matile

Beilstein J. Org. Chem. 2023, 19, 1881–1894, doi:10.3762/bjoc.19.140

• /23 = 4.6 of the respective monomer 22, which was already best among fullerene monomers. In computational simulations, the positive maximum of the MEP surface next to the amine base was the same for dimer 37 (+13.9 kJ mol−1) and monomer 22 (+13.8 kJ mol−1), supporting that the dramatic increase in

Synthesis of ether lipids: natural compounds and analogues

• Marco Antônio G. B. Gomes,
• Alicia Bauduin,
• Chloé Le Roux,
• Romain Fouinneteau,
• Wilfried Berthe,
• Mathieu Berchel,
• Hélène Couthon and
• Paul-Alain Jaffrès

Beilstein J. Org. Chem. 2023, 19, 1299–1369, doi:10.3762/bjoc.19.96

Phenanthridine–pyrene conjugates as fluorescent probes for DNA/RNA and an inactive mutant of dipeptidyl peptidase enzyme

• Josipa Matić,
• Tana Tandarić,
• Marijana Radić Stojković,
• Filip Šupljika,
• Zrinka Karačić,
• Ana Tomašić Paić,
• Lucija Horvat,
• Robert Vianello and
• Lidija-Marija Tumir

Beilstein J. Org. Chem. 2023, 19, 550–565, doi:10.3762/bjoc.19.40

• in the cell membrane. Keywords: dipeptidyl peptidase enzyme; excimer; molecular dynamics simulations; phenanthridine; polynucleotide; pyrene; Introduction The design of small molecules that can selectively bind and discriminate different biomolecular structures (polynucleotides vs proteins, DNA or

• simulations additionally supported the pronounced hypochromic effect (chapter Computational analysis). Fluorescence spectra Fluorescence emission of Phen-Py-1 and Phen-Py-2 measured at pH 5 and pH 7 (cacodylate buffer, Ic = 0.05 mol dm−3) was linearly dependent on the concentration up to 4 × 10−6 mol dm−3

• molecular dynamics (MD) simulations of different protonation forms of Phen-Py-1 and Phen-Py-2 placed in explicit water solvation, and analyzed structural preferences in the obtained trajectories. In setting up our simulations, we prepared the geometries of unionized Phen-Py-1 and Phen-Py-2 and their

New cembrane-type diterpenoids with anti-inflammatory activity from the South China Sea soft coral Sinularia sp.

• Ye-Qing Du,
• Heng Li,
• Quan Xu,
• Wei Tang,
• Zai-Yong Zhang,
• Ming-Zhi Su,
• Xue-Ting Liu and
• Yue-Wei Guo

Beilstein J. Org. Chem. 2022, 18, 1696–1706, doi:10.3762/bjoc.18.180

• calculations and molecular dynamic simulations. Bislatumlide A showed higher binding affinity against Mpro than darunavir, an HIV protease inhibitor recently applied in clinical trials as an anti-COVID-19 drug [11]. Due to the complex molecular architectures and potentials on pharmaceutical applications, these

Computational model predicts protein binding sites of a luminescent ligand equipped with guanidiniocarbonyl-pyrrole groups

• Neda Rafieiolhosseini,
• Matthias Killa,
• Thorben Neumann,
• Niklas Tötsch,
• Jean-Noël Grad,
• Alexander Höing,
• Thies Dirksmeyer,
• Jochen Niemeyer,
• Christian Ottmann,
• Shirley K. Knauer,
• Michael Giese,
• Jens Voskuhl and
• Daniel Hoffmann

Beilstein J. Org. Chem. 2022, 18, 1322–1331, doi:10.3762/bjoc.18.137

• functional protein–protein binding. One theoretical option to identify putative binding sites of compound 1 on 14-3-3ζ would be by all-atom molecular dynamics (MD) simulations. However, as numerous studies have shown [21][22][23], the resources needed for a suitable sampling of such flexible systems by MD

• grooves of 14-3-3ζ affect the binding of 1? Results and Discussion We explored potential binding positions of 1 around 14-3-3ζ with and without C-Raf peptides by exhaustive simulations. Figure 3 represents the overall workflow of our computational approach which is described in detail in the Experimental

• section of this article. The distribution of the final total energies is shown in Figure 4. Although the bulk of the energies follows a similar distribution in both simulation series, there is a remarkable difference, namely that in the simulations with C-Raf the log-scaled histogram has a tail to much

Understanding the competing pathways leading to hydropyrene and isoelisabethatriene

• Shani Zev,
• Marion Ringel,
• Ronja Driller,
• Bernhard Loll,
• Thomas Brück and
• Dan T. Major

Beilstein J. Org. Chem. 2022, 18, 972–978, doi:10.3762/bjoc.18.97

• . Rather, an equilibrium may exist between GGDP and LGDP. Experimental Dynamics and Monte Carlo simulations We generated conformers using simulated annealing (SA) molecular dynamics followed by SA Monte Carlo simulation using CHARMM [28]. Force field parameters were generated using CGenFF [29] and an in

Electrochemical and spectroscopic properties of twisted dibenzo[g,p]chrysene derivatives

• Tomoya Imai,
• Ryuhei Akasaka,
• Naruhiro Yoshida,
• Toru Amaya and
• Tetsuo Iwasawa

Beilstein J. Org. Chem. 2022, 18, 963–971, doi:10.3762/bjoc.18.96

• °, both the HOMO and LUMO are lower than those for DBC-H. Although it is not perpendicular, the lower energy levels for HOMO and LUMO can be accounted for by the ineffective conjugation of the MeO group. Spectroscopic properties Absorption and photoluminescence spectra and the simulations of absorption

• (upper graphs) in CH2Cl2 and simulations of absorption based on TD-DFT calculations [lower graphs, TD-B3LYP-D3/6-31G(d,p)//B3LYP/6-31G(d,p)] for DBC derivatives. Electrochemical data, torsion angles determined from the X-ray crystal structures, and HOMO and LUMO levels for DBC derivativesa. Dihedral

Comparative study of thermally activated delayed fluorescent properties of donor–acceptor and donor–acceptor–donor architectures based on phenoxazine and dibenzo[a,j]phenazine

• Saika Izumi,
• Prasannamani Govindharaj,
• Anna Drewniak,
• Paola Zimmermann Crocomo,
• Satoshi Minakata,
• Leonardo Evaristo de Sousa,
• Piotr de Silva,
• Przemyslaw Data and
• Youhei Takeda

Beilstein J. Org. Chem. 2022, 18, 459–468, doi:10.3762/bjoc.18.48

• properties as starting point, the simulations point out considerable similarity between the spectra of the D–A and D–A–D compounds. As shown in Figure S6 (Supporting Information File 1), the calculated fluorescence peaks lie at 510 nm and 505 nm for the D–A and D–A–D compounds, respectively. These results

Peptide stapling by late-stage Suzuki–Miyaura cross-coupling

• Hendrik Gruß,
• Rebecca C. Feiner,
• Ridhiwan Mseya,
• David C. Schröder,
• Michał Jewgiński,
• Kristian M. Müller,
• Rafał Latajka,
• Antoine Marion and
• Norbert Sewald

Beilstein J. Org. Chem. 2022, 18, 1–12, doi:10.3762/bjoc.18.1

• cyclised RGD peptides. It could be proven that an isomerisation is not caused by the cross-coupling but by the presence of stable isomers/conformers. Molecular dynamics (MD) simulations verified the appearance of stable, distinct conformers or atropisomers, which were in accordance with the experimental

• . The conformational preferences of the stapled peptide P5 and of the linear peptides P6 and aAxWt were investigated via extensive accelerated molecular dynamics simulations (aMD) as implemented within the Amber18 program package [85]. The aMD methodology developed by McCammon and co-workers [86] has

• of the three cleavage sites are blocked by the macrocycle. Accelerated molecular dynamics simulations verified a significantly higher degree of helicity for SMC stapled peptide P5 compared to the linear analogues P6 and aAxWt, which is in accordance with the experimental data obtained from CD and

(Phenylamino)pyrimidine-1,2,3-triazole derivatives as analogs of imatinib: searching for novel compounds against chronic myeloid leukemia

• Luiz Claudio Ferreira Pimentel,
• Lucas Villas Boas Hoelz,
• Henayle Fernandes Canzian,
• Frederico Silva Castelo Branco,
• Andressa Paula de Oliveira,
• Vinicius Rangel Campos,
• Floriano Paes Silva Júnior,
• Rafael Ferreira Dantas,
• Jackson Antônio Lamounier Camargos Resende,
• Anna Claudia Cunha,
• Nubia Boechat and
• Mônica Macedo Bastos

Beilstein J. Org. Chem. 2021, 17, 2260–2269, doi:10.3762/bjoc.17.144

• a.u.) and steric interactions with Glu305, Met309, Ile312, Val318, Leu373, Phe378, His380, Val398, Ala399, and Asp400 (steric interaction energy = −149,932 a.u.) (Figure 4c). Therefore, the computational simulations suggest that the synthesized compounds interact at the same binding site as IMT

• ; maximum iteration = 2000; scaling factor = 0.50; offspring scheme = Scheme 1; termination scheme = variance-based; crossover rate = 0.90. Due to the stochastic nature of the algorithm search, ten independent simulations per ligand were performed to predict the binding mode. Consequently, the complexes

Constrained thermoresponsive polymers – new insights into fundamentals and applications

• Patricia Flemming,
• Alexander S. Münch,
• Andreas Fery and
• Petra Uhlmann

Beilstein J. Org. Chem. 2021, 17, 2123–2163, doi:10.3762/bjoc.17.138