Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

Molecular Dynamics of Anthraquinone DNA Intercalators with Polyethylene Glycol Side Chains

Abstract

The interactions of a homologous series of four anthraquinone (AQ) intercalators with increasing lengths of polyethylene glycol (PEG) side chains with DNA have been studied via molecular dynamics (MD) simulations. The geometry, conformation, interactions, and hydration of the complexes were examined. The geometries of the four ligands were similar with parallel stacking to the long axis of the adjacent DNA base pairs. Hydrogen bonding between the AQ amide and DNA led to a preference for the trans-syn conformer. As the side chain lengthened, binding to DNA reduced the conformational space, resulting in an increase in unfavorable entropy. Increased localization of the PEG side chain in the DNA groove, indicating some interaction of the side chain with DNA, also contributed unfavorably to the entropy. The changes in free energy of binding due to entropic considerations (— 3.9 to—6.3 kcal/mol) of AQ I-IV were significant. The hydration of the PEG side chain decreased upon binding to DNA. Understanding of side chain conformations, interactions, and hydration changes that accompany the formation of a ligand-DNA complex may be important in the development of new applications of pegylated small molecules that target biological macromolecules.     

Molecular dynamics (MD) simulations and binding free energy calculation studies between inhibitors and type II dehydroquinase (DHQ2)

Type II dehydroquinase (DHQ2) is the third enzyme of the shikimic acid pathway, and it has been the effective target for tuberculosis (TB). So far, developing multiple potent inhibitors of the DHQ2 of Mycobacterium tuberculosis (DHQ2-Mt) has been considered to be the new therapy to TB. Molecular dynamics simulations followed by molecular mechanics-generalised Born surface area were carried out to calculate the free binding energy and to determine the affinity ability of the four chosen inhibitor molecules, L1, L2, L3 and L4. Energy decomposition analyses show the electrostatic interaction and van der Waals interaction of the ligands to every residue of the DHQ2-Mt. The results suggest that some important residues have different interactions with the four ligands, such as Arg19 and Tyr24. These interactions may have an effect on the ligand binding affinity. The binding affinity of monosubstituted inhibitor is higher than that of disubstituted inhibitor, due to some important interactions with the DHQ2-Mt residues. These computational works will be significant to the theoretical research in the future.     

Obesity genes and the regulation of body fat content

Physiological investigation has demonstrated that the central nervous system monitors body composition and adjusts energy intake and expenditure to stabilize total adipose tissue mass. Genetic variations in the signalling molecules involved in this regulatory system account for the heritable component of body fat content. The application of molecular techniques to rodent models of Mendelian obesity has resulted in the characterization of five loci at which mutations produce an abnormal accumulation of body fat. The genes at these loci include agouti, which encodes a molecule that antagonizes the binding of alpha melanocyte-stimulating hormone to its receptor; fat, which encodes carboxypeptidase E; tubby, which encodes a putative phosphodiesterase; obese, which encodes a circulating satiety protein; and diabetes, which encodes the receptor for the obese gene product. A more detailed understanding of the functional interrelationships of these genes should lead to important new insights into the causes and potential therapies for human obesity.     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Molecular Modelling of The Mechanism of Action of Organic Clay-Swelling Inhibitors

Abstract

It is well known that the sodium smectite class of clays swells macroscopically in contact with water, whereas under normal conditions the potassium form does not. In recent work using molecular simulation methods, we have provided a quantitative explanation both for the swelling behaviour of sodium smectite clays and the lack of swelling of potassium smectites [1]. In the present paper, we apply similar modelling methods to study the mechanism of inhibition of clay-swelling by a range of organic molecules.

Experimentally, it is known that polyalkylene glycols (polyethers) of intermediate to high relative molecular mass are effective inhibitors of smectite clay swelling. We use a range of atomistic simulation techniques, including Monte Carlo and molecular dynamics, to investigate the interactions between a selection of these compounds, water, and a model smectite clay mineral. These interactions occur by means of organised intercalation of water and organic molecules within the galleries between individual clay layers.

The atomic interaction potentials deployed in this work are not as highly optimised as those used in our clay-cation-water work [1]. Nevertheless, our simulations yield trends and results that are in qualitative and sometimes semi-quantitative agreement with experimental findings on similiar (but not identical) systems. The internal energy of adsorption of simple polyethers per unit mass on the model clay is not significantly different from that for water adsorption; our Monte Carlo studies indicate that entropy is the driving force for the sorption of the simpler organic molecules inside the clay layers: a single long chain polyethylene glycol can displace a large number of water molecules, each of whose translational entropy is greatly enhanced when outside the clay. Hydrophobically modified polyalkylene glycols also enjoy significant van der Waals interactions within the layers which they form within the clay galleries.

