A novel professional-use synergistic peel technology to reduce visible hyperpigmentation on face: Clinical evidence and mechanistic understanding by computational biology and optical biopsy.

Topicals and chemical peels are the standard of care for management of facial hyperpigmentation. However, traditional therapies have come under recent scrutiny, such as topical hydroquinone (HQ) has some regulatory restrictions, and high concentration trichloroacetic acid (TCA) peel pose a risk in patients with skin of colour. The objective of our research was to identify, investigate and elucidate the mechanism of action of a novel TCA- and HQ-free professional-use chemical peel to manage common types of facial hyperpigmentation. Using computational modelling and in vitro assays on tyrosinase, we identified proprietary multi-acid synergistic technology (MAST). After a single application on human skin explants, MAST peel was found to be more effective than a commercial HQ peel in inhibiting melanin (histochemical imaging and gene expression). All participants completed the case study (N = 9) without any adverse events. After administration of the MAST peel by a dermatologist, the scoring and VISIA photography reported improvements in hyperpigmentation, texture and erythema, which could be linked to underlying pathophysiological changes in skin after peeling, visualized by non-invasive optical biopsy of face. Using reflectance confocal microscopy (VivaScope®) and multiphoton tomography (MPTflex™), we observed reduction in melanin, increase in metabolic activity of keratinocytes, and no signs of inflammatory cells after peeling. Subsequent swabbing of the cheek skin found no microbiota dysbiosis resulting from the chemical peel. The strong efficacy with minimum downtime and no adverse events could be linked to the synergistic action of the ingredients in the novel HQ- and TCA-free professional peel technology.

Microbiome and lipidomic analysis reveal the interplay between skin bacteria and lipids in a cohort study.

Human skin acts as a protective barrier between the body and the external environment. Skin microbiome and intercellular lipids in the stratum corneum (SC) are essential for maintaining skin barrier function. However, the interplay between skin bacteria and the lipids is not fully understood. In this study, we characterized the skin microbiome and SC lipid profiles from the forearm and face in a cohort of 57 healthy participants. 16S rRNA gene sequencing showed the skin microbial composition is significantly different between body locations and genders. Female forearm samples have the highest microbial diversity. The relative abundance of Staphylococcus hominis, Micrococcus luteus, Corynebacterium tuberculostearicum, Finegoldia magna, and Moraxellaceae sp. are significantly higher in the forearm than the face. The predictive functional analysis of 16S rRNA gene sequencing by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) and ANCOM-BC showed different bacterial metabolic pathway profiles between body locations or genders, and identified 271 differential pathways, including arginine and polyamine biosynthesis, chorismate biosynthesis pathways, which are more abundant in the female forearm, and sulfur oxidation pathway, which is more abundant in the male face. The SC lipid profiles differ between the body locations as well. Total free fatty acids (FFA), cholesterol sulfate and sphingosine are more abundant in the face. Dihydro-/6-hydroxy/phyto-ceramides are more abundant in the forearm. The correlation analysis of 16S rRNA gene sequencing and lipids revealed novel interplay between the bacteria and skin lipids. Shannon entropy and S. hominis negatively correlated with FFA, cholesterol sulfate and sphingosine; while positively correlated with dihydro-/6-hydroxy/phyto-ceramides. The correlation of predictive pathway profiles and lipids identified pathways involved in amino acids metabolism, carbohydrates degradation, aromatic compounds metabolism and fatty acid degradation metabolism are positively correlated with dihydro-/6-hydroxy/phyto-ceramides and negatively correlated with FFA, cholesterol sulfate and sphingosine. This study provides insights on the potential correlation between skin microbiome and lipids.

Co-occurrence network analysis reveals the alterations of the skin microbiome and metabolome in adults with mild to moderate atopic dermatitis.

