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

In-layer inhomogeneity of molecular dynamics in quasi-liquid layers of ice

Quasi-liquid layers (QLLs) are present on the surface of ice and play a significant role in its distinctive chemical and physical properties. These layers exhibit considerable heterogeneity across different scales ranging from nanometers to millimeters. Although the formation of partially ice-like structures has been proposed, the molecular-level understanding of this heterogeneity remains unclear. Here, we examined the heterogeneity of molecular dynamics on QLLs based on molecular dynamics simulations and machine learning analysis of the simulation data. We demonstrated that the molecular dynamics of QLLs do not comprise a mixture of solid- and liquid water molecules. Rather, molecules having similar behaviors form dynamical domains that are associated with the dynamical heterogeneity of supercooled water. Nonetheless, molecules in the domains frequently switch their dynamical state. Furthermore, while there is no observable characteristic domain size, the long-range ordering strongly depends on the temperature and crystal face. Instead of a mixture of static solid- and liquid-like regions, our results indicate the presence of heterogeneous molecular dynamics in QLLs, which offers molecular-level insights into the surface properties of ice.

Graph-Neural-Network-Based Unsupervised Learning of the Temporal Similarity of Structural Features Observed in Molecular Dynamics Simulations

Classification of molecular structures is a crucial step in molecular dynamics (MD) simulations to detect various structures and phases within systems. Molecular structures, which are commonly identified using order parameters, were recently identified using machine learning (ML), that is, the ML models acquire structural features using labeled crystals or phases via supervised learning. However, these approaches may not identify unlabeled or unknown structures, such as the imperfect crystal structures observed in nonequilibrium systems and interfaces. In this study, we proposed the use of a novel unsupervised learning framework, denoted temporal self-supervised learning (TSSL), to learn structural features and design their parameters. In TSSL, the ML models learn that the structural similarity is learned via contrastive learning based on minor short-term variations caused by perturbations in MD simulations. This learning framework is applied to a sophisticated architecture of graph neural network models that use bond angle and length data of the neighboring atoms. TSSL successfully classifies water and ice crystals based on high local ordering, and furthermore, it detects imperfect structures typical of interfaces such as the water-ice and ice-vapor interfaces.

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.

Statistical Improvement of rGILCC 1 and rPOXA 1B Laccases Activity Assay Conditions Supported by Molecular Dynamics

Laccases (E.C. 1.10.3.2) are glycoproteins widely distributed in nature. Their structural conformation includes three copper sites in their catalytic center, which are responsible for facilitating substrate oxidation, leading to the generation of H2O instead of H2O2. The measurement of laccase activity (UL-1) results may vary depending on the type of laccase, buffer, redox mediators, and substrates employed. The aim was to select the best conditions for rGILCC 1 and rPOXA 1B laccases activity assay. After sequential statistical assays, the molecular dynamics proved to support this process, and we aimed to accumulate valuable insights into the potential application of these enzymes for the degradation of novel substrates with negative environmental implications. Citrate buffer treatment T2 (CB T2) (pH 3.0 ± 0.2; λ420nm, 2 mM ABTS) had the most favorable results, with 7.315 ± 0.131 UL-1 for rGILCC 1 and 5291.665 ± 45.83 UL-1 for rPOXA 1B. The use of citrate buffer increased the enzyme affinity for ABTS since lower Km values occurred for both enzymes (1.49 × 10-2 mM for rGILCC 1 and 3.72 × 10-2 mM for rPOXA 1B) compared to those obtained in acetate buffer (5.36 × 10-2 mM for rGILCC 1 and 1.72 mM for rPOXA 1B). The molecular dynamics of GILCC 1-ABTS and POXA 1B-ABTS showed stable behavior, with root mean square deviation (RMSD) values not exceeding 2.0 Å. Enzyme activities (rGILCC 1 and rPOXA 1B) and 3D model-ABTS interactions (GILCC 1-ABTS and POXA 1B-ABTS) were under the strong influence of pH, wavelength, ions, and ABTS concentration, supported by computational studies identifying the stabilizing residues and interactions. Integration of the experimental and computational approaches yielded a comprehensive understanding of enzyme-substrate interactions, offering potential applications in environmental substrate treatments.

