In silico thermodynamics stability change analysis involved in BH4 responsive mutations in phenylalanine hydroxylase: QM/MM and MD simulations analysis

Secondary metabolites of Trichoderma spp. as EGFR tyrosine kinase inhibitors: Evaluation of anticancer efficacy through computational approach

The present study explores the epidermal growth factor receptor (EGFR) tyrosine kinase inhibition efficacy of secondary metabolites in Trichoderma spp. through molecular docking, molecular dynamics (MD) simulation and MM-PBSA approach. The result of molecular docking confirmed that out of 200 metabolites screened, three metabolites such as Harzianelactone A, Pretrichodermamide G and Aspochalasin M, potentially bound with the active binding site of EGFR tyrosine kinase domain(PDB ID: 1M17) with a threshold docking score of ≤- 9.0 kcal/mol when compared with the standard EGFR inhibitor (Erlotinib). The MD simulation was run to investigate the potential for stable complex formation in EGFR tyrosine kinase domain-unbound/lead metabolite (Aspochalasin M)-bound/standard inhibitor (Erlotinib)-bound complex. The MD simulation analysis at 100 ns revealed that Aspochalasin M formed the stable complex with EGFR. Besides, the in silico predication of pharmacokinetic properties further confirmed that Aspochalasin M qualified the drug-likeness rules with no harmful side effects (viz., hERG toxicity, hepatotoxicity and skin sensitization), non-mutagenicity and favourable logBB value. Moreover, the BOILED-Egg model predicted that Aspochalasin M showed a higher gastrointestinal absorption with improved bioavailability when administered orally and removed from the central nervous system (CNS). The results of the computational studies concluded that Aspochalasin M possessed significant efficacy in binding EGFR's active sites compared to the known standard inhibitor (Erlotinib). Therefore, Aspochalasin M can be used as a possible anticancer drug candidate and further in vitro and in vivo experimental validation of Aspochalasin M of Trichoderma spp. are required to determine its anticancer potential.

Biological Evaluation and Binding Mechanism of 5-HT7 Specific Arylpiperazinyl-Alkyl Benzothiazolone: Radiobiology and Photo-physical Studies

Diversely substituted methoxy derivatives of arylpiperazinyl-alkyl benzothiazolone has been evaluated as specific probe for 5HT7. To determine the best methoxy derivative for 5HT7 receptor affinity, we synthesised a number of 2-benzothiazolone arylalkyl piperazine derivatives. In-vitro/vivo studies with C-2 substituted [11C]ABT showed 5HT7 specific binding. The radiochemical purity of [11C]ABT was found to be more than 99% with radiochemical stability persistence for more than 1.5 hr at 25 °C. The interaction of BSA and ABT has been analysed by photophysical studies for better understanding of properties such as adsortion, distribution, metabolism and elemination (ADME). The interaction between ABT and BSA was analyzed by using the UV-vis and fluorescence spectra. UV-vis spectra analyzed the changes in primary structure of BSA on its interaction with ABT. ABT showed quenched fluorescence emission intensity of tryptophan residues in BSA via static quenching mechanism. This study might help to understand how ABT binds to serum protein or subsequently to know the ADME of this drug candidate.

How living cells are affected during the cold atmospheric pressure plasma treatment

When the electric discharge process is limited by high voltage electrodes shielding, the ionization measure would be controlled to less than one percent and the temperature to less than 37 °C even at atmospheric pressure, so-called cold atmospheric pressure plasma (CAP). CAP has been found to have profound medical applications in association with its reactive oxygen and nitrogen species (ROS/RNS). In this way that during plasma exposure, the subjected medium (e.g. cell cytoplasmic membrane in plasma therapy) interacts with ROS/RNS. Accordingly, a precise study of the mentioned interactions and their consequences on the cells' behavior changes, is necessary. The results lead to the reduction of possible risks and provide the opportunity of optimizing the efficacy of CAP before the development of CAP applications in the field of plasma medicine. In this report molecular dynamic (MD) simulation is used to investigate the mentioned interactions and a proper and compatible comparison with the experimental results is presented. Based on this, the effects of H2O2, NO and O2 on the living cell's membrane are investigated in biological conditions. Our results show that: i) The hydration of phospholipid polar heads would be enhanced associated with the H2O2 presence. ii) A new definition of the surface area assigned to each phospholipid (APL), more reliable and compatible with the physical expectations, is introduced. iii) The long-term behavior of NO and O2 is their penetration into the lipid bilayer and sometimes passing through the membrane into the cell. The latter would be an indication of internal cells' pathways activation leading to modification of cells' function.

