Glass Transition Temperatures Studied by MD Simulation of Some Alkali Metasilicates

Abstract

Molecular dynamics simulation of some alkali metasilicates (M₂SiO₃, M = Li, Na, K) was performed to compare glass transition temperatures, T_g , defined in various ways. The potential parameters derived from ab initio MO calculations were used and found to reproduce the inflection of V-T relation on cooling the system. The T_g defined by the inflection point corresponds well to that defined by geometrical changes of coordination polyhedra found in previous work. The self-diffusion coefficients of the alkali ions in higher temperature regions were shown to be related to the amount of free volume in these systems.     

Molecular Dynamics Simulation Studies of the Density and Temperature Dependence of Self-diffusion in a Cylindrical Micropore

We report Molecular Dynamics calculations of radial density profiles and self-diffusion coefficients of Lennard-Jones fluids in a cylindrical pore of radius 2σ, for a wide range of temperatures and densities. At n _p σ³ = 0.825 the self-diffusion coefficient parallel to the pore walls D _∥*. follows a monotonic (nearly linear) increase with kT/ε and is very similar to that of the bulk self-diffusion coefficient D _b *. At n _p σ³ = 0.4 and kT/ε ≤ 1.0 the curve of D _∥* vs. kT/ε shows a distinct inflection in the region 0.7 ≤ kT/ε ≤ 0.9 and values of D _∥* are much less than D _b * decreasing to near solid state values at very low temperatures. At the highest temperature studied, kT/ε = 2.98, D _∥* is almost inversely proportional to density and in a fairly close agreement with that of D _b *. At KT/ε = 0.49, D _∥* is much smaller than D _b *. The motion of adsorbate particles normal to the walls is also discussed.     

A model of cell firing patterns during epileptic seizures

A simplified model is presented of the dynamics of excitatory and inhibitory neurons in the cerebral cortex. A key feature of the model is that neurons may cease to fire when strongly depolarized (spike inactivation). Computer simulations for different parameters reveal five classes of solutons: a) steady states in which neither excitatory nor inhibitory cells are active, b) steady states in which one or both types of cells fire repetitively, c) states in which one type of cell fluctuates rapidly between bursts of action potentials and inactivity due to strong depolarization, d) rhythmic activity in which both types of cells fire in unison followed by a period of spike inactivation and e) states similar to d but in which the inhibitory cells never produce action potentials. Solutions b, c, d, and e qualitatively resemble the different firing patterns observed during experimental seizures. It is shown that changes in those parameters that are functions of potassium concentration can induce changes in the type of solution. It is therefore proposed that the increase in extracellular potassium concentration during seizures may be responsible for the progressive changes observed in firing patterns and particularly for the transition from tonic to clonic patterns. A method is also outlined for testing the predictions of the model.     

Brownian Dynamics Simulations of Colloidal Liquids: Hydrodynamics and Stress Relaxation

Abstract

We describe the statistical mechanics background and additional algorithmic features of a recently proposed simple mean-field Brownian Dynamics algorithm formulated to include many-body hydrodynamics, using a local density approximation for the friction coefficient. We show that the equations of motion satisfy the incompressibility of phase space. We make further developments to the model, computing the hydrodynamic effects on the shear stress relaxation function. We show that stress relaxation takes place over two well-defined regimes, in both cases with and without mean field hydrodynamics, MFH. At short times ta ²/D ₀ < 10^?3, where a is the radius of the colloidal particle and D ₀ is the self-diffusion coefficient at infinite dilution, decay of the stress autocorrelation function, C_s(t) is essentially independent of volume fraction and does not fit to a simple analytic form. At longer times than ta ²/D ₀ < 10^?2 the decay has the fractional exponential form ～exp(-t ^β) with β ? 1. The transition between these two regimes coincides with a rapid fall in the time-dependent diffusion coefficient from the so-called short-time to long-time values. We do not find any evidence for power law decay in the C_s(t) as predicted by recent mode-coupling based analytical expansions.     

