CLEAPA: a framework for exploring the conformational landscape of cryo-EM using energy-aware pathfinding algorithm.

MOTIVATION: Cryo-electron microscopy (cryo-EM) is a powerful technique for studying macromolecules and holds the potential for identifying kinetically preferred transition sequences between conformational states. Typically, these sequences are explored within two-dimensional energy landscapes. However, due to the complexity of biomolecules, representing conformational changes in two dimensions can be challenging. Recent advancements in reconstruction models have successfully extracted structural heterogeneity from cryo-EM images using higher-dimension latent space. Nonetheless, creating high-dimensional conformational landscapes in the latent space and then searching for preferred paths continues to be a formidable task. RESULTS: This study introduces an innovative framework for identifying preferred trajectories within high-dimensional conformational landscapes. Our method encompasses the search for the minimum energy path in the graph, where edge weights are determined based on the energy estimation at each node using local density. The effectiveness of this approach is demonstrated by identifying accurate transition states in both synthetic and real-world datasets featuring continuous conformational changes. AVAILABILITY AND IMPLEMENTATION: The CLEAPA package is available at https://github.com/tengyulin/energy_aware_pathfinding/.

Three neo-Clerodane Diterpenoids from Tinospora cordifolia Stems and Their Arginase Inhibitory Activities.

Three neo-clerodane diterpenoids, including two new tinocordifoliols A (1) and B (2) and one known tinopanoid R (3), were isolated from the ethyl acetate-soluble fraction of the 70% ethanol extract of Tinospora cordifolia stems. The structures were elucidated by various spectroscopic methods, including one dimensional (1D) and 2D-NMR, high resolution-electrospray ionization (HR-ESI)-MS, and electronic circular dichroism (ECD) data. The T. cordifolia extract and all isolated compounds 1-3 possessed arginase I inhibitory activities. Among them, 3 exhibited moderate competitive inhibition of human arginase I (IC50 = 61.9 µM). Furthermore, docking studies revealed that the presence of a ß-substituted furan in 3 may play a key role in the arginase I inhibitory activities.

New tetrahydroisoquinoline-4-carbonitrile derivatives as potent agents against cyclin-dependent kinases, crystal structures, and computational studies.

The synthesis of two new hexahydroisoquinoline-4-carbonitrile derivatives (3a and 3b) is reported along with spectroscopic data and their crystal structures. In compound 3a, the intramolecular O-H···O hydrogen bond constraints the acetyl and hydroxyl groups to be syn. In the crystal, inversion dimers are generated by C-H···O hydrogen bonds and are connected into layers parallel to (10-1) by additional C-H···O hydrogen bonds. The layers are stacked with Cl···S contacts 0.17 Å less than the sum of the respective van der Waals radii. The conformation of the compound 3b is partially determined by the intramolecular O-H···O hydrogen bond. A puckering analysis of the tetrahydroisoquinoline unit was performed. In the crystal, O-H···O and C-H···O hydrogen bonds together with C-H···π(ring) interactions form layers parallel to (01-1) which pack with normal van der Waals interactions. To understand the binding efficiency and stability of the title molecules, molecular docking, and 100 ns dynamic simulation analyses were performed with CDK5A1. To rationalize their structure-activity relationship(s), a DFT study at the B3LYP/6-311++G** theoretical level was also done. The 3D Hirshfled surfaces were also taken to investigate the crystal packings of both compounds. In addition, their ADMET properties were explored.Communicated by Ramaswamy H. Sarma.

Elucidating the structural and conformational preferences of bioactive chromone and its monohydrate through advanced rotational spectroscopy.

Chromones are a class of naturally occurring compounds, renowned for their diverse biological activities with significant relevance in medicine and biochemistry. This study marks the first analysis of rotational spectra of both the chromone monomer and its monohydrate through Fourier transform microwave spectroscopy. The observation of nine mono-substituted 13C isotopologues facilitated a semi-experimental determination of the equilibrium structure of the chromone monomer. In the case of chromone monohydrate, two distinct isomers were identified, each characterized by a combination of O-H⋯O and C-H⋯O hydrogen bonds involving the chromone's carbonyl group. This study further delved into intermolecular non-covalent interactions, employing different theoretical approaches. The relative population ratio of the two identified isomers was estimated to be about 2:1 within the supersonic jet.

Lauinoids A-X: Labdane-type diterpenoids with anti-inflammatory activity from Croton laui.