In conjunction with experimental studies, our work furnishes valuable insights into the relative effectiveness of the compounds considered and reveals the generic features that high performance clay-swelling inhibitors should possess. For optimal inhibitory activity, these compounds should be reasonably long chain linear organic molecules with localised hydrophobic and hydrophilic regions along the chain. On intercalation of these molecules within the clay layers, the hydrophobic regions provide an effective seal against ingress of water, while the hydrophilic ones enhance the binding of the sodium cations to the clay surface, preventing their hydration and the ensuing clay swelling.     

An Extended Gibbs Ensemble

Abstract

A general extended Gibbs ensemble, obtained by augmenting the standard Gibbs ensemble by intermediate states in the spirit of the scaled particle method of Nezbeda and Kolafa [Molec. Simul., 5, 391 (1991)], is introduced. The intermediate states span the states with different number of particles in the simulation boxes and facilitate the transfer of particles even in such complex systems as e.g., mixtures of very different components, systems of flexible polymeric molecules, or systems at very high densities. A general formulation of the ensemble is given and two implementations are considered in detail. The method is very general and is exemplified by studying the fluid-fluid coexistence in a dense binary mixture of the hard-sphere and square-well fluids. It is found that its efficiency is about by factor three greater in comparison with the standard Gibbs ensemble simulations.     

Molecular Modelling of the Antagonist Binding Site on the Adenosine A,Receptor

Abstract

With the aid of molecular modelling both adenosine and adenosine A, receptor antagonists belonging to various chemical classes were compared. A model for the antagonist binding site was developed. As a consequence 1H-imidazo[4, 5-c]-quinolin-4-amines were synthesized, constituting a novel class of potent non-xanthine adenosine receptor antagonists.     

Tethered Molecular Sorbents: Enabling Metal‐Sulfur Battery Cathodes

A rechargeable battery that uses sulfur at the cathode and a metal (e.g., Li, Na, Mg, or Al) at the anode provides perhaps the most promising path to a solid‐state, rechargeable electrochemical storage device capable of high charge storage capacity. It is understood that solubilization in the electrolyte and loss of sulfur in the form of long‐chain lithium polysulfides (Li₂S_x, 2 < x < 8) has hindered development of the most studied of these devices, the rechargeable Li‐S battery. Beginning with density‐functional calculations of the structure and interactions of a generic lithium polysulfide species with nitrile containing molecules, it is shown that it is possible to design nitrile‐rich molecular sorbents that anchor to other components in a sulfur cathode and which exert high‐enough binding affinity to Li₂S_x to limit its loss to the electrolyte. It is found that sorbents based on amines and imidazolium chloride present barriers to dissolution of long‐chain Li₂S_x and that introduction of as little as 2 wt% of these molecules to a physical sulfur‐carbon blend leads to Li‐S battery cathodes that exhibit stable long‐term cycling behaviors at high and low charge/discharge rates.     

A CHARMm Based Force Field for Carbohydrates Using the CHEAT Approach: Carbohydrate Hydroxyl Groups Represented by Extended Atoms

Abstract

Computational studies of carbohydrates that do not consider explicit solvent molecules suffer from the strong tendency of the carbohydrate pendant hydroxyl groups to form intramolecular hydrogen bonds that are unlikely to be present in protic media. In this paper a novel approach towards molecular modelling of carbohydrates is described. The average effect of intra- and intermolecular hydrogen bonding is introduced into the potential energy function by adding a new (extended) atom type representing a carbohydrate hydroxyl group to the CHARMm force field; we coin the name CHEAT (Carbohydrate Hydroxyls represented by Extended AToms) for the resulting force field. As a training set for the parametrisation of CHEAT we used ethylene glycol, 10 cyclohexanols, 5 inositols, and 12 glycopyranoses for which in total 64 conformational energy differences were estimated using a set of steric interaction energies between hydroxyl and/or methyl groups on six-membered ring compounds as derived by Angyal (Angew. Chem., 8, 172-182, (1969)). The root-mean-square deviation between the estimated energy differences and the corresponding values obtained by our CHEAT approach amounts to 0.37 kcal/mol (n = 64). The CHEAT approach, which is claimed to calculate aqueous state compatible energetical and conformational properties of carbohydrates, is computationally very efficient and facilitates the calculation of nanosecond range MD trajectories as well as systematic conformational searches of oligosaccharides.     