Static and dynamic properties of mid-size liquid n- alkanes,C12∼C400: a molecular dynamics simulation study

ABSTRACT

In this paper, we have extended our previous study of the static and dynamic properties (self-diffusion coefficient D_self and friction coefficient ζ) of liquid n-alkane systems up C₄₀₀ at several temperatures (～2300?K) using molecular dynamics (MD) simulations in the canonical ensembles. For the small n-alkanes with n?≤?120 (n: the chain length), the chains are clearly ?R² _ee?/6?R² _g? ≥ 1 (1.06 ～ 1.44), which leads to the conclusion that the liquid n-alkanes are far away from the ideal chain regime. But for the n-alkanes of n?≥?160, the chains are ?R² _ee?/6?R² _g? ≈ 1, indicating that they are Gaussian. It is found that the long chains of these n-alkanes at high temperatures show abnormalities in density and friction coefficient. We observed a clear transition in the power law dependence of n-alkane self-diffusion coefficient on the molecular weight (M) of n-alkane, D_self ～ M^?γ, occurs in the range C₁₂₀～C₁₆₀ at temperatures of 318, and 618?K, corresponding to a crossover from the ‘oligomer’ to the ‘Rouse’ regime. The entanglement lengths (N_e) are calculated by the Z1 code and discussed shortly.     

Molecular modelling insights into a physiologically favourable approach to eicosanoid biosynthesis inhibition through novel thieno[2,3-b]pyridine derivatives

In this research, we exploited derivatives of thieno[2,3-b]pyridine as dual inhibitors of the key enzymes in eicosanoid biosynthesis, cyclooxygenase (COX, subtypes 1 and 2) and 5-lipoxygensase (5-LOX). Testing these compounds in a rat paw oedema model revealed potency higher than ibuprofen. The most active compounds 7a, 7b, 8b, and 8c were screened against COX-1/2 and 5-LOX enzymes. Compound 7a was the most powerful inhibitor of 5-LOX with IC₅₀?=?0.15?µM, while its p-chloro analogue 7b was more active against COX-2 (IC₅₀?=?7.5?µM). The less desirable target COX-1 was inhibited more potently by 8c with IC₅₀?=?7.7?µM. Surflex docking programme predicted that the more stable anti- conformer of compound (7a) formed a favourable complex with the active site of 5-LOX but not COX-1. This is in contrast to the binding mode of 8c, which resembles the syn-conformer of series 7 and binds favourably to COX-1.     

Synthesis of “Reversed” Methylenecyclopropane Analogues of Antiviral Phosphonates

Synthesis of “reversed” methylenecyclopropane analogues of nucleoside phosphonates 6a,7a, 6b, and 7b is described. 1-Bromo-1-bromomethylcyclopropane 8 was converted to the bromocyclopropyl phosphonate 9 by Michaelis-Arbuzov reaction with triisopropyl phosphite. Base-catalyzed β-elimination and deacetylation gave the key Z- and E-hydroxymethylcyclopropyl phosphonates 10 and 11 separated by chromatography. The Mitsunobu type of alkylation of 10 or 11 with adenine or 2-amino-6-chloropurine afforded phosphonates 12a, 12b, 13a, and 13b. Acid hydrolysis furnished the adenine and guanine analogues 6a, 7a, 6b, and 7b. The E and Z configuration was assigned on the basis of NOE experiments with phosphonates 6b and 7b. All Z- and E-isomers were also distinguished by different chemical shifts of CH₂O or CH₂N (H₄ or H_4′). Significant differences of the chemical shifts of the cyclopropane C_3(3’) carbons and coupling constants ³J_P,C2(2’) or ³J_P,C3(3’) selective for the Z- or E-isomers were also noted. Phosphonates 6a, 7a, 6b, and 7b are devoid of significant antiviral activity.     

Design,synthesis, in vitro and in vivo evaluation of benzylpiperidine-linked 1,3-dimethylbenzimidazolinones as cholinesterase inhibitors against Alzheimer’s disease

Abstract

