Effect of temperature on the structure and drug-release behaviour of inclusion complex of ß-cyclodextrin with cyclophosphamide: a molecular dynamics study.

Cyclodextrins (CDs) are suitable drug carriers because of their doughnut-shaped cavities with hydrophilic outer and hydrophobic inner surfaces. Temperature-responsive CD-based drug carriers are expected to be one of the most promising candidates for drug delivery systems. In this study, we performed molecular dynamics simulations of the inclusion complex of ß-CD with cyclophosphamide (CP) at temperatures from 300 K to 400 K to investigate the temperature dependency of the release behaviour of CP and structural changes of ß-CD in an aqueous solution. We analysed the distance between the centres of mass of ß-CD and CP and the radius of gyration of ß-CD. The CP molecule was released from the ß-CD cavity at 400 K, whereas two different inclusion complexes, partially and completely, were observed at T < 400 K. ß-CD encapsulating a CP molecule had a more spherical shape and rigidity than ß-CD without a CP, and the rigidity of their inclusion complex decreased with increasing temperature. Our findings provide fundamental insights into the behaviours of the ß-CD/CP complex and drug release at the molecular level and can facilitate the development of new temperature-responsive drug delivery systems with CD nanocarriers triggered by localised temperature increases using focused ultrasound.

In Silico Identification of Inhibitory Compounds for SARS-Cov-2 Papain-Like Protease.

Virtual screening with high-performance computers is a powerful and cost-effective technique in drug discovery. A chemical database is searched to find candidate compounds firmly bound to a target protein, judging from the binding poses and/or binding scores. The severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) infectious disease has spread worldwide for the last three years, causing severe slumps in economic and social activities. SARS-Cov-2 has two viral proteases: 3-chymotrypsin-like (3CL) and papain-like (PL) protease. While approved drugs have already been released for the 3CL protease, no approved agent is available for PL protease. In this work, we carried out in silico screening for the PL protease inhibitors, combining docking simulation and molecular mechanics calculation. Docking simulations were applied to 8,820 molecules in a chemical database of approved and investigational compounds. Based on the binding poses generated by the docking simulations, molecular mechanics calculations were performed to optimize the binding structures and to obtain the binding scores. Based on the binding scores, 57 compounds were selected for in vitro assay of the inhibitory activity. Five inhibitory compounds were identified from the in vitro measurement. The predicted binding structures of the identified five compounds were examined, and the significant interaction between the individual compound and the protease catalytic site was clarified. This work demonstrates that computational virtual screening by combining docking simulation with molecular mechanics calculation is effective for searching candidate compounds in drug discovery.

The force-sensing device region of α-catenin is an intrinsically disordered segment in the absence of intramolecular stabilization of the autoinhibitory form.

Mechanotransduction by α-catenin facilitates the force-dependent development of adherens junctions (AJs) by recruiting vinculin to reinforce actin anchoring of AJs. The α-catenin mechanotransducing action is facilitated by its force-sensing device region that autoinhibits the vinculin-binding site 1 (VBS1). Here, we report the high-resolution structure of the force-sensing device region of α-catenin, which shows the autoinhibited form comprised of helix bundles E, F and G. The cryptic VBS1 is embedded into helix bundle E stabilized by direct interactions with the autoinhibitory region forming helix bundles F and G. Our molecular dissection study showed that helix bundles F and G are stable in solution in each isolated form, whereas helix bundle E that contains VBS1 is unstable and intrinsically disordered in solution in the isolated form. We successfully identified key residues mediating the autoinhibition and produced mutated α-catenins that display variable force sensitivity and autoinhibition. Using these mutants, we demonstrate both in vitro and in vivo that, in the absence of this stabilization, the helix bundle containing VBS1 would adopt an unfolded form, thus exposing VBS for vinculin binding. We provide evidence for importance of mechanotransduction with the intrinsic force sensitivity for vinculin recruitment to adherens junctions of epithelial cell sheets with mutated α-catenins.

