Fast, efficient generation of high-quality atomic charges

The novel AM1-BCC charge model quickly and efficiently generates high-quality atomic charges for organic molecules suitable for computer simulations in the solution phase. The concept of the AM1-BCC charge model is to produce atomic charges that emulate the HF/6-31G* electrostatic potential (ESP). Underlying electronic structure features including formal charge and electron density delocalization are first captured by AM1 atomic charges; bond charge corrections (BCCs) are then simply added to these AM1 atomic charges to produce the AM1-BCC charges. The BCCs have been determined such that, when added to the AM1 atomic charges, the resulting AM1-BCC atomic charges emulate the HF/6-31G* ESP. The BCCs were parameterized against a training set of >2700 molecules sampling most organic functional groups and their combinations, as well as an extensive variety of cyclic and fused bicyclic heteroaryl systems. The resulting BCC parameters allow the AM1-BCC charging scheme to handle virtually all organic compounds in The Merck Index and the NCI Database. The AM1-BCC charge model reproduces ab initio dimer energies of a diverse set of molecules with an average error of 0.9 kcal/mol and experimental relative free energies of solvation of a diverse set of compounds with an average error 0.7 kcal/mol

Calculation of absolute free energy of binding for theophylline and its analogs to RNA aptamer using nonequilibrium work values

The massively parallel computation of absolute binding free energy with a well-equilibrated system (MP-CAFEE) has been developed [H. Fujitani, Y. Tanida, M. Ito, G. Jayachandran, C. D. Snow, M. R. Shirts, E. J. Sorin, and V. S. Pande, J. Chem. Phys. ${\bf 123}$, 084108 (2005)]. As an application, we perform the binding affinity calculations of six theophylline-related ligands with RNA aptamer. Basically, our method is applicable when using many compute nodes to accelerate simulations, thus a parallel computing system is also developed. To further reduce the computational cost, the adequate non-uniform intervals of coupling constant $\lambda$, connecting two equilibrium states, namely bound and unbound, are determined. The absolute binding energies $\Delta G$ thus obtained have effective linear relation between the computed and experimental values. If the results of two other different methods are compared, thermodynamic integration (TI) and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) by the paper of Gouda $et al$ [H. Gouda, I. D. Kuntz, D. A. Case, and P. A. Kollman, Biopolymers ${\bf 68}$, 16 (2003)], the predictive accuracy of the relative values $\Delta\Delta G$ is almost comparable to that of TI: the correlation coefficients (R) obtained are 0.99 (this work), 0.97 (TI), and 0.78 (MM-PBSA). On absolute binding energies meanwhile, a constant energy shift of $\sim$ -7 kcal/mol against the experimental values is evident. To solve this problem, several presumable reasons are investigated.Comment: 23 pages including 6 figure

Molecular simulation studies of cyanine-based chromonic mesogens: spontaneous symmetry breaking to form chiral aggregates and the formation of a novel lamellar structure

All‐atom molecular dynamics simulations are performed on two chromonic mesogens in aqueous solution: 5,5′‐dimethoxy‐bis‐(3,3′‐di‐sulphopropyl)‐thiacyanine triethylammonium salt (Dye A) and 5,5′‐dichloro‐bis‐(3,3′‐di‐sulphopropyl)‐thiacyanine triethylammonuim salt (Dye B). Simulations demonstrate the formation of self‐assembled chromonic aggregates with an interlayer distance of ≈0.35 nm, with neighboring molecules showing a predominantly head‐to‐tail antiparallel stacking arrangement to minimize electrostatic repulsion between hydrophilic groups. Strong overlap of the aromatic rings occurs within the self‐assembled columns, characteristic of H‐aggregation in aqueous solution. At low concentrations, aggregates of Dye A form chiral columns, despite the presence of strictly achiral species. Chirality arises out of the minimization of steric repulsion between methoxy groups, which would otherwise disrupt the stacking of aromatic molecular cores. At higher concentrations, simulations suggest the interaction of short columns leads to the formation of an achiral‐layered structure in which hydrophobic aromatic regions of the molecule are sandwiched between two layers of hydrophilic groups. This novel lamellar structure is suggested as a likely candidate for the structure of a J‐aggregate. The latter is known to exhibit intense red‐shifted absorption peaks in solution but their structure has not yet been characterized. Self‐organization of such structures provides a route to the formation of “smectic” chromonic mesophases

