Myeloperoxidase G-463A polymorphism and susceptibility to coronary artery disease: A meta-analysis

Published data on the association between the myeloperoxidase (MPO) G-463A polymorphism and coronary artery disease (CAD) are inconclusive. To derive a more precise estimation of the relationship, a meta-analysis on this topic was performed. PubMed, EMBASE and Chinese national knowledge infrastructure were searched for studies regarding the association between the MPO G-463A polymorphism and CAD. A logistic regression analysis was used to estimate the genetic effect and the possible genetic model of action. Summary odds ratios (ORs) with their corresponding 95% confidence intervals (CIs) were calculated. There was strong evidence for an association between the MPO G-463A polymorphism and CAD. The genetic model of action was most likely to be co-dominant. Overall, the data showed that AA and GA genotypes were significantly associated with reduced risk of CAD (AA vs. GG: OR = 0.37, 95% CI = 0.17–0.78; GA vs. GG: OR = 0.73, 95% CI = 0.57–0.92). In subgroup analyses by study population and sources of controls, statistically significant results were observed in the Chinese population (AA vs. GG: OR = 0.21, 95% CI = 0.10–0.43; GA vs. GG: OR = 0.57, 95% CI =0.44–0.74) and in hospital-based control studies (AA vs. GG: OR = 0.20, 95% CI = 0.10–0.39; GA vs. GG: OR = 0.61, 95% CI = 0.48–0.77). This meta-analysis suggests that the MPO G-463A variant genotypes may be associated with decreased risk of CAD. However, given the limited number of studies and the potential biases, the influence of this polymorphism on CAD risk needs further investigation.     

Scaling of species diversity and body mass in mammals: Cope’s rule and the evolutionary cost of large size

Abstract

Equations are constructed describing the inverse correlation of species diversity and body mass in extant and Cenozoic mammals. Cope’s rule, the tendency for many mammal clades to increase in body size through time, through phyletic change in single lineages or turnover within species groups, is interpreted as a probability function reducing diversity potential as a tradeoff for ecological/evolutionary gains. The inverse rule predicts that large species in clades will be less diverse than smaller species and, unless origination rates remain high among smaller clade members, clades conforming to Cope’s rule will decline in diversity, moving towards extinction. This proposition is evaluated in the Cenozoic histories of five North American mammal clades; cotton rats, felids, canids, hyaenodontids, and equids. Diversity potential of different size classes within the 3.75 million year phyletic history of the muskrat, Ondatra zibethicus, is also examined. A corollary prediction of the inverse rule, that large species should have longer durations (species lifespans) than small species, is unresolved. Successful clades maintain small size or a significant number of smaller species relative to clade average size. The potential loss of unique extant large mammal species justifies the conservation effort to protect them. The similarity of scaling exponents of species diversity to mass around a slope of -1.0 suggests that species diversity is correlated with home range size, the latter related to the probability of population fragmentation.     

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity

