Leveraging the Fragment Molecular Orbital and MM-GBSA Methods in Virtual Screening for the Discovery of Novel Non-Covalent Inhibitors Targeting the TEAD Lipid Binding Pocket.

The Hippo pathway controls organ size and homeostasis and is linked to numerous diseases, including cancer. The transcriptional enhanced associate domain (TEAD) family of transcription factors acts as a receptor for downstream effectors, namely yes-associated protein (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ), which binds to various transcription factors and is essential for stimulated gene transcription. YAP/TAZ-TEAD facilitates the upregulation of multiple genes involved in evolutionary cell proliferation and survival. TEAD1-4 overexpression has been observed in different cancers in various tissues, making TEAD an attractive target for drug development. The central drug-accessible pocket of TEAD is crucial because it undergoes a post-translational modification called auto-palmitoylation. Crystal structures of the C-terminal TEAD complex with small molecules are available in the Protein Data Bank, aiding structure-based drug design. In this study, we utilized the fragment molecular orbital (FMO) method, molecular dynamics (MD) simulations, shape-based screening, and molecular mechanics-generalized Born surface area (MM-GBSA) calculations for virtual screening, and we identified a novel non-covalent inhibitor-BC-001-with IC50 = 3.7 µM in a reporter assay. Subsequently, we optimized several analogs of BC-001 and found that the optimized compound BC-011 exhibited an IC50 of 72.43 nM. These findings can be used to design effective TEAD modulators with anticancer therapeutic implications.

Exploring the Conformational Ensembles of Protein-Protein Complex with Transformer-Based Generative Model.

Protein-protein interactions are the basis of many protein functions, and understanding the contact and conformational changes of protein-protein interactions is crucial for linking the protein structure to biological function. Although difficult to detect experimentally, molecular dynamics (MD) simulations are widely used to study the conformational ensembles and dynamics of protein-protein complexes, but there are significant limitations in sampling efficiency and computational costs. In this study, a generative neural network was trained on protein-protein complex conformations obtained from molecular simulations to directly generate novel conformations with physical realism. We demonstrated the use of a deep learning model based on the transformer architecture to explore the conformational ensembles of protein-protein complexes through MD simulations. The results showed that the learned latent space can be used to generate unsampled conformations of protein-protein complexes for obtaining new conformations complementing pre-existing ones, which can be used as an exploratory tool for the analysis and enhancement of molecular simulations of protein-protein complexes.

Transformer-Based Molecular Generative Model for Antiviral Drug Design.

Since the Simplified Molecular Input Line Entry System (SMILES) is oriented to the atomic-level representation of molecules and is not friendly in terms of human readability and editable, however, IUPAC is the closest to natural language and is very friendly in terms of human-oriented readability and performing molecular editing, we can manipulate IUPAC to generate corresponding new molecules and produce programming-friendly molecular forms of SMILES. In addition, antiviral drug design, especially analogue-based drug design, is also more appropriate to edit and design directly from the functional group level of IUPAC than from the atomic level of SMILES, since designing analogues involves altering the R group only, which is closer to the knowledge-based molecular design of a chemist. Herein, we present a novel data-driven self-supervised pretraining generative model called "TransAntivirus" to make select-and-replace edits and convert organic molecules into the desired properties for design of antiviral candidate analogues. The results indicated that TransAntivirus is significantly superior to the control models in terms of novelty, validity, uniqueness, and diversity. TransAntivirus showed excellent performance in the design and optimization of nucleoside and non-nucleoside analogues by chemical space analysis and property prediction analysis. Furthermore, to validate the applicability of TransAntivirus in the design of antiviral drugs, we conducted two case studies on the design of nucleoside analogues and non-nucleoside analogues and screened four candidate lead compounds against anticoronavirus disease (COVID-19). Finally, we recommend this framework for accelerating antiviral drug discovery.

Application of a deep generative model produces novel and diverse functional peptides against microbial resistance.

