[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.

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.

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.

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.

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.

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.

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study.

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Conformational Effects on the Absorption Spectra of Phytochromes.

Bacteriophytochrome is a photoreceptor protein that contains the biliverdin (BV) chromophore as its active component. The spectra of BV upon mutation remain remarkably unchanged, as far as spectral positions are concerned. This points toward the minimal effect of electrostatic effects on the electronic structure of the chromophore. However, the relative intensities of the Q and Soret bands of the chromophore change dramatically upon mutation. In this work, we delve into the molecular origin of this unusual intensity modulation. Using extensive classical MD and QM/MM calculations, we show that due to mutation, the conformational population of the chromophore changes significantly. The noncovalent interactions, especially the stacking interactions, lead to extra stabilization of the cyclic form in the D207H mutated species as opposed to the open form in the wild-type BV. Thus, unlike the commonly observed direct electrostatic effect on the spectral shift, in the case of BV the difference observed is in varying intensities, and this in turn is driven by a conformational shift due to enhanced stacking interaction.

Tertiary structure and conformational dynamics of the anti-amyloidogenic chaperone DNAJB6b at atomistic resolution.

DNAJB6b is a molecular chaperone of the heat shock protein network, shown to play a crucial role in preventing aggregation of several disease-related intrinsically disordered proteins. Using homology modeling and microsecond-long all-atom molecular dynamics (MD) simulations, we show that monomeric DNAJB6b is a transiently interconverting protein cycling between three states: a closed state, an open state (both abundant), and a less abundant extended state. Interestingly, the reported regulatory autoinhibitory anchor between helix V in the G/F1 region and helices II/III of the J-domain, which obstructs the access of Hsp70 to the J-domain remains present in all three states. This possibly suggests a mechanistically intriguing regulation in which DNAJB6b only becomes exposed when loaded with substrates that require Hsp70 processing. Our MD results of DNAJB6b carrying mutations in the G/F1 region that are linked to limb-girdle muscular dystrophy type D1 (LGMDD1) show that this G/F1 region becomes highly dynamic, pointing towards a spontaneous release of the autoinhibitory helix V from helices II/III. This would increase the probability of non-functional Hsp70 interactions to DNAJB6b without substrates. Our cellular data indeed confirm that non-substrate loaded LGMDD1 mutants have aberrant interactions with Hsp70.

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.

Maturation and Conformational Switching of a De Novo Designed Phase-Separating Polypeptide.

Cellular compartments formed by biomolecular condensation are widespread features of cell biology. These organelle-like assemblies compartmentalize macromolecules dynamically within the crowded intracellular environment. However, the intermolecular interactions that produce condensed droplets may also create arrested states and potentially pathological assemblies such as fibers, aggregates, and gels through droplet maturation. Protein liquid-liquid phase separation is a metastable process, so maturation may be an intrinsic property of phase-separating proteins, where nucleation of different phases or states arises in supersaturated condensates. Here, we describe the formation of both phase-separated droplets and proteinaceous fibers driven by a de novo designed polypeptide. We characterize the formation of supramolecular fibers in vitro and in bacterial cells. We show that client proteins can be targeted to the fibers in cells using a droplet-forming construct. Finally, we explore the interplay between phase separation and fiber formation of the de novo polypeptide, showing that the droplets mature with a post-translational switch to largely ß conformations, analogous to models of pathological phase separation.

Modeling Zinc Complexes Using Neural Networks.

Understanding the energetic landscapes of large molecules is necessary for the study of chemical and biological systems. Recently, deep learning has greatly accelerated the development of models based on quantum chemistry, making it possible to build potential energy surfaces and explore chemical space. However, most of this work has focused on organic molecules due to the simplicity of their electronic structures as well as the availability of data sets. In this work, we build a deep learning architecture to model the energetics of zinc organometallic complexes. To achieve this, we have compiled a configurationally and conformationally diverse data set of zinc complexes using metadynamics to overcome the limitations of traditional sampling methods. In terms of the neural network potentials, our results indicate that for zinc complexes, partial charges play an important role in modeling the long-range interactions with a neural network. Our developed model outperforms semiempirical methods in predicting the relative energy of zinc conformers, yielding a mean absolute error (MAE) of 1.32 kcal/mol with reference to the double-hybrid PWPB95 method.

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.

Expansive discovery of chemically diverse structured macrocyclic oligoamides.

Small macrocycles with four or fewer amino acids are among the most potent natural products known, but there is currently no way to systematically generate such compounds. We describe a computational method for identifying ordered macrocycles composed of alpha, beta, gamma, and 17 other amino acid backbone chemistries, which we used to predict 14.9 million closed cycles composed of >42,000 monomer combinations. We chemically synthesized 18 macrocycles predicted to adopt single low-energy states and determined their x-ray or nuclear magnetic resonance structures; 15 of these were very close to the design models. We illustrate the therapeutic potential of these macrocycle designs by developing selective inhibitors of three protein targets of current interest. By opening up a vast space of readily synthesizable drug-like macrocycles, our results should considerably enhance structure-based drug design.

Computational insights into the conformational transition of STING: Mechanistic, energetic considerations, and the influence of crucial mutations.

STING (stimulator of interferon genes) is a crucial protein in the innate immune system's response to viral and bacterial infections. In this study, we investigated the mechanistic and energetic mechanism of the conformational transition process of STING activated by cGAMP binding. We found that the STING connector region undergoes an energetically unfavorable rotation during this process, which is compensated by the favorable interaction between cGAMP and the STING ligand binding domain. We further studied several disease-causing mutations and found that the V155 M mutation facilitates a smoother transition in the STING connector region. However, the V147L mutation exhibits unfavorable conformational transition energy, suggesting it may hinder STING activation pathway that relies on connector region rotation. Despite being labeled as hyperactive, the widespread prevalence of V147L/V147I mutations across species implies a neutral character, indicating complexity in its role. Overall, our analysis deepens the understanding of STING activation within the connector region, and targeting this region with compounds may provide an alternative approach to interfering with STING's function.

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.

Structurally diverse 1,2-diarylpropanes from the fruit of Crataegus pinnatifida and the investigation on their mirror-image ECD spectra with the same absolute configurations.

1,2-diarylpropanes are a kind of abundant natural products formed by radical coupling. On account of molecular flexibility, it was challenged in the identifications of relative and absolute configurations of the 1,2-diarylpropanes. In this research, fourteen pairs of enantiomeric 1,2-diarylpropanes (1a/1b-14a/14b), comprising twelve previously undescribed pairs (1a/1b-4a/4b, 6a/6b-10a/10b, and 12a/12b-14a/14b), were isolated from the fruit of Crataegus pinnatifida. Their structures were determined through multiple NMR spectral analyses, empirical NMR rules, X-ray crystallography, and the comparison of experimental ECD spectra with calculated data. In addition, the analysis of ECD spectra revealed that substituent effects could generate an inverted chiroptical response, exhibiting in mirror-image ECD signals. This phenomenon was investigated by conformational analysis, molecular orbital analysis, the transition density matrix and hole/electron distributions. Moreover, a potential experimental rule was proposed for the rapid determination of the absolute configurations of the 1,2-diarylpropanes.

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.

Unraveling the Catalytic Mechanism of ß-Cyclodextrin in the Vitamin D Formation.

Previous experimental studies have shown that the isomerization reaction of previtamin D3 (PreD3) to vitamin D3 (VitD3) is accelerated 40-fold when it takes place within a ß-cyclodextrin dimer, in comparison to the reaction occurring in conventional isotropic solutions. In this study, we employ quantum mechanics-based molecular dynamics (MD) simulations and statistical multistructural variational transition state theory to unveil the origin of this acceleration. We find that the conformational landscape in the PreD3 isomerization is highly dependent on whether the system is encapsulated. In isotropic media, the triene moiety of the PreD3 exhibits a rich torsional flexibility. However, when encapsulated, such a flexibility is limited to a more confined conformational space. In both scenarios, our calculated rate constants are in close agreement with experimental results and allow us to identify the PreD3 flexibility restriction as the primary catalytic factor. These findings enhance our understanding of VitD3 isomerization and underscore the significance of MD and environmental factors in biochemical modeling.