In Silico Computational Studies of Bioactive Secondary Metabolites from Wedelia trilobata against Anti-Apoptotic B-Cell Lymphoma-2 (Bcl-2) Protein Associated with Cancer Cell Survival and Resistance

In the present study, the binding affinity of 52 bioactive secondary metabolites from Wedelia trilobata towards the anti-apoptotic B-cell lymphoma-2 (Bcl-2) protein (PDB: 2W3L) structure was identified by using in silico molecular docking and molecular dynamics simulation. The molecular docking results demonstrated that the binding energies of docked compounds with Bcl-2 protein ranged from -5.3 kcal/mol to -10.1 kcal/mol. However, the lowest binding energy (-10.1 kcal/mol) was offered by Friedelin against Bcl-2 protein when compared to other metabolites and the standard drug Obatoclax (-8.4 kcal/mol). The molecular dynamics simulations revealed that the Friedelin-Bcl-2 protein complex was found to be stable throughout the simulation period of 100 ns. Overall, the predicted Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties of Friedelin are relatively better than Obatoclax, with the most noticeable differences in many parameters where Friedelin has no AMES toxicity, hepatotoxicity, and skin sensitization. The ADMET profiling of selected compounds supported their in silico drug-likeness properties. Based on the computational analyses, the present study concluded that Friedelin of W. trilobata was found to be the potential inhibitor of the Bcl-2 protein, which merits attention for further in vitro and in vivo studies before clinical trials.

Phytoconstituents of Withania somnifera unveiled Ashwagandhanolide as a potential drug targeting breast cancer: Investigations through computational, molecular docking and conceptual DFT studies

Breast cancer is the second most common malignancy in females worldwide and poses a great challenge that necessitates the identification of novel therapeutic agents from several sources. This research aimed to study the molecular docking and molecular dynamics simulations of four proteins (such as PDB: 6CBZ, 1FDW, 5GWK and 2WTT) with the selected phytochemicals from Withania somnifera to identify the potential inhibitors for breast cancer. The molecular docking result showed that among 44 compounds, two of them, Ashwagandhanolide and Withanolide sulfoxide have the potential to inhibit estrogen receptor alpha (ERα), 17-beta-hydroxysteroid -dehydrogenase type 1 (17β-HSD1), topoisomerase II alpha (TOP2A) and p73 tetramerization domain that are expressed during breast cancer. The molecular dynamics (MD) simulations results suggested that Ashwagandhanolide remained inside the binding cavity of four targeted proteins and contributed favorably towards forming a stable protein-ligand complex throughout the simulation. Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties confirmed that Ashwagandhanolide is hydrophobic and has moderate intestinal permeability, good intestinal absorption, and poor skin permeability. The compound has a relatively low VDss value (-1.652) and can be transported across ABC transporter and good central nervous system (CNS) permeability but did not easily cross the blood-brain barrier (BBB). This compound does not possess any mutagenicity, hepatotoxicity and skin sensitization. Based on the results obtained, the present study highlights the anticancer potential of Ashwagandhanolide, a compound from W. somnifera. Furthermore, in vitro and in vivo studies are necessary to perform before clinical trials to prove the potentiality of Ashwagandhanolide.

Determination of Hybrid TSPO Ligands with Minimal Impact of SNP (rs6971) through Molecular Docking and MD Simulation Study

Background:

In an endeavor to ascertain high-affinity TSPO ligands with minimal single nucleotide polymorphism (SNP), six hybrid molecules have been identified as new leads for future inflammation PET imaging.

Objective:

Genesis for chemical design was encouraged from structural families of well-known ligands FEBMP and PBR28/ DAA1106 that have demonstrated remarkable TSPO binding characteristics.

