A novel set of ss-N-biaryl ether sulfonamide hydroxamates as potential MMPs inhibitors: Molecular dynamics simulations and molecular properties evaluation

Matrix metalloproteinases (MMPs) constitute a family of zinc-dependent proteases involved in the extracellular matrix degradation. MMP-2 and MMP9 are overexpressed in several human cancer types, including melanoma, thus the development of new compounds to inhibit MMPs' activity is desirable. Molecular dynamic simulation and molecular properties calculations were performed on a set of novel beta-N-biaryl ether sulfonamide-based hydroxamates, reported as MMP-2 and MMP-9 inhibitors, for providing data to develop an exploratory analysis. Thermodynamic, electronic, and steric descriptors have significantly discriminated highly active from moderately and less active inhibitors of MMP-2 whereas apparent partition coefficient at pH 1.5 was also significant for the MMP-9 data set. Compound 47 was considered an outlier in all analysis, indicating the presence of a bulky substituent group in R3 is crucial to this set of inhibitors for the establishment of molecular interactions with the S1 subsite of both enzymes, but there is a limit. (C) 2012 Wiley Periodicals, Inc.

Photodynamic Efficiency of Cationic meso-Porphyrins at Lipid Bilayers: Insights from Molecular Dynamics Simulations

Porphyrin derivatives have applications as photoactive drugs in photodynamic therapy. However, little is known about their interactions with phospholipid membranes at the molecular level. We employed molecular dynamics simulations to model the binding between a series of cationic meso-(N-methyl-4-pyridinium)phenylporphyrins and anionic phosphatidylglycerol lipid bilayers. This was done in the presence of molecular oxygen within the membrane. The ability of various porphyrins to cause photodamage was quantified in terms of their immersion depth and degree of exposition to a higher oxygen concentration inside the membrane. Simulations showed that the photodynamic efficiency could be improved as the number of hydrophobic phenyl substituents attached to the porphyrinic ring increased. In the specific case of porphyrins containing two hydrophobic and two charged substituents, the cis isomer was significantly more efficient than the trans. These results correlate well with previous experimental observations. They highlight the importance of both the total charge and amphiphilicity of the photosensitizer for its performance in photodynamic therapy.

Effect of the thermostat in the molecular dynamics simulation on the folding of the model protein chignolin

Molecular dynamics simulations of the model protein chignolin with explicit solvent were carried out, in order to analyze the influence of the Berendsen thermostat on the evolution and folding of the peptide. The dependence of the peptide behavior on temperature was tested with the commonly employed thermostat scheme consisting of one thermostat for the protein and another for the solvent. The thermostat coupling time of the protein was increased to infinity, when the protein is not in direct contact with the thermal bath, a situation known as minimally invasive thermostat. In agreement with other works, it was observed that only in the last situation the instantaneous temperature of the model protein obeys a canonical distribution. As for the folding studies, it was shown that, in the applications of the commonly utilized thermostat schemes, the systems are trapped in local minima regions from which it has difficulty escaping. With the minimally invasive thermostat the time that the protein needs to fold was reduced by two to three times. These results show that the obstacles to the evolution of the extended peptide to the folded structure can be overcome when the temperature of the peptide is not directly controlled.

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

Employing perichromism for probing the properties of carboxymethyl cellulose films: an expedient, accurate method for the determination of the degree of substitution of the biopolymer derivative

