Discovery of novel low-molecular-weight HIV-1 inhibitors interacting with cyclophilin A using in silico screening and biological evaluations

Cyclophilin A has attracted attention recently as a new target of anti-human immunodeficiency virus type 1 (HIV-1) drugs. However, so far no drug against HIV-1 infection exhibiting this mechanism of action has been approved. To identify new potent candidates for inhibitors, we performed in silico screening of a commercial database of more than 1,300 drug-like compounds by using receptor-based docking studies. The candidates selected from docking studies were subsequently tested using biological assays to assess anti-HIV activities. As a result, two compounds were identified as the most active. Specifically, both exhibited anti-HIV activity against viral replication at a low concentration and relatively low cytotoxicity at the effective concentration inhibiting viral growth by 50 %. Further modification of these molecules may lead to the elucidation of potent inhibitors of HIV-1.

Discovery of potent inhibitors for interleukin-2-inducible T-cell kinase: structure-based virtual screening and molecular dynamics simulation approaches

In our study, a structure-based virtual screening study was conducted to identify potent ITK inhibitors, as ITK is considered to play an important role in the treatment of inflammatory diseases. We developed a structure-based pharmacophore model using the crystal structure (PDB ID: 3MJ2) of ITK complexed with BMS-50944. The most predictive model, SB-Hypo1, consisted of six features: three hydrogen-bond acceptors (HBA), one hydrogen-bond donor (HBD), one ring aromatic (RA), and one hydrophobic (HY). The statistical significance of SB-Hypo1 was validated using wide range of test set molecules and a decoy set. The resulting well-validated model could then be confidently used as a 3D query to screen for drug-like molecules in a database, in order to retrieve new chemical scaffolds that may be potent ITK inhibitors. The hits retrieved from this search were filtered based on the maximum fit value, drug-likeness, and ADMET properties, and the hits that were retained were used in a molecular docking study to find the binding mode and molecular interactions with crucial residues at the active site of the protein. These hits were then fed into a molecular dynamics simulation to study the flexibility of the activation loop of ITK upon ligand binding. This combination of methodologies is a valuable tool for identifying structurally diverse molecules with desired biological activities, and for designing new classes of selective ITK inhibitors.

In vitro inhibitory profile of NDGA against AChE and its in silico structural modifications based on ADME profile

Acetylcholinesterase (AChE) inhibitors are currently in focus for the pharmacotherapy of Alzheimer’s disease (AD). These inhibitors increase the level of acetylcholine in the brain and facilitate cholinergic neurotransmission. AChE inhibitors such as rivastigmine, galantamine, physostigmine and huperzine are obtained from plants, indicating that plants can serve as a potential source for novel AChE inhibitors. We have performed a virtual screening of diverse natural products with distinct chemical structure against AChE. NDGA was one among the top scored compounds and was selected for enzyme kinetic studies. The IC₅₀ of NDGA on AChE was 46.2 μM. However, NDGA showed very poor central nervous system (CNS) activity and blood–brain barrier (BBB) penetration. In silico structural modification on NDGA was carried out in order to obtain derivatives with better CNS activity as well as BBB penetration. The studies revealed that some of the designed compounds can be used as lead molecules for the development of drugs against AD

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Exploration of various electronic properties along the reaction coordinate for hydration of Pt(II) and Ru(II) complexes; the CCSD, MPx, and DFT computational study

In the study behavior of molecular electrostatic potential, averaged local ionization energy, and reaction electronic flux along the reaction coordinate of hydration process of three representative Ru(II) and Pt(II) complexes were explored using both post-HF and DFT quantum chemical approximations. Previously determined reaction mechanisms were explored by more detailed insight into changes of electronic properties using ωB97XD functional and MP2 method with 6–311++G(2df,2pd) basis set and CCSD/6–31(+)G(d,p) approach. The dependences of all examined properties on reaction coordinate give more detailed understanding of the hydration process.

Stability and electronic properties of 3D covalent organic frameworks

Covalent organic frameworks (COFs) are a class of covalently linked crystalline nanoporous materials, versatile for nanoelectronic and storage applications. 3D COFs, in particular, have very large pores and low mass densities. Extensive theoretical studies of their energetic and mechanical stability, as well as their electronic properties, have been carried out for all known 3D COFs. COFs are energetically stable and their bulk modulus ranges from 3 to 20 GPa. Electronically, all COFs are semiconductors with band gaps corresponding to the HOMO–LUMO gaps of the building units.

Insight on the interaction of polychlorobiphenyl with nucleic acid–base

The interaction between one polychlorobiphenyl (3,3′,4,4′,-tetrachlorobiphenyl, coded PCB77) and the four DNA nucleic acid–base is studied by means of quantum mechanics calculations in stacked conformations. It is shown that even if the intermolecular dispersion energy is the largest component of the total interaction energy, some other contributions play a non negligible role. In particular the electrostatic dipole-dipole interaction and the charge transfer from the nucleobase to the PCB are responsible for the relative orientation of the monomers in the complexes. In addition, the charge transfer tends to flatten the PCB, which could therefore intercalate more easily between DNA base pairs. From these seminal results, we predict that PCB could intercalate completely between two base pairs, preferably between Guanine:Cytosine pairs.