A missense mutation I148M in PNPLA3 (patatin‐like phospholipase domain‐containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate‐access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate‐bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild‐type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A diterpenoid derivate compound targets selenocysteine of thioredoxin reductases and induces Bax/Bak-independent apoptosis

We have previously shown that the natural diterpenoid derivative S3 induced Bim upregulation and apoptosis in a Bax/Bak-independent manner. However, the exact molecular target(s) of S3 and the mechanism controlling Bim upregulation are still not clear. Here, we identify that S3 targets the selenoproteins TrxR1 and TrxR2 at the selenocysteine residue of the reactive center of the enzymes and inhibits their antioxidant activities. Consequently, cellular ROS is elevated, leading to the activation of FOXO3a, which contributes to Bim upregulation in Bax/Bak-deficient cells. Moreover, S3 retards tumor growth in subcutaneous xenograft tumors by inhibiting TrxR activity in vivo. Our studies delineate the signaling pathway controlling Bim upregulation, which results in Bax/Bak-independent apoptosis and provide evidence that the compounds can act as anticancer agents based on mammalian TrxRs inhibition.     

Rutin as promising drug for the treatment of Parkinson’s disease: an assessment of MAO-B inhibitory potential by docking,molecular dynamics and DFT studies

ABSTRACT

Inhibitors of monoamine oxidase (MAO)-B have been used for many years in the therapy of Parkinson’s disease (PD). Owing to the safety concerns of the currently used agents, the discovery of novel scaffolds is of considerable interest. MAO-B inhibitory potential of rutin, a flavonoid derived from natural sources, has been established in experimental findings. Hence, the current study seeks to examine the interactions between rutin and crystal structure of human MAO-B enzyme. Molecular docking calculations, as well as molecular dynamics simulations, were employed to investigate the binding mode and the stability of the rutin/MAO-B complex. Energies of highest occupied/lowest unoccupied molecular orbitals were computed through DFT studies and used to calculate electron affinity, hardness, chemical potential, electronegativity, and electrophilicity index in order to investigate the capability of these parameters to influence the ligand–receptor interactions. It was found that rutin traverses both the entrance cavity and the substrate cavity, forcing the Ile-199 ‘gate’ to rotate into its open conformation. It results in the fusion of the two cavities of the MAO-B binding site and directly leads to better binding interactions. Results of the current study can be used for lead modification and development of novel drugs for the treatment of PD.     

On the Application of Widom's Test Particle Method to Homogeneous and Inhomogeneous Fluids

Abstract

Chemical potentials of a homogeneous and an inhomogeneous Lennard-Jones fluid have been determined by molecular dynamics simulations on the vector computer CYBER 205 by applying essentially the fictitious test particle method of Widom. For the homogeneous fluid we find, contrary to the previous result of Guillot and Guissani, that the simulated chemical potential is independent of the particle number. The crucial point, however, is a sufficiently large cut-off radius in the evaluation of the Boltzmann factor. Comparing with our WCA-type perturbation theory, we get agreement in the chemical potentials within 0.1 kT up to the density n[sgrave]³ = 0.80 and a difference of 0.2 kT at n[sgrave]³ = 0.85. For the inhomogeneous case we consider a fluid in a cylindrical pore and integrate Widom's equation over a certain probe volume as suggested earlier by us. Chemical potentials are then calculated independently in five different probe volumes, which are cylindrical shells. The results agree well from the second to the fourth shell. Inaccuracies in the innermost cylinder can be easily explained by bad statistics. In the shell close to the wall the extremely high local density is responsible for the inaccuracies. Extending the probe volume over all cylindrical shells besides the one closest to the wall is thought to yield rather reliable results for the chemical potential. As a by-product of the simulations we also obtained diffusion coefficients, which are given in an appendix.     

Direct Evaluation of Vapour-Liquid Equilibria of Mixtures by Molecular Dynamics Using Gibbs-Duhem Integration

Abstract

We present an extension of the Gibbs-Duhem integration method that permits direct evaluation of vapour-liquid equilibria of mixtures by molecular dynamics. The Gibbs-Duhem integration combines the best elements of the Gibbs ensemble Monte Carlo technique and thermodynamic integration. Given conditions of coexistence of pure substances, simultaneous but independent molecular dynamics simulations of each phase at constant number of particles, constant pressure, constant temperature and constant fugacity fraction of species 2 are carried out in succession along coexistence lines. In each simulation, the coexistence pressure is adjusted to satisfy the Clapeyron-type equation. The Clapeyron-type equation is a first-order nonlinear differential equation that prescribes how the pressure must change with the fugacity fraction of species 2 to maintain coexistence at constant temperature. The Clapeyron-type equation is solved by the predictor-corrector method. Running averages of mole fraction and compressibility factor for the two phases are used to evaluate the right-hand side of the Clapeyron-type equation. The Gibbs-Duhem integration method is applied to three prototypes of binary mixtures of the two-centre Lennard-Jones fluid having various elongations. The starting points on the coexistence curve were taken from published data.     

Modelling of the mechanical and mass transport properties of auxetic molecular sieves: an idealised inorganic (zeolitic) host–guest system

Force field based simulations have been employed to model the structure, and mechanical and mass transport properties of the all-silica zeolite MFI (ZSM5—Si₉₆O₁₉₂). Undeformed and deformed MFI subject to uniaxial loading in each of the three principal directions were investigated. The mechanical properties are predicted to include negative on-axis Poisson's ratios (auxetic behaviour) in the x ₁–x ₃ plane of the undeformed structure, and are strain-dependent. Transformation from positive-to-negative Poisson's ratio behaviour, and vice versa, is predicted for most on-axis Poisson's ratios at critical loading strains. Simulations of the simultaneous sorption of neopentane and benzene guest molecules onto the undeformed host MFI framework indicate a low neopentane-to-benzene loading ratio, consistent with experimental observation. The sorption of these two molecular species onto deformed MFI is Poisson's ratio- and strain-dependent. Uniaxial tensile loading along a direction containing a negative on-axis Poisson's ratio leads to an increase in the loading of the larger neopentane molecules with respect to benzene, strongly correlated with the increase in volume associated with auxetic behaviour.     

Population density analysis for determining the protonation state of the catalytic dyad in BACE1-tertiary carbinamine-based inhibitor complex

BACE1 is an aspartyl protease with a very relevant role in medicinal chemistry related to Alzheimer Disease since it has demonstrated to be a promising therapeutic target for inhibition and possible control for the progress of the peptide accumulation characteristic of this pathology. The enzymatic activity of this protein is given by the aspartic dyad, Asp93 and Asp289, which can adopt several protonation states depending on the chemical nature of its inhibitors, this is, monoprotonated, diprotonated and di-deprotonated states. In the present study, the analysis of the population density, for a series of protein-inhibitor molecular dynamics simulations, was carried out to identify the most feasible protonation state adopted by the catalytic dyad in the presence of tertiary carbinamine (TC) transition state analog inhibitors. The results revealed that the monoprotonated Asp289i state, in which the Asp93 and Asp289 residue side chains are deprotonated and protonated on the inner oxygen, respectively, is the most preferred in the presence of TC family inhibitors. This result was obtained after evaluating, for all 9 possible protonation state configurations, the individual and combined population densities of a set of parameters sensitive to protonation state of the Aspartic dyad, using an X-ray experimental BACE1/TC crystallographic structure as reference. This case study demonstrates again the usefulness of the concept of population density as a quantitative tool to establish the most stable system settings, among all possible, by measuring the level of occurrence of simultaneous events obtained from a sampling over time. These results will help to clear the phenomena related to the TCs inhibitory pathway, as well as assist in the design of better TC inhibitors against Alzheimer’s protease.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.