Preferential Door for Ligand Binding and Unbinding Pathways in Inhibited Human Acetylcholinesterase.

Rising global population and increased food demands have resulted in the increased use of organophosphate pesticides (OPs), leading to toxin accumulation and transmission to humans. Pralidoxime (2-PAM), an FDA-approved drug, serves as an antidote for OP therapy. However, the atomic-level detoxification mechanisms regarding the design of novel antidotes remain unclear. This is the first study to examine the binding and unbinding pathways of 2-PAM to human acetylcholinesterase (HuAChE) through three identified doors using an enhanced sampling method called ligand-binding parallel cascade selection molecular dynamics (LB-PaCS-MD). Remarkably, LB-PaCS-MD could identify a predominant in-line binding mechanism through the acyl door at 63.79% ± 6.83%, also implicating it in a potential unbinding route (90.14% ± 4.22%). Interestingly, crucial conformational shifts in key residues, W86, Y341, and Y449, and the Ω loop significantly affect door dynamics and ligand binding modes. The LB-PaCS-MD technique can study ligand-binding pathways, thereby contributing to the design of antidotes and covalent drugs.

Site Identification and Next Choice Protocol for Hit-to-Lead Optimization.

Time efficiency and cost savings are major challenges in drug discovery and development. In this process, the hit-to-lead stage is expected to improve efficiency because it primarily exploits the trial-and-error approach of medicinal chemists. This study proposes a site identification and next choice (SINCHO) protocol to improve the hit-to-lead efficiency. This protocol selects an anchor atom and growth site pair, which is desirable for a hit-to-lead strategy starting from a 3D complex structure. We developed and fine-tuned the protocol using a training data set and assessed it using a test data set of the preceding hit-to-lead strategy. The protocol was tested for experimentally determined structures and molecular dynamics (MD) ensembles. The protocol had a high prediction accuracy for applying MD ensembles, owing to the consideration of protein flexibility. The SINCHO protocol enables medicinal chemists to visualize and modify functional groups in a hit-to-lead manner.

[Development of Integrated Database for Experiments and Simulations].

Researchers collect data and use various methods to organize it. Ensuring the reliability and reproducibility of data is crucial, and collaboration across different research fields is on the rise. However, when there is geographical distance, sharing data becomes a challenging task. Therefore, there is a need for the development of a mechanism for sharing data on the web. We have developed an integrated database to facilitate the sharing and management of research data, particularly focusing on small molecules. The integrated database serves as a platform for centralizing data related to small molecules, including their chemical structures, wet lab experimental data, simulation data, and more. It has been constructed as a web application, offering features such as library management for small molecules, registration and viewing of wet lab experiment results, generation of initial conformations for simulations, and data visualization. This enables researchers to efficiently share their research data and collaborate seamlessly, whether within their research group or via cloud-based access that allows project and team members to connect from anywhere. This integrated database plays a critical role in connecting wet lab experiments and simulations, enabling researchers to cross-reference and analyze experimental data comprehensively. It serves as an essential tool to advance research and foster idea generation.

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

Cocrystalline Matrices for Hyperpolarization at Room Temperature Using Photoexcited Electrons.

We propose using cocrystals as effective polarization matrices for triplet dynamic nuclear polarization (DNP) at room temperature. The polarization source can be uniformly doped into cocrystals formed through acid-acid, amide-amide, and acid-amide synthons. The dense-packing crystal structures, facilitated by multiple hydrogen bonding and π-π interactions, result in extended T1 relaxation times, enabling efficient polarization diffusion within the crystals. Our study demonstrates the successful polarization of a DNP-magnetic resonance imaging molecular probe, such as urea, within a cocrystal matrix at room temperature using triplet-DNP.

Development of an FeII Complex Exhibiting Intermolecular Proton Shifting Coupled Spin Transition.

In this study, we investigated reversible intermolecular proton shifting (IPS) coupled with spin transition (ST) in a novel FeII complex. The host FeII complex and the guest carboxylic acid anion were connected by intermolecular hydrogen bonds (IHBs). We extended the intramolecular proton transfer coupled ST phenomenon to the intermolecular system. The dynamic phenomenon was confirmed by variable-temperature single-crystal X-ray diffraction, neutron crystallography, and infrared spectroscopy. The mechanism of IPS was further validated using density functional theory calculations. The discovery of IPS-coupled ST in crystalline molecular materials provides good insights into fundamental processes and promotes the design of novel multifunctional materials with tunable properties for various applications, such as optoelectronics, information storage, and molecular devices.

