Machine learning-based QSAR and LB-PaCS-MD guided design of SARS-CoV-2 main protease inhibitors

The global outbreak of the COVID-19 pandemic caused by the SARS-CoV-2 virus had led to profound respiratory health implications. This study focused on designing organoselenium-based inhibitors targeting the SARS-CoV-2 main protease (Mpro). The ligand-binding pathway sampling method based on parallel cascade selection molecular dynamics (LB-PaCS-MD) simulations was employed to elucidate plausible paths and conformations of ebselen, a synthetic organoselenium drug, within the Mpro catalytic site. Ebselen effectively engaged the active site, adopting proximity to H41 and interacting through the benzoisoselenazole ring in a π-π T-shaped arrangement, with an additional π-sulfur interaction with C145. In addition, the ligand-based drug design using the QSAR with GFA-MLR, RF, and ANN models were employed for biological activity prediction. The QSAR-ANN model showed robust statistical performance, with an r2training exceeding 0.98 and an RMSEtest of 0.21, indicating its suitability for predicting biological activities. Integration the ANN model with the LB-PaCS-MD insights enabled the rational design of novel compounds anchored in the ebselen core structure, identifying promising candidates with favorable predicted IC50 values. The designed compounds exhibited suitable drug-like characteristics and adopted an active conformation similar to ebselen, inhibiting Mpro function. These findings represent a synergistic approach merging ligand and structure-based drug design; with the potential to guide experimental synthesis and enzyme assay testing.

Free-Energy Profiles for Membrane Permeation of Compounds Calculated Using Rare-Event Sampling Methods

The free-energy profile of a compound is an essential measurement in evaluating the membrane permeation process by means of theoretical methods. Computationally, molecular dynamics (MD) simulation allows the free-energy profile calculation. However, MD simulations frequently fail to sample membrane permeation because they are rare events induced in longer timescales than the accessible timescale of MD, leading to an insufficient conformational search to calculate an incorrect free-energy profile. To achieve a sufficient conformational search, several enhanced sampling methods have been developed and elucidated the membrane permeation process. In addition to these enhanced sampling methods, we proposed a simple yet powerful free-energy calculation of a compound for the membrane permeation process based on originally rare-event sampling methods developed by us. Our methods have a weak dependency on external biases and their optimizations to promote the membrane permeation process. Based on distributed computing, our methods only require the selection of initial structures and their conformational resampling, whereas the enhanced sampling methods may be required to adjust external biases. Furthermore, our methods efficiently search membrane permeation processes with simple scripts without modifying any MD program. As demonstrations, we calculated the free-energy profiles of seven linear compounds for their membrane permeation based on a hybrid conformational search using two rare-event sampling methods, that is, (1) parallel cascade selection MD (PaCS-MD) and (2) outlier flooding method (OFLOOD), combined with a Markov state model (MSM) construction. In the first step, PaCS-MD generated initial membrane permeation paths of a compound. In the second step, OFLOOD expanded the unsearched conformational area around the initial paths, allowing for a broad conformational search. Finally, the trajectories were employed to construct reliable MSMs, enabling correct free-energy profile calculations. Furthermore, we estimated the membrane permeability coefficients of all compounds by constructing the reliable MSMs for their membrane permeation. In conclusion, the calculated coefficients were qualitatively correlated with the experimental measurements (correlation coefficient (R²) = 0.8689), indicating that the hybrid conformational search successfully calculated the free-energy profiles and membrane permeability coefficients of the seven compounds.

Simulation toolkits at the molecular scale for trans-scale thermal signaling

Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.

Ligand Binding Path Sampling Based on Parallel Cascade Selection Molecular Dynamics: LB-PaCS-MD

Parallel cascade selection molecular dynamics (PaCS-MD) is a rare-event sampling method that generates transition pathways between a reactant and product. To sample the transition pathways, PaCS-MD repeats short-time MD simulations from important configurations as conformational resampling cycles. In this study, PaCS-MD was extended to sample ligand binding pathways toward a target protein, which is referred to as LB-PaCS-MD. In a ligand-concentrated environment, where multiple ligand copies are randomly arranged around the target protein, LB-PaCS-MD allows for the frequent sampling of ligand binding pathways. To select the important configurations, we specified the center of mass (COM) distance between each ligand and the relevant binding site of the target protein, where snapshots generated by the short-time MD simulations were ranked by their COM distance values. From each cycle, snapshots with smaller COM distance values were selected as the important configurations to be resampled using the short-time MD simulations. By repeating conformational resampling cycles, the COM distance values gradually decreased and converged to constants, meaning that a set of ligand binding pathways toward the target protein was sampled by LB-PaCS-MD. To demonstrate relative efficiency, LB-PaCS-MD was applied to several proteins, and their ligand binding pathways were sampled more frequently than conventional MD simulations.