The properties of films of carboxymethyl cellulose, CMC, of different degree of substitution, DS, have been examined by the use of perichromic indicators (probes). The film properties that have been determined are: empirical polarity, E-T(33); "acidity", alpha; "basicity", beta; and dipolarity/polarizability, pi*. This has been achieved by employing the following perichromic probes: 4-nitroaniline, 4-nitroanisole, 4-nitro-N,N-dimethylaniline, and 2,6-dichloro-4-(2,4,6-triphenyl-pyridinium-1-yl)phenolate, WB. The correlations between both E-T(33)- or pi* and DS were found to be linear; that between beta and DS is a second order polynomial; no obvious correlation was found between alpha and DS. The polarities of CMC films are in the range of those of butyl alcohols. As models for CMC, we have employed cellulose plus CMC of high DS; oxidized cellulose with degree of oxidation = 0.5; sodium glucuronate. The former model behaved akin to CMC, but the plots of the perichromic properties versus DS showed different slopes/intercepts. FTIR data and molecular dynamics simulations on the solvation of WB have shown that this difference can be traced to more efficient hydrogen bonding between the film of the model and the probe. This affects the intra-molecular charge-transfer energy of the latter, leading to different responses to the variation of DS. Based on the excellent linear correlation between E-T(33) and DS, for CMC from different origins, we suggest that perichromism is a simple, accurate, and expedient alternative for the determination of DS of the biopolymer derivative.

Hb S-Sao Paulo: A new sickling hemoglobin with stable polymers and decreased oxygen affinity

Hb S-Sao Paulo (SP) [HBB:c.20A > T p.Glu6Val: c.196A > G p.Lys65Glu] is a new double-mutant hemoglobin that was found in heterozygosis in an 18-month-old Brazilian male with moderate anemia. It behaves like Hb S in acid electrophoresis, isoelectric focusing and solubility testing but shows different behavior in alkaline electrophoresis, cation-exchange HPLC and RP-HPLC. The variant is slightly unstable, showed reduced oxygen affinity and also appeared to form polymers more stable than the Hb S. Molecular dynamics simulation suggests that the polymerization is favored by interfacial electrostatic interactions. This provides a plausible explanation for some of the reported experimental observations. (C) 2012 Elsevier Inc. All rights reserved.

Extensive structural change of the envelope protein of dengue virus induced by a tuned ionic strength: conformational and energetic analyses

The Dengue has become a global public health threat, with over 100 million infections annually; to date there is no specific vaccine or any antiviral drug. The structures of the envelope (E) proteins of the four known serotype of the dengue virus (DENV) are already known, but there are insufficient molecular details of their structural behavior in solution in the distinct environmental conditions in which the DENVs are submitted, from the digestive tract of the mosquito up to its replication inside the host cell. Such detailed knowledge becomes important because of the multifunctional character of the E protein: it mediates the early events in cell entry, via receptor endocytosis and, as a class II protein, participates determinately in the process of membrane fusion. The proposed infection mechanism asserts that once in the endosome, at low pH, the E homodimers dissociate and insert into the endosomal lipid membrane, after an extensive conformational change, mainly on the relative arrangement of its three domains. In this work we employ all-atom explicit solvent Molecular Dynamics simulations to specify the thermodynamic conditions in that the E proteins are induced to experience extensive structural changes, such as during the process of reducing pH. We study the structural behavior of the E protein monomer at acid pH solution of distinct ionic strength. Extensive simulations are carried out with all the histidine residues in its full protonated form at four distinct ionic strengths. The results are analyzed in detail from structural and energetic perspectives, and the virtual protein movements are described by means of the principal component analyses. As the main result, we found that at acid pH and physiological ionic strength, the E protein suffers a major structural change; for lower or higher ionic strengths, the crystal structure is essentially maintained along of all extensive simulations. On the other hand, at basic pH, when all histidine residues are in the unprotonated form, the protein structure is very stable for ionic strengths ranging from 0 to 225 mM. Therefore, our findings support the hypothesis that the histidines constitute the hot points that induce configurational changes of E protein in acid pH, and give extra motivation to the development of new ideas for antivirus compound design.

Interactions between 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and gamma-Al2O3 (100) surface calculated by density functional theory

The main aim of this work is to investigate the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4C1Im]+[BF4]-) ionic liquid (IL) adsorption on the gamma-Al2O3 (100) by density functional theory calculations to try to rationalize the adsorption as an electrostatic phenomenon. Optimized geometries and interaction energies of IL one-monolayer on the gamma-Al2O3 were obtained on high surface coverage (one cationanion pair per 94.96 nm2). A study of dispersion force was made to estimate its contribution to the adsorption. Overall, the process is ruled by electrostatic interaction between ions and surface. Adsorption of the anion [BF4]- and cation [C4C1Im]+ was also studied by Bader charge analysis and charge density difference for supported and unsupported situations. It is suggested that the IL ions have their charges maintained with significant anion cloud polarization inward to the acid aluminum sites. (c) 2012 Wiley Periodicals, Inc.