Cholinesterase inhibitor plays an important role in the treatment of patients with Alzheimer’s disease (AD). Herein, we report the medicinal chemistry efforts leading to a new series of 1,3-dimethylbenzimidazolinone derivatives. Among the synthesised compounds, 15b and 15j showed submicromolar IC₅₀ values (15b, eeAChE IC₅₀?=?0.39?±?0.11?µM; 15j, eqBChE IC₅₀?=?0.16?±?0.04?µM) towards acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Kinetic and molecular modelling studies revealed that 15b and 15j act in a competitive manner. 15b and 15j showed neuroprotective effect against H₂O₂-induced oxidative damage on PC12 cells. This effect was further supported by their antioxidant activity determined in a DPPH assay in vitro. Morris water maze test confirmed the memory amelioration effect of the two compounds in a scopolamine-induced mouse model. Moreover, the hepatotoxicity of 15b and 15j was lower than tacrine. In summary, these data suggest 15b and 15j are promising multifunctional agents against AD.     

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (μ _s P _s)/(μ _v P _v), where μ _s and μ _v denote mutation rate of transitions and transversions, respectively, andP _s andP _v denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions,P _s/P_v ≈ 1 at fourfold degenerate sites, so that thes/v ratio at fourfold degenerate sites is approximately μ _s /μ _v , which is a measure of mutational contribution to transition bias. Similarly, thes/v ratio at nondegenerate sites is also an estimate of μ _s /μ _v if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochromeb (cyt-b) in pocket gophers, thes/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, thes/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the μ _s /μ _v ratio is about 2, theP _s/P_v ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content ons/v ratio of protein-coding genes. Correspondence to: X. Xia     

Synthesis,antitumour activities and molecular docking of thiocarboxylic acid ester-based NSAID scaffolds: COX-2 inhibition and mechanistic studies

A new series of NSAID thioesters were synthesized and evaluated for their in vitro antitumor effects against a panel of four human tumor cell lines, namely: HepG2, MCF-7, HCT-116 and Caco-2, using the MTT assay. Compared to the reference drugs 5-FU, afatinib and celecoxib, compounds 2b, 3b, 6a, 7a, 7b and 8a showed potent broad-spectrum antitumor activity against the selected tumour cell lines. Accordingly, these compounds were selected for mechanistic studies about COX inhibition and kinase assays. In vitro COX-1/COX-2 enzyme inhibition assay results indicated that compounds 2b, 3b, 6a, 7a, 7b, 8a and 8?b selectively inhibited the COX-2 enzyme (IC₅₀?=?～0.20–0.69?μM), with SI values of (>72.5–250) compared with celecoxib (IC₅₀?=?0.16?μM, COX-2 SI:?>?312.5); however, all the tested compounds did not inhibit the COX-1 enzyme (IC₅₀?>?50?μM). On the other hand, EGFR, HER2, HER4 and cSrc kinase inhibition assays were evaluated at a 10?μM concentration. The selected candidates displayed limited activities against the various tested kinases; the compounds 2a, 3b, 6a, 7a, 7b and 8a showed no activity to weak activity (% inhibition?=?～0–10%). The molecular docking study revealed the importance of the thioester moiety for the interaction of the drugs with the amino acids in the active sites of COX-2. The aforementioned results indicated that thioester based on NSAID scaffolds derivatives may serve as new antitumor compounds.     

Synthesis,antitumour and antioxidant activities of novel α,β-unsaturated ketones and related heterocyclic analogues: EGFR inhibition and molecular modelling study