Influence of an Interfacial Effect on the Laser Performance of a Rhodamine 6G/Cellulose Acetate Waveguide on a Vinylidene Fluoride Copolymer Layer.

The fluorescent properties of dyes and fluorophores in condensed matter significantly affect the laser performance of organic dye lasers and fluorescent polymer lasers. Concentration quenching of fluorescence is commonly observed in condensed matter. Several approaches have been presented to suppress such quenching, such as the use of a dendrimer and the use of effective energy transfer in a guest-host system. The enhanced fluorescence of rhodamine 6G (R6G) dye on a vinylidene fluoride polymer is an alternative method for enhancing laser performance because of the roughness of the P(VDF-TrFE) surface and the interaction between polar ß-crystals of P(VDF-TrFE) and R6G dye. In this paper, a significant improvement in slope efficiency (SE) is demonstrated without a significant depression in the lasing threshold for distributed feedback (DFB) and distributed Bragg reflector (DBR) lasers fabricated using an R6G-dispersed cellulose acetate (CA) matrix spin-coated onto a copolymer of vinylidene fluoride and trifluoroethylene P(VDF-TrFE) thin film. SEs of 3.4 and 1.3% were measured for DBR and DFB laser devices with CA/R6G on a P(VDF-TrFE) thin film, respectively, whereas an SE of less than 1.0% was measured for both corresponding laser devices without a P(VDF-TrFE) thin film. From the aspect of simple fabrication procedures, repeatability in device fabrication and performance, stability of the device, time for device fabrication, the present approach is the most preferable way for industrial applications, requiring only the additional step of spin-coating of a P(VDF-TrFE) thin film.

Efficient expression of SRK intracellular domain by a modeling-based protein engineering.

S-locus protein kinase (SRK) is a receptor kinase that plays a critical role in self-recognition in the Brassicaceae self-incompatibility (SI) response. SRK is activated by binding of its ligand S-locus protein 11 (SP11) and subsequently induced phosphorylation of the intracellular kinase domain. However, a detailed activation mechanism of SRK is still largely unknown because of the difficulty in stably expressing SRK recombinant proteins. Here, we performed modeling-based protein engineering of the SRK kinase domain for stable expression in Escherichia coli. The engineered SRK intracellular domain was expressed about 54-fold higher production than wild type SRK, without loss of the kinase activity, suggesting it could be useful for further biochemical and structural studies.

Fluorogenic Enhancement of an in Vitro-Selected Peptide Ligand by Replacement of a Fluorescent Group.

To prepare a fluorogenic peptide ligand which binds to an arbitrary target, we previously succeeded in seeking a fluorogenic ligand to calmodulin using in vitro selection. In this study the environment-sensitive fluorescent group in the selected peptide ligand was replaced with other fluorescent groups to find the possibility to increase the fluorogenic activity. Surface plasmon resonance measurement exhibited that the binding affinity was held even after the replacement. However, the replacement significantly affected the fluorogenic activity. It depended on the kind of incorporated fluorophors and linker length. As a result, the incorporation of 4-N,N-dimethylamino-1,8-naphthalimide enhanced the fluorescence intensity over 100-fold in the presence of target calcium-bound calmodulin. This study demonstrated that the functionality of in vitro selected peptide can be tuned with keeping the binding affinity.

Structural and functional analysis of the yeast N-acetyltransferase Mpr1 involved in oxidative stress tolerance via proline metabolism.

