Understanding Ring Puckering in Small Molecules and Cyclic Peptides

Accelerated Lignocellulosic Molecule Adsorption Structure Determination

Here, we present a study combining Bayesian optimization structural inference with the machine learning interatomic potential Neural Equivariant Interatomic Potential (NequIP) to accelerate and enable the study of the adsorption of the conformationally flexible lignocellulosic molecules β-d-xylose and 1,4-β-d-xylotetraose on a copper surface. The number of structure evaluations needed to map out the relevant potential energy surfaces are reduced by Bayesian optimization, while NequIP minimizes the time spent on each evaluation, ultimately resulting in cost-efficient and reliable sampling of large systems and configurational spaces. Although the applicability of Bayesian optimization for the conformational analysis of the more flexible xylotetraose molecule is restricted by the sample complexity bottleneck, the latter can be effectively bypassed with external conformer search tools, such as the Conformer-Rotamer Ensemble Sampling Tool, facilitating the subsequent lower-dimensional global minimum adsorption structure determination. Finally, we demonstrate the applicability of the described approach to find adsorption structures practically equivalent to the density functional theory counterparts at a fraction of the computational cost.

Systematic Comparison of Experimental Crystallographic Geometries and Gas-Phase Computed Conformers for Torsion Preferences

We performed exhaustive torsion sampling on more than 3 million compounds using the GFN2-xTB method and performed a comparison of experimental crystallographic and gas-phase conformers. Many conformer sampling methods derive torsional angle distributions from experimental crystallographic data, limiting the torsion preferences to molecules that must be stable, synthetically accessible, and able to be crystallized. In this work, we evaluate the differences in torsional preferences of experimental crystallographic geometries and gas-phase computed conformers from a broad selection of compounds to determine whether torsional angle distributions obtained from semiempirical methods are suitable priors for conformer sampling. We find that differences in torsion preferences can be mostly attributed to a lack of available experimental crystallographic data with small deviations derived from gas-phase geometry differences. GFN2 demonstrates the ability to provide accurate and reliable torsional preferences that can provide a basis for new methods free from the limitations of experimental data collection. We provide Gaussian-based fits and sampling distributions suitable for torsion sampling and propose an alternative to the widely used "experimental torsion and knowledge distance geometry" (ETKDG) method using quantum torsion-derived distance geometry (QTDG) methods.

Fluorination Influences the Bioisostery of Myo-Inositol Pyrophosphate Analogs

Inositol pyrophosphates (PP-IPs) are densely phosphorylated messenger molecules involved in numerous biological processes. PP-IPs contain one or two pyrophosphate group(s) attached to a phosphorylated myo-inositol ring. 5PP-IP5 is the most abundant PP-IP in human cells. To investigate the function and regulation by PP-IPs in biological contexts, metabolically stable analogs have been developed. Here, we report the synthesis of a new fluorinated phosphoramidite reagent and its application for the synthesis of a difluoromethylene bisphosphonate analog of 5PP-IP5 . Subsequently, the properties of all currently reported analogs were benchmarked using a number of biophysical and biochemical methods, including co-crystallization, ITC, kinase activity assays and chromatography. Together, the results showcase how small structural alterations of the analogs can have notable effects on their properties in a biochemical setting and will guide in the choice of the most suitable analog(s) for future investigations.

DDD-028: A potent, neuroprotective, non-opioid compound for the treatment of diabetic neuropathy

Diabetic neuropathy (DN) is a painful, chronic ailment that affects a large segment of diabetic population worldwide. Current medications such as pregabalin or duloxetine treat only the pain symptom associated with DN, but not the underlying nerve damage. DDD-028 (1) is a small molecule that displays potent pain-relieving activity in streptozotocin (STZ)-induced rodent model of DN. Combined with other studies indicating that DDD-028 suppresses astrogliosis and nerve damage induced by the anti-cancer drug, paclitaxel, the present study suggests that DDD-028 would be useful as a disease modifying therapeutic in the treatment of DN. The 3-dimensional structure of DDD-028 was confirmed by single crystal X-ray crystallography.

Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R²) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.

