We present an all-atom additive empirical force field for the hexopyranose monosaccharide form of glucose and its diastereomers allose, altrose, galactose, gulose, idose, mannose, and talose. The model is developed to be consistent with the CHARMM all-atom biomolecular force fields, and the same parameters are used for all diastereomers, including both the α-and β-anomers of each monosaccharide. The force field is developed in a hierarchical manner and reproduces the gas-phase and condensed-phase properties of small-molecule model compounds corresponding to fragments of pyranose monosaccharides. The resultant parameters are transferred to the full pyranose monosaccharides and additional parameter development is done to achieve a complete hexopyranose monosaccharide force field. Parametrization target data include vibrational frequencies, crystal geometries, solute -water interaction energies, molecular volumes, heats of vaporization, and conformational energies, including those for over 1800 monosaccharide conformations at the MP2/ cc-pVTZ//MP2/6-31G(d) level of theory. Though not targeted during parametrization, free energies of aqueous solvation for the model compounds compare favorably with experimental values. Also well-reproduced are monosaccharide crystal unit cell dimensions and ring pucker, densities of concentrated aqueous glucose systems, and the thermodynamic and dynamic properties of the exocyclic torsion in dilute aqueous systems. The new parameter set expands the CHARMM additive force field to allow for simulation of heterogeneous systems that include hexopyranose monosaccharides in addition to proteins, nucleic acids, and lipids.

Empirical force field parameters consistent with the CHARMM additive and classical Drude based polarizable force fields are presented for linear and cyclic ethers. Initiation of the optimization process involved validation of the aliphatic parameters based on linear alkanes and cyclic alkanes. Results showed the transfer to cyclohexane to yield satisfactory agreement with target data; however, in the case of cyclopentane direct transfer of the Lennard-Jones parameters was not sufficient due to ring strain, requiring additional optimization of these parameters for this molecule. Parameters for the ethers were then developed starting with the available aliphatic parameters, with the nonbond parameters for the oxygens optimized to reproduce both gas- and condensed-phase properties. Nonbond parameters for the polarizable model include the use of an anisotropic electrostatic model on the oxygens. Parameter optimization emphasized the development of transferable parameters between the ethers of a given class. The ether models are shown to be in satisfactory agreement with both pure solvent and aqueous solvation properties, and the resulting parameters are transferable to test molecules. The presented force field will allow for simulation studies of ethers in condensed phase and provides a basis for ongoing developments in both additive and polarizable force fields for biological molecules.

Toll-like receptor (TLR) signaling is initiated by dimerization of intracellular Toll/IL-1 receptor resistance (TIR) domains. For all TLRs except TLR3, recruitment of the adapter, myeloid differentiation primary response gene 88 (MyD88), to TLR TIR domains results in downstream signaling culminating in proinflammatory cytokine production. Therefore, blocking TLR TIR dimerization may ameliorate TLR2-mediated hyperinflammatory states. The BB loop within the TLR TIR domain is critical for mediating certain protein-protein interactions. Examination of the human TLR2 TIR domain crystal structure revealed a pocket adjacent to the highly conserved P681 and G682 BB loop residues. Using computer-aided drug design (CADD), we sought to identify a small molecule inhibitor(s) that would fit within this pocket and potentially disrupt TLR2 signaling. In silico screening identified 149 compounds and 20 US Food and Drug Administration-approved drugs based on their predicted ability to bind in the BB loop pocket. These compounds were screened in HEK293T-TLR2 transfectants for the ability to inhibit TLR2-mediated IL-8 mRNA. C 16 H 15 NO 4 (C29) was identified as a potential TLR2 inhibitor. C29, and its derivative, ortho-vanillin (o-vanillin), inhibited TLR2/1 and TLR2/6 signaling induced by synthetic and bacterial TLR2 agonists in human HEK-TLR2 and THP-1 cells, but only TLR2/1 signaling in murine macrophages. C29 failed to inhibit signaling induced by other TLR agonists and TNF-α. Mutagenesis of BB loop pocket residues revealed an indispensable role for TLR2/1, but not TLR2/6, signaling, suggesting divergent roles. Mice treated with o-vanillin exhibited reduced TLR2-induced inflammation. Our data provide proof of principle that targeting the BB loop pocket is an effective approach for identification of TLR2 signaling inhibitors.small molecule inhibitor | BB loop | TLR2 pocket | CADD T oll-like receptors (TLRs) are type I transmembrane receptors that detect conserved "pathogen-associated molecular patterns" from microbes, as well as host-derived "danger-associated molecular patterns" (1). TLR2 heterodimerizes with TLR6 or TLR1 to recognize diacyl lipopeptides or triacyl lipopeptides, respectively (2, 3), present in gram-positive and gram-negative bacteria (4-9).Ligand engagement of TLR2/1 or TLR2/6 activates the myeloid differentiation primary response gene 88 (MyD88)-dependent pathway (i.e., nuclear translocation of NF-κB, activation of MAPKs), resulting in production of proinflammatory cytokines (10). Dysregulated TLR2 signaling has been implicated in numerous diseases (e.g., sepsis, atherosclerosis, tumor metastasis, ischemia/reperfusion injury) (11)(12)(13)(14). Several inhibitors of TLR2 signaling have been developed (15-18), yet none is licensed for human use. A better understanding of the Toll/IL-1 receptor resistance (TIR) domain interactions involved in TLR2 signaling could lead to novel therapeutic agents.Both TLRs and adapter proteins contain a cytoplasmic TIR domain that mediates homotypic and heterotypic interactions ...