Mpr1 (sigma1278b gene for proline-analog resistance 1), which was originally isolated as N-acetyltransferase detoxifying the proline analog L-azetidine-2-carboxylate, protects yeast cells from various oxidative stresses. Mpr1 mediates the L-proline and L-arginine metabolism by acetylating L-Δ(1)-pyrroline-5-carboxylate, leading to the L-arginine-dependent production of nitric oxide, which confers oxidative stress tolerance. Mpr1 belongs to the Gcn5-related N-acetyltransferase (GNAT) superfamily, but exhibits poor sequence homology with the GNAT enzymes and unique substrate specificity. Here, we present the X-ray crystal structure of Mpr1 and its complex with the substrate cis-4-hydroxy-L-proline at 1.9 and 2.3 Å resolution, respectively. Mpr1 is folded into α/ß-structure with eight-stranded mixed ß-sheets and six α-helices. The substrate binds to Asn135 and the backbone amide of Asn172 and Leu173, and the predicted acetyl-CoA-binding site is located near the backbone amide of Phe138 and the side chain of Asn178. Alanine substitution of Asn178, which can interact with the sulfur of acetyl-CoA, caused a large reduction in the apparent kcat value. The replacement of Asn135 led to a remarkable increase in the apparent Km value. These results indicate that Asn178 and Asn135 play an important role in catalysis and substrate recognition, respectively. Such a catalytic mechanism has not been reported in the GNAT proteins. Importantly, the amino acid substitutions in these residues increased the L-Δ(1)-pyrroline-5-carboxylate level in yeast cells exposed to heat stress, indicating that these residues are also crucial for its physiological functions. These studies provide some benefits of Mpr1 applications, such as the breeding of industrial yeasts and the development of antifungal drugs.

Structural basis of cargo recognition by the myosin-X MyTH4-FERM domain.

Myosin-X is an important unconventional myosin that is critical for cargo transportation to filopodia tips and is also utilized in spindle assembly by interacting with microtubules. We present a series of structural and biochemical studies of the myosin-X tail domain cassette, consisting of myosin tail homology 4 (MyTH4) and FERM domains in complex with its specific cargo, a netrin receptor DCC (deleted in colorectal cancer). The MyTH4 domain is folded into a helical VHS-like structure and is associated with the FERM domain. We found an unexpected binding mode of the DCC peptide to the subdomain C groove of the FERM domain, which is distinct from previously reported ß-ß associations found in radixin-adhesion molecule complexes. We also revealed direct interactions between the MyTH4-FERM cassette and tubulin C-terminal acidic tails, and identified a positively charged patch of the MyTH4 domain, which is involved in tubulin binding. We demonstrated that both DCC and integrin bindings interfere with microtubule binding and that DCC binding interferes with integrin binding. Our results provide the molecular basis by which myosin-X facilitates alternative dual binding to cargos and microtubules.

Mechanism of high-affinity abscisic acid binding to PYL9/RCAR1.

Arabidopsis receptors of abscisic acid (ABA), the key plant hormone for adaptation to water stress, comprise 14 PYR/PYLs/RCARs proteins classified into three subfamilies I, II, and III, which suggests functional differentiation. Although their monomer-dimer equilibria may be correlated with differences in their ABA-binding affinities, how the dimerization decreases the affinity is unclear. Comparative structural and binding studies between PYL9, which is a representative of high-affinity subfamily I, and low-affinity members of subfamily III reveals that the nonpolar triplet (Ile110, Val162, and Leu165) and Pro64 contribute to enhance ABA-binding affinity by inducing a shift of the ABA carboxyl group to form additional direct hydrogen bonds with conserved Asn169. Our mutation studies of PYL1 successfully produced a monomeric mutant PYL1 exhibiting low ABA affinity and also a dimeric mutant PYL1 exhibiting high ABA-binding affinity, suggesting that dimer formation of ABA receptors is not essential for their low ABA-binding affinity. Our study contributes toward establishing the structural basis for the higher ABA-binding affinity of the subfamily receptors and provides a clue for understanding the broad spectrum of hormone actions in plants manifested by the different hormone-binding affinity of multiple receptors.

Structures of D14 and D14L in the strigolactone and karrikin signaling pathways.