Skin microbiome can be altered in patients with atopic dermatitis (AD). An understanding of the changes from healthy to atopic skin can help develop new targets for treatment by identifying microbial and molecular biomarkers. This study investigates the skin microbiome and metabolome of healthy adult subjects and lesion (ADL) and non-lesion (ADNL) of AD patients by 16S rRNA gene sequencing and mass spectrometry, respectively. Samples from AD patients showed alterations in the diversity and composition of the skin microbiome, with ADL skin having the greatest divergence. Staphylococcus species, especially S. aureus, were significantly increased in AD patients. Metabolomic profiles were also different between the groups. Dipeptide derivatives are more abundant in ADL, which may be related to skin inflammation. Co-occurrence network analysis of the microbiome and metabolomics data revealed higher co-occurrence of metabolites and bacteria in healthy ADNL compared to ADL. S. aureus co-occurred with dipeptide derivatives in ADL, while phytosphingosine-derived compounds showed co-occurrences with commensal bacteria, for example, Paracoccus sp., Pseudomonas sp., Prevotella bivia, Lactobacillus iners, Anaerococcus sp., Micrococcus sp., Corynebacterium ureicelerivorans, Corynebacterium massiliense, Streptococcus thermophilus, and Roseomonas mucosa, in healthy and ADNL groups. Therefore, these findings provide valuable insights into how AD affects the human skin metabolome and microbiome.IMPORTANCEThis study provides valuable insight into changes in the skin microbiome and associated metabolomic profiles in an adult population with mild to moderate atopic dermatitis. It also identifies new therapeutic targets that may be useful for developing personalized treatments for individuals with atopic dermatitis based on their unique skin microbiome and metabolic profiles.

MALDI Mass Spectrometry Imaging and Semi-Quantification of Topically Delivered Lactic Acid.

BACKGROUND: Lactic acid is a common active ingredient in many topical skincare products; however, measuring its delivery into the skin is challenging due to the presence of a large level of endogenous lactic acid. In this study, matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) was used to quantitatively and qualitatively measure the delivery of lactic acid into the skin from a range of topical skincare products. MATERIALS AND METHODS: Porcine skin samples were treated with various skincare products containing lactic acid. After 24 h, skin samples were sectioned and treated via H&E staining or prepared for MALDI-MSI using chemical derivatization. Samples were then analyzed by MALDI-MSI imaging to obtain lactic acid distribution in the entire skin section. RESULTS: Due to the high level of endogenous lactic acid in the skin, a "triple isotope" of lactic acid (L-Lactic acid-13 C3 ), was needed to provide full resolution from the skin's background signal with MALDI-MSI. With this approach, the topically delivered lactic acid could be quantitatively and qualitatively analyzed from a variety of skincare products. CONCLUSIONS: The combination of L-Lactic acid-13 C3 and MALDI-MSI was successfully used to quantitatively and qualitatively measure the topical delivery of lactic acid from a variety of skincare products. This approach could be used in future work to better understand the mode of action of lactic acid as an active ingredient in skincare products.

Multi-omic approach to decipher the impact of skincare products with pre/postbiotics on skin microbiome and metabolome.

Introduction: Although pre/pro/postbiotics have become more prevalent in dermatologic and cosmetic fields, the mode of action when topically applied is largely unknown. A multi-omic approach was applied to decipher the impact of the skincare products with pre/postbiotics on skin microbiome and metabolome. Methods: Subjects with dry skin applied a body wash and body lotion with or without pre/postbiotics for 6 weeks. Skin hydration was measured at baseline, 3 and 6 weeks. Skin swabs were collected for 16S rRNA gene sequencing, metagenomics and metabolomics analysis. Results: Skin hydration significantly increased in both groups. The prebiotic group significantly reduced opportunistic pathogens, e.g., Pseudomonas stutzeri and Sphingomonas anadarae, and increased the commensals, e.g., Staphylococcus equorum, Streptococcus mitis, Halomonas desiderata. Bacterial sugar degradation pathways were enriched in the prebiotic group, while fatty acid biosynthesis pathways were reduced in control. The changes on skin metabolome profiles by the products were more prominent. The prebiotic group performed greater modulation on many clinically-relevant metabolites compared to control. Correlation analysis showed H. desiderata and S. mitis positively correlated with skin hydration, P. stutzeri and S. anadarae negatively correlated with the metabolites that are positively associated with skin hydration improvement. Conclusion: This holistic study supported a hypothesis that the pre/postbiotics increased skin hydration through the modulation of skin microbiome, metabolic pathways and metabolome.