Fe-type nitrile hydratase (NHase) is a non-heme Fe(III)-dependent enzyme that catalyzes the hydration of nitriles to the corresponding amides. Despite experimental studies of the enzyme and model Fe(III)-containing complexes, many questions concerning the electronic structure and spectroscopic transitions of the metal center remain unanswered. In addition, the catalytic mechanism of nitrile hydration has not yet been determined. We now report density functional theory (B3LYP/6-31G) calculations on three models of the Fe(III) center in the active site of NHase corresponding to hypothetical intermediates in the enzyme-catalyzed hydration of acetonitrile. Together with natural bond orbital (NBO) analysis of the chemical bonding in these active-site models and INDO/S CIS calculations of their electronic spectra, this theoretical investigation gives new insight into the molecular origin of the unusual low-spin preference and spectroscopic properties of the Fe(III) center. In addition, the low-energy electronic transition observed for the active form of NHase is assigned to a dd transition that is coupled with charge-transfer transitions involving the metal and its sulfur ligands. Calculations of isodesmic ligand-exchange reaction energies provide support for coordination of the Fe(III) center in free NHase by a water molecule rather than a hydroxide ion and suggest that the activation of the nitrile substrate by binding to the metal in the sixth coordination site during catalytic turnover cannot yet be definitively ruled out.

The unusual metal coordination and spin-state of the Fe(III) center in nitrile hydratase (NHase) has stimulated the synthesis of numerous model complexes in efforts to understand the reactivity and spectroscopic properties of the enzyme. A particular problem has been the development of model Fe(III) complexes that exhibit reversible, photolabile binding to nitric oxide (NO) in a manner similar to that observed for the NHase metal center. We now report a detailed NBO analysis of the ground-state chemical bonding in three [Fe-NO](6) complexes that exhibit different responses to irradiation, together with investigations of their spectroscopic properties using semiempirical INDO/S CI singles calculations. Our computational studies reveal a correlation between the photolability of these complexes and the existence of low-energy transitions that promote an electron into the Fe-NO pi(*) antibonding molecular orbital. In addition to providing detailed insights into how the ligand field influences the spectroscopy of these mononuclear complexes, these studies strengthen our previous conclusions regarding the role of post-translational cysteine modification in modulating the photoreactivity of the inactive, NO complex of NHase.

No abstract

The inactive, nitrosyl bound form of Fe-type nitrile hydratase (NHase) contains two active site cysteine residues that are post-translationally modified to sulfenate (SO-) and sulfinate (SO2-) ligands. DFT and INDO/S calculations support the hypothesis that these unusual modifications play a key role in modulating the electronic absorption spectra and photoreactivity of the Fe(III) centre in the enzyme.

We report the calculated visible spectrum of [FeIII(PyPepS)2]- in aqueous solution. From all-classical molecular dynamics simulations on the solute and 200 water molecules with a polarizable force field, 25 solute/solvent configurations were chosen at random from a 50 ps production run and subjected the systems to calculations using time-dependent density functional theory (TD-DFT) for the solute, combined with a solvation model in which the water molecules carry charges and polarizabilities. In each calculation the first 60 excited states were collected in order to span the experimental spectrum. Since the solute has a doublet ground state several excitations to states are of type "three electrons in three orbitals," each of which gives rise to a manifold of a quartet and two doublet states which cannot properly be represented by single Slater determinants. We applied a tentative scheme to analyze this type of spin contamination in terms of Delta and Delta transitions between the same orbital pairs. Assuming the associated states as pure single determinants obtained from restricted calculations, we construct conformation state functions (CFSs), i.e., eigenfunctions of the Hamiltonian Sz and S2, for the two doublets and the quartet for each Delta,Delta pair, the necessary parameters coming from regular and spin-flip calculations. It appears that the lower final states remain where they were originally calculated, while the higher states move up by some tenths of an eV. In this case filtering out these higher states gives a spectrum that compares very well with experiment, but nevertheless we suggest investigating a possible (re)formulation of TD-DFT in terms of CFSs rather than determinants.