Strigolactones (SLs) are plant hormones that inhibit shoot branching. DWARF14 (D14) inhibits rice tillering and is an SL receptor candidate in the branching inhibition pathway, whereas the close homologue DWARF14-LIKE (D14L) participates in the signaling pathway of karrikins (KARs), which are derived from burnt vegetation as smoke stimulants of seed germination. We provide the first evidence for direct binding of the bioactive SL analogue GR24 to D14. Isothermal titration calorimetry measurements show a D14-GR24 binding affinity in the sub-micromolar range. Similarly, bioactive KAR1 directly binds D14L in the micromolar range. The crystal structure of rice D14 shows a compact α-/ß-fold hydrolase domain forming a deep ligand-binding pocket capable of accommodating GR24. Insertion of four α-helices between ß6 strand and αD helix forms the helical cap of the pocket, although the pocket is open to the solvent. The pocket contains the conserved catalytic triad Ser-His-Asp aligned with the oxyanion hole, suggesting hydrolase activity. Although these structural characteristics are conserved in D14L, the D14L pocket is smaller than that of D14. The KAR-insensitive mutation kai2-1 is located at the prominent long ß6-αD1 loop, which is characteristic in D14 and D14L, but not in related α-/ß-fold hydrolases.

Gibberellin-induced DELLA recognition by the gibberellin receptor GID1.

Gibberellins control a range of growth and developmental processes in higher plants and have been widely used in the agricultural industry. By binding to a nuclear receptor, GIBBERELLIN INSENSITIVE DWARF1 (GID1), gibberellins regulate gene expression by promoting degradation of the transcriptional regulator DELLA proteins, including GIBBERELLIN INSENSITIVE (GAI). The precise manner in which GID1 discriminates and becomes activated by bioactive gibberellins for specific binding to DELLA proteins remains unclear. Here we present the crystal structure of a ternary complex of Arabidopsis thaliana GID1A, a bioactive gibberellin and the amino-terminal DELLA domain of GAI. In this complex, GID1A occludes gibberellin in a deep binding pocket covered by its N-terminal helical switch region, which in turn interacts with the DELLA domain containing DELLA, VHYNP and LExLE motifs. Our results establish a structural model of a plant hormone receptor that is distinct from the mechanism of the hormone perception and effector recognition of the known auxin receptors.

Cryo-EM analysis reveals human SID-1 transmembrane family member 1 dynamics underlying lipid hydrolytic activity.

Two mammalian homologs of systemic RNA interference defective protein 1 (SID-1) (SIDT1/2) are suggested to function as double-stranded RNA (dsRNA) transporters for extracellular dsRNA uptake or for release of incorporated dsRNA from lysosome to cytoplasm. SIDT1/2 is also suggested to be involved in cholesterol transport and lipid metabolism. Here, we determine the cryo-electron microscopy structures of human SIDT1, homodimer in a side-by-side arrangement, with two distinct conformations, the cholesterol-bound form and the unbound form. Our structures reveal that the membrane-spanning region of SIDT1 harbors conserved histidine and aspartate residues coordinating to putative zinc ion, in a structurally similar manner to alkaline ceramidases or adiponectin receptors that require zinc for ceramidase activity. We identify that SIDT1 has a ceramidase activity that is attenuated by cholesterol binding. Observations from two structures suggest that cholesterol molecules serve as allosteric regulator that binds the transmembrane region of SIDT1 and induces the conformation change and the reorientation of the catalytic residues. This study represents a contribution to the elucidation of the cholesterol-mediated mechanisms of lipid hydrolytic activity and RNA transport in the SID-1 family proteins.

Structural features of aquaporin 4 supporting the formation of arrays and junctions in biomembranes.