Antimicrobial resistance could threaten millions of lives in the immediate future. Antimicrobial peptides (AMPs) are an alternative to conventional antibiotics practice against infectious diseases. Despite the potential contribution of AMPs to the antibiotic's world, their development and optimization have encountered serious challenges. Cutting-edge methods with novel and improved selectivity toward resistant targets must be established to create AMPs-driven treatments. Here, we present AMPTrans-lstm, a deep generative network-based approach for the rational design of AMPs. The AMPTrans-lstm pipeline involves pre-training, transfer learning, and module identification. The AMPTrans-lstm model has two sub-models, namely, (long short-term memory) LSTM sampler and Transformer converter, which can be connected in series to make full use of the stability of LSTM and the novelty of Transformer model. These elements could generate AMPs candidates, which can then be tailored for specific applications. By analyzing the generated sequence and trained AMPs, we prove that AMPTrans-lstm can expand the design space of the trained AMPs and produce reasonable and brand-new AMPs sequences. AMPTrans-lstm can generate functional peptides for antimicrobial resistance with good novelty and diversity, so it is an efficient AMPs design tool.

Explore drug-like space with deep generative models.

The process of design/discovery of drugs involves the identification and design of novel molecules that have the desired properties and bind well to a given disease-relevant target. One of the main challenges to effectively identify potential drug candidates is to explore the vast drug-like chemical space to find novel chemical structures with desired physicochemical properties and biological characteristics. Moreover, the chemical space of currently available molecular libraries is only a small fraction of the total possible drug-like chemical space. Deep molecular generative models have received much attention and provide an alternative approach to the design and discovery of molecules. To efficiently explore the drug-like space, we first constructed the drug-like dataset and then performed the generative design of drug-like molecules using a Conditional Randomized Transformer approach with the molecular access system (MACCS) fingerprint as a condition and compared it with previously published molecular generative models. The results show that the deep molecular generative model explores the wider drug-like chemical space. The generated drug-like molecules share the chemical space with known drugs, and the drug-like space captured by the combination of quantitative estimation of drug-likeness (QED) and quantitative estimate of protein-protein interaction targeting drug-likeness (QEPPI) can cover a larger drug-like space. Finally, we show the potential application of the model in design of inhibitors of MDM2-p53 protein-protein interaction. Our results demonstrate the potential application of deep molecular generative models for guided exploration in drug-like chemical space and molecular design.

Development of an innovative data-driven system to generate descriptive prediction equation of dielectric constant on small sample sets.

Dielectric constant (DC, ε) is a fundamental parameter in material sciences to measure polarizability of the system. In industrial processes, its value is an imperative indicator, which demonstrates the dielectric property of material and compiles information including separation information, chemical equilibrium, chemical reactivity analysis, and solubility modeling. Since, the available ε-prediction models are fairly primitive and frequently suffer from serious failures especially when deals with strong polar compounds. Therefore, we have developed a novel data-driven system to improve the efficiency and wide-range applicability of ε using in material sciences. This innovative scheme adopts the correlation distance and genetic algorithm to discriminate features' combination and avoid overfitting. Herein, the prediction output of the single ML model as a coding to estimate the target value by simulating the layer-by-layer extraction in deep learning, and enabling instant search for the optimal combination of features is recruited. Our model established an improved correlation value of 0.956 with target as compared to the previously available best traditional ML result of 0.877. Our framework established a profound improvement, especially for material systems possessing ε value >50. In terms of interpretability, we have derived a conceptual computational equation from a minimum generating tree. Our innovative data-driven system is preferentially superior over other methods due to its application for the prediction of dielectric constants as well as for the prediction of overall micro and macro-properties of any multi-components complex.

De novo molecular design with deep molecular generative models for PPI inhibitors.

We construct a protein-protein interaction (PPI) targeted drug-likeness dataset and propose a deep molecular generative framework to generate novel drug-likeness molecules from the features of the seed compounds. This framework gains inspiration from published molecular generative models, uses the key features associated with PPI inhibitors as input and develops deep molecular generative models for de novo molecular design of PPI inhibitors. For the first time, quantitative estimation index for compounds targeting PPI was applied to the evaluation of the molecular generation model for de novo design of PPI-targeted compounds. Our results estimated that the generated molecules had better PPI-targeted drug-likeness and drug-likeness. Additionally, our model also exhibits comparable performance to other several state-of-the-art molecule generation models. The generated molecules share chemical space with iPPI-DB inhibitors as demonstrated by chemical space analysis. The peptide characterization-oriented design of PPI inhibitors and the ligand-based design of PPI inhibitors are explored. Finally, we recommend that this framework will be an important step forward for the de novo design of PPI-targeted therapeutics.