GQ-16, a Novel Peroxisome Proliferator-activated Receptor gamma (PPAR gamma) Ligand, Promotes Insulin Sensitization without Weight Gain

The recent discovery that peroxisome proliferator-activated receptor gamma (PPAR gamma) targeted anti-diabetic drugs function by inhibiting Cdk5-mediated phosphorylation of the receptor has provided a new viewpoint to evaluate and perhaps develop improved insulin-sensitizing agents. Herein we report the development of a novel thiazolidinedione that retains similar anti-diabetic efficacy as rosiglitazone in mice yet does not elicit weight gain or edema, common side effects associated with full PPAR gamma activation. Further characterization of this compound shows GQ-16 to be an effective inhibitor of Cdk5-mediated phosphorylation of PPAR gamma. The structure of GQ-16 bound to PPAR gamma demonstrates that the compound utilizes a binding mode distinct from other reported PPAR gamma ligands, although it does share some structural features with other partial agonists, such as MRL-24 and PA-082, that have similarly been reported to dissociate insulin sensitization from weight gain. Hydrogen/deuterium exchange studies reveal that GQ-16 strongly stabilizes the beta-sheet region of the receptor, presumably explaining the compound's efficacy in inhibiting Cdk5-mediated phosphorylation of Ser-273. Molecular dynamics simulations suggest that the partial agonist activity of GQ-16 results from the compound's weak ability to stabilize helix 12 in its active conformation. Our results suggest that the emerging model, whereby "ideal" PPAR gamma-based therapeutics stabilize the beta-sheet/Ser-273 region and inhibit Cdk5-mediated phosphorylation while minimally invoking adipogenesis and classical agonism, is indeed a valid framework to develop improved PPAR gamma modulators that retain antidiabetic actions while minimizing untoward effects.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

The Role of Molecular Conformation and Polarizable Embedding for One- and Two-Photon Absorption of Disperse Orange 3 in Solution

Solvent effects on the one- and two-photon absorption (IPA and 2PA) of disperse orange 3 (DO3) in dimethyl sulfoxide (DMSO) are studied using a discrete polarizable embedding (PE) response theory. The scheme comprises a quantum region containing the chromophore and an atomically granulated classical region for the solvent accounting for full interactions within and between the two regions. Either classical molecular dynamics (MD) or hybrid Car-Parrinello (CP) quantum/classical (QM/MM) molecular dynamics simulations are employed to describe the solvation of DO3 in DMSO, allowing for an analysis of the effect of the intermolecular short-range repulsion, long-range attraction, and electrostatic interactions on the conformational changes of the chromophore and also the effect of the solute-solvent polarization. PE linear response calculations are performed to verify the character, solvatochromic shift, and overlap of the two lowest energy transitions responsible for the linear absorption spectrum of DO3 in DMSO in the visible spectral region. Results of the PE linear and quadratic response calculations, performed using uncorrelated solute-solvent configurations sampled from either the classical or hybrid CP QM/MM MD simulations, are used to estimate the width of the line shape function of the two electronic lowest energy excited states, which allow a prediction of the 2PA cross-sections without the use of empirical parameters. Appropriate exchange-correlation functionals have been employed in order to describe the charge-transfer process following the electronic transitions of the chromophore in solution.

Mode of Peroxisome Proliferator-Activated Receptor gamma Activation by Luteolin

