Temperature-dependent spin crossover in neuronal nitric oxide synthase bound with the heme-coordinating thioether inhibitors

A series of L-arginine analogue nitric oxide synthase inhibitors with a thioether tail have been shown to form an Fe-S thioether interaction as evidenced by continuous electron density between the Fe and S atoms. Even so, the Fe-S thioether interaction was found to be far less important for inhibitor binding than the hydrophobic interactions between the alkyl group in the thioether tail and surrounding protein (Martell et al. J. Am. Chem. Soc. 2010 , 132 , 798). However, among the few thioether inhibitors that showed Fe-S thioether interaction in crystal structures, variations in spin state (high-spin or low-spin) were observed dependent upon the heme iron oxidation state and temperature. Since modern synchrotron X-ray data collection is typically carried out at cryogenic temperatures, we reasoned that some of the discrepancies between cryo-crystal structures and room-temperature UV-visible spectroscopy could be the result of temperature-dependent spin-state changes. We, therefore, have characterized some of these neuronal nitric oxide synthase (nNOS)-thioether inhibitor complexes in both crystal and solution using EPR and UV-visible absorption spectrometry as a function of temperature and the heme iron redox state. We found that some thioether inhibitors switch from high to low spin at lower temperatures similar to the "spin crossover" phenomenon observed in many transition metal complexes.     

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.     

Probing the neutral graphene-ionic liquid interface: insights from molecular dynamics simulations

We study basic mechanisms of the interfacial layer formation at the neutral graphite monolayer (graphene)-ionic liquid (1,3-dimethylimidazolium chloride, [dmim][Cl]) interface by fully atomistic molecular dynamics simulations. We probe the interface area by a spherical probe varying the charge (-1e, 0, +1e) as well as the size of the probe (diameter 0.50 nm and 0.38 nm). The molecular modelling results suggest that: there is a significant enrichment of ionic liquid cations at the surface. This cationic layer attracts Cl(-) anions that leads to the formation of several distinct ionic liquid layers at the surface. There is strong asymmetry in cationic/anionic probe interactions with the graphene wall due to the preferential adsorption of the ionic liquid cations at the graphene surface. The high density of ionic liquid cations at the interface adds an additional high energy barrier for the cationic probe to come to the wall compared to the anionic probe. Qualitatively the results from probes with diameter 0.50 nm and 0.38 nm are similar although the smaller probe can approach closer to the wall. We discuss the simulation results in light of available experimental data on the interfacial structure in ionic liquids.     

Minimal pharmacophoric elements and fragment hopping, an approach directed at molecular diversity and isozyme selectivity. Design of selective neuronal nitric oxide synthase inhibitors

Fragment hopping, a new fragment-based approach for de novo inhibitor design focusing on ligand diversity and isozyme selectivity, is described. The core of this approach is the derivation of the minimal pharmacophoric element for each pharmacophore. Sites for both ligand binding and isozyme selectivity are considered in deriving the minimal pharmacophoric elements. Five general-purpose libraries are established: the basic fragment library, the bioisostere library, the rules for metabolic stability, the toxicophore library, and the side chain library. These libraries are employed to generate focused fragment libraries to match the minimal pharmacophoric elements for each pharmacophore and then to link the fragment to the desired molecule. This method was successfully applied to neuronal nitric oxide synthase (nNOS), which is implicated in stroke and neurodegenerative diseases. Starting with the nitroarginine-containing dipeptide inhibitors we developed previously, a small organic molecule with a totally different chemical structure was designed, which showed nanomolar nNOS inhibitory potency and more than 1000-fold nNOS selectivity. The crystallographic analysis confirms that the small organic molecule with a constrained conformation can exactly mimic the mode of action of the dipeptide nNOS inhibitors. Therefore, a new peptidomimetic strategy, referred to as fragment hopping, which creates small organic molecules that mimic the biological function of peptides by a pharmacophore-driven strategy for fragment-based de novo design, has been established as a new type of fragment-based inhibitor design. As an open system, the newly established approach efficiently incorporates the concept of early "ADME/Tox" considerations and provides a basic platform for medicinal chemistry-driven efforts.     

Unfolding the conformational behavior of peptide dendrimers: insights from molecular dynamics simulations