Grinding-Induced Water Solubility Exhibited by Mechanochromic Luminescent Supramolecular Fibers.

Most mechanochromic luminescent compounds are crystalline and highly hydrophobic; however, mechanochromic luminescent molecular assemblies comprising amphiphilic molecules have rarely been explored. This study investigated mechanochromic luminescent supramolecular fibers composed of dumbbell-shaped 9,10-bis(phenylethynyl)anthracene-based amphiphiles without any tetraethylene glycol (TEG) substituents or with two TEG substituents. Both amphiphiles formed water-insoluble supramolecular fibers via linear hydrogen bond formation. Both compounds acquired water solubility when solid samples composed of supramolecular fibers are ground. Grinding induces the conversion of 1D supramolecular fibers into micellar assemblies where fluorophores can form excimers, thereby resulting in a large redshift in the fluorescence spectra. Excimer emission from the ground amphiphile without TEG chains is retained after dissolution in water. The micelles are stable in water because hydrophilic dendrons surround the hydrophobic luminophores. By contrast, when water is added to a ground amphiphile having TEG substituents, fragmented supramolecular fibers with the same molecular arrangement as the initial supramolecular fibers are observed, because fragmented fibers are thermodynamically preferable to micelles as the hydrophobic arrays of fluorophores are covered with hydrophilic TEG chains. This leads to the recovery of the initial fluorescent properties for the latter amphiphile. These supramolecular fibers can be used as practical mechanosensors to detect forces at the mesoscale.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors.

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Machine learning-guided design of potent darunavir analogs targeting HIV-1 proteases: A computational approach for antiretroviral drug discovery.

In the pursuit of novel antiretroviral therapies for human immunodeficiency virus type-1 (HIV-1) proteases (PRs), recent improvements in drug discovery have embraced machine learning (ML) techniques to guide the design process. This study employs ensemble learning models to identify crucial substructures as significant features for drug development. Using molecular docking techniques, a collection of 160 darunavir (DRV) analogs was designed based on these key substructures and subsequently screened using molecular docking techniques. Chemical structures with high fitness scores were selected, combined, and one-dimensional (1D) screening based on beyond Lipinski's rule of five (bRo5) and ADME (absorption, distribution, metabolism, and excretion) prediction implemented in the Combined Analog generator Tool (CAT) program. A total of 473 screened analogs were subjected to docking analysis through convolutional neural networks scoring function against both the wild-type (WT) and 12 major mutated PRs. DRV analogs with negative changes in binding free energy ( ΔΔ G bind ) compared to DRV could be categorized into four attractive groups based on their interactions with the majority of vital PRs. The analysis of interaction profiles revealed that potent designed analogs, targeting both WT and mutant PRs, exhibited interactions with common key amino acid residues. This observation further confirms that the ML model-guided approach effectively identified the substructures that play a crucial role in potent analogs. It is expected to function as a powerful computational tool, offering valuable guidance in the identification of chemical substructures for synthesis and subsequent experimental testing.

Manipulation of photosynthetic energy transfer by vibrational strong coupling.

Uncovering the mystery of efficient and directional energy transfer in photosynthetic organisms remains a critical challenge in quantum biology. Recent experimental evidence and quantum theory developments indicate the significance of quantum features of molecular vibrations in assisting photosynthetic energy transfer, which provides the possibility of manipulating the process by controlling molecular vibrations. Here, we propose and theoretically demonstrate efficient manipulation of photosynthetic energy transfer by using vibrational strong coupling between the vibrational state of a Fenna-Matthews-Olson (FMO) complex and the vacuum state of an optical cavity. Specifically, based on a full-quantum analytical model to describe the strong coupling effect between the optical cavity and molecular vibration, we realize efficient manipulation of energy transfer efficiency (from 58% to 92%) and energy transfer time (from 20 to 500 ps) in one branch of FMO complex by actively controlling the coupling strength and the quality factor of the optical cavity under both near-resonant and off-resonant conditions, respectively. Our work provides a practical scenario to manipulate photosynthetic energy transfer by externally interfering molecular vibrations via an optical cavity and a comprehensible conceptual framework for researching other similar systems.

Latrunculin resistance mechanism of non-conventional actin NAP1 uncovered by molecular dynamics simulations.