Croton laui (Euphorbiaceae) is a traditional medicinal plant used by the Li ethnic group in China to treat headaches, stomachaches, and diphtheria. To understand the pharmacological basis of its medicinal use, an extensive investigation of the ethanolic extract of the bark of C. laui was performed. After repeated chromatography, twenty-four undescribed labdane-type diterpenoids, lauinoids A-X (1-24), and five known analogs (25-29) were isolated. Their structures and absolute configurations were established using a combination of spectroscopic analyses, electronic circular dichroism, nuclear magnetic resonance calculations, and single-crystal X-ray diffraction. Among them, compounds 1-3 exhibited an 11(12 â†’ 13)-abeo-16-nor-labdane skeleton, which originated putatively from 9 through a plausible pathway that involves a semipinacol rearrangement process. Compounds 11 and 12 belong to the rare class of 14,15-dinor-labdane diterpenoids. Compounds 18 and 28 exhibited substantial inhibitory effects by suppressing lipopolysaccharide-induced NO production in RAW 264.7 macrophages, with IC50 values of 3.37 ± 0.23 and 5.82 ± 0.28 µM, respectively. This study has greatly expanded the chemical diversity of labdane diterpenoids from C. laui and will guide future research on this ethnomedicinal plant.

Accurate and Efficient Conformer Sampling of Cyclic Drug-Like Molecules with Inverse Kinematics.

Identification of all of the influential conformers of biomolecules is a crucial step in many tasks of computational biochemistry. Specifically, molecular docking, a key component of in silico drug development, requires a comprehensive set of conformations for potential candidates in order to generate the optimal ligand-receptor poses and, ultimately, find the best drug candidates. However, the presence of flexible cycles in a molecule complicates the initial search for conformers since exhaustive sampling algorithms via torsional random and systematic searches become very inefficient. The devised inverse-kinematics-based Monte Carlo with refinement (MCR) algorithm identifies independently rotatable dihedral angles in (poly)cyclic molecules and uses them to perform global conformational sampling, outperforming popular alternatives (MacroModel, CREST, and RDKit) in terms of speed and diversity of the resulting conformer ensembles. Moreover, MCR quickly and accurately recovers naturally occurring macrocycle conformations for most of the considered molecules.

Exploring the Conformational Landscape of Poly(l-lysine) Dendrimers Using Ion Mobility Mass Spectrometry.

Ion mobility mass spectrometry (IM-MS) measures the mass, size, and shape of ions in the same experiment, and structural information is provided via collision cross-section (CCS) values. The majority of commercially available IM-MS instrumentation relies on the use of CCS calibrants, and here, we present data from a family of poly(l-lysine) dendrimers and explore their suitability for this purpose. In order to test these compounds, we employed three different IM-MS platforms (Agilent 6560 IM-QToF, Waters Synapt G2, and a home-built variable temperature drift tube IM-MS) and used them to investigate six different generations of dendrimers in two buffer gases (helium and nitrogen). Each molecule gives a highly discrete CCS distribution suggestive of single conformers for each m/z value. The DTCCSN2 values of this series of molecules (molecular weight: 330-16,214 Da) range from 182 to 2941 Å2, which spans the CCS range that would be found by many synthetic molecules including supramolecular compounds and many biopolymers. The CCS values for each charge state were highly reproducible in day-to-day analysis on each instrument, although we found small variations in the absolute CCS values between instruments. The rigidity of each dendrimer was probed using collisionally activated and high-temperature IM-MS experiments, where no evidence for a significant CCS change ensued. Taken together, this data indicates that these polymers are candidates for CCS calibration and could also help to reconcile differences found in CCS measurements on different instrument geometries.

Cevanine-type steroidal alkaloids from the bulbs of Fritillaria cirrhosa D. Don.

Eight previously undescribed cevanine-type steroidal alkaloids, cirrhosinones I-N and cirrhosinols A-B, along with five known analogs, were isolated from the bulbs of Fritillaria cirrhosa D. Don. Their structures were elucidated on the basis of comprehensive analysis of HRESIMS, 1D and 2D NMR spectroscopic data, and single-crystal X-ray diffraction analyses. All compounds revealed weak NO inhibitory activities in the LPS-stimulated NR8383 cells at the concentration of 20 µM, with inhibition ratios ranging from 5.1% to 14.3%.

[Development of Membrane Permeability Coefficient by Means of Novel Molecular Dynamics Methods].