Naïve B cells with low differentiation improve the immune reconstitution of HIV-infected patients.

Incomplete immune reconstitution happens in some HIV-infected patients who have achieved persistent viral suppression under antiretroviral therapy (ART). We performed single-cell RNA sequencing for peripheral blood mononuclear cells to analyze B cell receptor (BCR) repertoire and B cell subtypes in health controls (non-HIV-infected, HCs), HIV-infected immunological responders (IRs), and immunological nonresponders (INRs). We found that the dominant usage of IGHV gene segments of naïve B cells and memory B cells were IGHV3 and IGHV4, and the diversity of BCR repertoire was decreased in INRs. Differentiation trajectory analysis showed that the low differentiation of naïve B cells was related to satisfactory immune status. The cell cycle of B cells with immune-specific genes of IgD+ B cells was degraded in INRs, which was mediated by the anaphase-promoting complex/cyclosome pathway in the phase of G2/M checkpoints. These findings provide significant insights to understand the function of B cell-mediated immune response in immune reconstitution after HIV infection.

Messenger RNA expression profiles and bioinformatics analysis of mouse hippocampi during exercise alleviates methamphetamine dependence via mRNA profile change in hippocampi.

Background: Methamphetamine (METH) is a highly addictive, psychoactive drug which can harm individual health and lead to great social problems. Various approaches have been adopted to address the problems arising from METH addiction, but relapse rates remain high. Recently, it has been found that comprehensive treatment combined with scientific and appropriate exercise interventions can improve the mental state and physical fitness of drug addicts and promote their physical and mental rehabilitation. Long-term, regular exercise improves the symptoms of METH withdrawal and reduces METH relapse. This study aimed to investigate the effects and regulated gene expression related to running exercise in METH-addicted mice. Methods: Male C57BL/6J mice were used to construct a METH addiction model. We performed a running exercise intervention and used conditioned place preference (CPP) to measure the effects of the running intervention on the METH-addicted mice. We also performed RNA sequencing (RNA-seq) and transcriptome analysis on the mice hippocampi, and the functions and differentially expressed genes (DEGs) that were significantly regulated by exercise intervention in the METH-addicted mice were analyzed and noted. Results: The results showed that days of CPP were shortened to 3 days in METH-addicted mice that underwent moderate exercise intervention, compared to 6 days in METH-addicted mice that went without exercise intervention. In addition, hippocampal transcriptome analysis revealed 12 DEGs significantly regulated by exercise intervention. By performing Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, it was revealed that the function of immune responses was significantly enriched in the METH-addicted mice undertaking exercise. The expression of 12 DEGs was verified by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), which showed that relative messenger RNA (mRNA) expression of DEGs was consistent with the RNA-seq results. Conclusions: A running intervention can promote the recovery of METH addiction in mice, and the 12 candidate DEGs from the mouse hippocampus can be used for further research on the regulatory mechanisms of exercise in METH-addicted mice.

Transcriptomic Profiling Reveals Underlying Immunoregulation Mechanisms of Resistant Hypertension in Injection Drug Users.