Monomeric G-actin polymerizes into F-actin to perform various cellular functions. Actin depolymerization drugs, such as latrunculin-A (Lat-A), inhibit filament formation and disrupt the cytoskeleton. Interestingly, the green algae Chlamydomonas alternatively produces a non-conventional actin, NAP1, that responds to inhibition by latrunculin. However, the molecular mechanism underlying latrunculin resistance of NAP1 remains unclear because of the difficulty due to its low in vitro polymerizability. Instead of biochemical experiments, we performed molecular dynamics (MD) simulations to investigate whether NAP1 has a lower affinity for Lat-A than the conventional actins. Our phylogenetic comparison of the binding free energies shows that Lat-A is evolutionarily optimized for skeletal muscles. By decomposing the binding free energy into each amino acid residue, we found that some residues in NAP1 play an important role in latrunculin resistance, suggesting that the primary mechanism of latrunculin resistance is the loss of affinity for Lat-A due to substitutions. In conclusion, our binding-free-energy calculations using MD simulations provide the critical insight that loss of affinity is the direct mechanism of latrunculin resistance.

Pathological mutations promote proteolysis of mitochondrial tRNA-specific 2-thiouridylase 1 (MTU1) via mitochondrial caseinolytic peptidase (CLPP).

MTU1 controls intramitochondrial protein synthesis by catalyzing the 2-thiouridine modification of mitochondrial transfer RNAs (mt-tRNAs). Missense mutations in the MTU1 gene are associated with life-threatening reversible infantile hepatic failure. However, the molecular pathogenesis is not well understood. Here, we investigated 17 mutations associated with this disease, and our results showed that most disease-related mutations are partial loss-of-function mutations, with three mutations being particularly severe. Mutant MTU1 is rapidly degraded by mitochondrial caseinolytic peptidase (CLPP) through a direct interaction with its chaperone protein CLPX. Notably, knockdown of CLPP significantly increased mutant MTU1 protein expression and mt-tRNA 2-thiolation, suggesting that accelerated proteolysis of mutant MTU1 plays a role in disease pathogenesis. In addition, molecular dynamics simulations demonstrated that disease-associated mutations may lead to abnormal intermolecular interactions, thereby impairing MTU1 enzyme activity. Finally, clinical data analysis underscores a significant correlation between patient prognosis and residual 2-thiolation levels, which is partially consistent with the AlphaMissense predictions. These findings provide a comprehensive understanding of MTU1-related diseases, offering prospects for modification-based diagnostics and novel therapeutic strategies centered on targeting CLPP.

Synthesis and biological evaluation of 2'-hydroxychalcone derivatives as AMPK activators.

A series of 2'-hydroxychalcone derivatives with various substituents on B-ring were synthesized and evaluated for AMP-activated protein kinase (AMPK) activation activity in podocyte cells. The results displayed that hydroxy, methoxy and methylenedioxy groups on B-ring could enhance the activitiy better than O-saturated alkyl, O-unsaturated alkyl or other alkoxy groups. Compounds 27 and 29 possess the highest fold change of 2.48 and 2.73, respectively, which were higher than those of reference compound (8) (1.28) and metformin (1.88). Compounds 27 and 29 were then subjected to a concentration-response study to obtain the EC50 values of 2.0 and 4.8 µM, respectively and MTT assays also showed that cell viability was not influenced by the exposure of podocytes to compounds 27 and 29 at concentrations up to 50 µM. In addition, compound 27 was proved to activate AMPK via calcium/calmodulin-dependent protein kinase kinase ß (CaMKKß)-dependent pathway without affecting intracellular calcium levels. The computational study showed that the potent compounds exhibited stronger ligand-binding strength to CaMKKß, particularly compounds 27 (-8.4 kcal/mol) and 29 (-8.0 kcal/mol), compared to compound 8 (-7.5 kcal/mol). Fragment molecular orbital (FMO) calculation demonstrated that compound 27 was superior to compound 29 due to the presence of methyl group, which amplified the binding by hydrophobic interactions. Therefore, compound 27 would represent a promising AMPK activator for further investigation of the treatment of diabetes and diabetic nephropathy.

Structural and thermodynamic insights into antibody light chain tetramer formation through 3D domain swapping.