We present here the first comprehensive structural characterization of peptide dendrimers using molecular simulation methods. Multiple long molecular dynamics simulations are used to extensively sample the conformational preferences of five third-generation peptide dendrimers, including some known to bind aquacobalamine. We start by analyzing the compactness of the conformations thus sampled using their radius of gyration profiles. A more detailed analysis is then performed using dissimilarity measures, principal coordinate analysis, and free energy landscapes, with the aim of identifying groups of similar conformations. The results point to a high conformational flexibility of these molecules, with no clear "folded state", although two markedly distinct behaviors were found: one of the dendrimers displayed mostly compact conformations clustered into distinct basins (rough landscape), while the remaining dendrimers displayed mainly noncompact conformations with no significant clustering (downhill landscape). This study brings new insight into the conformational behavior of peptide dendrimers and may provide better routes for their functional design. In particular, we propose a yet unsynthesized peptide dendrimer that might exhibit enhanced ability to coordinate aquocobalamin.     

Phenanthrenoids from Juncus acutus L., new natural lipopolysaccharide-inducible nitric oxide synthase inhibitors

The novel natural product juncutol (1), 1,4,7-trimethyl-8,9-dihydro-4H-cyclopenta[def]phenanthrene-2,6-diol, along with the three related metabolites juncusol (2), dehydrojuncusol (3), and 6-hydroxymethyl-1-methyl-5-vinyl-9,10-dihydrophenanthrene-2-ol (4), were isolated from the rhizomes of Juncus acutus L. (Juncaceae) growing in Egypt. The structural identity of 1 was determined on the basis of spectroscopic analyses, including 2D NMR spectroscopy. The inhibitory effect of these natural products on the expression of inducible nitric oxide synthase (iNOS) in lipopolysaccharide-stimulated RAW264.7 macrophage cells was determined for the first time. The unprecedented symmetrical compound juncutol (1) was found to be the most potent inhibitor against the induction of the proinflammatory iNOS protein.     

How natural materials remove heavy metals from water: mechanistic insights from molecular dynamics simulations

Water pollution by heavy metals is of increasing concern due to its devastating effects on the environment and on human health. For the removal of heavy metals from water sources, natural materials, such as spent-coffee-grains or orange/banana/chestnut peels, appear to offer a potential cheap alternative to more sophisticated and costly technologies currently in use. However, in order to employ them effectively, it is necessary to gain a deeper understanding – at the molecular level – of the heavy metals-bioorganic-water system and exploit the power of computer simulations. As a step in this direction, we investigate via atomistic simulations the capture of lead ions from water by hemicellulose – the latter being representative of the polysaccharides that are common components of vegetables and fruit peels − as well as the reverse process. A series of independent molecular dynamics simulations, both classical and ab initio, reveals a coherent scenario which is consistent with what one would expect of an efficient capture, i.e. that it be fast and irreversible: (i) binding of the metal ions via adsorption is found to happen spontaneously on both carboxylate and hydroxide functional groups; (ii) in contrast, metal ion desorption, leading to solvation in water, involves sizable free-energy barriers.

We investigate via atomistic simulations the capture of lead ions from water by hemicellulose – as representative of the polysaccharides that are common components of vegetables and fruit peels – and the reverse process.     

Vesicle shapes from molecular dynamics simulations

Lipid bilayer membranes are known to form various structures such as large sheets or vesicles. When the two leaflets of the bilayer have an equal composition, the membrane preferentially forms a flat sheet or a spherical vesicle. However, a difference in the composition of the two leaflets may result in a curved bilayer or in a wide variety of vesicle shapes. Vesicles with different shapes have already been shown in experiments and diverse vesicle shapes have been predicted theoretically from energy minimization of continuous curves. Here we present a molecular dynamics study of the effect of small changes in the phospholipid headgroups on the spontaneous curvature of the bilayer and on the resulting vesicle shape transformations. Small asymmetries in the bilayers already result in high spontaneous curvature and large vesicle deformations. Vesicle shapes that are formed include ellipsoids, discoids, pear-shaped vesicles, cup-shaped vesicles, as well as budded vesicles. Comparison of these vesicles with theoretically derived vesicle shapes shows both resemblances and differences.     

Estimating entropies from molecular dynamics simulations

While the determination of free-energy differences by MD simulation has become a standard procedure for which many techniques have been developed, total entropies and entropy differences are still hardly ever computed. An overview of techniques to determine entropy differences is given, and the accuracy and convergence behavior of five methods based on thermodynamic integration and perturbation techniques was evaluated using liquid water as a test system. Reasonably accurate entropy differences are obtained through thermodynamic integration in which many copies of a solute are desolvated. When only one solute molecule is involved, only two methods seem to yield useful results, the calculation of solute-solvent entropy through thermodynamic integration, and the calculation of solvation entropy through the temperature derivative of the corresponding free-energy difference. One-step perturbation methods seem unsuitable to obtain entropy estimates.     

Gating NO release from nitric oxide synthase

