Dataset of AMBER force field parameters of drugs, natural products and steroids for simulations using GROMACS

Identification of FTY720 and COH29 as novel topoisomerase I catalytic inhibitors by experimental and computational studies

The development of novel topoisomerase I (TOP1) inhibitors is crucial for overcoming the drawbacks and limitations of current TOP1 poisons. Here, we identified two potential TOP1 inhibitors, namely, FTY720 (a sphingosine 1-phosphate antagonist) and COH29 (a ribonucleotide reductase inhibitor), through experimental screening of known active compounds. Biological experiments verified that FTY720 and COH29 were nonintercalative TOP1 catalytic inhibitors that did not induce the formation of DNA-TOP1 covalent complexes. Molecular docking revealed that FTY720 and COH29 interacted favorably with TOP1. Molecular dynamics simulations revealed that FTY720 and COH29 could affect the catalytic domain of TOP1, thus resulting in altered DNA-binding cavity size. The alanine scanning and interaction entropy identified Arg536 as a hotspot residue. In addition, the bioinformatics analysis predicted that FTY720 and COH29 could be effective in treating malignant breast tumors. Biological experiments verified their antitumor activities using MCF-7 breast cancer cells. Their combinatory effects with TOP1 poisons were also investigated. Further, FTY720 and COH29 were found to cause less DNA damage compared with TOP1 poisons. The findings provide reliable lead compounds for the development of novel TOP1 catalytic inhibitors and offer new insights into the potential clinical applications of FTY720 and COH29 in targeting TOP1.

Differential molecular mechanisms of substrate recognition by selenium methyltransferases, INMT and TPMT, in selenium detoxification and excretion

It is known that the recommended dietary allowance of selenium (Se) is dangerously close to its tolerable upper intake level. Se is detoxified and excreted in urine as trimethylselenonium ion (TMSe) when the amount ingested exceeds the nutritional level. Recently, we demonstrated that the production of TMSe requires two methyltransferases: thiopurine S-methyltransferase (TPMT) and indolethylamine N-methyltransferase (INMT). In this study, we investigated the substrate recognition mechanisms of INMT and TPMT in the Se-methylation reaction. Examination of the Se-methyltransferase activities of two paralogs of INMT, namely, nicotinamide N-methyltransferase and phenylethanolamine N-methyltransferase, revealed that only INMT exhibited Se-methyltransferase activity. Consistently, molecular dynamics simulations demonstrated that dimethylselenide was preferentially associated with the active center of INMT. Using the fragment molecular orbital method, we identified hydrophobic residues involved in the binding of dimethylselenide to the active center of INMT. The INMT-L164R mutation resulted in a deficiency in Se- and N-methyltransferase activities. Similarly, TPMT-R152, which occupies the same position as INMT-L164, played a crucial role in the Se-methyltransferase activity of TPMT. Our findings suggest that TPMT recognizes negatively charged substrates, whereas INMT recognizes electrically neutral substrates in the hydrophobic active center embedded within the protein. These observations explain the sequential requirement of the two methyltransferases in producing TMSe.

Molecular Dynamics Investigations of Human DNA-Topoisomerase I Interacting with Novel Dewar Valence Photo-Adducts: Insights into Inhibitory Activity

Chronic exposure to ultraviolet (UV) radiation is known to induce the formation of DNA photo-adducts, including cyclobutane pyrimidine dimers (CPDs) and Dewar valence derivatives (DVs). While CPDs usually occur at higher frequency than DVs, recent studies have shown that the latter display superior selectivity and significant stability in interaction with the human DNA/topoisomerase 1 complex (TOP1). With the aim to deeply investigate the mechanism of interaction of DVs with TOP1, we report here four all-atom molecular dynamic simulations spanning one microsecond. These simulations are focused on the stability and conformational changes of two DNA/TOP1-DV complexes in solution, the data being compared with the biomimetic thymine dimer counterparts. Results from root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analyses unequivocally confirmed increased stability of the DNA/TOP1-DV complexes throughout the simulation duration. Detailed interaction analyses, uncovering the presence of salt bridges, hydrogen bonds, water-mediated interactions, and hydrophobic interactions, as well as pinpointing the non-covalent interactions within the complexes, enabled the identification of specific TOP1 residues involved in the interactions over time and suggested a potential TOP1 inhibition mechanism in action.