New α,β-unsaturated ketones 4a,b; 5a–c; and 6a,b; as well as 4-H pyran 7; pyrazoline 8a,b; isoxazoline 9; pyridine 10–11; and quinoline-4-carboxylic acid 12a,b derivatives were synthesized and evaluated for in vitro antitumour activity against HepG2, MCF-7, HeLa, and PC-3 cancer cell lines. Antioxidant activity was investigated by the ability of these compounds to scavenge the 2,2′-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) radical cation (ABTS^?+). Compounds 6a, 6b, 7, and 8b exhibited potent antitumour activities against all tested cell lines with [IC₅₀] ?5.5–18.1 µΜ), in addition to significantly high ABTS^?+ scavenging activities. In vitro EGFR kinase assay for 6a, 6b, 7, and 8b as the most potent antitumour compounds showed that; compounds 6b, and 7 exhibited worthy EGFR inhibition activity with IC₅₀ values of 0.56 and 1.6?µM, respectively, while compounds 6a and 8b showed good inhibition activity with IC₅₀ values of 4.66 and 2.16?µM, respectively, compared with sorafenib reference drug (IC₅₀?=?1.28?µM). Molecular modelling studies for compounds 6b, 7, and 8b were conducted to exhibit the binding mode towards EGFR kinase, which showed similar interaction with erlotinib.     

Synthesis and potential antitumor activity of 7-(4-substituted piperazin-1-yl)-4-oxoquinolines based on ciprofloxacin and norfloxacin scaffolds: in silico studies

The potential antitumor activities of a series of 7-(4-substituted piperazin-1-yl)fluoroquinolone derivatives (1–14a,b) using ciprofloxacin and norfloxacin as scaffolds are described. These compounds exhibit potent and broad spectrum antitumor activities using 60 human cell lines in addition to the inherent antibacterial activity. Compounds 1a, 2a, 3b, 6b and 7a were found to be the most potent, while 2b, 5b, and 6a were found to have an average activity. The results of this study demonstrated that compounds 1a, 2a, 3b, 6b and 7a (mean GI₅₀; 2.63–3.09?µM) are nearly 7-fold more potent compared with the positive control 5-fluorouracil (mean GI₅₀; 22.60?µM). More interestingly, compounds 1a, 2a, 3b, 6b and 7a have an almost antitumor activity similar to gefitinib (mean GI₅₀; 3.24?µM) and are nearly 2-fold more potent compared to erlotinib (mean GI₅₀; 7.29?µM). In silico study and ADME-Tox prediction methodology were used to study the antitumor activity of the most active compounds and to identify the structural features required for antitumor activity.     

The Codon-Degeneracy Model of Molecular Evolution

Mitochondrial genetic codons can be categorized by four patterns of nucleotide-site degeneracy based on varying combinations of twofold- or nondegenerate sites at first codon positions and twofold- or fourfold-degenerate sites at third codon positions. Herein, a model of molecular evolution is introduced that uses these patterns to calculate expected substitution frequencies for each codon position and substitution type relative to overall number of synonymous or nonsynonymous substitutions. Regions of the pocket gopher cytochrome oxidase subunit I (COI) and cytochrome b (cyt-b) genes are analyzed using this model. Chi-square distributions are used to produce relative goodness-of-fit (GF) scores for measuring the difference between substitution frequencies predicted by the codon-degeneracy model (CDM), and frequencies inferred using a well-supported phylogenetic tree of closely related species. The GF scores for expected and observed synonymous (GF_syn= 0.429, p= 0.807) and nonsynonymous (GF_ns= 2.309, p= 0.679) substitution frequencies resulted in a failure to reject the CDM as a null hypothesis for the molecular evolution of COI and cyt-b in pocket gophers. Alternative tree topologies and calculations of transition bias for these data result in higher GF scores. Received: 25 March 1999 / Accepted: 17 September 1999     

An NMR self-diffusion study of the interactions between spermidine and oligonucleotides

Self-diffusion coefficients have been determined by pulsed field gradient nmr methods for spermidine in solutions of the oligonucleotides d(GC)₄ and d(GGAATTCC). The self-diffusion behavior of spermidine in solution of d(GC)₄ is very similar to that observed previously for methylspermidine (completely N-methylated spermidine). Moreover, the self-diffusion behaviors of spermidine in solutions of d(GC)₄ and d(GGAATTCC) are also quite similar, indicating that there is no significant influenceon on self-diffusion of oligonucleotide base composition. Furthermore, self-diffusion coefficients of the oligonucleotide d(GC)₈ show only a small dependence on oligonucleotide concentration, and no measurable dependence on sodium ion or magnesium ion concentration. © 1996 John Wiley & Sons, Inc.     

