Journal of Computer-Aided Molecular Design最新文献

This study aimed to reveal interactions of the stress response sigma subunits (factors) σ^D and σ^H of RNA polymerase and promoters in Gram-positive bacterium Corynebacterium glutamicum by combining wet-lab obtained data and in silico modeling. Computer modeling-guided point mutagenesis of C. glutamicum σ^H subunit led to the creation of a panel of σ^H variants. Their ability to initiate transcription from naturally occurring hybrid σ^D/σ^H-dependent promoter Pcg0441 and two control canonical promoters (σ^D-dependent PrsdA and σ^H-dependent PuvrD3) was measured and interpreted using molecular dynamics simulations of homology models of all complexes. The results led us to design the artificial hybrid promoter PD₃₅H₁₀ combining the −10 element of the PuvrD3 promoter and the −35 element of the PrsdA promoter. This artificial hybrid promoter PD_35-rsdAH_10-uvrD3 showed almost optimal properties needed for the bio-orthogonal transcription (not interfering with the native biological processes).

The influence of polar, water-miscible organic solvents (POS) on protein structure, stability, and functional activity is a subject of significant interest and complexity. This study examines the effects of acetonitrile (ACN), a semipolar, aprotic solvent, on the solvation properties of blocked Ace-Gly-X-Gly-Nme tripeptides (where Ace and Nme stands for acetyl and N-methyl amide groups respectively and X is any amino acid) through extensive molecular dynamics simulations. Individual simulations were conducted for each peptide, encompassing five different ACN concentrations within the range of χ_ACN = 0.1–0.9. The preferential solvation parameter (Γ) calculated using the Kirkwood-Buff integral method was used for the assessment of peptide interactions with water/ACN. Additionally, weighted Voronoi tessellation was applied to obtain a three-way data set containing four time-averaged contact surface area types between peptide atoms and water/ACN atoms. A mathematical technique known as N-way Partial Least Squares (NPLS) was utilized to anticipate the preferential interactions between peptides and water/ACN from the contact surface areas. Furthermore, the temperature dependency of peptide-solvent interactions was investigated using a subset of 10 amino acids representing a range of hydrophobicities. MD simulations were conducted at five temperatures, spanning from 283 to 343 K, with subsequent analysis of data focusing on both preferential solvation and peptide-solvent contact surface areas. The results demonstrate the efficacy of utilizing contact surface areas between the peptide and solvent constituents for successfully predicting preferential interactions in water/ACN mixtures across various ACN concentrations and temperatures.

PIK3CA gene encoding PI3K p110α is one of the most frequently mutated and overexpressed in majority of human cancers. Development of potent and selective novel inhibitors targeting PI3Kα was considered as the most promising approaches for cancer treatment. In this investigation, a virtual screening platform for PI3Kα inhibitors was established by employing machine learning methods, pharmacophore modeling, and molecular docking approaches. 28 potential PI3Kα inhibitors with different scaffolds were selected from the databases with 295,024 compounds. Among the 28 hits, hit15 exhibited the best inhibitory effect against PI3Kα with IC₅₀ value less than 1.0 µM. The molecular dynamics simulation indicated that hit15 could stably bind to the active site of PI3Kα, interact with some residues by hydrophobic, electrostatic and hydrogen bonding interactions, and finally induced PI3Kα active pocket substantial conformation changes. Stable H-bond interactions were formed between hit15 and residues of Lys776, Asp810 and Asp933. The binding free energy of PI3Kα-hit15 was − 65.3 kJ/mol. The free energy decomposition indicated that key residues of Asp805, Ile848 and Ile932 contributed stronger energies to the binding free energy. The above results indicated that hit15 with novel scaffold was a potent PI3Kα inhibitor and considered as a promising candidate for further drug development to treat various cancers with PI3Kα over activated.

Graphical Abstract

Relative solubilities, i.e. whether a given molecule is more soluble in one solvent compared to others, is a critical parameter for pharmaceutical and agricultural formulation development and chemical synthesis, material science, and environmental chemistry. In silico predictions of this crucial variable can help reducing experiments, waste of solvents and synthesis optimization. In this study, we evaluate the performance of different physics-based methods for predicting relative solubilities. Our assessment involves quantum mechanics-based COSMO-RS and molecular dynamics-based free energy methods using OPLS4, the open-source OpenFF Sage, and GAFF force fields, spanning over 200 solvent–solute combinations. Our investigation highlights the important role of compound multimerization, an effect which must be accounted for to obtain accurate relative solubility predictions. The performance landscape of these methods is varied, with significant differences in precision depending on both the method used and the solute considered, thereby offering an improved understanding of the predictive power of physics-based methods in chemical research.

G-quadruplexes refer to a large group of nucleic acid–based structures. In recent years, they have been attracting attention due to their biological roles in the telomeres and promoter regions. These structures show wide diversity in topology, however, development of methods for structural classification of G-quadruplexes has been evaded for a long time. There has been a limited number of studies aiming to bring forth a secondary structure classification method. The situation was even more complex than imagined, since the discovery of bulged and mismatched G-quadruplexes while most of the available tools fail to distinguish these non-canonical G-quadruplex motifs. Moreover, the interpretation of their analysis output still requires expert knowledge. In this study, we propose a new method for identification of unimolecular G-Quadruplexes and classification by secondary structures based on three-dimensional structural data. Briefly, coordinates of guanines are processed to identify tetrads, loops and bulges. Then, we present the secondary structure in the form of a depiction which shows the loop types, bulges, and guanines that participate in each tetrad. Moreover, CIIS-GQ identifies non-guanine nucleotides that joins the G-tetrads and forms multiplets. Finally, the results of our study are compared with DSSR and ElTetrado classification methods, and the advantages of the proposed depiction method for representing secondary structures were discussed. The source code of the method can be accessed via https://github.com/TugayDirek/CIIS-GQ.