Regulation of Tumor Apoptosis of Poriae cutis-Derived Lanostane Triterpenes by AKT/PI3K and MAPK Signaling Pathways In Vitro

Poria cocos is traditionally used as both food and medicine. Triterpenoids in Poria cocos have a wide range of pharmacological activities, such as diuretic, sedative and tonic properties. In this study, the anti-tumor activities of poricoic acid A (PAA) and poricoic acid B (PAB), purified by high-speed counter-current chromatography, as well as their mechanisms and signaling pathways, were investigated using a HepG2 cell model. After treatment with PAA and PAB on HepG2 cells, the apoptosis was obviously increased (p < 0.05), and the cell cycle arrested in the G2/M phase. Studies showed that PAA and PAB can also inhibit the occurrence and development of tumor cells by stimulating the generation of ROS in tumor cells and inhibiting tumor migration and invasion. Combined Polymerase Chain Reaction and computer simulation of molecular docking were employed to explore the mechanism of tumor proliferation inhibition by PAA and PAB. By interfering with phosphatidylinositol-3-kinase/protein kinase B, Mitogen-activated protein kinases and p53 signaling pathways; and further affecting the expression of downstream caspases; matrix metalloproteinase family, cyclin-dependent kinase -cyclin, Intercellular adhesion molecules-1, Vascular Cell Adhesion Molecule-1 and Cyclooxygenase -2, may be responsible for their anti-tumor activity. Overall, the results suggested that PAA and PAB induced apoptosis, halted the cell cycle, and inhibited tumor migration and invasion through multi-pathway interactions, which may serve as a potential therapeutic agent against cancer.

Rational design of small-molecule responsive protein switches

Small-molecule responsive protein switches are powerful tools for controlling cellular processes. These switches are designed to respond rapidly and specifically to their inducer. They have been used in numerous applications, including the regulation of gene expression, post-translational protein modification, and signal transduction. Typically, small-molecule responsive protein switches consist of two proteins that interact with each other in the presence or absence of a small molecule. Recent advances in computational protein design already contributed to the development of protein switches with an expanded range of small-molecule inducers and increasingly sophisticated switch mechanisms. Further progress in the engineering of small-molecule responsive switches is fueled by cutting-edge computational design approaches, which will enable more complex and precise control over cellular processes and advance synthetic biology applications in biotechnology and medicine. Here, we discuss recent milestones and how technological advances are impacting the development of chemical switches.

Prediction of the Peptide-TIM3 Binding Site in Inhibiting TIM3-Galectin 9 Binding Pathways

T-cell immunoglobulin and mucin domain-containing protein-3 (TIM3) is an important receptor protein that modulates the immune system. The binding of TIM3 with Galectin 9 (GAL9) triggers immune system suppression, but the TIM3-GAL9 binding can be inhibited by binding of the peptide P26 to TIM3. A fast and accurate prediction of the P26-TIM3 binding site is crucial and a prerequisite for the investigation of P26-TIM3 interactions and TIM3-GAL9 binding pathways. Here, we present a machine learning approach, which considers protein conformational changes, to quickly identify the ligand-binding site on TIM3. Our results show that the P26 binding site is located near the C″-D loop of TIM3. Further simulations show that the binding pose is stabilized by a variety of electrostatic and hydrophobic interactions. Binding of P26 can alter the conformations of nearby glycan side chains on TIM3, providing possible mechanisms of how P26 inhibits TIM3-GAL9 binding pathways. The insights from this work will facilitate the identification of other peptides or antibodies that may also inhibit the TIM3-GAL9 pathways and eventually lead to improved attempts in the modulation of the TIM3-GAL9 immunosuppression pathways. The strategies and machine learning method can be generalized to study ligand-receptor binding when the conformational changes during the binding are important.

DataDTA: a multi-feature and dual-interaction aggregation framework for drug-target binding affinity prediction

MOTIVATION