A limited class of aquaporins has been described to form regular arrays and junctions in membranes. The biological significance of these structures, however, remains uncertain. Here we analyze the underlying physical principles with the help of a computational procedure that takes into account protein-protein as well as protein-membrane interactions. Experimentally observed array/junction structures are systematically (dis)assembled and major driving forces identified. Aquaporin 4 was found to be markedly different from the non-junction forming aquaporin 1. The environmental stabilization resulting from embedding into the biomembrane was identified as the main driving force. This highlights the role of protein-membrane interactions in aquaporin 4. Analysis of the type presented here can help to decipher the biological role of membrane arrays and junctions formed by aquaporin.

Common force field thermodynamics of cholesterol.

Four different force fields are examined for dynamic characteristics using cholesterol as a case study. The extent to which various types of internal degrees of freedom become thermodynamically relevant is evaluated by means of principal component analysis. More complex degrees of freedom (angle bending, dihedral rotations) show a trend towards force field independence. Moreover, charge assignments for membrane-embedded compounds are revealed to be critical with significant impact on biological reasoning.

Wetting hysteresis induces effective unidirectional water transport through a fluctuating nanochannel.

We propose a water pump that actively transports water molecules through nanochannels. Spatially asymmetric noise fluctuations imposed on the channel radius cause unidirectional water flow without osmotic pressure, which can be attributed to hysteresis in the cyclic transition between the wetting/drying states. We show that the water transport depends on fluctuations, such as white, Brownian, and pink noises. Because of the high-frequency components in white noise, fast switching of open and closed states inhibits channel wetting. Conversely, pink and Brownian noises generate high-pass filtered net flow. Brownian fluctuation leads to a faster water transport rate, whereas pink noise has a higher capability to overcome pressure differences in the opposite direction. A trade-off relationship exists between the resonant frequency of the fluctuation and the flow amplification. The proposed pump can be considered as an analogy for the reversed Carnot cycle, which is the upper limit of the energy conversion efficiency.

Unsupervised deep learning for molecular dynamics simulations: a novel analysis of protein-ligand interactions in SARS-CoV-2 Mpro.

Molecular dynamics (MD) simulations, which are central to drug discovery, offer detailed insights into protein-ligand interactions. However, analyzing large MD datasets remains a challenge. Current machine-learning solutions are predominantly supervised and have data labelling and standardisation issues. In this study, we adopted an unsupervised deep-learning framework, previously benchmarked for rigid proteins, to study the more flexible SARS-CoV-2 main protease (Mpro). We ran MD simulations of Mpro with various ligands and refined the data by focusing on binding-site residues and time frames in stable protein conformations. The optimal descriptor chosen was the distance between the residues and the center of the binding pocket. Using this approach, a local dynamic ensemble was generated and fed into our neural network to compute Wasserstein distances across system pairs, revealing ligand-induced conformational differences in Mpro. Dimensionality reduction yielded an embedding map that correlated ligand-induced dynamics and binding affinity. Notably, the high-affinity compounds showed pronounced effects on the protein's conformations. We also identified the key residues that contributed to these differences. Our findings emphasize the potential of combining unsupervised deep learning with MD simulations to extract valuable information and accelerate drug discovery.

Bisabosqual A: A novel asparagine synthetase inhibitor suppressing the proliferation and migration of human non-small cell lung cancer A549 cells.

Asparagine synthetase (ASNS) is a crucial enzyme for the de novo biosynthesis of endogenous asparagine (Asn), and ASNS shows the positive relationship with the growth of several solid tumors. Most of ASNS inhibitors are analogs of transition-state in ASNS reaction, but their low cell permeability hinders their anticancer activity. Therefore, novel ASNS inhibitors with a new pharmacophore urgently need to be developed. In this study, we established and applied a system for in vitro screening of ASNS inhibitors, and found a promising unique bisabolane-type meroterpenoid molecule, bisabosqual A (Bis A), able to covalently modify K556 site of ASNS protein. Bis A targeted ASNS to suppress cell proliferation of human non-small cell lung cancer A549 cells and exhibited a synergistic effect with L-asparaginase (L-ASNase). Mechanistically, Bis A promoted oxidative stress and apoptosis, while inhibiting autophagy, cell migration and epithelial-mesenchymal transition (EMT), impeding cancer cell development. Moreover, Bis A induced negative feedback pathways containing the GCN2-eIF2α-ATF4, PI3K-AKT-mTORC1 and RAF-MEK-ERK axes, but combination treatment of Bis A and rapamycin/torin-1 overcame the potential drug resistance triggered by mTOR pathways. Our study demonstrates that ASNS inhibition is promising for cancer chemotherapy, and Bis A is a potential lead ASNS inhibitor for anticancer development.