The membrane permeability, and its evaluation, is crucial factor in the process of uptake of compounds from outside to inside the cell and in the inhibition of the activity of disease-causing target proteins. Although molecular dynamics (MD) simulations have been shown to be able to reproduce the conformational changes of compounds occurring during membrane permeation, it is still challenging to extract the membrane permeability at an affordable computational workload solely by conventional MD. Indeed, the time scale accessible by MD is far below the one characterizing the actual permeation process. Phenomena occurring in living organisms escaping the reach of standard MD are generally referred to as biological rare events, and the membrane permeation process is one of them. To overcome this time-scale problem, several enhanced sampling methods have been proposed over the years to improve conformational sampling. In this review, a hybrid sampling method that combines the parallel cascade selection MD (PaCS-MD) and the outlier flooding method (OFLOOD), introduced and developed by our group, is proposed as a tool to study the membrane permeation from structural sampling (rare-event sampling). The obtained trajectories are used to estimate the free energy profiles for the membrane permeation and to compute the membrane permeation coefficients. Moreover, we present an example of application of the free energy reaction network method as a versatile way for incorporating explicitly into reaction coordinates the degrees of freedom related to internal motion.

Conformational analysis and spectroscopic properties of antichagasic nifurtimox.

Considering the health relevance of Chagas' disease, recent research efforts have focused on developing more efficient drug delivery systems containing nifurtimox (NFX). This paper comprehensively investigates NFX through conformational analysis and spectroscopic characterization. Using a conformer-rotamer ensemble sampling tool (CREST-xtb), five distinct conformers of NFX were sampled within a 3.0 kcal mol-1 relative energy window. Subsequently, such structures were used as inputs for geometry optimization by density functional theory (DFT) at B3LYP-def2-TZVP level of theory. Notably, harmonic vibrational frequencies were calculated to establish an in-depth comparison with experimental results and existing literature for the NFX or similar molecules and functional groups, thereby achieving a widely reasoned assignment of the mid-infrared band absorptions for the first time. Moreover, UV-VIS spectra of NFX were obtained in several solvents, enabling the determination of the molar absorptivity coefficient for the two electronic transitions observed for NFX. Among the aprotic solvents, a bathochromic effect was observed in the function of the dielectric constants. Furthermore, a hypochromic effect was observed when the drug was dissolved in protic solvents. These findings offer crucial support for new drug delivery systems containing NFX while demonstrating the potential of spectrophotometric studies in establishing quality control assays for NFX drug products.

Relationships between conformational and vibrational features of tryptophan characteristic Raman markers.

Tryptophan (Trp) residue provides characteristic vibrational markers to the middle wavenumber spectral region of the Raman spectra recorded from peptides and proteins. In this report, we were particularly interested in eight Trp Raman markers, referred to as Wi (i = 1,…,8). All responsible for pronounced Raman lines, these markers originate from indole moiety, a bicyclic conjugated segment involved in the Trp structure. Numerous investigations have previously attempted to relate the variations observed in the spectral features of these markers to the environmental changes of Trp residues. To emphasize the most important points we can mention (i) the variations in the Raman profile of W4 (∼1360 cm-1) and W5 (∼1340 cm-1), frequently observed as a doublet with variable intensity ratio. These two markers were thought to result from a Fermi-resonance effect between certain planar and nonplanar modes; (ii) the changes observed in the wavenumbers and relative intensities of W4, W7 (∼880 cm-1) and W8 (∼760 cm-1) were supposed to be related to the accessibility of Trp to surrounding water molecules; and (iii) the wavenumber fluctuations of W3 (∼1550 cm-1), taken as a Trp side chain orientational marker. However, some ambiguities still exist regarding the interpretation of these markers, needing further clarification. Herein, upon a joint experimental and theoretical analysis based on a multiconformational approach, attention was paid to the relationships between structural and vibrational features of three indole-containing compounds with increasing structural complexity, i.e., skatole (3-methylindole), tryptophan, and tripeptide Gly-Trp-Gly. This study clearly shows that the existing assignments given to certain Trp Raman markers should be reconsidered, especially those based on the Fermi-resonance origin of W4-W5 (∼1360-1340 cm-1) doublet, as well as the purely environmental dependence of W7 and W8 markers.

Isolation and Structure Elucidation of JBIR-157, a Skeletally Novel Aromatic Polyketide Produced by the Heterologous Expression of a Cryptic Gene Cluster.