Background: Hypertension is a common complication in injection drug users (IDU), especially a high proportion of resistant hypertension occurs among them. However, the involving mechanisms remain largely unknown. Methods: We here investigated the key signaling moieties in resistant hypertension in drug users. Analyses were performed with high-throughput transcriptomic sequencing data of peripheral blood from individuals with drug-sensitive hypertension (Ctrl-DS), IDU with resistant hypertension (IDU-DR), and IDU with sensitive hypertension (IDU-DS). Results: We showed that 17 and 1 genes in IDU-DS, 48 and 4 genes in IDU-DR were upregulated and downregulated compared Ctrl-DS, and 2 and 4 genes were upregulated and downregulated in IDU-DR compared with IDU-DS, respectively (p ≤ 0.01 and |log2(FC)| ≥ 1). Differentially expressed genes (DEGs) between Ctrl-DS and IDU-DS were mainly involved in Gene ontology terms of immunoglobulin complex and blood microparticle. DEGs between IDU-DS and IDU-DR were mainly involved in immune system process and immunoglobulin complex. DEGs between Ctrl-DS and IDU-DR were mainly involved in immunoglobulin complex, blood microparticle and cytoplasmic vesicle lumen. We identified 2 gene clusters (brown modules, MEbrown; turquoise module, MEturquoise) correlated with IDU-DR and a gene cluster (magenta module, MEmagenta) correlated with IDU-DS by weighted gene co-expression network analysis (WGCNA). Functional analysis demonstrated that pathways of focal adhesion and focalin-1-rich granule lumen were involved in the development of IDU-DR, and the cytosolic large ribosomal subunit may relate to IDU-DR. Further, immune cell infiltration analysis demonstrated that the abundance of dendritic cells (DCs), natural Treg cells (nTreg), and exhausted T cells (Tex) in IDU-DR and IDU-DS, naïve CD8+ T cells in IDU-DS was significantly different compared with that in Ctrl-DS. The abundance of cytotoxic T cells (Tc) was significantly different between IDU-DS and IDU-DR. Conclusion: Our findings indicated a potential function of immunoregulation mechanisms for resistant hypertension.

Effects and associated transcriptomic landscape changes of methamphetamine on immune cells.

BACKGROUND: Methamphetamine (METH) abuse causes serious health problems, including injury to the immune system, leading to increased incidence of infections and even making withdrawal more difficult. Of course, immune cells, an important part of the immune system, are also injured in methamphetamine abuse. However, due to different research models and the lack of bioinformatics, the mechanism of METH injury to immune cells has not been clarified. METHODS: We examined the response of three common immune cell lines, namely Jurkat, NK-92 and THP-1 cell lines, to methamphetamine by cell viability and apoptosis assay in vitro, and examined their response patterns at the mRNA level by RNA-sequencing. Differential expression analysis of two conditions (control and METH treatment) in three types of immune cells was performed using the DESeq2 R package (1.20.0). And some of the differentially expressed genes were verified by qPCR. We performed Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis of differentially expressed genes by the clusterProfiler R package (3.14.3). And gene enrichment analysis was also performed using MetaScape ( www.metascape.org ). RESULTS: The viability of the three immune cells was differentially affected by methamphetamine, and the rate of NK-cell apoptosis was significantly increased. At the mRNA level, we found disorders of cholesterol metabolism in Jurkat cells, activation of ERK1 and ERK2 cascade in NK-92 cells, and disruption of calcium transport channels in THP-1 cells. In addition, all three cells showed changes in the phospholipid metabolic process. CONCLUSIONS: The results suggest that both innate and adaptive immune cells are affected by METH abuse, and there may be commonalities between different immune cells at the transcriptome level. These results provide new insights into the potential effects by which METH injures the immune cells.

Holistic approach to visualize and quantify collagen organization at macro, micro, and nano-scale.

BACKGROUND: There is scarcity of imaging and image processing techniques for accurate discrimination and quantitation of the dermal extracellular matrix (ECM), primarily collagen. The aim of this study was to develop and demonstrate a holistic imaging and image processing approach to visualize and quantify collagen remodeling at the macro-, micro- and nano-scale using histochemical imaging, Reflectance Confocal Microscopy (RCM), and Atomic Force Microscopy (AFM), respectively. MATERIAL AND METHODS: For proof-of-concept, a commercial anti-aging product known to induce collagen neo-synthesis and re-organization was tested ex vivo on human skin biopsies from two aged females. RESULTS: Relative to untreated skin, collagen fibers (RCM) and fibrils (AFM) were longer and aligned after treatment. The content of collagen and elastin (histochemical imaging and ELISA) statistically improved after treatment. CONCLUSION: Based on our findings, we can conclude: (1) AFM, RCM, and histochemical imaging can accurately discriminate collagen from other ECM components in the skin and (2) the image processing methods can enable quantitation and hence capture small improvements in collagen remodeling after treatment (commercial cosmetic product with collagen organizer technology as proof-of-concept). The reported holistic imaging approach has direct clinical implications for scientists and dermatologists to make quick, real-time, and accurate decisions in skin research and diagnostics.