Molecular dynamics simulation of temperature induced unfolding of animal prion protein

To elucidate the structural stability and the unfolding dynamics of the animal prion protein, the temperature induced structural evolution of turtle prion protein (tPrPc) and bank vole prion protein (bvPrPc) have been performed with molecular dynamics (MD) simulation. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet. Extension and disruption of β-sheet commonly appeared in the temperature induced unfolding process. The conversion of α-helix to π-helix occurred more readily at the elevating temperature. Furthermore, it was suggested in this work that the unfolding of prion protein could be regulated by the temperature.

Perspectives on the reaction force constant

A synchronous, concerted chemical process is rigorously divided by the reaction force F(R), the negative gradient of V(R), into “reactant” and “product” regions which are dominated by structural changes and an intervening “transition” region which is electronically intensive. The reaction force constant κ(R), the second derivative of V(R), is negative throughout the transition region, not just at the nominal transition state, at which κ(R) has a minimum. This is consistent with experimental evidence that there is a transition region, not simply a specific point. We show graphically that significant nonsynchronicity in the process is associated with the development of a maximum of κ(R) in the transition region, which increases as the process becomes more nonsynchronous. (We speculate that for a nonconcerted process this maximum is actually positive.) Thus, κ(R) can serve as an indicator of the level of nonsynchronicity.

DFT study of the Lewis acid mediated synthesis of 3-acyltetramic acids

The recent emergence and re-emergence of alphaviruses, in particular the chikungunya virus (CHIKV), in numerous countries has invoked a worldwide threat to human health, while simultaneously generating an economic burden on affected countries. There are currently no vaccines or effective drugs available for the treatment of the CHIKV, and with few lead compounds reported, the vital medicinal chemistry is significantly more challenging. This study reports on the discovery of potential inhibitors for the nsP3 macro domain of CHIKV using molecular docking, virtual screening, and molecular dynamics simulations, as well as work done to evaluate and confirm the active site of nsP3. Virtual screening was carried out based on blind docking as well as focused docking, using the database of 1541 compounds from NCI Diversity Set II, to identify hit compounds for nsP3. The top hit compounds were further subjected to molecular dynamic simulations, yielding a greater understanding of the dynamic behavior of nsP3 and its complexes with various ligands, concurrently confirming the outcomes of docking, and establishing in silico lead compounds which target the CHIKV nsP3 enzyme.

A theoretical study on the hydrogen adducts of diamidocarbenes and diaminocarbenes

The hybrid-meta GGA DFT functional M06-2X was used to examine the potential of N,N′-diamidocarbenes for use as hydrogen storage materials. We previously discovered that borylene, which is isoelectronic with an Arduengo-type carbene, was a suitable candidate for a hydrogen storage material. We compared the capabilities of N,N′-diamidocarbenes and N-heterocyclic carbenes as hydrogen storage materials. The results indicate that diamidocarbenes are not suitable hydrogen storage materials because the removal of H₂ is more endothermic for diamidocarbenes than for diaminocarbenes.

Regioselectivity in Sonogashira synthesis of 6-(4-nitrobenzyl)-2-phenylthiazolo[3,2-b]1,2,4-triazole: a quantum chemistry study

In the present research, the experimentally observed regioselectivity in Sonogashira synthesis of 6-(4-nitrobenzyl)-2-phenylthiazolo[3,2-b]1,2,4triazole has been modeled by means of density functional theory (DFT) employed to investigate the structural and thermochemical aspects of this synthesis in the gas and solution phases. Comparison of our calculated structural parameters of the title compound with the available X-ray crystallographical data demonstrate a reliable agreement. Then, the effect of two different solvents, DMF and ethanol, are examined via polarized continuum model calculations, showing a significant decrease in the computed values of the reaction enthalpy and free energy changes compared with the gas phase results. We have also considered two tautomeric structures of the intermediate species that it seems the mode of its intermolecular cyclization has an important role in regioselectivity of the final products. Moreover, all obtained results in the gas and solution phases also confirm that the synthesis of the title compound is thermodynamically more favorable than the other regioisomeric product. We also discuss the thermodynamical feasibility of this reaction at higher temperatures. Finally, we concentrate on the survey of substituent effect by choosing electron-withdrawing and electron-donating groups on the aryl iodide. Our calculated thermochemical data in the gas and solution phases indicate that the use of electron-withdrawing moieties is more favorable thermodynamically than electron-donating ones which has been previously concluded via the experimental elucidations.

Conformational entropy of a polymer chain grafted to rough surfaces