Molecular Mechanics Calculations of Proteins. Comparison of Different Energy Minimization Strategies

Abstract

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant ? = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.     

Sphingomyelins of Rat Liver: Biliary Enrichment with Molecular Species Containing 16:0 Fatty Acids as Compared to Canalicular-Enriched Plasma Membranes

We harvested canalicular-enriched plasma membranes of hepatocytes and collected fistula bile from male rats and isolated the sphingomyelins. Following sphingomyelinase hydrolysis, we identified the sphingomyelin molecular species on the basis of their benzoylated ceramide derivatives employing high performance liquid chromatography. Sphingomyelin constitutes ≤3% of total biliary phospholipids (which are mostly sn-1 16:0 long-chain phosphatidylcholines) and approximately 30% of canalicular-enriched membranes. In both cases, the principal molecular species were composed of 16:0, 18:0, 20:0, 22:0, 23:0, 24:0, 24:1 and 24:2 fatty acid classes. However, the 16:0 fatty acid species was enriched in biliary sphingomyelin to a significantly greater degree than in sphingomyelins of canalicular-enriched plasma membranes (46% vs. 25% of total). We argue a physical-chemical case for laterally separated domains of very long chain sphingomyelins on the exoplasmic leaflet of the canalicular membrane. We bolster our hypothesis by the likelihood that the least hydrophobic, e.g., 16:0 sphingomyelin molecular species, are miscible with biliary phosphatidylcholines, and are secreted into bile. Laterally separated domains of very long chain sphingomyelins on the exoplasmic leaflet of the canalicular membrane could provide a means of sequestering cholesterol molecules prior to secretion into bile. Received: 19 March 1998/Revised: 8 October 1998     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Composition of Woody Species in a Dynamicforest–woodland–savannah Mosaic in Uganda: Implications for Conservation and Management

Forest–woodland–savannah mosaics are a common feature in the East African landscape. For the conservation of the woody species that occur in such landscapes, the species patterns and the factors that maintain it need to be understood. We studied the woody species distribution in a forest–woodland–savannah mosaic in Budongo Forest Reserve, Uganda. The existing vegetation gradients were analyzed using data from a total of 591 plots of 400 or 500 m² each. Remotely sensed data was used to explore current vegetation cover and the gradients there in for the whole area. A clear species gradient exists in the study area ranging from forest, where there is least disturbance, to wooded grassland, where frequent fire disturbance occurs. Most species are not limited to a specific part of the gradient although many show a maximum abundance at some point along the gradient. Fire and accessibility to the protected area were closely related to variation in species composition along the ordination axis with species like Cynometra alexandri and Uvariopsis congensis occurring at one end of the gradient and Combretum guenzi and Lonchocarpus laxiflorus at the other. The vegetation cover classes identified in the area differed in diversity, density and, especially, basal area. All vegetation cover classes, except open woodland, had indicator species. Diospyros abyssinica, Uvariopsis congensis, Holoptelea grandis and all Celtis species were the indicator species for the forest class, Terminalia velutina and Albizia grandbracteata for closed woodland, Grewia mollis and Combretum mole for very open woodland and Lonchocarpus laxiflorus, Grewia bicolor and Combretum guenzi for the wooded grassland class. Eleven of the species occurred in all cover classes and most of the species that occurred in more than one vegetation cover class showed peak abundance in a specific cover class. Species composition in the study area changes gradually from forest to savannah. Along the gradient, the cover classes are distinguishable in terms of species composition and vegetation structure. These classes are, however, interrelated in species composition. For conservation of the full range of the species within this East African landscape, the mosaic has to be managed as an integrated whole. Burning should be varied over the area with the forest not being burnt at all and the wooded grassland burnt regularly. The different vegetation types that occur between these two extremes should be maintained using a varied fire regime.     

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Insights into the Substrate Specificity of a Thioesterase Rv0098 of Mycobacterium Tuberculosis through X-ray Crystallographic and Molecular Dynamics Studies

Abstract