Methods:

All proposed hybrid ligands 1-6 are subjected to molecular docking and simulation studies with wild and mutant protein to study their interactions, binding, consistency of active conformations and are correlated with well-established TSPO ligands.

Results:

Each hybrid ligand demonstrate better docking score > -11.00 kcal/mol with TSPO with respect to gold standard PK11195, i.e., -11.00 kcal/mol for 4UC3 and -12.94 kcal/mol for 4UC1. On comparison with FEBMP and GE-180 (-12.57, -7.24 kcal/mol for 4UC3 and -14.10, -11.32 kcal/mol for 4UC1), ligand 3 demonstrates maximum docking energy (> -15.50 kcal/mol) with minimum SNP (0.26 kcal/mol).

Discussion:

Presence of strong hydrogen bond Arg148-3.27Å (4UC1) and Trp50-2.43Å, Asp28- 2.57Å (4UC3) apart from short-range interactions, including π–π interactions with the aromatic residues, such as (Trp39, Phe46, Trp135) and (Trp39, Trp108), attributes towards its strong binding.

Conclusion:

Utilizing the results of binding energy, we concluded stable complex formation of these hybrid ligands that could bind to TSPO with the least effect of SNP with similar interactions to known ligands. Overall, ligand 3 stands out as the best ligand having insignificant deviations per residue of protein that can be further explored and assessed in detail for future inflammation PET application after subsequent detailed biological evaluation.

Synthesis, molecular docking, molecular dynamics and evaluation of Drug-Likeness properties of the fused N-Formyl pyrazoline substituted new dehydroepiandrosterone derivatives

The hybrid molecules bearing heterocyclic structures in the A or D rings of steroids have significant biological activity. 16 (E)-Hetereoarylidene steroids were synthesized from the reaction of different heteroaromatic carbaldehydes and trans-Dehydroepiandrosterone (DHEA) in a basic medium. Then, synthesis of the N-formyl pyrazoline substituted new DHEA derivatives were carried out from the reaction of hydrazine hydrate and 16 (E)-hetereoarylidene steroids. The structures of the synthesized compounds were elucidated by elemental analysis, FT-IR, ¹H NMR, and ¹³C NMR spectroscopy. To investigate the activation pathway of synthesized N-formyl pyrazoline substituted steroid derivatives, a molecular docking study was performed on human cytochrome P450-(CYP17A1: PDB ID 5IRQ) with the help of the free AutoDock Vina. 100 ns molecular dynamic simulation process was performed to monitor the behavior of the complex structure formed by CYP17A1 and to calculate the stability over time of 2a and 2d (-9.8 kcal/mol), which gave the lowest value according to the results obtained in the molecular docking study with AutoDock Vina. Accordingly, RMSD, RMSF, Rg, and SASA analyzes of 2a and 2d were performed, and MMPBSA was calculated. Lastly, the ADMET (absorption, distribution, metabolism, excretion, and toxicity) analyses of the novel steroid derivatives were investigated.Communicated by Ramaswamy H. Sarma.

Lanthanide (Ln3+) complexes of bifunctional chelate: Synthesis, physicochemical study and interaction with human serum albumin (HSA)

Bifunctional chelate EDTA-bis amide (N,N'-bis (tyramide)ethylenediamine-N,N'-diacetic acid) that has ability to mimic natural amino acids was synthesized and analyzed by various spectroscopic techniques. The physicochemical studies were performed to calculate the various thermodynamic and kinetic parameters for the synthesized poly-amino carboxylate ligand. The two protonation constant (pka's = 3.460 and 6.722) of the prepared ligand and stability constants (log K_ML's = 15.8, 18.1, 16.2, 18.4, 17.5, 18.9, 13.6 and 12.8) of the complexes formed with Ce³⁺, Sm³⁺, Eu³⁺, Gd³⁺, Tb³⁺, Lu³⁺, Zn²⁺ and Cu²⁺ were determined by potentiometric titration using 0.1 M Me₄NOH as non-aqueous base. The formation kinetics of [EuEDTA-TA₂]⁺ and [CeEDTA-TA₂]⁺ was studied and the rate constants were found to be 2.95 × 10-5 s^-1 and 4.414 × 10-5 s^-1respectively including the exchange reaction of [EuEDTA-TA₂]⁺ with Zn²⁺ and Cu²⁺ spectrophotometrically. The Eu(III) complex of EDTA(TA)₂ gives three emission bands at 480 nm, 540 nm and 610 nm (λ_max = 270 nm, excitation) which shows efficacy of the ligand as an optical imaging agent. Molecular docking studies with Human Serum Albumin (HSA: PDB 1E78) showed binding pattern with the residues Arg218, Arg222, Lys195 and Lys444 in sub domain II A of site I via hydrogen bond and identifies the ligand-HSA interaction and specific insight for transportation to the target sites. Subsequently, fluorescence spectroscopy was performed at λ_ex = 350 nm binding constant for HSA was 5.847 × 10⁴ M^-1 which showed effective quenching effect.