Expanding the Chemical Space of 1,2-Difunctionalized Cyclobutanes

An efficient approach to the synthesis of previously unavailable or hardly accessible 1,2-difunctionalized cyclobutanes (mostly with NH₂/NHBoc, OH, SH, or SO₂F groups attached to the carbocycle either directly or via a CH₂ unit) relying on the divergent strategy is described. This class of compounds provides sp³-enriched and conformationally restricted building blocks that are of special demand for medicinal chemistry. The target compounds were prepared not only as pure racemic (±)-cis- and (±)-trans-diastereomers but in some cases also as single enantiomers. The developed procedures are readily scaled up and allow obtaining the target compounds on an up to hundred-gram scale. On the basis of the results of 20 X-ray diffraction experiments, structural characterization of the 1,2-difunctionalized cyclobutane core was performed using the extended Cremer-Pople puckering parameters and exit vector (EVP) plots.

A Dynamic Proton Bond: MH+·H2O ⇌ M·H3O+ Interconversion in Loosely Coordinated Environments

The interaction of organic molecules with oxonium cations within their solvation shell may lead to the emergence of dynamic supramolecular structures with recurrently changing host-guest chemical identity. We illustrate this phenomenon in benchmark proton-bonded complexes of water with polyether macrocyles. Despite the smaller proton affinity of water versus the ether group, water in fact retains the proton in the form of H₃O⁺, with increasing stability as the coordination number increases. Hindrance in many-fold coordination induces dynamic reversible (ether)·H₃O⁺ ⇌ (etherH⁺)·H₂O interconversion. We perform infrared action ion spectroscopy over a broad spectral range to expose the vibrational signatures of the loose proton bonding in these systems. Remarkably, characteristic bands for the two limiting proton bonding configurations are observed in the experimental vibrational spectra, superimposed onto diffuse bands associated with proton delocalization. These features cannot be described by static equilibrium structures but are accurately modeled within the framework of ab initio molecular dynamics.

The crystal structure of trans-dichlorido-(ethylenediamine-κ 2 N,N′)-bis(triphenylphosphine-κ 1 P)ruthenium(II), C38H38Cl2N2P2Ru

C38H38Cl2N2P2Ru, monoclinic, P21/c (no. 14), a = 10.8827(3) Å, b = 30.6226(7) Å, c = 20.2089(4) Å, β = 93.8560(10)°, V = 6719.5(3) Å3, Z = 8, Rgt (F) = 0.0306, wRref (F 2) = 0.0636, T = 173(2) K.

Influence of Sugar Modifications on the Nucleoside Conformation and Oligonucleotide Stability: A Critical Review

Ribofuranose sugar conformation plays an important role in the structure and dynamics of functional nucleic acids such as siRNAs, AONs, aptamers, miRNAs, etc. To improve their therapeutic potential, several chemical modifications have been introduced into the sugar moiety over the years. The stability of the oligonucleotide duplexes as well as the formation of stable and functional protein-oligonucleotide complexes are dictated by the conformation and dynamics of the sugar moiety. In this review, we systematically categorise various ribofuranose sugar modifications employed in DNAs and RNAs so far. We discuss different stereoelectronic effects imparted by different substituents on the sugar ring and how these effects control sugar puckering. Using this data, it would be possible to predict the precise use of chemical modifications and design novel sugar-modified nucleosides for therapeutic oligonucleotides that can improve their physicochemical properties.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Crystal structure of propane-1,3-diylbis(diphenylphosphine sulfide) ethanol solvate, C27H26P2S2

C27H26P2S2, tetragonal, I4 (no. 79), a = 16.9196(12) Å, c = 9.4187(8) Å, V = 2696.3(4) Å3, Z = 4, R gt(F) = 0.0341, wR ref(F 2) = 0.0726, T = 173(2) K.

Zero-Dimensional (Piperidinium)2MnBr4: Ring Puckering-Induced Isostructural Transition and Strong Electron-Phonon Coupling-Mediated Self-Trapped Exciton Emission