Independent Nontargeted Parallel Cascade Selection Molecular Dynamics (Ino-PaCS-MD) to Enhance the Conformational Sampling of Proteins

Biological functions are related to long-time protein dynamics (rare events) that are induced over microseconds. Such protein dynamics can be investigated using molecular dynamics (MD) simulations. However, the detection of rare events remains challenging using conventional MD (cMD) since the accessible timescales of cMD are shorter than those of the biological functions. Recently, the parallel cascade selection MD (PaCS-MD) has been proposed to detect such rare events, wherein transition paths are generated between a given reactant and product. As an extension, the nontargeted PaCS-MD (nt-PaCS-MD) has been proposed to predict the transition paths without requiring reference to any product. Thus, as a further extension, we herein propose independent nt-PaCS-MD, namely, Ino-PaCS-MD, wherein multiple walkers are launched from a set of different starting configurations. Each walker repeats a cycle of restarting short-time MD simulations from configurations with high potentials for making transitions to neighboring metastable states. To further enhance the sampling ability, Ino-PaCS-MD temporarily stops the conformational search and periodically resets the starting configurations so that they are uniformly distributed in a conformational subspace, thereby preventing a given protein from being trapped in one of the metastable states. As a demonstration, Ino-PaCS-MD successfully detects rare events of a maltose-binding protein as open-close transitions with a nanosecond-order simulation time, although a microsecond-order cMD simulation failed to detect these rare events, showing the high sampling efficiency of Ino-PaCS-MD.

Activation Pathways and Free Energy Landscapes of the SARS-CoV-2 Spike Protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses a spike protein (S-protein) to recognize the receptor protein ACE2 of human cells and initiate infection, during which the conformational transition of the S-protein from inactive (down) state to active (up) state is one of the key molecular events determining the infectivity but the underlying mechanism remains poorly understood. In this work, we investigated the activation pathways and free energy landscape of the S-protein of SARS-CoV-2 and compared with those of the closely related counterpart SARS-CoV using molecular dynamics simulations. Our results revealed a large difference between the activation pathways of the two S-proteins. The transition from inactive to an active state for the S-protein of SARS-CoV-2 is more cooperative, involving simultaneous disruptions of several key interfacial hydrogen bonds, and the transition encounters a much higher free energy barrier. In addition, the conformational equilibrium of the SARS-CoV-2 S-protein is more biased to the inactive state compared to that of the SARS-CoV S-protein, suggesting the transient feature of the active state before binding to the receptor protein of the host cell. The key interactions contributing to the difference of the activation pathways and free energy landscapes were discussed. The results provide insights into the molecular mechanism involved in the initial stage of the SARS-CoV-2 infection.

The Effects of Quercetin on the Apoptosis of Human Breast Cancer Cell Lines MCF-7 and MDA-MB-231: A Systematic Review

This systematic review was performed with a focus on the effects of quercetin (QT) on the human breast cancer cell lines MCF-7 and MDA-MB-231. PubMed, Scopus, Science Direct, and Google Scholar databases were searched up to May 2020 using relevant keywords. All articles written in English evaluating the effects of QT on the human breast cancer cell lines MCF-7 and/or MDA-MB-231 were eligible for the review. Totally, 31 articles were included in this review. Out of them, 23 studies investigated the effects of QT on MCF-7 cells and indicated that QT induces apoptosis in the cells. Of 15 studies that examined the effects of QT on MDA-MB-231 cells, 14 reports showed successful apoptosis. It is concluded that QT might be beneficial in the eliminating of breast cancer cells. However, further clinical trials are warranted to further verify these outcomes.

Enhanced Conformational Sampling Method Based on Anomaly Detection Parallel Cascade Selection Molecular Dynamics: ad-PaCS-MD