Accurate prediction of drug-target binding affinity (DTA) is crucial for drug discovery. The increase in the publication of large-scale DTA datasets enables the development of various computational methods for DTA prediction. Numerous deep learning-based methods have been proposed to predict affinities, some of which only utilize original sequence information or complex structures, but the effective combination of various information and protein-binding pockets have not been fully mined. Therefore, a new method that integrates available key information is urgently needed to predict DTA and accelerate the drug discovery process.

RESULTS

In this study, we propose a novel deep learning-based predictor termed DataDTA to estimate the affinities of drug-target pairs. DataDTA utilizes descriptors of predicted pockets and sequences of proteins, as well as low-dimensional molecular features and SMILES strings of compounds as inputs. Specifically, the pockets were predicted from the three-dimensional structure of proteins and their descriptors were extracted as the partial input features for DTA prediction. The molecular representation of compounds based on algebraic graph features was collected to supplement the input information of targets. Furthermore, to ensure effective learning of multiscale interaction features, a dual-interaction aggregation neural network strategy was developed. DataDTA was compared with state-of-the-art methods on different datasets, and the results showed that DataDTA is a reliable prediction tool for affinities estimation. Specifically, the concordance index (CI) of DataDTA is 0.806 and the Pearson correlation coefficient (R) value is 0.814 on the test dataset, which is higher than other methods.

AVAILABILITY AND IMPLEMENTATION

The codes and datasets of DataDTA are available at https://github.com/YanZhu06/DataDTA.

Molecular Dynamics Simulation of the Complex of PDE5 and Evodiamine

Alzheimer's disease is an irreversible neurological disorder for which there are no effective small molecule therapeutics. A phosphodiesterase 5 (PDE5) inhibitor is a candidate medicine for the treatment of Alzheimer's disease. Rutaecarpine, an indole alkaloid found in Euodiae Fructus, has inhibitory activity for PDE5. Euodiae Fructus contains more evodiamine than rutaecarpine. Therefore, we performed molecular dynamics simulations of the complex of PDE5 and evodiamine. The results showed that the PDE5 and (-)-evodiamine complexes were placed inside the reaction center compared to the case of PDE5 and (+)-evodiamine complex. The binding of (-)-evodiamine to PDE5 increased the root-mean-square deviation and radius of gyration of PDE5. In the PDE5 with (-)-evodiamine complex, the value of the root-mean-square fluctuation of the M-loop, which is thought to be important for activity, increased. This result suggests that (-)-evodiamine may have inhibitory activity.

Combining Molecular Dynamics and Machine Learning to Analyze Shear Thinning for Alkane and Globular Lubricants in the Low Shear Regime

Lubricants with desirable frictional properties are important in achieving an energy-saving society. Lubricants at the interfaces of mechanical components are confined under high shear rates and pressures and behave quite differently from the bulk material. Computational approaches such as nonequilibrium molecular dynamics (NEMD) simulations have been performed to probe the molecular behavior of lubricants. However, the low-shear-velocity regions of the materials have rarely been simulated owing to the expensive calculations necessary to do so, and the molecular dynamics under shear velocities comparable with that in the experiments are not clearly understood. In this study, we performed NEMD simulations of extremely confined lubricants, i.e., two molecular layers for four types of lubricants confined in mica walls, under shear velocities from 0.001 to 1 m/s. While we confirmed shear thinning, the velocity profiles could not show the flow behavior when the shear velocity was much slower than thermal fluctuations. Therefore, we used an unsupervised machine learning approach to detect molecular movements that contribute to shear thinning. First, we extracted the simple features of molecular movements from large amounts of MD data, which were found to correlate with the effective viscosity. Subsequently, the extracted features were interpreted by examining the trajectories contributing to these features. The magnitude of diffusion corresponded to the viscosity, and the location of slips that varied depending on the spherical and chain lubricants was irrelevant. Finally, we attempted to apply a modified Stokes-Einstein relation at equilibrium to the nonequilibrium and confined systems. While systems with low shear rates obeyed the relation sufficiently, large deviations were observed under large shear rates.