We report on the synthesis, structure, and photophysical properties of a lead-free organic-inorganic hybrid halide, (Piperidinium)₂MnBr₄ (PipMBr). It crystallizes in a monoclinic P2₁/n structure, with isolated MnBr₄ tetrahedra representing a zero-dimensional compound. It undergoes a reversible isostructural transition at 422/417 K in the heating/cooling cycle owing to the hydrogen-bonding rearrangement mediated by ring puckering of piperidinium cations. This compound exhibits green emission with a photoluminescence quantum yield of 51%. Interestingly, strong electron-longitudinal optical phonon coupling with γ_LO of 237 meV is evidenced from the broadening of the temperature-dependent emission linewidth and the Raman spectrum. Such strong electron-phonon coupling and a relatively low Debye temperature (137 K) suggest the self-trapped exciton emission in this compound.

Towards understanding solvation effects on the conformational entropy of non-rigid molecules

The absolute molecular entropy is a fundamental quantity for the accurate description of thermodynamic properties. For non-rigid molecules, a substantial part of the entropy can be attributed to a conformational contribution. Systems and properties where this is relevant, e.g., protein-ligand binding affinities or pK_a values refer usually to the liquid phase. In this work, the influence of solvation on the conformational entropy is investigated. A recently introduced state-of-the-art and automated computational protocol for the computation of conformational entropies [Pracht et al., Chem. Sci., 2021, 12, 6551-6568.] is applied in combination with fast and accurate semiempirical quantum-chemical methods and implicit solvation models for a set of 25 commercially available drug molecules and five transition metal compounds. Computed gas-phase conformational entropies are compared with values obtained in implicit n-hexane and water. It is found that implicit solvation can have a substantial effect of several cal mol^-1 K^-1 on the entropy as a result of large conformational changes in the different phases. We conclude that for flexible molecules chemical accuracy for free energies in solution can only be achieved if solvation effects on the conformational ensemble are considered.

Self-Consistent Parameterization of DNA Residues for the Non-Polarizable AMBER Force Fields

Fixed-charge (non-polarizable) forcefields are accurate and computationally efficient tools for modeling the molecular dynamics of nucleic acid polymers, particularly DNA, well into the µs timescale. The continued utility of these forcefields depends in part on expanding the residue set in step with advancing nucleic acid chemistry and biology. A key step in parameterizing new residues is charge derivation which is self-consistent with the existing residues. As atomic charges are derived by fitting against molecular electrostatic potentials, appropriate structural models are critical. Benchmarking against the existing charge set used in current AMBER nucleic acid forcefields, we report that quantum mechanical models of deoxynucleosides, even at a high level of theory, are not optimal structures for charge derivation. Instead, structures from molecular mechanics minimization yield charges with up to 6-fold lower RMS deviation from the published values, due to the choice of such an approach in the derivation of the original charge set. We present a contemporary protocol for rendering self-consistent charges as well as optimized charges for a panel of nine non-canonical residues that will permit comparison with literature as well as studying the dynamics of novel DNA polymers.

Design of Protein Segments and Peptides for Binding to Protein Targets

Recent years have witnessed a rise in methods for accurate prediction of structure and design of novel functional proteins. Design of functional protein fragments and peptides occupy a small, albeit unique, space within the general field of protein design. While the smaller size of these peptides allows for more exhaustive computational methods, flexibility in their structure and sparsity of data compared to proteins, as well as presence of noncanonical building blocks, add additional challenges to their design. This review summarizes the current advances in the design of protein fragments and peptides for binding to targets and discusses the challenges in the field, with an eye toward future directions.

Impact of Deuteration and Temperature on Furan Ring Dynamics

Despite significant progress in conformational analysis of cyclic molecules, the number of computational studies is still limited while most of that available in the literature data have been obtained long time ago with outdated methods. In present research, we have studied temperature driven conformational changes of the furan ring at three different temperatures. Additionally, the effect of deuteration on the ring dynamics is discussed; in addition, the aromaticity indices following the Bird and HOMA schemes are computed along all trajectories. Our ab initio molecular dynamic simulations revealed that deuteration has changed the furan ring dynamics and the obvious consequences; in addition, the shape and size of molecule are expected to be different.