A small molecule iCDM-34 identified by in silico screening suppresses HBV DNA through activation of aryl hydrocarbon receptor.

IFN-alpha have been reported to suppress hepatitis B virus (HBV) cccDNA via APOBEC3 cytidine deaminase activity through interferon signaling. To develop a novel anti-HBV drug for a functional cure, we performed in silico screening of the binding compounds fitting the steric structure of the IFN-alpha-binding pocket in IFNAR2. We identified 37 compounds and named them in silico cccDNA modulator (iCDM)-1-37. We found that iCDM-34, a new small molecule with a pyrazole moiety, showed anti-HCV and anti-HBV activities. We measured the anti-HBV activity of iCDM-34 dependent on or independent of entecavir (ETV). iCDM-34 suppressed HBV DNA, pgRNA, HBsAg, and HBeAg, and also clearly exhibited additive inhibitory effects on the suppression of HBV DNA with ETV. We confirmed metabolic stability of iCDM-34 was stable in human liver microsomal fraction. Furthermore, anti-HBV activity in human hepatocyte-chimeric mice revealed that iCDM-34 was not effective as a single reagent, but when combined with ETV, it suppressed HBV DNA compared to ETV alone. Phosphoproteome and Western blotting analysis showed that iCDM-34 did not activate IFN-signaling. The transcriptome analysis of interferon-stimulated genes revealed no increase in expression, whereas downstream factors of aryl hydrocarbon receptor (AhR) showed increased levels of the expression. CDK1/2 and phospho-SAMHD1 levels decreased under iCDM-34 treatment. In addition, AhR knockdown inhibited anti-HCV activity of iCDM-34 in HCV replicon cells. These results suggest that iCDM-34 decreases the phosphorylation of SAMHD1 through CDK1/2, and suppresses HCV replicon RNA, HBV DNA, and pgRNA formation.

Differences in ligand-induced protein dynamics extracted from an unsupervised deep learning approach correlate with protein-ligand binding affinities.

Prediction of protein-ligand binding affinity is a major goal in drug discovery. Generally, free energy gap is calculated between two states (e.g., ligand binding and unbinding). The energy gap implicitly includes the effects of changes in protein dynamics induced by ligand binding. However, the relationship between protein dynamics and binding affinity remains unclear. Here, we propose a method that represents ligand-binding-induced protein behavioral change with a simple feature that can be used to predict protein-ligand affinity. From unbiased molecular simulation data, an unsupervised deep learning method measures the differences in protein dynamics at a ligand-binding site depending on the bound ligands. A dimension reduction method extracts a dynamic feature that strongly correlates to the binding affinities. Moreover, the residues that play important roles in protein-ligand interactions are specified based on their contribution to the differences. These results indicate the potential for binding dynamics-based drug discovery.

Non-Gaussian fluctuations resulting from power-law trapping in a lipid bilayer.

Anomalous diffusion in lipid bilayers is usually attributed to viscoelastic behavior. We compute the scaling exponent of relative fluctuations of the time-averaged mean square displacement in a lipid bilayer, by using a molecular dynamics simulation. According to the continuous time random walk theory, this exponent indicates non-Gaussian behavior caused by a power-law trapping time. Our results provide the first evidence that a lipid bilayer has not only viscoelastic properties but also trapping times distributed according to a power law.