We have investigated the kinetics of NO escape from Geobacillus stearothermophilus nitric oxide synthase (gsNOS). Previous work indicated that NO release was gated at position 223 in mammalian enzymes; our kinetics experiments include mutants at that position along with measurements on the wild type enzyme. Employing stopped-flow UV-vis methods, reactions were triggered by mixing a reduced enzyme/N-hydroxy-l-arginine complex with an aerated buffer solution. NO release kinetics were obtained for wt NOS and three mutants (H134S, I223V, H134S/I223V). We have confirmed that wt gsNOS has the lowest NO release rate of known NOS enzymes, whether bacterial or mammalian. We also have found that steric clashes at positions 223 and 134 hinder NO escape, as judged by enhanced rates in the single mutants. The empirical rate of NO release from the gsNOS double mutant (H134/I223V) is nearly as rapid as that of the fastest mammalian enzymes, demonstrating that both positions 223 and 134 function as gates for escape of the product diatomic molecule.     

Binding of novel azole-bridged dinuclear platinum(II) anticancer drugs to DNA: insights from hybrid QM/MM molecular dynamics simulations

Dinuclear Pt-containing compounds might be used to overcome the intrinsic and acquired cell resistance of widely used anticancer drugs such as cisplatin. Recently, the complexes [[cis-Pt(NH3)2]2(mu-OH)(mu-pz)](NO3)2 (with pz = pyrazolate) (1), [[cis-Pt(NH3)2]2(mu-OH)(mu-1,2,3-ta-N(1),N(2))](NO3)2 (with ta = 1,2,3-triazolate) (2), and the binding of 1 to d(CpTpCpTpG*pG*pTpCpTpCp) have been characterized. Here we provide the structural and electronic properties of the free drugs, of the intermediates of binding to guanine bases, and of the products, by performing DFT calculations. Our results show that in 2 an isomerization of the Pt-coordination sphere from N(2) to N(3) of the triazolate unit determines a thermodynamic stabilization of approximately 20 kcal/mol as a consequence of the formation of an allylic structure. In addition, hybrid quantum-classical molecular dynamics simulations of 1 and 2 DNA adducts have shed light on the structural distortions that the drugs induce to the DNA duplex. Our calculations show that the rise and the tilt of the two adjacent guanines are identical in the presence of 1 and 2, but they markedly increase when 2 binds in the N(1),N(3) fashion. In addition, the drugs do not provoke any kink upon binding to the double-stranded DNA, suggesting that they may act with a mechanism different than that of cisplatin. The accuracy of our calculations is established by a comparison with the NMR data for the corresponding complex with 1.     

Properties of inhibitors of methane hydrate formation via molecular dynamics simulations

Within the framework of a proposed two-step mechanism for hydrate inhibition, the energy of binding of four inhibitor molecules (PEO, PVP, PVCap, and VIMA) to a hydrate surface is estimated with molecular dynamic simulations. One key feature of this proposed mechanism is that the binding of an inhibitor molecule to the surface of an ensuing hydrate crystal disrupts growth and therein crystallization. It is found through the molecular dynamic simulations that inhibitor molecules that experimentally exhibit better inhibition strength also have higher free energies of binding, an indirect confirmation of our proposed mechanism. Inhibitors increasing in effectiveness, PEO < PVP < PVCap < VIMA, have increasingly negative (exothermic) binding energies of -0.2 < -20.6 < -37.5 < -45.8 kcal/mol and binding free energies of increasing favorability (+0.4 approximately = +0.5 < -9.4 < -15.1 kcal/mol). Furthermore, the effect of an inhibitor molecule on the local liquid water structure under hydrate-forming conditions was examined and correlated to the experimental effectiveness of the inhibitors. Two molecular characteristics that lead to strongly binding inhibitors were found: (1) a charge distribution on the edge of the inhibitor that mimics the charge separation in the water molecules on the surface of the hydrate and (2) the congruence of the size of the inhibitor with respect to the available space at the hydrate-surface binding site. Equipped with this molecular-level understanding of the process of hydrate inhibition via low-dosage kinetic hydrate inhibitors we can design new, more effective inhibitor molecules.     

Selective ablation of Xe from silicon surfaces: molecular dynamics simulations and experimental laser patterning

The mechanism of laser-induced removal of Xe overlayers from a Si substrate has been investigated employing MD simulations and evaluated by buffer layer assisted laser patterning experiments. Two distinct regimes of overlayer removal are identified in the simulations of a uniform heating of the Si substrate by a 5 ns laser pulse: The intensive evaporation from the surface of the Xe overlayer and the detachment of the entire Xe overlayer driven by explosive boiling in the vicinity of the hot substrate. Simulations of selective heating of only a fraction of the silicon substrate suggest that the lateral heat transfer and bonding to the unheated, colder regions of the Xe overlayer is very efficient and suppresses the separation of a fraction of the overlayer from the substrate. Interaction with surrounding cold Xe is responsible for significant increase in the substrate temperature required for achieving the spatially selective ablation of the overlayer. The predictions of the MD simulations are found to be in a qualitative agreement with the results of experimental measurements of the threshold laser power required for the removal of Xe overlayers of different thickness and the shapes of metallic stripes generated by buffer-assisted laser patterning.     