The peroxisome proliferator-activated receptor gamma (PPAR gamma) is a target for treatment of type II diabetes and other conditions. PPAR gamma full agonists, such as thiazolidinediones (TZDs), are effective insulin sensitizers and anti-inflammatory agents, but their use is limited by adverse side effects. Luteolin is a flavonoid with anti-inflammatory actions that binds PPAR gamma but, unlike TZDs, does not promote adipocyte differentiation. However, previous reports suggested variously that luteolin is a PPAR gamma agonist or an antagonist. We show that luteolin exhibits weak partial agonist/antagonist activity in transfections, inhibits several PPAR gamma target genes in 3T3-L1 cells (LPL, ORL1, and CEBP alpha) and PPAR gamma-dependent adipogenesis, but activates GLUT4 to a similar degree as rosiglitazone, implying gene-specific partial agonism. The crystal structure of the PPAR gamma ligand-binding domain (LBD) reveals that luteolin occupies a buried ligand-binding pocket (LBP) but binds an inactive PPAR gamma LBD conformer and occupies a space near the beta-sheet region far from the activation helix (H12), consistent with partial agonist/antagonist actions. A single myristic acid molecule simultaneously binds the LBP, suggesting that luteolin may cooperate with other ligands to bind PPAR gamma, and molecular dynamics simulations show that luteolin and myristic acid cooperate to stabilize the Omega-loop among H2', H3, and the beta-sheet region. It is noteworthy that luteolin strongly suppresses hypertonicity-induced release of the pro-inflammatory interleukin-8 from human corneal epithelial cells and reverses reductions in transepithelial electrical resistance. This effect is PPAR gamma-dependent. We propose that activities of luteolin are related to its singular binding mode, that anti-inflammatory activity does not require H12 stabilization, and that our structure can be useful in developing safe selective PPAR gamma modulators.

Enzyme Microheterogeneous Hydration and Stabilization in Supercritical Carbon Dioxide

Supercritical carbon dioxide is a promising green-chemistry solvent for many enzyme-catalyzed chemical reactions, yet the striking stability of some enzymes in such unconventional environments is not well understood. Here, we investigate the stabilization of the Candida antarctica Lipase B (CALB) in supercritical carbon dioxide-water biphasic systems using molecular dynamics simulations. The preservation of the enzyme structure and optimal activity depend on the presence of small amounts of water in the supercritical dispersing medium. When the protein is at least partially hydrated, water molecules bind to specific sites on the enzyme surface and prevent carbon dioxide from penetrating its catalytic core. Strikingly, water and supercritical carbon dioxide cover the protein surface quite heterogeneously. In the first solvation layer, the hydrophilic residues at the surface of the protein are able to pin down patches of water, whereas carbon dioxide solvates preferentially hydrophobic surface residues. In the outer solvation shells, water molecules tend to cluster predominantly on top of the larger water patches of the first solvation layer instead of spreading evenly around the remainder of the protein surface. For CALB, this exposes the substrate-binding region of the enzyme to carbon dioxide, possibly facilitating diffusion of nonpolar substrates into the catalytic funnel. Therefore, by means of microheterogeneous solvation, enhanced accessibility of hydrophobic substrates to the active site can be achieved, while preserving the functional structure of the enzyme. Our results provide a molecular picture on the nature of the stability of proteins in nonaqueous media.

Recoupled separated-local-field experiments and applications to study intermediate-regime molecular motions

A specific separated-local-field NMR experiment, dubbed Dipolar-Chemical-Shift Correlation (DIPSHIFT) is frequently used to study molecular motions by probing reorientations through the changes in XH dipolar coupling and T-2. In systems where the coupling is weak or the reorientation angle is small, a recoupled variant of the DIPSHIFT experiment is applied, where the effective dipolar coupling is amplified by a REDOR-like pi-pulse train. However, a previously described constant-time variant of this experiment is not sensitive to the motion-induced T-2 effect, which precludes the observation of motions over a large range of rates ranging from hundreds of Hz to around a MHz. We present a DIPSHIFT implementation which amplifies the dipolar couplings and is still sensitive to T-2 effects. Spin dynamics simulations, analytical calculations and experiments demonstrate the sensitivity of the technique to molecular motions, and suggest the best experimental conditions to avoid imperfections. Furthermore, an in-depth theoretical analysis of the interplay of REDOR-like recoupling and proton decoupling based on Average-Hamiltonian Theory was performed, which allowed explaining the origin of many artifacts found in literature data. (C) 2012 Elsevier Inc. All rights reserved.