Molecular dynamic simulation and DFT study on the Drug-DNA interaction; Crocetin as an anti-cancer and DNA nanostructure model

Computer-based identification of olive oil components as a potential inhibitor of neirisaral adhesion a regulatory protein

In silico methods such as molecular docking and molecular dynamic (MD) simulations have significant interest due to their ability to identify the protein-ligand interactions at the atomic level. In this work, different computational methods were used to elucidate the ability of some olive oil components to act as Neisseria adhesion A Regulatory protein (NadR) inhibitors. The frontier molecular orbitals (FMOs) and the global properties such as global hardness, electronegativity, and global softness of ten olive oil components (α-Tocopherol, Erythrodiol, Hydroxytyrosol, Linoleic acid, Apigenin, Luteolin, Oleic acid, Oleocanthal, Palmitic acid, and Tyrosol) were reported using Density Functional Theory (DFT) methods. Among all investigated compounds, Erythrodiol, Apigenin, and Luteolin demonstrated the highest binding affinities (-8.72, -7.12, and -8.24 kcal/mol, respectively) against NadR, compared to -8.21 kcal/mol of the native ligand based on molecular docking calculations. ADMET properties and physicochemical features showed that Erythrodiol, Apigenin, and Luteolin have good physicochemical features and can act as drugs candidate. Molecular dynamics (MD) simulations demonstrated that Erythrodiol, Apigenin, and Luteolin show stable binding affinity and molecular interaction with NadR. Further Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analyses using the MD trajectories also demonstrated the higher binding affinity of Erythrodiol, Apigenin and Luteolin inside NadR protein. The overall study provides a rationale to use Erythrodiol, Apigenin, and Luteolin in the drug development as anti-adhesive drugs lead. Communicated by Ramaswamy H. Sarma.

Crocetin: A Systematic Review

Crocetin is an aglycone of crocin naturally occurring in saffron and produced in biological systems by hydrolysis of crocin as a bioactive metabolite. It is known to exist in several medicinal plants, the desiccative ripe fruit of the cape jasmine belonging to the Rubiaceae family, and stigmas of the saffron plant of the Iridaceae family. According to modern pharmacological investigations, crocetin possesses cardioprotective, hepatoprotective, neuroprotective, antidepressant, antiviral, anticancer, atherosclerotic, antidiabetic, and memory-enhancing properties. Although poor bioavailability hinders therapeutic applications, derivatization and formulation preparation technologies have broadened the application prospects for crocetin. To promote the research and development of crocetin, we summarized the distribution, preparation and production, total synthesis and derivatization technology, pharmacological activity, pharmacokinetics, drug safety, drug formulations, and preparation of crocetin.

Bis-indole alkaloid caulerpin from a new source Sargassum platycarpum: isolation, characterization, in vitro anticancer activity, binding with nucleobases by DFT calculations and MD simulation

Caulerpin, a bis-indole alkaloid is isolated from a new source Sargassum platycarpum, brown alga (family Sargassaceae) for the first time. The structure of caulerpin was characterized by IR, H1NMR, C13 NMR, HSQC, HMBC, EI-MS spectroscopy. Antifungal results suggest that caulerpin has been inhibited Cryptococcus neoformas (12 mm) and Candida albicans (7 mm) than other microbes. In vitro anticancer activity of caulerpin has been explored by cell viability assay against new human cancer cell line (liver-HepG2). The results show that caulerpin has low IC50 value (24.6 ± 2.1 µg/mL) against HepG-2. Based on the least toxic activity of caulerpin, these results encourage for future in vivo anticancer study. The binding of caulerpin molecule with the two nucleobases (T/U) bases has been studied by DFT methods. According to the AIM analysis, there are two types of interactions between caulerpin and T/U bases partially covalent partially electrostatic and electrostatic in gas and water phases. Based on NBO analysis, the charges were transferred from the lone-pair (n) in orbitals of O atoms of caulerpin to the σ* orbitals of T/U bases atoms. ΔEbin in the state of caulerpin-T bases complexes are lower than those in the caulerpin-U bases complexes in both gas and water phase. MD simulation supported that caulerpin-T/U bases complexes are stable in presence of explicit water phase. Thus, the findings of our study will be useful for giving an insight into the caulerpin/bases complexes that could be helpful in future experimental studies to develop the performance of caulerpin molecules as natural candidate drug. Communicated by Ramaswamy H. Sarma.

Cytotoxic activity, molecular docking, pharmacokinetic properties and quantum mechanics calculations of the brown macroalga Cystoseira trinodis compounds