Insights into uranyl chemistry from molecular dynamics simulations

First-principles and purely classical molecular dynamics (MD) simulations for complexes of the uranyl ion (UO(2)(2+)) are reviewed. Validation of Car-Parrinello MD simulations for small uranyl complexes in aqueous solution is discussed. Special attention is called to the mechanism of ligand-exchange reactions at the uranyl centre and to effects of solvation and hydration on coordination and structural properties. Large-scale classical MD simulations are surveyed in the context of liquid-liquid extraction, with uranyl complexes ranging from simple hydrates to calixarenes, and nonaqueous phases from simple organic solvents and supercritical CO(2) to ionic liquids.     

Micellar crystals in solution from molecular dynamics simulations

Polymers with both soluble and insoluble blocks typically self-assemble into micelles, which are aggregates of a finite number of polymers where the soluble blocks shield the insoluble ones from contact with the solvent. Upon increasing concentration, these micelles often form gels that exhibit crystalline order in many systems. In this paper, we present a study of both the dynamics and the equilibrium properties of micellar crystals of triblock polymers using molecular dynamics simulations. Our results show that equilibration of single micelle degrees of freedom and crystal formation occur by polymer transfer between micelles, a process that is described by transition state theory. Near the disordered (or melting) transition, bcc lattices are favored for all triblocks studied. Lattices with fcc ordering are also found but only at lower kinetic temperatures and for triblocks with short hydrophilic blocks. Our results lead to a number of theoretical considerations and suggest a range of implications to experimental systems with a particular emphasis on Pluronic polymers.     

Nonequilibrium dynamics simulations of nitric oxide release: comparative study of nitrophorin and myoglobin

Nitrophorin 4 (NP4) is a heme protein that reversibly binds nitric oxide (NO), with release rates modulated by pH change. High-resolution structures of NP4 revealed that pH changes and NO binding induce a large conformational rearrangement in two loops that serve to protect the heme-bound NO molecule from solvent. We used extended (110 ns) molecular dynamics simulations of NP4 at pH 5 and pH 7, modeled by selective deprotonation of acidic groups. Conformational and dynamic changes were observed, consistent with those found in the crystal. Further, major solvent movement and NO escape were observed at pH 7, while the ligand remained in the heme binding pocket at pH 5. As a control, we also performed molecular dynamics (MD) simulations of sperm whale myoglobin, where NO migration into the interior cavities of the protein was observed, consistent with previous reports. We constructed a kinetic model of ligand escape to quantitatively relate the microscopic rate constants to the observed rates, and tested the predictions against the experimental data. The results suggest that release rates of diatomic molecules from heme proteins can be varied by several orders of magnitude through modest adjustments in geminate rebinding and gating behavior.     

Kinetic insights into the role of the solvent in the polymorphism of 5-fluorouracil from molecular dynamics simulations

We investigate the fundamental factors controlling polymorphism in 5-fluorouracil by performing molecular dynamics simulations of solutions of the compound in water, nitromethane, and wet nitromethane. Analysis of the effect of solvent on the initial aggregation of 5-fluorouracil molecules shows that the strong binding of water to the 5-fluorouracil molecule hinders the formation of the doubly hydrogen-bonded dimer and, by default, promotes close hydrophobic F...F interactions that are a feature of the unusual (Z' = 4) structure of form I. In contrast, doubly hydrogen-bonded dimers are observed to form readily in solution in dry nitromethane, consistent with the crystallization of the doubly hydrogen-bonded ribbon structure of form II from this solvent. When nitromethane is doped with water, the water forms hydrogen bonds to the solute, interfering with the formation of the doubly hydrogen-bonded dimers, which is consistent with the crystallization of form I from this hygroscopic solvent when it is not dried. Overall, the molecular dynamics simulations provide an atomistic picture of how solvent-solute interactions can significantly affect the initial association of 5-fluorouracil molecules to the extent that they determine the polymorphic outcome of the crystallization.     

Nucleation rate isotherms of argon from molecular dynamics simulations

We report six nucleation rate isotherms of vapor-liquid nucleation of Lennard-Jones argon from molecular dynamics simulations. The isotherms span three orders of magnitude in nucleation rates, 10(23)相似文献