Professional-Grade TCA-Lactic Acid Chemical Peel: Elucidating Mode of Action to Treat Photoaging and Hyperpigmentation.

Chemical peeling is usually performed by dermatologists, plastic surgeons, and aestheticians for the treatment of photo-aged skin, dyspigmented skin, skin prone to acne eruption, and pre-cancerous skin lesions, etc. In this research paper, we report our investigative findings to understand the mode of action of a commercial professional chemical peel to treat hyperpigmented and photoaged skin. In the in-vitro experiments, we found that the peel inhibits enzymes that are responsible for degradation of collagen and elastin, and the production of melanin pigment. It was surprising to observe that trichloroacetic acid (TCA), which is considered a workhorse of chemical peels for its cauterant action, could synergistically promote the inhibitory action of lactic acid. The rationale behind this synergistic effect could be the conformational change in TCA from linear structure to ring-like structure, which was elucidated through sequential docking using Rosetta software. The in-vitro results on collagen and elastin were corroborated by up-regulation of COL1A, COL3B, fibronectin, and elastin gene expression from 3D human skin equivalents treated with the peel. The findings were further validated through ex-vivo testing on human skin biopsy. The peel significantly inhibits the production of total melanin, and ameliorates photo-damage that was evident through repair of the collagen in the skin exposed to a biological effective dose of UV daily light (6 J/cm2). These research findings have implications for product developers and users (dermatologists, plastic surgeons, and aestheticians) in improving safety and efficacy of chemical peels/peeling.

A novel skin brightening topical technology.

BACKGROUND: Effective skin lightening remains an unmet need in over-the-counter formulations. AIMS: This research examined a topical facial formulation containing hexylresorcinol, silymarin, 20% vitamin C, and 5% vitamin E in a proprietary anhydrous vehicle in skin explants for UVB photoprotective effects and clinical benefits. PATIENTS/METHOD: In vitro investigation examined 12 skin explants to assess the test product and vehicle. Six skin explants received 10 µL of the study product, and six skin explants received the 10 µL of the vehicle. After 96 hours, half the skin samples were exposed to 250 mJ/cm2 of UVB radiation while the other half unexposed. Clinically, 42 female subjects with normal or dry skin 35-55 years with skin types I-VI were enrolled possessing discoloration, uneven skin tone, and fine lines. The dermatologist investigator evaluated brightening, evenness, fine lines, wrinkles, and global appearance. RESULTS: Explants treated with the study product experienced no significant change in gene marker expression of pro-collagen and pro-inflammatory gene markers upon UVB exposure. In contrast, skin explants treated with the vehicle experienced significant decreases in pro-collagen expression and significant increases in pro-inflammatory gene marker expression. Clinically, the greatest improvement as compared to baseline was seen at week 12 (P < .001) with 45% improvement in brightening, 27% improvement in evenness, 25% improvement in lines, and 25% improvement in global facial appearance. CONCLUSION: Hexylresorcinol, silymarin, 20% vitamin C, and 5% vitamin E in a proprietary anhydrous vehicle are effective in decreasing UVB-induced photodamage in skin explants while clinically producing skin brightness improvement.

Segmentation of Online Ferrograph Images with Strong Interference Based on Uniform Discrete Curvelet Transformation.

Through real-time acquisition of the visual characteristics of wear debris in lube oil, an on-line visual ferrograph (OLVF) achieves online monitoring of equipment wear in practice. However, since a large number of bubbles can exist in lube oil and appear as a dynamically changing interference shadow in OLVF ferrograms, traditional algorithms may easily misidentify the interference shadow as wear debris, resulting in a large error in the extracted wear debris characteristic. Based on this possibility, a jam-proof uniform discrete curvelet transformation (UDCT)-based method for the binarization of wear debris images was proposed. Through multiscale analysis of the OLVF ferrograms using UDCT and nonlinear transformation of UDCT coefficients, low-frequency suppression and high-frequency denoising of wear debris images were conducted. Then, the Otsu algorithm was used to achieve binarization of wear debris images under strong interference influence.