Bioinformatic investigation of Nipah virus surface protein mutations: Molecular docking with Ephrin B2 receptor, molecular dynamics simulation, and structural impact analysis

The SARS-CoV-2 outbreak resulted in significant challenges and loss of life. The Nipah virus, known for its high infectivity and severity, was designated an emergency concern by the World Health Organization. To understand its mutations, the Nipah virus proteins were analyzed extensively, with a focus on the essential G and F proteins responsible for viral entry into host cells. Our bioinformatics analysis unveiled multiple mutations, including simultaneous mutations within a single sequence. Notably, the G273S mutation in the F protein was identified as a potential cause of structural damage, which carries significant implications for vaccine development. Comparing the docking scores of G and F proteins with the Ephrin B2 receptor, it was found that the Y228H mutation in the G protein and the D252G mutation in the F protein likely affect virus entry into host cells. Moreover, our investigation into stability and deformability highlighted the impact of the Y228H mutation in the G protein complex. Molecular dynamics simulations revealed increased flexibility and conformational changes in the G protein complex with the Y228H mutation compared with the known complex. Furthermore, evaluating the root mean square deviation variation demonstrated greater dynamic behavior in the G protein complex and the Ephrin B2 receptor complex. This comprehensive study provides valuable insights into Nipah virus mutations, their significance for vaccine development, and the importance of understanding protein complex behavior in drug discovery. The identified mutations, especially G273S and Y228H, hold crucial implications for future research and potential interventions against the Nipah virus.

A Theobromine Derivative with Anticancer Properties Targeting VEGFR-2: Semisynthesis, in silico and in vitro Studies

A computer-assisted drug design (CADD) approach was utilized to design a new acetamido-N-(para-fluorophenyl)benzamide) derivative of the naturally occurring alkaloid, theobromine, (T-1-APFPB), following the pharmacophoric features of VEGFR-2 inhibitors. The stability and reactivity of T-1-AFPB were assessed through density functional theory (DFT) calculations. Molecular docking assessments showed T-1-AFPB's potential to bind with and inhibit VEGFR-2. The precise binding of T-1-AFPB against VEGFR-2 with optimal energy was further confirmed through several molecular dynamics (MD) simulations, PLIP, MM-GBSA, and PCA studies. Then, T-1-AFPB (4-(2-(3,7-Dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-1-yl)acetamido)-N-(4-fluorophenyl)benzamide) was semi-synthesized and the in vitro assays showed its potential to inhibit VEGFR-2 with an IC50 value of 69 nM (sorafenib's IC50 was 56 nM) and to inhibit the growth of HepG2 and MCF-7 cancer cell lines with IC50 values of 2.24±0.02 and 3.26±0.02 μM, respectively. Moreover, T-1-AFPB displayed very high selectivity indices against normal Vero cell lines. Furthermore, T-1-AFPB induced early (from 0.72 to 19.12) and late (from 0.13 to 6.37) apoptosis in HepG2 cell lines. In conclusion, the combined computational and experimental approaches demonstrated the efficacy and safety of T-1-APFPB providing it as a promising lead VEGFR-2 inhibitor for further development aiming at cancer therapy.

Structure-based virtual screening, molecular docking, molecular dynamics simulation and MM/PBSA calculations towards identification of steroidal and non-steroidal selective glucocorticoid receptor modulators