Ferrous Iron Binding Key to Mms6 Magnetite Biomineralisation: A Mechanistic Study to Understand Magnetite Formation Using pH Titration and NMR Spectroscopy

Formation of magnetite nanocrystals by magn eto-tactic bacteria is controlled by specificproteins whichregu-late the particles’ nucleation and growth.One such proteinis Mms6. This small, amphiph ilic protein can self-assembleand bind ferric ions to aid in magnetite formation. To under-stand the role of Mms6 duringinvitro iron oxide precipita-tion we have performed in situ pH titrations. We find Mms6has little effect duringferric salt precipitation, but exertsgreatest influence during the incorporation of ferrous ionsand conversion of this salt to mixed-valence iron mineral s,suggesting Mms6 has ahitherto unrecorded ferrous iron in-teracting property which promotes the formation of mag-netiteinferrous-rich solutions.Weshow ferrous binding tothe DEEVE motif within the C-term inal region of Mms6 byNMR spectroscopy,and model these binding events usingmolecular simulations. We conclude that Mms6 functions asamagnetite nucleating protein under cond itions where fer-rous ions predominate

Farnesyl pyrophosphate regulates adipocyte functions as an endogenous PPARγ agonist

The cholesterol biosynthetic pathway produces not only sterols but also non-sterol mevalonate metabolites involved in isoprenoid synthesis. Mevalonate metabolites affect transcriptional and post-transcriptional events that in turn affect various biological processes including energy metabolism. In the present study, we examine whether mevalonate metabolites activate PPARγ (peroxisome-proliferator-activated receptor γ), a ligand-dependent transcription factor playing a central role in adipocyte differentiation. In the luciferase reporter assay using both GAL4 chimaera and full-length PPARγ systems, a mevalonate metabolite, FPP (farnesyl pyrophosphate), which is the precursor of almost all isoprenoids and is positioned at branch points leading to the synthesis of other longer-chain isoprenoids, activated PPARγ in a dose-dependent manner. FPP induced the in vitro binding of a co-activator, SRC-1 (steroid receptor co-activator-1), to GST (glutathione transferase)–PPARγ. Direct binding of FPP to PPARγ was also indicated by docking simulation studies. Moreover, the addition of FPP up-regulated the mRNA expression levels of PPARγ target genes during adipocyte differentiation induction. In the presence of lovastatin, an HMG-CoA (3-hydroxy-3-methylglutaryl-CoA) reductase inhibitor, both intracellular FPP levels and PPARγ-target gene expressions were decreased. In contrast, the increase in intracellular FPP level after the addition of zaragozic acid, a squalene synthase inhibitor, induced PPARγ-target gene expression. The addition of FPP and zaragozic acid promotes lipid accumulation during adipocyte differentiation. These findings indicated that FPP might function as an endogenous PPARγ agonist and regulate gene expression in adipocytes

Sphingolipids contribute to acetic acid resistance in Zygosaccharomyces bailii