A polymer molecule (represented by a statistical chain) end-grafted to a topologically rough surface was studied by static MC simulations. A modified self-avoiding walk on a cubic lattice was used to model the polymer in an athermal solution. Different statistical models of surface roughness were applied. Conformational entropies of chains attached to uncorrelated Gaussian, Brownian, and fractional Brownian surfaces were calculated. Results were compared with the predictions of a simple analytical model of a macromolecule end-grafted to a fractal surface.

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies

We describe a procedure for performing quantitative analyses of fields f(r) on molecular surfaces, including statistical quantities and locating and evaluating their local extrema. Our approach avoids the need for explicit mathematical representation of the surface and can be implemented easily in existing graphical software, as it is based on the very popular representation of a surface as collection of polygons. We discuss applications involving the volumes, surface areas and molecular surface electrostatic potentials, and local ionization energies of a group of 11 molecules.

Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif

Protein farnesyltransferase (FTase) is an important target in many research fields, more markedly so in cancer investigation since several proteins known to be involved in human cancer development are thought to serve as substrates for FTase and to require farnesylation for proper biological activity. Several FTase inhibitors (FTIs) have advanced into clinical testing. Nevertheless, despite the progress in the field several functional and mechanistic doubts on the FTase catalytic activity have persisted. This work provides some crucial information on this important enzyme by describing the application of molecular dynamics simulations using specifically designed molecular mechanical parameters for a variety of 22 CaaX peptides known to work as natural substrates or inhibitors for this enzyme. The study involves a comparative analysis of several important molecular aspects, at the mechanistic level, of the behavior of substrates and inhibitors at the dynamic level, including the behavior of the enzyme and peptides, as well as their interaction, together with the effect of the solvent. Properties evaluated include the radial distribution function of the water molecules around the catalytically important zinc metal atom and cysteine sulfur of CaaX, the conformations of the substrate and inhibitor and the corresponding RMSF values, critical hydrogen bonds, and several catalytically relevant distances. These results are discussed in light of recent experimental and computational evidence that provides new insights into the activity of this enzyme.

Coarse-grained simulations for organic molecular liquids based on Gay-Berne and electric multipole potentials

Coarse-grained studies of CH₃SH, CH₃CHO and CHCl₃ liquids, based on anisotropic Gay-Berne (GB) and electric multipole potentials (EMP), demonstrate that the coarse-grained model is able to qualitatively reproduce the results obtained from the atomistic model (AMOEBA polarizable force field) and allows for significant saving in computation time. It should be pointed out that the accuracy of the coarse-grained model is very sensitive to how well the anisotropic GB particle is defined and how satisfactorily the EMP sites are chosen.

Mitochondrial complex I inhibitors, acetogenins, induce HepG2 cell death through the induction of the complete apoptotic mitochondrial pathway

The development of new anti-neoplastic drugs is a key issue for cancer chemotherapy due to the reality that, most likely, certain cancer cells are resistant to current chemotherapy. The past two decades have witnessed tremendous advances in our understanding of the pathogenesis of cancer. These advances have allowed identification new targets as oncogenes, tumor supressor genes and the possible implication of the mitochondria (Fulda et al. Nat Rev Drug Discov 9:447–464, 2010). Annonaceous Acetogenins (ACGs) have been described as the most potent inhibitors of the respiratory chain because of their interaction with mitochondrial Complex I (Degli Esposti and Ghelli Biochim Biophys Acta 1187:116–120, 1994; Zafra-Polo et al. Phytochemistry 42:253–271, 1996; Miyoshi et al. Biochim Biophys Acta 1365:443–452, 1998; Tormo et al. Arch Biochem Biophys 369:119–126, 1999; Motoyama et al. Bioorg Med Chem Lett 12:2089–2092, 2002). To explore a possible application of natural products from Annonaceous plants to cancer treatment, we have selected four bis-tetrahydrofuranic ACGs, three from Annona cherimolia (cherimolin-1, motrilin and laherradurin) and one from Rollinia mucosa (rollinianstatin-1) in order to fully describe their mechanisms responsible within the cell (Fig. 1). In this study, using a hepato-carcinoma cell line (HepG2) as a model, we showed that the bis-THF ACGs caused cell death through the induction of the apoptotic mitochondrial pathway. Their potency and behavior were compared with the classical mitochondrial respiratory chain Complex I inhibitor rotenone in every apoptotic pathway step.

Cubitane: a rare diterpenoid skeleton

There are only 15 diterpenoids known sharing the rare cubitane skeleton (1,3-diisopropyl-6,10-dimethylcyclododecane, Fig. 1) named after the natural product (+)-cubitene ((+)-1) from the termite Cubitermes umbratus. Cubitane-type diterpenoids (“cubitanoids”) have since then been isolated from gorgonian corals and there are reports that cubitene also occurs in plants. The twelve-membered ring of the cubitanoids is an interesting feature for which we have developed a novel synthesis involving bicyclic precursors. In this review, we discuss the state-of-research regarding isolation, structure elucidation, biological activity, biosynthesis, and total synthesis of cubitane-type diterpenoids.

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0–3,500 cm^?1 of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.