Stability of the analytical solution of Penna model of biological aging

There are some analytical solutions of the Penna model of biological aging; here, we discuss the approach by Coe et al. (Phys. Rev. Lett. 89, 288103, 2002), based on the concept of self-consistent solution of a master equation representing the Penna model. The equation describes transition of the population distribution at time t to next time step (t + 1). For the steady state, the population n(a, l, t) at age a and for given genome length l becomes time-independent. In this paper we discuss the stability of the analytical solution at various ranges of the model parameters—the birth rate b or mutation rate m. The map for the transition from n(a, l, t) to the next time step population distribution n(a + 1, l, t + 1) is constructed. Then the fix point (the steady state solution) brings recovery of Coe et al. results. From the analysis of the stability matrix, the Lyapunov coefficients, indicative of the stability of the solutions, are extracted. The results lead to phase diagram of the stable solutions in the space of model parameters (b, m, h), where h is the hunt rate. With increasing birth rate b, we observe critical b ₀ below which population is extinct, followed by non-zero stable single solution. Further increase in b leads to typical series of bifurcations with the cycle doubling until the chaos is reached at some b _c. Limiting cases such as those leading to the logistic model are also discussed.     

Measurements of the gas temperature in an oxygen plasma by spectroscopy of the O2( b 1Σ g+ ) → O2( X 3Σ g? ) transition

A method for measuring the gas temperature in an oxygen plasma by spectroscopy of the electronic transition from the O₂(b ¹Σ _g ⁺ , v = 0) metastable state of molecular oxygen into the O₂(X ³Σ _g ⁻ , v = 0) ground state is considered in detail. The method is verified experimentally for the plasma of dc glow discharge in pure oxygen. It is shown that the gas temperature can be determined by analyzing high-resolution spectra of the P branch of this transition, no matter whether its fine structure (^P P and ^P Q branches) is resolved or masked, provided that the rotational structure of the spectrum is resolved. The feasibility of the method proposed in 1999 by P. Maco and P. Veis for determining the gas temperature from the ratio between the intensity maxima of the R and P branches of the O₂(b ¹Σ _g ⁺ , v = 0) → O₂(X ³Σ _g ⁻ , v = 0) transition in a poorly resolved spectrum was studied experimentally. It is shown that, in order to use this method, it is necessary to know the spectrograph instrumental function. The effect of the spatial inhomogeneity of the temperature and concentration of O₂(b ¹Σ _g ⁺ ) molecules on the accuracy of integral (over the plasma volume) measurements of the gas temperature is investigated using spatially resolved spectroscopy of the O₂(b ¹Σ _g ⁺ , v = 0) → O₂(X ³Σ _g ⁻ , v = 0) transition. It is shown that precise measurements of the temperature require that the optical measurement system be thoroughly adjusted in order for the temperature and concentration of the emitting particles to vary insignificantly over the optically selected volume. Original Russian Text ? S.M. Zyryanov, D.V. Lopaev, 2007, published in Fizika Plazmy, 2007, Vol. 33, No. 6, pp. 563–574.     

Novel 1,5-diphenyl-6-substituted 1H-pyrazolo[3,4-d]pyrimidin-4(5H)-ones induced apoptosis in RKO colon cancer cells