Therapeutic monoclonal antibodies are an effective method of treating acute infectious diseases. However, knowing which of the produced antibodies in the vast number of human antibodies can cure the disease requires a long time and advanced technology. The previously introduced iBRAB method relies on studied antibodies to design a broad-spectrum antibody capable of neutralizing antigens of many different Influenza A viral strains. To evaluate the antigen-binding fragment as an applicable drug, the therapeutic antibody profiles providing guidelines collected from clinically staged therapeutic antibodies were used to access different measurements. Although the evaluated values were within an accepted range, the modification in the amino acid sequence is required for better properties. Thus, using the steered molecular dynamics (SMD) simulation to determine the binding capacity of amino acids in the functional region, the profile of interacted amino acids of Fab with the antigen was established for modified reference. As a result, the model was modified with amino acids elimination at positions 96–97 in the heavy chain and 26–27, 91, 96–97, and 102–103 in the light chain, which has better Therapeutic Antibody Profiler evaluations than the original designation. Thus again, SMD simulation is a promising computational approach for post-modification in rational drug design.

So-called “cross-docking” is the prediction of the bound configuration of small-molecule ligands that differ from the cognate ligand of a protein co-crystal structure. This is a much more challenging problem than re-docking the cognate ligand, particularly when the new ligand is structurally dissimilar from prior known ones. We have updated the previously introduced PINC (“PINC Is Not Cognate”) benchmark which introduced the idea of temporal segregation to measure cross-docking performance. The temporal set encompasses 846 future ligands for ten targets based on information from the earliest 25% of X-ray co-crystal structures known for each target. Here, we extend the benchmark to include thirteen targets where the bound poses of 128 macrocyclic ligands are to be predicted based on knowledge from structures of bound non-macrocyclic ligands. Performance was roughly equivalent for both the temporally-split non-macrocyclic ligand set and the macrocycle prediction set. Using standard and fully automatic protocols for the Surflex-Dock and ForceGen methods, across the combined 974 non-macrocyclic and macrocyclic ligands, the top-scoring pose family was correct 68% of the time, with the top-two pose families achieving a 79% success rate. Correct poses among all those predicted were identified 92% of the time. These success rates far exceeded those observed for the alternative methods AutoDock Vina and Gnina on both sets.

Over the last decade, automatic chemical design frameworks for discovering molecules with drug-like properties have significantly progressed. Among them, the variational autoencoder (VAE) is a cutting-edge approach that models the tractable latent space of the molecular space. In particular, the usage of a VAE along with a property estimator has attracted considerable interest because it enables gradient-based optimization of a given molecule. However, although successful results have been achieved experimentally, the theoretical background and prerequisites for the correct operation of this method have not yet been clarified. In view of the above, we theoretically analyze and rigorously reconstruct the entire framework. From the perspective of parameterized distribution and the information theory, we first describe how the previous model overcomes the limitations of the beta VAE in discovering molecules with the desired properties. Furthermore, we describe the prerequisites for training the above model. Next, from the log-likelihood perspective of each term, we reformulate the objectives for exploring latent space to generate drug-like molecules. The distributional constraints are defined in this study, which will break away from the invalid molecular search. We demonstrated that our model could discover a novel chemical compound for targeting BCL-2 family proteins in de novo approach. Through the theoretical analysis and practical implementation, the importance of the aforementioned prerequisites and constraints to operate the model was verified.

Human Hippo signaling pathway is an evolutionarily conserved regulator network that controls organ development and has been implicated in various cancers. Transcriptional enhanced associate domain-4 (TEAD4) is the final nuclear effector of Hippo pathway, which is activated by Yes-associated protein (YAP) through binding to two separated YAP regions of α1-helix and Ω-loop. Previous efforts have all been addressed on deriving peptide inhibitors from the YAP to target TEAD4. Instead, we herein attempted to rationally design a so-called ‘YAP helix^α1-trap’ based on the TEAD4 to target YAP by using dynamics simulation and energetics analysis as well as experimental assays at molecular and cellular levels. The trap represents a native double-stranded helical hairpin covering a specific YAP-binding site on TEAD4 surface, which is expected to form a three-helix bundle with the α1-helical region of YAP, thus competitively disrupting TEAD4–YAP interaction. The hairpin was further stapled by a disulfide bridge across its two helical arms. Circular dichroism characterized that the stapling can effectively constrain the trap into a native-like structured conformation in free state, thus largely minimizing the entropy penalty upon its binding to YAP. Affinity assays revealed that the stapling can considerably improve the trap binding potency to YAP α1-helix by up to 8.5-fold at molecular level, which also exhibited a good tumor-suppressing effect at cellular level if fused with TAT cell permeation sequence. In this respect, it is considered that the YAP helix^α1-trap-mediated blockade of Hippo pathway may be a new and promising therapeutic strategy against cancers.