Overexpression of antibody light chains in small plasma cell clones can lead to misfolding and aggregation. On the other hand, the formation of amyloid fibrils from antibody light chains is related to amyloidosis. Although aggregation of antibody light chain is an important issue, atomic-level structural examinations of antibody light chain aggregates are sparse. In this study, we present an antibody light chain that maintains an equilibrium between its monomeric and tetrameric states. According to data from X-ray crystallography, thermodynamic and kinetic measurements, as well as theoretical studies, this antibody light chain engages in 3D domain swapping within its variable region. Here, a pair of domain-swapped dimers creates a tetramer through hydrophobic interactions, facilitating the revelation of the domain-swapped structure. The negative cotton effect linked to the ß-sheet structure, observed around 215 nm in the circular dichroism (CD) spectrum of the tetrameric variable region, is more pronounced than that of the monomer. This suggests that the monomer contains less ß-sheet structures and exhibits greater flexibility than the tetramer in solution. These findings not only clarify the domain-swapped structure of the antibody light chain but also contribute to controlling antibody quality and advancing the development of future molecular recognition agents and drugs.

A sulfonamide chalcone inhibited dengue virus with a potential target at the SAM-binding site of viral methyltransferase.

Dengue infection is a global health problem as climate change facilitates the spread of mosquito vectors. Infected patients could progress to severe plasma leakage and hemorrhagic shock, where current standard treatment remains supportive. Previous reports suggested that several flavonoid derivatives inhibited mosquito-borne flaviviruses. This work aimed to explore sulfonamide chalcone derivatives as dengue inhibitors and to identify molecular targets. We initially screened 27 sulfonamide chalcones using cell-based antiviral and cytotoxic screenings. Two potential compounds, SC22 and SC27, were identified with DENV1-4 EC50s in the range of 0.71-0.94 and 3.15-4.46 µM, and CC50s at 14.63 and 31.02 µM, respectively. The compounds did not show any elevation in ALT or Cr in C57BL/6 mice on the 1st, 3rd, and 7th days after being administered intraperitoneally with 50 mg/kg SC22 or SC27 in a single dose. Moreover, the SAM-binding site of NS5 methyltransferase was a potential target of SC27 identified by computational and enzyme-based assays. The main target of SC22 was in a late stage of viral replication, but the exact target molecule had yet to be identified. In summary, a sulfonamide chalcone, SC27, was a potential DENV inhibitor that targeted viral methyltransferase. Further investigation should be the study of the structure-activity relationship of SC27 derivatives for higher potency and lower toxicity.

Design of electron-donating group substituted 2-PAM analogs as antidotes for organophosphate insecticide poisoning.