Heterologous expression of natural compound biosynthetic gene clusters (BGCs) is a robust approach for not only revealing the biosynthetic mechanisms leading to the compounds, but also for discovering new products from uncharacterized BGCs. We established a heterologous expression technique applicable to huge biosynthetic gene clusters for generating large molecular secondary metabolites such as type-I polyketides. As an example, we targeted concanamycin BGC from Streptomyces neyagawaensis IFO13477 (the cluster size of 99 kbp), and obtained a bacterial artificial chromosome (BAC) clone with an insert size of 211 kbp that contains the entire concanamycin BGC. Interestingly, heterologous expression for this BAC clone resulted in two additional aromatic polyketides, ent-gephyromycin, and a new compound designated as JBIR-157, together with the expected concanamycin. Bioinformatic and biochemical analyses revealed that a cryptic biosynthetic gene cluster in this BAC clone was responsible for the production of these type-II polyketide synthases (PKS) compounds. Here, we describe the production, isolation, and structure elucidation of JBIR-157, determined primarily by a series of NMR spectral analyses.

COSMIC: Molecular Conformation Space Modeling in Internal Coordinates with an Adversarial Framework.

The fast and accurate conformation space modeling is an essential part of computational approaches for solving ligand and structure-based drug discovery problems. Recent state-of-the-art diffusion models for molecular conformation generation show promising distribution coverage and physical plausibility metrics but suffer from a slow sampling procedure. We propose a novel adversarial generative framework, COSMIC, that shows comparable generative performance but provides a time-efficient sampling and training procedure. Given a molecular graph and random noise, the generator produces a conformation in two stages. First, it constructs a conformation in a rotation and translation invariant representation─internal coordinates. In the second step, the model predicts the distances between neighboring atoms and performs a few fast optimization steps to refine the initial conformation. The proposed model considers conformation energy, achieving comparable space coverage, and diversity metrics results.

Exploring long-range fluorine-carbon J-coupling for conformational analysis of deoxyfluorinated disaccharides: A combined computational and NMR study.

In this study, we investigated the potential of long-range fluorine-carbon J-coupling for determining the structures of deoxyfluorinated disaccharides. Three disaccharides, previously synthesized as potential galectin inhibitors, exhibited through-space fluorine-carbon J-couplings. In our independent conformational analysis of these disaccharide derivatives, we employed a combination of density functional theory (DFT) calculations and nuclear magnetic resonance (NMR) experiments. By comparing the calculated nuclear shieldings with the experimental carbon chemical shifts, we were able to identify the most probable conformers for each compound. A model comprising fluoromethane and methane molecules was used to study the relationship between molecular arrangements and intermolecular through-space J-coupling. Our study demonstrates the important effect of internuclear distance and molecular orientation on the magnitude of fluorine-carbon coupling. The experimental values for the fluorine-carbon through-space couplings (TSCs) of the disaccharides corresponded with values calculated for the most probable conformers identified by the conformational analysis. These results unlock the broader application of fluorine-carbon TSCs as powerful tools for conformational analysis of flexible molecules, offering valuable insights for future structural investigations.

Fagopyrins from Buckwheat Flowers: Structural and Stereochemical Characterization Through Combined NMR/CD Spectroscopy and Theoretical Calculations.

Fagopyrins are phenantroperylenequinones present in the flowers of Fagopyrum esculentum (buckwheat) endowed with photodynamic activity. It has been reported that fagopyrin extracts actually contain a complex mixture of closely related compounds, differing only on the nature of the perylenequinone substituents. We report our systematic and detailed study on the chemical composition of fagopyrin extracts by a combination of preparative and analytical techniques. The combined use of 1H-NMR and CD spectroscopy was found to be particularly suited to fully characterize all stereochemical aspects of the extracted fagopyrins. For the first time nine isomers have been structurally characterized and their stereochemistry fully elucidated. The presence of two different heterocyclic ring substituents, two stereogenic centers and the inherent axial chirality of the aromatic system provides a complex stereochemical relationships among isomers, thus giving account of the high level of molecular multiplicity found in the extract.

Daphmacrimines A-K, Daphniphyllum alkaloids from Daphniphyllum macropodum Miq.

Daphmacrimines A-K (1-11) were isolated from the leaves and stems of Daphniphyllum macropodum Miq. Their structures and stereochemistries were determined by extensive techniques, including HRESIMS, NMR, ECD, IR, and single-crystal X-ray crystallography. Daphmacrimines A-D (1-4) are unprecedented Daphniphyllum alkaloids with a 2-oxazolidinone ring. Daphmacrimine I (9) contains a nitrile group, which is relatively rare in naturally occurring alkaloids. The abilities of daphmacrimines A-D and daphmacrimines G-K to enhance lysosomal biogenesis were evaluated through LysoTracker Red staining. Daphmacrimine K (11) can induce lysosomal biogenesis and promote autophagic flux.