Novel 1,5-diphenyl-6-substituted-1H-pyrazolo[3,4-d]pyrimidin-4(5H)-ones were synthesized and characterized. All compounds were screened for their anti-proliferative activities in five different cancer cell lines. The results showed that compounds 7a and 7b comprising aminoguanidino or guanidino moiety at position 6 inhibited proliferation of RKO colon cancer cells with IC₅₀ of 8 and 4?μM, respectively. Compounds 7a and 7b induced apoptosis in RKO cells, which was confirmed by TUNEL and annexin V-FITC assays. Flow cytometric analysis indicated that compounds 7a and 7b arrested RKO cells in the G1 phase and the most active compound 7b increased levels of p53, p21, Bax, ERK1/2 and reduced levels of Bcl2 and Akt. Compound 7b also activates release of cytochrome c, which is consistent with activation of caspase-9. Additionally, compound 7b increased caspase-3 activity and cleaved PARP-1 in RKO cells. Collectively, these findings could establish a molecular basis for the development of new anti-cancer agents.     

Fixation of beneficial mutations in the presence of epistatic interactions

We investigate the effect of deleterious mutations on the process of fixation of new advantageous mutants in an asexual population. In particular we wish to study the dependence of the process on the strength of the deleterious mutations. We suppose the existence of epistatic interaction between the genes. We study the model by means of branching process theory and also by numerical simulations. Our results show the occurrence of two distinct regimes of behavior for the probability of fixation of these variants. The occurrence of either regime depends on the ratio between the selective advantage of the beneficial mutation s _b and on the selective parameter for deleterious mutations s _b . In the former, which takes place for s _b /s _d ≲ 1, the probability of fixation increases with the epistasis parameter α, whereas for s _b /s _d ≫ 1 the probability of fixation is a complex function of α and the mutation rate U. Surprisingly, we find that for the multiplicative landscape (α = 1) the probability of fixation P _fix is given by where π (s _b ) is the probability of fixation for the two-allele model in the absence of mutations as calculated by Haldane (1927, Proc. Camb. Phil. Soc., 26, 220–230) and Kimura (1962, Genetics, 47, 713–719).     

3D pharmacophore models for thromboxane A2 receptor antagonists

Thromboxane A₂ (TXA₂) is an endogenous arachidonic acid derivative closely correlated to thrombosis and other cardiovascular diseases. The action of TXA₂ can be effectively inhibited with TXA₂ receptor antagonists (TXRAs). Previous studies have attempted to describe the interactions between the TXA₂ receptor and its ligands, but their conclusions are still controversial. In this study, ligand-based computational drug design is used as a new and effective way to investigate the structure-activity relationship of TXRAs. Three-dimensional pharmacophore models of TXRAs were built with HypoGenRefine and HipHop modules in CATALYST software. The optimal HypoGenRefine model was developed on the basis of 25 TXRAs. It consists of two hydrophobic groups, one aromatic ring, one hydrogen-bond acceptor and four excluded volumes. The optimal HipHop model contains two hydrophobic groups and two hydrogen-bond acceptors. These models describe the key structure-activity relationship of TXRAs, can predict their activities, and can thus be used to design novel antagonists. Figure Optimal three-dimensional pharmacophore models of TXA₂ receptor antagonists (TXRAs) built with HypoGenRefine (a) and HipHop (b) modules. a Hypo-1 model mapped with compounds 11 (purple), and 20 (green). b Hypo-2 model mapped with compounds 31 (green) and 64 (yellow). Spheres: Green Hydrogen bond acceptors (HBA), cyanhydrophobic (H), orange aromatic rings (RA), black excluded volumes (EV) Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Crystal and molecular structure of V-anhydrous amylose

The conformation and crystalline packing of V-anhydrous amylose has been investigated by a combination of linked atom model building and X-ray diffraction analysis. The unit cell, the P2₁2₁2₁ space group, the left-handed sixfold helical conformation with all O(6) in gt rotational positions, and the intrahelical O(2)---O(3) and O(2)---O(6) hydrogen bonds are substantially in agreement with previous studies. A new model for packing of the chains in the unit cell and the presence of crystallographic water is proposed. Packing appears to be stabilized by corner-to-center chain O(2)---O(2) hydrogen bonds. The nature of the transition from the amylose–DMSO complex to V_a-amylose was considered and it is shown that the transition involves translation of the amylose chains parallel to the a and b unit cell axes with only slight changes in the orientation of the helix. No significant conformational changes result from the transition.