A Wear Debris Segmentation Method for Direct Reflection Online Visual Ferrography.

Wear debris in lube oil was observed using a direct reflection online visual ferrograph (OLVF) to monitor the machine running condition and judge wear failure online. The existing research has mainly concentrated on extraction of wear debris concentration and size according to ferrograms under transmitted light. Reports on the segmentation algorithm of the wear debris ferrograms under reflected light are lacking. In this paper, a wear debris segmentation algorithm based on edge detection and contour classification is proposed. The optimal segmentation threshold is obtained by an adaptive canny algorithm, and the contour classification filling method is applied to overcome the problems of excessive brightness or darkness of some wear debris that is often neglected by traditional segmentation algorithms such as the Otsu and Kittler algorithms.

Solid-state NMR, crystallographic, and computational investigation of bisphosphonates and farnesyl diphosphate synthase-bisphosphonate complexes.

Bisphosphonates are a class of molecules in widespread use in treating bone resorption diseases and are also of interest as immunomodulators and anti-infectives. They function by inhibiting the enzyme farnesyl diphosphate synthase (FPPS), but the details of how these molecules bind are not fully understood. Here, we report the results of a solid-state (13)C, (15)N, and (31)P magic-angle sample spinning (MAS) NMR and quantum chemical investigation of several bisphosphonates, both as pure compounds and when bound to FPPS, to provide information about side-chain and phosphonate backbone protonation states when bound to the enzyme. We then used computational docking methods (with the charges assigned by NMR) to predict how several bisphosphonates bind to FPPS. Finally, we used X-ray crystallography to determine the structures of two potent bisphosphonate inhibitors, finding good agreement with the computational results, opening up the possibility of using the combination of NMR, quantum chemistry and molecular docking to facilitate the design of other, novel prenytransferase inhibitors.

Crystallization and preliminary X-ray diffraction study of the farnesyl diphosphate synthase from Trypanosoma brucei.

Farnesyl diphosphate synthase (FPPS) catalyses the formation of farnesyl diphosphate from dimethylallyl diphosphate and isopentenyl diphosphate and is an RNAi-validated drug target in Trypanosoma brucei, the causative agent of African sleeping sickness. A T. brucei FPPS (390 amino acids) has been expressed in Escherichia coli and the recombinant protein has been crystallized in the absence and presence of the bisphosphonate inhibitor minodronate. Diffraction data were collected at 100 K using synchrotron radiation from both crystal types. Crystals obtained in the absence of minodronate belong to space group I222, with unit-cell parameters a = 61.43, b = 118.12, c = 120.04 A, while crystals grown in the presence of minodronate belong to space group C2, with unit-cell parameters a = 131.98, b = 118.10, c = 63.25 A, beta = 112.48 degrees. An initial model of the drug-free protein has been built using a homology model with the molecular-replacement method and refined to 3.3 A resolution. It shows mostly helical structure and resembles the structure of avian farnesyl diphosphate synthase, but with the addition of two loop regions.

3-D QSAR investigations of the inhibition of Leishmania major farnesyl pyrophosphate synthase by bisphosphonates.