Glucocorticoids have been used in the treatment of many diseases including inflammatory and autoimmune diseases. Despite the wide therapeutic effects of synthetic glucocorticoids, the use of these compounds has been limited due to side effects such as osteoporosis, immunodeficiency, and hyperglycaemia. To this end, extensive studies have been performed to discover new glucocorticoid modulators with the aim of increasing affinity for the receptor and thus less side effects. In the present work, structure-based virtual screening was used for the identification of novel potent compounds with glucocorticoid effects. The molecules derived from ZINC database were screened on account of structural similarity with some glucocorticoid agonists as the template. Subsequently, molecular docking was performed on 200 selected compounds to obtain the best steroidal and non-steroidal conformations. Three compounds, namely ZINC_000002083318, ZINC_000253697499 and ZINC_000003845653, were selected with the binding energies of -11.5, -10.5, and -9.5 kcal/mol, respectively. Molecular dynamic simulations on superior structures were accomplished with the glucocorticoid receptor. Additionally, root mean square deviations, root mean square fluctuation, radius of gyration, hydrogen bonds, and binding-free energy analysis showed the binding stability of the proposed compounds compared to budesonide as an approved drug. The results demonstrated that all the compounds had suitable binding stability compared to budesonide, while ZINC_000002083318 showed a tighter binding energy compared to the other compounds.Communicated by Ramaswamy H. Sarma.

Aggregation patterns of curcumin and piperine mixtures in different polar media

This work reports a thorough molecular dynamics investigation on the aggregation patterns of curcumin and piperine in water, ethanol and a mixture of both solvents. The low solubility of curcumin in water results in a rapid formation of very stable dimers for both keto and enol tautomers. In agreement with a higher solubility, piperine molecules move closer and farther apart several times during the simulation, which indicates the formation of a less stable dimer in water. In contrast, both curcumin and piperine are soluble in ethanol and, thus, dimers can hardly be formed in this media. In comparison with a pure-water solvent, a 30 : 70 mixture of ethanol and water significantly reduces the probability of formation of most dimers of curcumin and piperine molecules. The simulations show that larger clusters may be complex structures, but the formation of stacks (in the case of piperine and enol tautomer of curcumin) and cages (when the keto tautomer of curcumin is involved) are not rare. Furthermore, it is shown that each single molecule presents a certain degree of mobility in the cluster, especially on the surface, but without leading to dissociation.

MD-Ligand-Receptor: A High-Performance Computing Tool for Characterizing Ligand-Receptor Binding Interactions in Molecular Dynamics Trajectories

Molecular dynamics simulation is a widely employed computational technique for studying the dynamic behavior of molecular systems over time. By simulating macromolecular biological systems consisting of a drug, a receptor and a solvated environment with thousands of water molecules, MD allows for realistic ligand-receptor binding interactions (lrbi) to be studied. In this study, we present MD-ligand-receptor (MDLR), a state-of-the-art software designed to explore the intricate interactions between ligands and receptors over time using molecular dynamics trajectories. Unlike traditional static analysis tools, MDLR goes beyond simply taking a snapshot of ligand-receptor binding interactions (lrbi), uncovering long-lasting molecular interactions and predicting the time-dependent inhibitory activity of specific drugs. With MDLR, researchers can gain insights into the dynamic behavior of complex ligand-receptor systems. Our pipeline is optimized for high-performance computing, capable of efficiently processing vast molecular dynamics trajectories on multicore Linux servers or even multinode HPC clusters. In the latter case, MDLR allows the user to analyze large trajectories in a very short time. To facilitate the exploration and visualization of lrbi, we provide an intuitive Python notebook (Jupyter), which allows users to examine and interpret the results through various graphical representations.

Ligand-Based Virtual Screening, Molecular Docking, and Molecular Dynamic Simulations of New β-Estrogen Receptor Activators with Potential for Pharmacological Obesity Treatment