Lignocellulosic raw material plays a crucial role in the development of sustainable processes for the production of fuels and chemicals. Weak acids such as acetic acid and formic acid are troublesome inhibitors restricting efficient microbial conversion of the biomass to desired products. To improve our understanding of weak acid inhibition, and to identify engineering strategies to reduce acetic acid toxicity, the highly acetic-acid-tolerant yeast Zygosaccharomyces bailii was studied. The impact of acetic acid membrane permeability on acetic acid tolerance in Z. bailii was investigated with particular focus on how the previously demonstrated high sphingolipid content in the plasma membrane influences acetic acid tolerance and membrane permeability. Through molecular dynamics simulations we concluded that membranes with a high content of sphingolipids are thicker and more dense, increasing the free energy barrier for the permeation of acetic acid through the membrane. Z. bailii cultured with the drug myriocin, known to decrease cellular sphingolipid levels, exhibited significant growth inhibition in the presence of acetic acid, while growth in medium without acetic acid was unaffected by the myriocin addition. Furthermore, following an acetic acid pulse, the intracellular pH decreased more in myriocin-treated cells than in control cells. This indicates a higher inflow rate of acetic acid, and confirms that the reduction in growth of cells cultured with myriocin in the medium with acetic acid, was due to an increase in membrane permeability, thereby demonstrating the importance of a high fraction of sphingolipids in the membrane of Z. bailii to facilitate acetic acid resistance; a property potentially transferable to desired production organisms suffering from weak acid stres

Identification of a novel polyfluorinated compound as a lead to inhibit human enzymes aldose reductase and AKR1B10 : structure determination of both ternary complexes and implications for drug design

Aldo-keto reductases (AKRs) are mostly monomeric enzymes which fold into a highly conserved ([alpha]/[beta])8 barrel, while their substrate specificity and inhibitor selectivity are determined by interaction with residues located in three highly variable external loops. The closely related human enzymes aldose reductase (AR or AKR1B1) and AKR1B10 are of biomedical interest because of their involvement in secondary diabetic complications (AR) and in cancer, e.g. hepatocellular carcinoma and smoking-related lung cancer (AKR1B10). After characterization of the IC50 values of both AKRs with a series of polyhalogenated compounds, 2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol (JF0064) was identified as a lead inhibitor of both enzymes with a new scaffold (a 1,1'-biphenyl-4,4'-diol). An ultrahigh-resolution X-ray structure of the AR-­NADP+-JF0064 complex has been determined at 0.85 Å resolution, allowing it to be observed that JF0064 interacts with the catalytic residue Tyr48 through a negatively charged hydroxyl group (i.e. the acidic phenol). The non-competitive inhibition pattern observed for JF0064 with both enzymes suggests that this acidic hydroxyl group is also present in the case of AKR1B10. Moreover, the combination of surface lysine methylation and the introduction of K125R and V301L mutations enabled the determination of the X-ray crystallo­graphic structure of the corresponding AKR1B10-NADP+-JF0064 complex. Comparison of the two structures has unveiled some important hints for subsequent structure-based drug-design efforts

The Impact of Small Molecule Binding on the Energy Landscape of the Intrinsically Disordered Protein C-Myc

Intrinsically disordered proteins are attractive therapeutic targets owing to their prevalence in several diseases. Yet their lack of well-defined structure renders ligand discovery a challenging task. An intriguing example is provided by the oncoprotein c-Myc, a transcription factor that is over expressed in a broad range of cancers. Transcriptional activity of c-Myc is dependent on heterodimerization with partner protein Max. This protein-protein interaction is disrupted by the small molecule 10058-F4 (1), that binds to monomeric and disordered c-Myc. To rationalize the mechanism of inhibition, structural ensembles for the segment of the c-Myc domain that binds to 1 were computed in the absence and presence of the ligand using classical force fields and explicit solvent metadynamics molecular simulations. The accuracy of the computed structural ensembles was assessed by comparison of predicted and measured NMR chemical shifts. The small molecule 1 was found to perturb the composition of the apo equilibrium ensemble and to bind weakly to multiple distinct c-Myc conformations. Comparison of the apo and holo equilibrium ensembles reveals that the c-Myc conformations binding 1 are already partially formed in the apo ensemble, suggesting that 1 binds to c-Myc through an extended conformational selection mechanism. The present results have important implications for rational ligand design efforts targeting intrinsically disordered proteins