Structure-based screening and validation of bioactive compounds as Zika virus methyltransferase (MTase) inhibitors through first-principle density functional theory, classical molecular simulation and QM/MM affinity estimation

Recent Zika virus (ZIKV) outbreak and association with human diseases such as neurological disorders have raised global health concerns. However, in the absence of an approved anti-ZIKV drug has generated urgency for the drug development against ZIKV infection. Here, structure-based virtual screening of 8589 bioactive compounds, screened at the substrate-binding site of ZIKV nonstructural 5 (NS5)-based structure N-terminal methyltransferase (MTase) domain followed by ADMET (absorption, distribution, metabolism, excretion and toxicity) profiling concluded the four potential lead inhibitors, i.e. (4-acetylamino-benzenesulfonylamino)-acetic acid (F3342-0450), 3-(5-methylfuran-2-yl)-N-(4-sulfamoylphenyl)propanamide (F1736-0142), 8-(2-hydroxy-ethylamino)-1,3-dimethyl-7-(3-methyl-benzyl)-3,7-dihydro-purine-2,6-dione (F0886-0080) and N-[4-(aminosulfonyl)phenyl]-2,3-dihydro-1,4-benzodioxine-2-carboxamide (F0451-2187). Collectively, extra precision docking and Density Functional Theory(DFT) calculations studies identified the F3342-0450 molecule, having strong interactions on the active site of MTase, further supported by molecular dynamics simulation, binding affinity and hybrid QM/MM calculations, suggest a new drug molecule for the antiviral drug development against ZIKV infection. Communicated by Ramaswamy H. Sarma.

Computational prediction of interaction and pharmacokinetics profile study for polyamino-polycarboxylic ligands on binding with human serum albumin

Human serum albumin (HSA) is one of the most abundant plasma proteins available in blood and responsible for transport of fatty acids, drugs and metabolites at its binding sites which are very important for the assessment of pharmacokinetics profile of the polyamino-polycarboxylic ligands.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.

Evaluation of potential flavonoid inhibitors of glyoxalase-I based on virtual screening and in vitro studies

Glyoxalase-I (GLO-I) is a component of the ubiquitous detoxification system involved in the conversion of methylglyoxal (MG) to d-lactate in the glycolytic pathway. MG toxicity arises from its ability to form advanced glycation end products. GLO-I has been reported to be frequently overexpressed in various types of cancer cells. In this study, we performed structure-based virtual screening of focused flavonoids commercial library to identify potential and specific inhibitors of GLO-I. The compounds were ranked based on Glide extra precision docking score and five hits (curcumin, quercetin, morin, naringin and silibinin) were selected on the basis of their interaction with active site amino acid residues of GLO-I. Mixed mode QM/MM calculation was performed on the top-scoring hit to ascertain the role of zinc ion in ligand binding. In addition, the identified hits were subjected to MM/GBSA binding energy prediction, ADME prediction and similarity studies. The hits were tested in vitro for cell viability, and GLO-I inhibition. Naringin (ST072162) was found to be most potent inhibitor of GLO-I among the identified hits with highest glide XP dock score of -14.906. These findings suggest that naringin could be a new scaffold for designing inhibitors against GLO-I with potential application as anticancer agents.