Obesity is a pandemic and a serious health problem in developed and undeveloped countries. Activation of estrogen receptor beta (ERβ) has been shown to promote weight loss without modifying caloric intake, making it an attractive target for developing new drugs against obesity. This work aimed to predict new small molecules as potential ERβ activators. A ligand-based virtual screening of the ZINC15, PubChem, and Molport databases by substructure and similarity was carried out using the three-dimensional organization of known ligands as a reference. A molecular docking screening of FDA-approved drugs was also conducted as a repositioning strategy. Finally, selected compounds were evaluated by molecular dynamic simulations. Compounds 1 (-24.27 ± 0.34 kcal/mol), 2 (-23.33 ± 0.3 kcal/mol), and 6 (-29.55 ± 0.51 kcal/mol) showed the best stability on the active site in complex with ERβ with an RMSD < 3.3 Å. RMSF analysis showed that these compounds do not affect the fluctuation of the Cα of ERβ nor the compactness according to the radius of gyration. Finally, an in silico evaluation of ADMET showed they are safe molecules. These results suggest that new ERβ ligands could be promising molecules for obesity control.

Nanoscale Iron-Based Metal-Organic Frameworks: Incorporation of Functionalized Drugs and Degradation in Biological Media

Metal-organic frameworks (MOFs) attract growing interest in biomedical applications. Among thousands of MOF structures, the mesoporous iron(III) carboxylate MIL-100(Fe) (MIL stands for the Materials of Lavoisier Institute) is among the most studied MOF nanocarrier, owing to its high porosity, biodegradability, and lack of toxicity. Nanosized MIL-100(Fe) particles (nanoMOFs) readily coordinate with drugs leading to unprecedented payloads and controlled release. Here, we show how the functional groups of the challenging anticancer drug prednisolone influence their interactions with the nanoMOFs and their release in various media. Molecular modeling enabled predicting the strength of interactions between prednisolone-bearing or not phosphate or sulfate moieties (PP and PS, respectively) and the oxo-trimer of MIL-100(Fe) as well as understanding the pore filling of MIL-100(Fe). Noticeably, PP showed the strongest interactions (drug loading up to 30 wt %, encapsulation efficiency > 98%) and slowed down the nanoMOFs' degradation in simulated body fluid. This drug was shown to bind to the iron Lewis acid sites and was not displaced by other ions in the suspension media. On the contrary, PS was entrapped with lower efficiencies and was easily displaced by phosphates in the release media. Noticeably, the nanoMOFs maintained their size and faceted structures after drug loading and even after degradation in blood or serum after losing almost the totality of the constitutive trimesate ligands. Scanning electron microscopy with high annular dark field (STEM-HAADF) in conjunction with X-Ray energy-dispersive spectrometry (XEDS) was a powerful tool enabling the unraveling of the main elements to gain insights on the MOF structural evolution after drug loading and/or upon degradation.

Molecular docking, 3D-QSAR and simulation studies for identifying pharmacophoric features of indole derivatives as 17β-hydroxysteroid dehydrogenase type 5 (17β-HSD5) inhibitors

Excess of androgens leads to various diseases such as Poly-Cystic Ovarian Syndrome, Prostate Cancer, Hirsutism, Obesity and Acne. 17β-Hydroxysteroid Dehydrogenase type 5 (17β-HSD5) converts androstenedione into testosterone peripherally, thereby significantly contributing to the development of these diseases. Indole-bearing scaffolds are reported as potential 17β-HSD5 inhibitors for the manifestation of diseases arising due to androgen excess. In the present work, we have extensively performed a combination of molecular docking, Gaussian field-based 3D-QSAR, Pharmacophore mapping and MD-simulation studies (100 ns) to identify the pharmacophoric features of indole-based compounds as potent 17β-HSD5 inhibitors. Molecular simulation studies of the most potent compound in the binding pocket of enzyme revealed that the compound 11 was stable in the binding pocket and showed good binding affinity through interactions with various residues of active site pocket. The Molecular mechanics Generalized Born surface area continuum solvation (MM/GBSA) and Molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations revealed that the compound 11 possessed a free binding energy of -36.36 kcal/mol and -7.00 kcal/mol, respectively, which was better as compared to reference compound Desmethyl indomethacin (DES). The developed pharmacophore will be helpful to design novel indole-based molecules as potent 17β-HSD5 inhibitors for the treatment of various androgenic disorders.Communicated by Ramaswamy H. Sarma.