We report the activities of 62 bisphosphonates as inhibitors of the Leishmania major mevalonate/isoprene biosynthesis pathway enzyme, farnesyl pyrophosphate synthase. The compounds investigated exhibit activities (IC(50) values) ranging from approximately 100 nM to approximately 80 microM (corresponding to K(i) values as low as 10 nM). The most active compounds were found to be zoledronate (whose single-crystal X-ray structure is reported), pyridinyl-ethane-1-hydroxy-1,1-bisphosphonates or picolyl aminomethylene bisphosphonates. However, N-alicyclic aminomethylene bisphosphonates, such as incadronate (N-cycloheptyl aminomethylene bisphosphonate), as well as aliphatic aminomethylene bisphosphonates containing short (n = 4, 5) alkyl chains, were also active, with IC(50) values in the 200-1700 nM range (corresponding to K(i) values of approximately 20-170 nM). Bisphosphonates containing longer or multiple (N,N-) alkyl substitutions were inactive, as were aromatic species lacking an o- or m-nitrogen atom in the ring, or possessing multiple halogen substitutions or a p-amino group. To put these observations on a more quantitative structural basis, we used three-dimensional quantitative structure-activity relationship techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), to investigate which structural features correlated with high activity. Training set results (N = 62 compounds) yielded good correlations with each technique (R(2) = 0.87 and 0.88, respectively), and were further validated by using a training/test set approach. Test set results (N = 24 compounds) indicated that IC(50) values could be predicted within factors of 2.9 and 2.7 for the CoMFA and CoMSIA methods, respectively. The CoMSIA fields indicated that a positive charge in the bisphosphonate side chain and a hydrophobic feature contributed significantly to activity. Overall, these results are of general interest since they represent the first detailed quantitative structure-activity relationship study of the inhibition of an expressed farnesyl pyrophosphate synthase enzyme by bisphosphonate inhibitors and that the activity of these inhibitors can be predicted within about a factor of 3 by using 3D-QSAR techniques.

Nuclear magnetic resonance shifts in paramagnetic metalloporphyrins and metalloproteins.

We report the first detailed investigation of the (1)H, (13)C, (15)N, and (19)F nuclear magnetic resonance (NMR) spectroscopic shifts in paramagnetic metalloprotein and metalloporphyrin systems. The >3500 ppm range in experimentally observed hyperfine shifts can be well predicted by using density functional theory (DFT) methods. Using spin-unrestricted methods together with large, locally dense basis sets, we obtain very good correlations between experimental and theoretical results: R(2) = 0.941 (N = 37, p < 0.0001) when using the pure BPW91 functional and R(2) = 0.981 (N = 37, p < 0.0001) when using the hybrid functional, B3LYP. The correlations are even better for C(alpha) and C(beta) shifts alone: C(alpha), R(2) = 0.996 (N = 8, p < 0.0001, B3LYP); C(beta), R(2) = 0.995 (N = 8, p < 0.0001, B3LYP), but are worse for C(meso), in part because of the small range in C(meso) shifts. The results of these theoretical calculations also lead to a revision of previous heme and proximal histidine residue (13)C NMR assignments in deoxymyoglobin which are confirmed by new quantitative NMR measurements. Molecular orbital (MO) analyses of the resulting wave functions provide a graphical representation of the spin density distribution in the [Fe(TPP)(CN)(2)](-) (TPP = 5,10,15,20-tetraphenylporphyrinato) system (S = (1)/(2)), where the spin density is shown to be localized primarily in the d(xz) (or d(yz)) orbital, together with an analysis of the frontier MOs in Fe(TPP)Cl (S = (5)/(2)), Mn(TPP)Cl (S = 2), and a deoxymyoglobin model (S = 2). The ability to now begin to predict essentially all heavy atom NMR hyperfine shifts in paramagnetic metalloporphyrins and metalloproteins using quantum chemical methods should open up new areas of research aimed at structure prediction and refinement in paramagnetic systems in much the same way that DFT methods have been used successfully in the past to predict/refine elements of diamagnetic heme protein structures.

Mössbauer quadrupole splittings and electronic structure in heme proteins and model systems: a density functional theory investigation.