In the present study, we propose a rare-event sampling method called anomaly detection parallel cascade selection molecular dynamics (ad-PaCS-MD). The original PaCS-MD was designed to generate conformational transition pathways from a given reactant to a product when the latter is known a priori. As an extension of the original method, ad-PaCS-MD has been designed to efficiently search transition pathways from a given reactant without referring to a given product. In ad-PaCS-MD, rarely occurring but essential states (configurations) of proteins for the transitions are identified based on the degrees of an anomaly. In more detail, ad-PaCS-MD adopts an algorithm called an anomaly detection generative adversarial network (anoGAN) as a measure for detecting rarely occurring states to be resampled. Here, the essential configurations with higher degrees of the anomaly are selected with anoGAN and intensively resampled by restarting short-time MD simulations from the selected configurations. By repeating the detections and resampling of configurations with the higher degrees of the anomaly, ad-PaCS-MD automatically and efficiently promotes the rare events and gives a wide range of the free energy landscape by combining with the Markov state model construction. As demonstrations, open-closed transitions of two globular proteins (T4 lysozyme and maltose-binding protein) were promoted with ad-PaCS-MD by referring only to the given starting configurations. In each demonstration, ad-PaCS-MD promoted the large-amplitude open-closed transitions with nanosecond-order simulation times. In conclusion, our demonstrations showed a higher conformational sampling efficiency for ad-PaCS-MD than conventional MD (CMD) because CMD required computational costs of more than microsecond-order simulation times to promote the rare events.

Nontargeted Parallel Cascade Selection Molecular Dynamics Based on a Nonredundant Selection Rule for Initial Structures Enhances Conformational Sampling of Proteins

Nontargeted parallel cascade selection molecular dynamics (nt-PaCS-MD) is a method for enhanced conformational sampling of proteins. To search a broad conformational subspace, nt-PaCS-MD repeats cycles of conformational resampling from relevant initial structures. Generally, the conformational sampling efficiency of nt-PaCS-MD depends on a selection rule for the initial structures. In the original nt-PaCS-MD, the initial structures were selected by referring to structural distributions of protein configurations generated by conformational resampling (multiple short-time MD simulations). However, their structural redundancy among the initial structures was neglected for the cycles of conformational resampling, indicating that similar protein configurations might be frequently specified and resampled in every cycle in the original nt-PaCS-MD. To reduce the possibility of resampling from redundant initial structures, we propose an alternative selection rule that accounts for structural similarity among the initial structures. Specifically, a pairwise root-mean-square deviation (RMSD) is defined for all of the initial structures selected for all of the past cycles. Then a set of protein configurations with a larger pairwise RMSD is sequentially specified and resampled in the next cycle, which is regarded to as a history-dependent selection of initial structures by considering a profile of the past specified initial structures. The present scheme, termed extended nt-PaCS-MD, prevents us from resampling a set of redundant protein configurations. To check the conformational sampling efficiency of the extended nt-PaCS-MD, we used a middle-sized protein, T4 lysozyme, in explicit water. Through the assessment, this extended nt-PaCS-MD identified the open-closed transitions of T4 lysozyme more efficiently than the original nt-PaCS-MD.

Selection Rules for Outliers in Outlier Flooding Method Regulate Its Conformational Sampling Efficiency

The outlier flooding method (OFLOOD) has been proposed as an enhanced conformational sampling method of proteins. In OFLOOD, rarely occurring states of proteins are detected as sparse conformational distributions (outliers) with a clustering algorithm. The detected outliers are intensively resampled with short-time molecular dynamics (MD) simulations. As a set of cycles, OFLOOD repeats selections of outliers and their conformational resampling. Herein, as an essential issue to be tackled to perform OFLOOD efficiently, a selection rule for outliers should be carefully specified. Generally, many outliers are detected from distributions on conformational subspaces with the clustering. Judging from its computational costs, it is unreasonable to select all the detected outliers upon the conformational resampling. Therefore, it is important to consider which outliers should be selected from the sparse distributions when restarting their short-time MD simulations with limited computational costs. In this sense, we investigated the conformational sampling efficiency of OFLOOD by changing the selection rules for outliers. To address the conformational sampling efficiency of OFLOOD depending on its selection rules, outliers to be resampled were selected by focusing their probability occurrences (populations of outliers). As a comparison, a random selection rule for outliers was also considered. Through the present assessment, the random selection of outliers showed the most efficient conformational sampling efficiency compared to the other OFLOOD trials using the biased selection rules, indicating that a variety of outliers should be selected and resampled during the OFLOOD cycles. In conclusion, the random outlier selection rule is the best strategy to perform OFLOOD efficiently.