Novel Lennard-Jones Parameters for Cysteine and Selenocysteine in the AMBER Force Field

Cysteine is a common amino acid with a thiol group that plays a pivotal role in a variety of scenarios in redox biochemistry. In contrast, selenocysteine, the 21st amino acid, is only present in 25 human proteins. Classical force-field parameters for cysteine and selenocysteine are still scarce. In this context, we present a methodology to obtain Lennard-Jones parameters for cysteine and selenocysteine in different physiologically relevant oxidation and protonation states. The new force field parameters obtained in this work are available at https://github.com/MALBECC/AMBER-parameters-database. The parameters were adjusted to reproduce water radial distribution functions obtained by density functional theory ab initio molecular dynamics. We validated the results by evaluating the impact of the choice of parameters on the structure and dynamics in classical molecular dynamics simulations of representative proteins containing catalytic cysteine/selenocysteine residues. There are significant changes in protein structure and dynamics depending on the parameters choice, specifically affecting the residues close to the catalytic sites.

Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of New Drugs against Human Topoisomerase I Receptor

Human Topoisomerase I (hTop1p) is a ubiquitous enzyme that relaxes supercoiled DNA through a conserved mechanism involving transient breakage, rotation, and binding. Htop1p is the molecular target of the chemotherapeutic drug camptothecin (CPT). It causes the hTop1p-DNA complex to slow down the binding process and clash with the replicative machinery during the S phase of the cell cycle, forcing cells to activate the apoptotic response. This gives hTop1p a central role in cancer therapy. Recently, two artesunic acid derivatives (compounds c6 and c7) have been proposed as promising inhibitors of hTop1p with possible antitumor activity. We used several computational approaches to obtain in silico confirmations of the experimental data and to form a comprehensive dynamic description of the ligand-receptor system. We performed molecular docking analyses to verify the ability of the two new derivatives to access the enzyme-DNA interface, and a classical molecular dynamics simulation was performed to assess the capacity of the two compounds to maintain a stable binding pose over time. Finally, we calculated the noncovalent interactions between the two new derivatives and the hTop1p receptor in order to propose a possible inhibitory mechanism like that adopted by CPT.

Novel inhibitors of the main protease enzyme of SARS-CoV-2 identified via molecular dynamics simulation-guided in vitro assay

For the COVID-19 pandemic caused by SARS-CoV-2, there are currently no effective drugs or vaccines to treat this coronavirus infection. In this study, we focus on the main protease enzyme of SARS-CoV-2, 3CL^pro, which is critical for viral replication. We employ explicit solvent molecular dynamics simulations of about 150 compounds docked into 3CL^pro’s binding site and that had emerged as good main protease ligands from our previous in silico screening of over 1.2 million compounds. By incoporating protein dynamics and applying a range of structural descriptors, such as the ability to form specific contacts with the catalytic dyad residues of 3CL^pro and the structural fluctuations of the ligands in the binding site, we are able to further refine our compound selection. Fourteen compounds including estradiol shown to be the most promising based on our calculations were procured and screened against recombinant 3CL^pro in a fluorescence assay. Eight of these compounds have significant activity in inhibiting the SARS-CoV-2 main protease. Among these are corilagin, a gallotannin, and lurasidone, an antipsychotic drug, which emerged as the most promising natural product and drug, respectively, and might thus be candidates for drug repurposing for the treatment of COVID-19. In addition, we also tested the inhibitory activity of testosterone, and our results reveal testosterone as possessing moderate inhibitory potency against the 3CL^pro enzyme, which may thus provide an explanation why older men are more severely affected by COVID-19.