4,4-Difluoroproline as a Unique 19F NMR Probe of Proline Conformation.

Despite the importance of proline conformational equilibria (trans versus cis amide and exo versus endo ring pucker) on protein structure and function, there is a lack of convenient ways to probe proline conformation. 4,4-Difluoroproline (Dfp) was identified to be a sensitive 19F NMR-based probe of proline conformational biases and cis-trans isomerism. Within model compounds and disordered peptides, the diastereotopic fluorines of Dfp exhibit similar chemical shifts (ΔδFF = 0-3 ppm) when a trans X-Dfp amide bond is present. In contrast, the diastereotopic fluorines exhibit a large (ΔδFF = 5-12 ppm) difference in chemical shift in a cis X-Dfp prolyl amide bond. DFT calculations, X-ray crystallography, and solid-state NMR spectroscopy indicated that ΔδFF directly reports on the relative preference of one proline ring pucker over the other: a fluorine which is pseudo-axial (i.e., the pro-4R-F in an exo ring pucker, or the pro-4S-F in an endo ring pucker) is downfield, while a fluorine which is pseudo-equatorial (i.e., pro-4S-F when exo, or pro-4R-F when endo) is upfield. Thus, when a proline is disordered (a mixture of exo and endo ring puckers, as at trans-Pro in peptides in water), it exhibits a small Δδ. In contrast, when the Pro is ordered (i.e., when one ring pucker is strongly preferred, as in cis-Pro amide bonds, where the endo ring pucker is strongly favored), a large Δδ is observed. Dfp can be used to identify inherent induced order in peptides and to quantify proline cis-trans isomerism. Using Dfp, we discovered that the stable polyproline II helix (PPII) formed in the denatured state (8 M urea) exhibits essentially equal populations of the exo and endo proline ring puckers. In addition, the data with Dfp suggested the specific stabilization of PPII by water over other polar solvents. These data strongly support the importance of carbonyl solvation and n → π* interactions for the stabilization of PPII. Dfp was also employed to quantify proline cis-trans isomerism as a function of phosphorylation and the R406W mutation in peptides derived from the intrinsically disordered protein tau. Dfp is minimally sterically disruptive and can be incorporated in expressed proteins, suggesting its broad application in understanding proline cis-trans isomerization, protein folding, and local order in intrinsically disordered proteins.

Conformational and Electronic Variations in 1,2- and 1,5a-Cyclophellitols and their Impact on Retaining α-Glucosidase Inhibition.

Glycoside hydrolases (glycosidases) take part in myriad biological processes and are important therapeutic targets. Competitive and mechanism-based inhibitors are useful tools to dissect their biological role and comprise a good starting point for drug discovery. The natural product, cyclophellitol, a mechanism-based, covalent and irreversible retaining ß-glucosidase inhibitor has inspired the design of diverse α- and ß-glycosidase inhibitor and activity-based probe scaffolds. Here, we sought to deepen our understanding of the structural and functional requirements of cyclophellitol-type compounds for effective human α-glucosidase inhibition. We synthesized a comprehensive set of α-configured 1,2- and 1,5a-cyclophellitol analogues bearing a variety of electrophilic traps. The inhibitory potency of these compounds was assessed towards both lysosomal and ER retaining α-glucosidases. These studies revealed the 1,5a-cyclophellitols to be the most potent retaining α-glucosidase inhibitors, with the nature of the electrophile determining inhibitory mode of action (covalent or non-covalent). DFT calculations support the ability of the 1,5a-cyclophellitols, but not the 1,2-congeners, to adopt conformations that mimic either the Michaelis complex or transition state of α-glucosidases.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Using Three-Dimensional Information to Predict and Interpret the Facial Selectivities of Nucleophilic Additions to Cyclic Ketones.

In this study, we devised a new method to predict facial selectivity by quantifying steric and orbital factors for the nucleophile approaching both π-plane faces. Using this method, we quantified the total electron density and frontier orbital distributions of 163 cyclic ketones with various structures and quantitatively explained the surface selectivity of 323 reactions with eight nucleophiles (BH3, LiAlH4, NaBH4, LiAl(OMe)3H, MeLi, MeMgI, PhLi, and PnMgI). Importance analysis showed a large orbital effect for BH3, LiAlH4, and NaBH4 and the dominance of the steric effect for LiAl(OMe)3H, MeLi, MeMgI, PhLi, and PhMgI. Our method analyzes three-dimensional features based on Gaussian cube files, which can be easily obtained using mainstream computational chemistry software packages, and this approach should prove useful for predicting the rates and facial selectivity of other reactions.