The use of organophosphate (OPs) pesticides is widespread in agriculture and horticulture, but these chemicals can be lethal to humans, causing fatalities and deaths each year. The inhibition of acetylcholinesterase (AChE) by OPs leads to the overstimulation of cholinergic receptors, ultimately resulting in respiratory arrest, seizures, and death. Although 2-pralidoxime (2-PAM) is the FDA-approved drug for treating OP poisoning, there is difficulty in blood-brain barrier permeation. To address this issue, we designed and evaluated a series of 2-PAM analogs by substituting electron-donating groups on the para and/or ortho positions of the pyridinium core using in silico techniques. Our PCM-ONIOM2 (MP2/6-31G*:PM7//B3LYP/6-31G*:UFF) binding energy results demonstrated that 13 compounds exhibited higher binding energy than 2-PAM. The analog with phenyl and methyl groups substituted on the para and ortho positions, respectively, showed the most favorable binding characteristics, with aromatic residues in the active site (Y124, W286, F297, W338, and Y341) and the catalytic residue S203 covalently bonding with paraoxon. The results of DS-MD simulation revealed a highly favorable apical conformation of the potent analog, which has the potential to enhance reactivation of AChE. Importantly, newly designed compound demonstrated appropriate drug-likeness properties and blood-brain barrier penetration. These results provide a rational guide for developing new antidotes to treat organophosphate insecticide toxicity.

Structure-yeast α-glucosidase inhibitory activity relationship of 9-O-berberrubine carboxylates.

Thirty-five 9-O-berberrubine carboxylate derivatives were synthesized and evaluated for yeast α-glucosidase inhibitory activity. All compounds demonstrated better inhibitory activities than the parent compounds berberine (BBR) and berberrubine (BBRB), and a positive control, acarbose. The structure-activity correlation study indicated that most of the substituents on the benzoate moiety such as methoxy, hydroxy, methylenedioxy, benzyloxy, halogen, trifluoromethyl, nitro and alkyl can contribute to the activities except multi-methoxy, fluoro and cyano. In addition, replacing benzoate with naphthoate, cinnamate, piperate or diphenylacetate also led to an increase in inhibitory activities except with phenyl acetate. 9, 26, 27, 28 and 33 exhibited the most potent α-glucosidase inhibitory activities with the IC50 values in the range of 1.61-2.67 µM. Kinetic study revealed that 9, 26, 28 and 33 interacted with the enzyme via competitive mode. These four compounds were also proved to be not cytotoxic at their IC50 values. The competitive inhibition mechanism of these four compounds against yeast α-glucosidase was investigated using molecular docking and molecular dynamics simulations. The binding free energy calculations suggest that 26 exhibited the strongest binding affinity, and its binding stability is supported by hydrophobic interactions with D68, F157, F158 and F177. Therefore, 9, 26, 28 and 33 would be promising candidates for further studies of antidiabetic activity.

Deciphering the Potential of Multidimensional Carbon Materials for Surface-Enhanced Raman Spectroscopy through Density Functional Theory.

Surface-enhanced Raman spectroscopy (SERS) is a potent analytical tool, particularly for molecular identification and structural analysis. Conventional metallic SERS substrates, however, suffer from low reproducibility and compatibility with biological molecules. Recently, metal-free SERS substrates based on chemical enhancement have emerged as a promising alternative with carbon-based materials offering excellent reproducibility and compatibility. Nevertheless, our understanding of carbon materials in SERS remains limited, which hinders their rational design. Here we systematically explore multidimensional carbon materials, including zero-dimensional fullerenes (C60), one-dimensional carbon nanotubes, two-dimensional graphene, and their B-, N-, and O-doped derivatives, for SERS applications. Using density functional theory, we elucidate the nonresonant polarizability-enhanced and resonant charge-transfer-based chemical enhancement mechanisms of these materials by evaluating their static/dynamic polarizability and electron excitation properties. This work provides a critical reference for the future design of carbon-based SERS substrates, opening a new avenue in this field.

The 8-bromobaicalein alleviated chikungunya-induced musculoskeletal inflammation and reduced the viral load in healthy adult mice.

Chikungunya virus is a re-emerging arbovirus that has caused epidemic outbreaks in recent decades. Patients in older age groups with high viral load and severe immunologic response during acute infection are likely to develop chronic arthritis and severe joint pain. Currently, no antiviral drug is available. Previous studies suggested that a flavone derivative, 8-bromobaicalein, was a potential dengue and Zika replication inhibitor in a cell-based system targeting flaviviral polymerase. Here we characterized that 8-bromobaicalein inhibited chikungunya virus replication with EC50 of 0.49 ± 0.11 µM in Vero cells. The molecular target predicted at viral nsP1 methyltransferase using molecular binding and fragment molecular orbital calculation. Additionally, oral administration of 250 mg/kg twice daily treatment alleviated chikungunya-induced musculoskeletal inflammation and reduced viral load in healthy adult mice. Pharmacokinetic analysis indicated that the 250 mg/kg administration maintained the compound level above EC99.9 for 12 h. Therefore, 8-bromobaicalein should be a potential candidate for further development as a pan-arboviral drug.

Origin of Homochirality in Amino Acids Induced by Lyman-α Irradiation in the Early Stage of the Milky Way.

The enantiomeric excess (ee) of l-form amino acids found in the Murchison meteorite poses some issues about the cosmic origin of their chirality. Circular dichroism (CD) spectra of amino acids in the far-ultraviolet (FUV) at around 6.8 eV (182 nm) indicate that the circularly polarized light can induce ee through photochemical reactions. Here, we resort to ab initio calculations to extract the CD spectra up to the vacuum-ultraviolet (VUV) region (∼11 eV), and we propose a novel equation to compute the ee applicable to a wider range of light frequency than what is available to date. This allows us to show that the strength of the induced ee (|ee|) in the 10 eV VUV region is comparable to the one in the 6.8 eV FUV region. This feature is common for some key amino acids (alanine, 2-aminobutyric acid, and valine). In space, intense Lyman-α (Lyα) light of 10.2 eV is emitted from star forming regions. This study provides a theoretical basis that Lyα emitter from an early starburst in the Milky Way plays a crucial role in initiating the ee of amino acids.