We report the results of a series of density functional theory (DFT) calculations aimed at predicting the (57)Fe Mössbauer electric field gradient (EFG) tensors (quadrupole splittings and asymmetry parameters) and their orientations in S = 0, (1)/(2), 1, (3)/(2), 2, and (5)/(2) metalloproteins and/or model systems. Excellent results were found by using a Wachter's all electron basis set for iron, 6-311G for other heavy atoms, and 6-31G for hydrogen atoms, BPW91 and B3LYP exchange-correlation functionals, and spin-unrestricted methods for the paramagnetic systems. For the theory versus experiment correlation, we found R(2) = 0.975, slope = 0.99, intercept = -0.08 mm sec(-)(1), rmsd = 0.30 mm sec(-)(1) (N = 23 points) covering a DeltaE(Q) range of 5.63 mm s(-)(1) when using the BPW91 functional and R(2) = 0.978, slope = 1.12, intercept = -0.26 mm sec(-)(1), rmsd = 0.31 mm sec(-)(1) when using the B3LYP functional. DeltaE(Q) values in the following systems were successfully predicted: (1) ferric low-spin (S = (1)/(2)) systems, including one iron porphyrin with the usual (d(xy))(2)(d(xz)d(yz))(3) electronic configuration and two iron porphyrins with the more unusual (d(xz)d(yz))(4)(d(xy))(1) electronic configuration; (2) ferrous NO-heme model compounds (S = (1)/(2)); (3) ferrous intermediate spin (S = 1) tetraphenylporphinato iron(II); (4) a ferric intermediate spin (S = (3)/(2)) iron porphyrin; (5) ferrous high-spin (S = 2) deoxymyoglobin and deoxyhemoglobin; and (6) ferric high spin (S = (5)/(2)) metmyoglobin plus two five-coordinate and one six-coordinate iron porphyrins. In addition, seven diamagnetic (S = 0, d(6) and d(8)) systems studied previously were reinvestigated using the same functionals and basis set scheme as used for the paramagnetic systems. All computed asymmetry parameters were found to be in good agreement with the available experimental data as were the electric field gradient tensor orientations. In addition, we investigated the electronic structures of several systems, including the (d(xy))(2)(d(xz),d(yz))(3) and (d(xz),d(yz))(4)(d(xy))(1) [Fe(III)/porphyrinate](+) cations as well as the NO adduct of Fe(II)(octaethylporphinate), where interesting information on the spin density distributions can be readily obtained from the computed wave functions.

(57)Fe Mössbauer isomer shifts of heme protein model systems: electronic structure calculations.

We report the results of density functional theory (DFT) calculations of the (57)Fe Mössbauer isomer shifts (delta(Fe)) for a series of 24 inorganic, organometallic, and metalloprotein/metalloporphyrin model systems in S = 0, (1)/(2), 1, (3)/(2), 2, and (5)/(2) spin states. We find an excellent correlation between calculation and experiment over the entire 2.34 mm s(-1) range of isomer shifts: a 0.07-0.08 mm s(-1) rms deviation between calculation and experiment (corresponding to 3-4% of the total delta(Fe) range, depending on the functionals used) with R(2) values of 0.973 and 0.981 (p < 0.0001). The best results are obtained by using the hybrid exchange-correlation functional B3LYP, used previously for (57)Fe Mössbauer quadrupole splittings and (57)Fe NMR chemical shifts and chemical shielding anisotropies. The relativistically corrected value of alpha, alpha(rel), converges with the large basis set used in this work, but the exact values vary somewhat with the methods used: -0.253 a(0)(3) mm s(-1) (Hartree-Fock; HF); -0.316 a(0)(3) mm s(-1) (hybrid HF-DFT; B3LYP), or -0.367 a(0)(3) mm s(-1) (pure DFT; BPW91). Both normal and intermediate spin state isomer shifts are well reproduced by the calculations, as is the broad range of delta(Fe) values: from [Fe(VI)O(4)](2-) (-0.90 mm s(-1) expt; -1.01 mm s(-1) calc) to KFe(II)F(3) (1.44 mm s(-1) expt; 1.46 mm s(-1) calc). Molecular orbital analyses of all inorganic solids as well as all organometallic and metalloporphyrin systems studied reveal that there are three major core MO contributions to rho(tot)(0), the total charge density at the iron nucleus (and hence delta(Fe)), that do not vary with changes in chemistry, while the valence MO contributions are highly correlated with delta(Fe) (R(2) = 0.915-0.938, depending on the functionals used), and the correlation between the valence MO contributions and the total MO contribution is even better (R(2) = 0.965-0.976, depending on the functionals used). These results are of general interest since they demonstrate that DFT methods now enable the accurate prediction of delta(Fe) values in inorganic, organometallic, and metalloporphyrin systems in all spin states and over a very wide range of delta(Fe) values with a very small rms error.