The crystal structure of Rv0098, a long-chain fatty acyl-CoA thioesterase from Mycobacterium tuberculosis with bound dodecanoic acid at the active site provided insights into the mode of substrate binding but did not reveal the structural basis of substrate specificities of varying chain length. Molecular dynamics studies demonstrated that certain residues of the substrate binding tunnel are flexible and thus modulate the length of the tunnel. The flexibility of the loop at the base of the tunnel was also found to be important for determining the length of the tunnel for accommodating appropriate substrates. A combination of crystallographic and molecular dynamics studies thus explained the structural basis of accommodating long chain substrates by Rv0098 of M. tuberculosis.     

Inherent Anti-amyloidogenic Activity of Human Immunoglobulin �� Heavy Chains

We have previously shown that a subpopulation of naturally occurring human IgGs were cross-reactive against conformational epitopes on pathologic aggregates of Aβ, a peptide that forms amyloid fibrils in the brains of patients with Alzheimer disease, inhibited amyloid fibril growth, and dissociated amyloid in vivo. Here, we describe similar anti-amyloidogenic activity that is a general property of free human Ig γ heavy chains. A γ₁ heavy chain, F1, had nanomolar binding to an amyloid fibril-related conformational epitope on synthetic oligomers and fibrils as well as on amyloid-laden tissue sections. F1 did not bind to native Aβ monomers, further indicating the conformational nature of its binding site. The inherent anti-amyloidogenic activity of Ig γ heavy chains was demonstrated by nanomolar amyloid fibril and oligomer binding by polyclonal and monoclonal human heavy chains that were isolated from inert or weakly reactive antibodies. Most importantly, the F1 heavy chain prevented in vitro fibril growth and reduced in vivo soluble Aβ oligomer-induced impairment of rodent hippocampal long term potentiation, a cellular mechanism of learning and memory. These findings demonstrate that free human Ig γ heavy chains comprise a novel class of molecules for developing potential therapeutics for Alzheimer disease and other amyloid disorders. Moreover, establishing the molecular basis for heavy chain-amyloidogenic conformer interactions should advance understanding on the types of interactions that these pathologic assemblies have with biological molecules.     

Impermeant potential-sensitive oxonol dyes: II. The dependence of the absorption signal on the length of alkyl substituents attached to the dye

Summary We have measured the potential-dependent light absorption changes of 43 impermeant oxonol dyes with an oxidized cholesterol bilayer lipid membrane system. The size of the signal is strongly dependent on the chain length of alkyl groups attached to the chromophore. Dye molecules with intermediate chain lengths give the largest signals. To better understand the dependence of the absorbance signal on alkyl chain length, a simple equilibrium thermodynamic analysis has been derived. The analysis uses the free energy of dye binding to the membrane and the on-off model (E.B. George et al.,J. Membrane Biol.,103:245–253, 1988a) for the potential-sensing mechanism. In this model, a population of dye molecules in nonpolar membrane binding sites is in a potential-dependent equilibrium with a second population of dye that resides in an unstirred layer adjacent to the membrane. Dye in the unstirred layer is in a separate equilibrium with dye in the bulk bathing solution. The equilibrium binding theory predicts a sigmoidally shaped increase in signal with increasing alkyl chain length, even for very nonpolar dyes. We suggest that aggregation of the more hydrophobic dyes in the membrane bathing solution may be responsible for their low signals, which are not predicted by the theory.     

<i>Insilico</i> Anticancer Peptide Prediction from <i>Curcuma longa</i>

Cancer is the second biggest cause of death globally, and the use of therapeutic peptides to specifically target and destroy cancer cells has gotten much interest. Cancer peptides or vaccinations are utilized to treat cancer nowadays, apart from chemotherapy, which has significant discomfort, side effects and costly. It is time demanding to identify and predict potential anticancer peptides using computational biology approaches. Thus, 3-D molecular modeling is being used to find possible ACP candidates. In this research, Curcuma longa has predicted peptide sequences were docked on breast cancer receptors and used a molecular docking technique to assess the anticipated peptides’ binding affinities to MHC molecules. A similar approach was utilized to simulate the interactions of the chosen peptide with the TCR. Additionally, the Pep10 LIRQHVASNIGIAKSKIREPIV was examined, and our findings indicated interaction with MHC classes I and II. However, the maximum binding energy was obtained with TCR at 695.61, giving strength through eight hydrogen bonds. Similarly, the Pep20, GAIIGNRKIKLQPHIIIRID, the projected, has the most significant overall binding energy with MHC classes I and II but a lower global E total value with TCR, namely −600.97 kj/Mol, and also four hydrogen bonds. This research could lead to the development of novel anticancer drugs based on the anticancer activity of the Curcuma longa medicinal plant.