An explicit-solvent conformation search method using open software

Widely-targeted in silico and in vitro evaluation of veratrum alkaloid analogs as FAK inhibitors and dual targeting of FAK and Hh/SMO pathways for cancer therapy: A critical analysis

Focal Adhesive Kinase (FAK), a key player in aggressive cancers, mediates signals crucial for progression, invasion, and metastasis. Despite advances in targeted therapies, drug resistance is still a challenge, and survival rates remain low, particularly for late-stage patients, emphasizing the need for innovative cancer therapeutics. Cyclopamine, a veratrum alkaloid, has shown promising anti-tumor properties, but the search for more potent analogs with enhanced affinity for the biological target continues. This study employs a hybrid virtual screening approach combining pharmacophore model-based virtual screening (PB-VS) and docking-based virtual screening (DB-VS) to identify potential inhibitors of the FAK catalytic domain. PB-VS on the PubChem database yielded a set of hits, which were then docked with the FAK catalytic domain in two stages (1st and 2nd DB-VS). Hits were ranked based on docking scores and interactions with the active site. The top three compounds underwent molecular dynamics simulations, alongside two control compounds (SMO inhibitor(s) and FAK inhibitor(s)), to assess stability through RMSD, RMSF, Rg, and SASA analyses. ADMET properties were evaluated, and compounds were filtered based on drug-likeness criteria. Molecular dynamics simulations demonstrated the stability of compounds when complexed with the FAK catalytic domain. Compounds 16 (-25 kcal/mol), 88 (-27.47 kcal/mol), and 87 (-18.94 kcal/mol) exhibited comparable docking scores, interaction profiles, stability, and binding energies, indicating their potential as lead candidates. However, further validation and optimization through quantitative structure-activity relationship (QSAR) studies are essential to refine their efficacy and therapeutic potential. The in vitro cell-based assay demonstrated that compound 101PF, a FAK inhibitor, significantly inhibited the proliferation and migration of A549 cells. However, the results regarding the combined effects of FAK and SMO inhibitors were inconclusive, highlighting the need for further investigation. This study contributes to developing more effective anti-cancer drugs by improving the understanding of potential cyclopamine-based veratrum alkaloid analogs with enhanced interactions with the FAK catalytic domain.

PyConSolv: A Python Package for Conformer Generation of (Metal-Containing) Systems in Explicit Solvent

We introduce PyConSolv, a freely available Python package that automates the generation of conformers of metal- and nonmetal-containing complexes in explicit solvent, through classical molecular dynamics simulations. Using a streamlined workflow and interfacing with widely used computational chemistry software, PyConSolv is an all-in-one tool for the generation of conformers in any solvent. Input requirements are minimal; only the geometry of the structure and the desired solvent in xyz (XMOL) format are needed. The package can also account for charged systems, by including arbitrary counterions in the simulation. A bonded model parametrization is performed automatically, utilizing AmberTools, ORCA, and Multiwfn software packages. PyConSolv provides a selection of preparametrized solvents and counterions for use in classical molecular dynamics simulations. We show the applicability of our package on a number of (transition-metal-containing) systems. The software is provided open source and free of charge.

Stereochemistry of Simple Molecules inside Nanotubes and Fullerenes: Unusual Behavior of Usual Systems

Over the past three decades, carbon nanotubes and fullerenes have become remarkable objects for starting the implementation of new models and technologies in different branches of science. To a great extent, this is defined by the unique electronic and spatial properties of nanocavities due to the ramified π-electron systems. This provides an opportunity for the formation of endohedral complexes containing non-covalently bonded atoms or molecules inside fullerenes and nanotubes. The guest species are exposed to the force field of the nanocavity, which can be described as a combination of electronic and steric requirements. Its action significantly changes conformational properties of even relatively simple molecules, including ethane and its analogs, as well as compounds with C-O, C-S, B-B, B-O, B-N, N-N, Al-Al, Si-Si and Ge-Ge bonds. Besides that, the cavity of the host molecule dramatically alters the stereochemical characteristics of cyclic and heterocyclic systems, affects the energy of pyramidal nitrogen inversion in amines, changes the relative stability of cis and trans isomers and, in the case of chiral nanotubes, strongly influences the properties of R- and S- enantiomers. The present review aims at primary compilation of such unusual stereochemical effects and initial evaluation of the nature of the force field inside nanotubes and fullerenes.

Quality Threshold Clustering of Molecular Dynamics: A Word of Caution

Clustering Molecular Dynamics trajectories is a common analysis that allows grouping together similar conformations. Several algorithms have been designed and optimized to perform this routine task, and among them, Quality Threshold stands as a very attractive option. This algorithm guarantees that in retrieved clusters no pair of frames will have a similarity value greater than a specified threshold, and hence, a set of strongly correlated frames are obtained for each cluster. In this work, it is shown that various commonly used software implementations are flawed by confusing Quality Threshold with another simplistic well-known clustering algorithm published by Daura et al. (Daura, X.; van Gunsteren, W. F.; Jaun, B.; Mark, A. E.; Gademann, K.; Seebach, D. Peptide Folding: When Simulation Meets Experiment. Angew. Chemie Int. Ed. 1999, 38 (1/2), 236-240). Daura's algorithm does not impose any quality threshold for the frames contained in retrieved clusters, bringing unrelated structural configurations altogether. The advantages of using Quality Threshold whenever possible to explore Molecular Dynamic trajectories is exemplified. An in-house implementation of the original Quality Threshold algorithm has been developed in order to illustrate our comments, and its code is freely available for further use by the scientific community.

Modeling covalent-modifier drugs

In this review, we present a summary of how computer modeling has been used in the development of covalent-modifier drugs. Covalent-modifier drugs bind by forming a chemical bond with their target. This covalent binding can improve the selectivity of the drug for a target with complementary reactivity and result in increased binding affinities due to the strength of the covalent bond formed. In some cases, this results in irreversible inhibition of the target, but some targeted covalent inhibitor (TCI) drugs bind covalently but reversibly. Computer modeling is widely used in drug discovery, but different computational methods must be used to model covalent modifiers because of the chemical bonds formed. Structural and bioinformatic analysis has identified sites of modification that could yield selectivity for a chosen target. Docking methods, which are used to rank binding poses of large sets of inhibitors, have been augmented to support the formation of protein-ligand bonds and are now capable of predicting the binding pose of covalent modifiers accurately. The pK_a's of amino acids can be calculated in order to assess their reactivity towards electrophiles. QM/MM methods have been used to model the reaction mechanisms of covalent modification. The continued development of these tools will allow computation to aid in the development of new covalent-modifier drugs. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.