In this study, nine compounds were isolated, eight of them were isolated for the first time from Cystoseira trinodis. The biological activity of the extract, fractions and pure compounds was evaluated. The antimicrobial activity was investigated against 3 fungi species, 3 gram + ve and 3 gram -ve bacteria. The crude extract and fractions showed moderate inhibition against some of the tested microorganisms, especially the butanol fraction exhibited the maximum inhibition zone against Salmonella typhimurium (16 ± 0.60 mm). Cytotoxicity was evaluated against HepG-2 and MCF-7 cell lines. Hexane fraction exhibited the highest cytotoxic effect against HepG-2 and MCF-7 cell lines with an IC₅₀ value of 14.3 ± 0.8 and 19.2 ± 0.7 µg/ml, respectively with compared to other fractions. The isolates were identified as octacosanoic acid (1), glyceryl trilinoleate (2), oleic acid (3), and the epimeric mixture of saringosterols (4, 5), β-sitosterol (6), glycoglycerolipid (7) and a mixture of kjellmanianone and loliolide (8, 9) by spectroscopic analysis. Among the all tested compounds kjellmanianone and loliolide mixture exhibited significant cytotoxic activity with an IC₅₀ value of 7.27 µg/ml against HepG-2 cells. The major and minor constituents of the extract and fractions were identified using GC-MS analysis. Molecular docking analysis confirmed that most of the studied compounds especially compounds 8 and 9 strongly interact with TPK and VEGFR-2 with highest binding energies supported that the high cytotoxicity of these compounds against human hepatocellular cancer in the experimental part. The energetic, geometric and topological properties of compounds 8 and 9 binding with cytosine base were computed by DFT methods. Molecular properties descriptors, bioactivity score and ADMET analysis confirmed that most of the studied compounds especially compounds 8 and 9 exhibit significant biological activities and have a better chance to be developed as drug leads. Communicated by Ramaswamy H. Sarma.

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches

Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E _HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E _HB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇²ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E _HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Structural insights into the anti-cancer activity of quercetin on G-tetrad, mixed G-tetrad, and G-quadruplex DNA using quantum chemical and molecular dynamics simulations

Human telomerase referred as 'terminal transferase' is a nucleoprotein enzyme which inhibits the disintegration of telomere length and act as a drug target for the anticancer therapy. The tandem repeating structure of telomere sequence forms the guanine-rich quadruplex structures that stabilize stacked tetrads. In our present work, we have investigated the interaction of quercetin with DNA tetrads using DFT. Geometrical analysis revealed that the influence of quercetin drug induces the structural changes into the DNA tetrads. Among DNA tetrads, the quercetin stacked with GCGC tetrad has the highest interaction energy of -88.08 kcal/mol. The binding mode and the structural stability are verified by the absorption spectroscopy method. The longer wavelength was found at 380 nm and it exhibits bathochromic shift. The findings help us to understand the binding nature of quercetin drug with DNA tetrads and it also inhibits the telomerase activity. Further, the quercetin drug interacted with G-quadruplex DNA by using molecular dynamics (MD) simulation studies for 100 ns simulation at different temperatures and different pH levels (T = 298 K, 320 K and pH = 7.4, 5.4). The structural stability of the quercetin with G-quadruplex structure is confirmed by RMSD. For the acidic condition (pH = 5.4), the binding affinity is higher toward G-quadruplex DNA, this result resembles that the quercetin drug is well interacted with G-quadruplex DNA at acidic condition (pH = 7.4) than the neutral condition. The obtained results show that quercetin drug stabilizes the G-quadruplex DNA, which regulates telomerase enzyme and it potentially acts as a novel anti-cancer agent.Communicated by Ramaswamy H. Sarma.

Drug–DNA interaction, a joint DFT-D3/MD study on safranal as an anticancer and DNA nanostructure model

In this research, using a combination of quantum mechanics and molecular dynamic (MD) simulations, the interaction of safranal (2,6,6-trimethylcyclohexa-1,3-dien-1-carboxaldehyde) as an anti-cancer drug and Dickerson B-DNA was studied. MD simulations were executed for 35 ns in water. Binding energy analysis in three definite parts of the B-DNA and comparison between different contributions of the binding energy shows that the van der Waals energy part of the interaction is impressive among the standard molecular mechanic energy terms. On the basis of Gibbs energies, it is confirmed that the most important interactions in the safranal complex are related to the A–T and C–G rich regions, which is in agreement with the experimental data. Quantum theory of atoms in molecules and natural bond orbital analyses were applied. A diminution in the electronic chemical potential of the safranal–DNA complex in comparison with the isolated DNA, 0.026 and 0.022 au for the S1 region and 0.012 and 0.017 au for the S2 region, was obtained in the gas phase and water, respectively, which increases the complex stability. An enhancement in the electrophilicity character, during the complexation process, shows the electron charge flux between the safranal and DNA, especially in water. The strengths of the CH⋯O bonds at the center of safranal–DNA interaction were also evaluated. A mean value of 0.06 au for the electron density of the bond critical point of the H⋯O in the complex confirms the H-bond formation during the complexation.

Identification of protegrin-1 as a stable and nontoxic scaffold among protegrin family - a computational approach

Achieving both, nontoxicity and stability in antimicrobial peptides (AMP) is a challenge. This study predicts a structurally stable, nontoxic scaffold among the protegrin family, for future therapeutic peptide analogs. Protegrins (PG) are a class of pharmaceutically approved, in demand AMPs, which require further improvement in terms of nontoxicity and stability. Out of five protegrins viz., PG1, PG2, PG3, PG4 and PG5, PG1 has been predicted as best scaffold. Prediction was based upon sequential elimination of other protegrins, using computational methods to assess the extracellular bacterial membrane penetrability, nontoxicity and structural stability by geometric observables. Initially, PG2 and PG4 showing the lowest membrane penetrability and highest toxicity respectively, were screened out. Among the remaining three protegrins, PG1 displayed both lowest root mean square deviation and radius of gyration, with a considerable occupancy of seven H-bonds and established uniform secondary structure profile throughout its ensembles. Therefore, the authors claim the superiority of PG1 as a nontoxic stable scaffold among its family. Communicated by Ramaswamy H. Sarma.

Capecitabine as a minor groove binder of DNA: molecular docking, molecular dynamics, and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^-1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.