Comprehensive strategies of machine-learning-based quantitative structure-activity relationship models.

Early quantitative structure-activity relationship (QSAR) technologies have unsatisfactory versatility and accuracy in fields such as drug discovery because they are based on traditional machine learning and interpretive expert features. The development of Big Data and deep learning technologies significantly improve the processing of unstructured data and unleash the great potential of QSAR. Here we discuss the integration of wet experiments (which provide experimental data and reliable verification), molecular dynamics simulation (which provides mechanistic interpretation at the atomic/molecular levels), and machine learning (including deep learning) techniques to improve QSAR models. We first review the history of traditional QSAR and point out its problems. We then propose a better QSAR model characterized by a new iterative framework to integrate machine learning with disparate data input. Finally, we discuss the application of QSAR and machine learning to many practical research fields, including drug development and clinical trials.

A new class of 1,3,5-triazine-based selective estrogen receptor degraders (SERDs): Lead optimization, molecular docking and dynamic simulation.

Selective estrogen receptor degrader (SERD) that acts as not only ER antagonist, but also ER degrader, would be useful for the treatment for drug-resistance ER+ breast cancer. However, most of currently available SERD candidates involve very limited molecular scaffolds and are still in clinical trials. In this study, we introduced a 1,3,5-triazine ring into a homobibenzyl motif extracted from amounts of ER ligands and synthesized sixteen SERDs bearing acrylic acid or acrylic amide side chains that possess both ERα antagonism and degradation properties. And all compounds were screened for their anti-proliferative activity against ER+ MCF-7 and Ishikawa cell lines. Among them, compound XHA1614 displayed potent growth inhibition activity against MCF-7 and Ishikawa cells with IC50 values of 3.15 µM and 3.11 µM, respectively. Moreover, XHA1614 could dramatically degrade ER level at 1 nM in a Western blotting assay and afforded an outstanding antagonistic activity via suppressing the expression of progesterone receptor messenger RNA in MCF-7 cells in a RT-PCR assay. Further molecular docking and dynamic simulation on properly selected derivative furnished insights into its binding profile within ERα. Our findings suggest that the 1,3,5-triazine core was a feasible alternative to currently reported SERD scaffold, and provide information that will be useful for further development of promising SERDs candidates for breast cancer therapies.

Seeking mTORC1 Inhibitors Through Molecular Dynamics Simulation of Arginine Analogs Inhibiting CASTOR1.

BACKGROUND: Hyperactivity of the mechanistic target of rapamycin complex 1 (mTORC1) is implicated in a variety of diseases such as cancer and diabetes. Treatment may benefit from effective mTORC1 inhibition, which can be achieved by preventing arginine from disrupting the cytosolic arginine sensor for mTORC1 subunit 1 (CASTOR1)-GTPase-activating proteins toward RAGS subcomplex 2 (GATOR2) complex through binding with CASTOR1. An attractive idea is to determine analogues of arginine that are as competent as arginine in binding with CASTOR1, but without disrupting the CASTOR1-GATOR2 interaction. MATERIALS AND METHODS: Molecular dynamics simulations were performed for binding of arginine analogues with CASTOR1 and binding free energy, hydrogen bond formation, and root mean squared deviation and root mean square fluctuation kinetics were then calculated. RESULTS: The binding free energy calculations revealed that Nα-acetyl-arginine, citrulline, and norarginine have sufficient binding affinity with CASTOR1 to compete with arginine. The hydrogen bond analysis revealed that norarginine, Nα-acetyl-arginine and D-arginine have proficient H-bonds that can facilitate their entering the narrow binding pocket. CONCLUSION: Norarginine and Nα-acetyl-arginine are the top drug candidates for mTORC1 inhibition, with Nα-acetyl-arginine being the best choice.