A Humanized Mouse Model Coupled with Computational Analysis Identifies Potent Glycolipid Agonist of Invariant NKT Cells

Invariant natural killer T (iNKT) cells play an important role in many innate and adaptive immune responses, with potential applications in cancer immunotherapy. The glycolipid KRN7000, an α-galactosylceramide, potently activates iNKT cells but has shown limited anticancer effects in human clinical trials conducted so far. In spite of almost three decades of structure-activity relationship studies, no alternative glycolipid has yet emerged as a superior clinical candidate. One reason for the slow progress in this area is that standard mouse models do not accurately reflect the specific ligand recognition by human iNKT cells and their requirements for activation. Here we evaluated a series of KRN7000 analogues using a recently developed humanized mouse model that expresses a human αTCR chain sequence and human CD1d. In this process, a more stimulatory, previously reported but largely overlooked glycolipid was identified, and its activity was probed and rationalized via molecular simulations.

Re-Examining Cannabidiol: Conversion to Tetrahydrocannabinol Using Only Heat

Introduction: In the last decade, the market for Cannabidiol (CBD) has grown to become a near $2 billion dollar industry in the United States alone. This growth can be attributed to a growing social acceptance of marijuana, a more detailed understanding of many health benefits attributed to cannabinoids, and the low cost and wide availibility of hemp-derived cannabinoids. Due to the complex legal histories of marijuana and cannabinoids, the stability and safety of CBD is still an area of interest as research has been restricted globally. Conversion of CBD to its psychoactive isomers, most notably delta-9-Tetrahydrocannabinol (Δ9-THC), presents a significant safety issue for consumers and producers of CBD products. Methods: Previous studies investigating the stability of CBD have focused mainly on replicating conditions experienced during long-term storage at room temperature or lower. Here, we report the thermal stability of CBD at 175°C. Dynamic 1H-NMR experiments and computational electronic structure calculations were used to characterize possible reaction paths from CBD to THC. Results: After 30 minutes of heating, Δ9-THC was produced in detectable amounts in aerobic and anaerobic conditions without an acid catalyst. Conclusions: Our findings support an energetically feasible reaction route that is favorable due to both an increase in phenol acidity at high temperatures and the presence of intramolecular OH-π hydrogen bonding.

Photoactivation of the orange carotenoid protein requires two light-driven reactions mediated by a metastable monomeric intermediate

The orange carotenoid protein (OCP) functions as a sensor of the ambient light intensity and as a quencher of bilin excitons when it binds to the core of the cyanobacterial phycobilisome. We show herein that the photoactivation mechanism that converts the resting, orange-colored state, OCPO, to the active red-colored state, OCPR, requires a sequence of two reactions, each requiring absorption of a single photon by an intrinsic ketocarotenoid chromophore. Global analysis of absorption spectra recorded during continuous illumination of OCPO preparations from Synechocystis sp. PCC 6803 detects the reversible formation of a metastable intermediate, OCPI, in which the ketocarotenoid canthaxanthin exhibits an absorption spectrum with a partial red shift and a broadened vibronic structure compared to that of the OCPO state. While the dark recovery from OCPR to OCPI is a first-order, unimolecular reaction, the subsequent conversion of OCPI to the resting OCPO state is bimolecular, involving association of two OCPO monomers to form the dark-stable OCPO dimer aggregate. These results indicate that photodissociation of the OCPO dimer to form the monomeric OCPO intermediate is the first step in the photoactivation mechanism. Formation of the OCPO monomer from the dimer increases the mean value and broadens the distribution of the solvent-accessible surface area of the canthaxanthin chromophore measured in molecular dynamics trajectories at 300 K. The second step in the photoactivation mechanism is initiated by absorption of a second photon, by canthaxanthin in the OCPO monomer, which obtains the fully red-shifted and broadened absorption spectrum detected in the OCPR product state owing to displacement of the C-terminal domain and the translocation of canthaxanthin more than 12 Å into the N-terminal domain. Both steps in the photoactivation reaction of OCP are likely to involve changes in the structure of the C-terminal domain elicited by excited-state conformational motions of the ketocarotenoid.

Syntheses and Aromaticity Parameters of Hexahydroxypyrrocorphin, Porphotrilactones, and Their Oxidation State Intermediates

Triple OsO₄-mediated dihydroxylation of meso-tetrakis(pentafluorophenyl)porphyrin formed a non-aromatic hexahydroxypyrrocorphin as a single stereo-isomer. A one-step oxidative conversion of all three diol functionalities to lactone moieties generated three out of the four possible porphotrilactone regioisomers that were spectroscopically and structurally characterized. This conversion recovered most of the porphyrinic macrocycle aromatic ring current, as seen in their ¹H NMR spectra and modeled using DFT computations. Stepwise OsO₄-mediated dihydroxylations of porpho-mono- and -di-lactones generated intermediate oxidation state compounds between the pyrrole-three pyrroline macrocycle of the pyrrocorphin and the pyrrole-three oxazolone chromophore of the trilactones. The aromaticity of these chromophores was reduced with increasing number of oxazolone to pyrroline replacements, showing the importance for the presence of three lactone moieties for the retention of the macrocycle aromaticity in the tris-β,β'-modified macrocycles. This work first describes hexahydoxypyrrocorphins, porphotrislactones, and the oxidation state intermediates between them; furthers the understanding of the roles of β-lactone moieties in the expression of porphyrinic macrocycle aromaticity; and generally broadens access to chemically stable pyrrocorphins and pyrrocorphin analogues.

Microsolvation Effects in the Spectral Tuning of Heliorhodopsin

Heliorhodopsins (HeR) are a new category of heptahelical transmembrane photoactive proteins with a covalently linked all-trans retinal. The protonated Schiff base (PSB) nitrogen in the retinal is stabilized by a negatively charged counterion. It is well-known that stronger or weaker electrostatic interactions with the counterion cause a significant spectral blue- or red-shift, respectively, in both microbial and animal rhodopsins. In HeR, however, while Glu107 acts as the counterion, mutations of this residue are not directly correlated with a spectral shift. A molecular dynamics analysis revealed that a water cluster pocket produces a microsolvation effect on the Schiff base, compensating to various extents the replacement of the native counterion. Using a combination of molecular dynamics and quantum mechanical/molecular mechanics (QM/MM), we study this microsolvation effect on the electronic absorption of the retinylidene Schiff base chromophore of HeR.

A favorable path to domain separation in the orange carotenoid protein

The orange carotenoid protein (OCP) is responsible for nonphotochemical quenching (NPQ) in cyanobacteria, a defense mechanism against potentially damaging effects of excess light conditions. This soluble two-domain protein undergoes profound conformational changes upon photoactivation, involving translocation of the ketocarotenoid inside the cavity followed by domain separation. Domain separation is a critical step in the photocycle of OCP because it exposes the N-terminal domain (NTD) to perform quenching of the phycobilisomes. Many details regarding the mechanism and energetics of OCP domain separation remain unknown. In this work, we apply metadynamics to elucidate the protein rearrangements that lead to the active, domain-separated, form of OCP. We find that translocation of the ketocarotenoid canthaxanthin has a profound effect on the energetic landscape and that domain separation only becomes favorable following translocation. We further explore, characterize, and validate the free energy surface (FES) using equilibrium simulations initiated from different states on the FES. Through pathway optimization methods, we characterize the most probable path to domain separation and reveal the barriers along that pathway. We find that the free energy barriers are relatively small (<5 kcal/mol), but the overall estimated kinetic rate is consistent with experimental measurements (>1 ms). Overall, our results provide detailed information on the requirement for canthaxanthin translocation to precede domain separation and an energetically feasible pathway to dissociation.

Electronic energy levels of porphyrins are influenced by the local chemical environment

Self-assembled islands of 5,10,15,20-tetrakis(pentafluoro-phenyl)porphyrin (2HTFPP) on Au(111) contain two bistable molecular species that differ by shifted electronic energy levels. Interactions with the underlying gold herringbone reconstruction and neighboring 2HTFPP molecules cause approximately 60% of molecules to have shifted electronic energy levels. We observed the packing density decrease from 0.64 ± 0.04 molecules per nm² to 0.38 ± 0.03 molecules per nm² after annealing to 200 °C. The molecules with shifted electronic energy levels show longer-range hexagonal packing or are adjacent to molecular vacancies, indicating that molecule–molecule and molecule–substrate interactions contribute to the shifted energies. Multilayers of porphyrins do not exhibit the same shifting of electronic energy levels which strongly suggests that molecule–substrate interactions play a critical role in stabilization of two electronic species of 2HTFPP on Au(111).

Self-assembled islands of 5,10,15,20-tetrakis(pentafluoro-phenyl)porphyrin (2HTFPP) on Au(111) contain two bistable molecular species that differ by shifted electronic energy levels.

Easily accessible non-aromatic heterocycles with handles: 4-bromo-2,3-dihydrofurans from 1,2-dibromohomoallylic alcohols

The first general preparation of 4-bromo-2,3-dihydrofurans is reported. These non-aromatic heterocycles containing a useful coupling handle are accessed via Cu-catalyzed intramolecular cyclization of 1,2-dibromohomoallylic alcohols, which are themselves available in just two steps from aromatic and aliphatic aldehydes and ketones. Molecular dynamics simulations using the simple substrates and key geometric parameters provide a rationale for the selectivities observed. The synthetic utility of the 4-bromodihydrofurans is also demonstrated.

High-Throughput Non-targeted Chemical Structure Identification Using Gas-Phase Infrared Spectra

The high-throughput identification of unknown metabolites in biological samples remains challenging. Most current non-targeted metabolomics studies rely on mass spectrometry, followed by computational methods that rank thousands of candidate structures based on how closely their predicted mass spectra match the experimental mass spectrum of an unknown. We reasoned that the infrared (IR) spectra could be used in an analogous manner and could add orthologous structure discrimination; however, this has never been evaluated on large data sets. Here, we present results of a high-throughput computational method for predicting IR spectra of candidate compounds obtained from the PubChem database. Predicted spectra were ranked based on their similarity to gas-phase experimental IR spectra of test compounds obtained from the NIST. Our computational workflow (IRdentify) consists of a fast semiempirical quantum mechanical method for initial IR spectra prediction, ranking, and triaging, followed by a final IR spectra prediction and ranking using density functional theory. This approach resulted in the correct identification of 47% of 258 test compounds. On average, there were 2152 candidate structures evaluated for each test compound, giving a total of approximately 555,200 candidate structures evaluated. We discuss several variables that influenced the identification accuracy and then demonstrate the potential application of this approach in three areas: (1) combining IR and mass spectra rankings into a single composite rank score, (2) identifying the precursor and fragment ions using cryogenic ion vibrational spectroscopy, and (3) the incorporation of a trimethylsilyl derivatization step to extend the method compatibility to less-volatile compounds. Overall, our results suggest that matching computational with experimental IR spectra is a potentially powerful orthogonal option for adding significant high-throughput chemical structure discrimination when used with other non-targeted chemical structure identification methods.

Amide-Linked C4″-Saccharide Modification of KRN7000 Provides Potent Stimulation of Human Invariant NKT Cells and Anti-Tumor Immunity in a Humanized Mouse Model

Activation of invariant natural killer T (iNKT) cells by α-galactosylceramides (α-GalCers) stimulates strong immune responses and potent anti-tumor immunity. Numerous modifications of the glycolipid structure have been assessed to derive activating ligands for these T cells with altered and potentially advantageous properties in the induction of immune responses. Here, we synthesized variants of the prototypical α-GalCer, KRN7000, with amide-linked phenyl alkane substitutions on the C4″-position of the galactose ring. We show that these variants have weak iNKT cell stimulating activity in mouse models but substantially greater activity for human iNKT cells. The most active of the C4″-amides in our study showed strong anti-tumor effects in a partially humanized mouse model for iNKT cell responses. In silico analysis suggested that the tether length and degree of flexibility of the amide substituent affected the recognition by iNKT cell antigen receptors of the C4″-amide substituted glycolipids in complex with their antigen presenting molecule CD1d. Our findings establish the use of stable C4″-amide linked additions to the sugar moiety for further exploration of the immunological effects of structural modifications of iNKT cell activating glycolipids and highlight the critical need for more accurate animal models to assess these compounds for immunotherapeutic potential in humans.

Spectral Signatures of Canthaxanthin Translocation in the Orange Carotenoid Protein

The orange carotenoid protein (OCP) is involved in the photoprotective processes in cyanobacteria via nonphotochemical quenching. Triggered by blue-green light absorption, the carotenoid chromophore undergoes translocation, displacing around 12 Å from the C-terminal domain (CTD) to the N-terminal domain (NTD). The detailed molecular rearrangements that occur within the carotenoid and the protein during this process remain largely elusive. By using a combination of molecular dynamics, well-tempered metadynamics, and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations, we were able to mimic the translocation of the carotenoid from the inactive OCP^O and obtain metastable red-shifted states in the photoactivation mechanism, replicating the λ_max values of reference experimental spectra. In addition, our simulations give insight into the structure of the red-shifted form of the inactive state of OCP.

Fluorescence Anisotropy Detection of Barrier Crossing and Ultrafast Conformational Dynamics in the S2 State of β-Carotene

Carotenoids are usually only weakly fluorescent despite being very strong absorbers in the mid-visible region because their first two excited singlet states, S₁ and S₂, have very short lifetimes. To probe the structural mechanisms that promote the nonradiative decay of the S₂ state to the S₁ state, we have carried out a series of fluorescence lineshape and anisotropy measurements with a prototype carotenoid, β-carotene, in four aprotic solvents. The anisotropy values observed in the fluorescence emission bands originating from the S₂ and S₁ states reveal that the large internal rotations of the emission transition dipole moment, as much as 50° relative to that of the absorption transition dipole moment, are initiated during ultrafast evolution on the S₂ state potential energy surface and persist upon nonradiative decay to the S₁ state. Electronic structure calculations of the orientation of the transition dipole moment account for the anisotropy results in terms of torsional and pyramidal distortions near the center of the isoprenoid backbone. The excitation wavelength dependence of the fluorescence anisotropy indicates that these out-of-plane conformational motions are initiated by passage over a low-activation energy barrier from the Franck-Condon S₂ structure. This conclusion is consistent with detection over the 80-200 K range of a broad, red-shifted fluorescence band from a dynamic intermediate evolving on a steep gradient of the S₂ state potential energy surface after crossing the activation barrier. The temperature dependence of the oscillator strength and anisotropy indicate that nonadiabatic passage from S₂ through a conical intersection seam to S₁ is promoted by the out-of-plane motions of the isoprenoid backbone with strong hindrance by solvent friction.

Aluminosilicate Mineralogy and the Sorption of Organic Cations: Interplay between Electrostatic Barriers and Compound Structural Features

Current predictive models of organic cation sorption assume that sorbates interact with all sites on aluminosilicate minerals in the same manner. To examine whether differences in aluminosilicate structure and the resultant changes in electrostatic potential influence the sorption of organic cations, seven smectites were chosen with different proportions of isomorphic substitutions (origin of clay charge) located in octahedral versus tetrahedral layers and with the presence or absence of aluminosilicate interlayers. Sorption coefficients for 14 benzylamine derivatives with systematic differences in compound structures were collected to understand the possible influence of aluminosilicate mineralogy. Benzylamine compounds with methyl group substitution on the charged amine or with electron-donating or -withdrawing ring substituents displayed decreases in cation exchange-normalized sorption coefficients (K_CEC), by up to one order of magnitude, between hectorite (100% isomorphic substitution in the octahedral layer) and nontronite (100% isomorphic substitution in the tetrahedral layer). To understand this difference across aluminosilicates, stochastic molecular models of the various aluminosilicate minerals with interlayers were performed. These models showed that negative charge density associated with tetrahedral sites results in high positive electrostatic energy barriers within the interlayer, creating a penalty for compounds with positive charge spread over a larger compound surface area as occurs from primary to quaternary amines. Conversely, clays with charge originating from octahedral sites produce low electrostatic potential barriers within the interlayer, decreasing the penalty for quaternary amine sorption. Trends for nine cationic pharmaceutical compounds, which varied in size, group alkylation, and/or polar substituents, demonstrated similar decreases in K_CEC values to aluminosilicate minerals with high electrostatic energy barriers. Overall, aluminosilicate mineralogy was found to exert a large influence (0.5-1 order of magnitude in sorption coefficients) on organic cation sorption. The application of atomistic electrostatic potential mapping of both sorbent and sorbate structures provided insights to explain trends in sorption coefficients that could not be described by the basic electrostatic potential theory or by assuming that sorbate structure moieties yielded additive sorption contributions.

Origins of the Electronic Modulations of Bacterio- and Isobacteriodilactone Regioisomers

Advances in the utilization of porphyrinoids for photomedicine, catalysis, and artificial photosynthesis require a fundamental understanding of the relationships between their molecular connectivity and resulting electronic structures. Herein, we analyze how the replacement of two pyrrolic C_β═C_β bonds of a porphyrin by two lactone (O═C-O) moieties modulates the ground-state thermodynamic stability and electronic structure of the resulting five possible pyrrole-modified porphyrin isomers. We made these determinations based on density functional theory (DFT) and time-dependent DFT computations of the optical spectra of all regioisomers. We also analyzed the computed magnetically induced currents of their aromatic π-systems. All regioisomers adopt the tautomeric state that maximizes aromaticity, whether or not transannular steric strains are incurred. In all isomers, the O═C_β-O_β bonds were found to support a macrocycle diatropic ring current. We attributed this to the delocalization of nonbonding electrons from the ring oxa- and oxo-atoms into the macrocycle. As a consequence of this delocalization, the dilactone regioisomers are as-or even more-aromatic than their hydroporphyrin congeners. The electronic structures follow different trends for the bacteriochlorin- and isobacteriochlorin-type isomers. The presence of either oxo- or oxa-oxygens conjugated with the macrocyclic π-system was found to be the minimal structural requirement for the regioisomers to exhibit distinct electronic properties. Our computational methods and mechanistic insights provide a basis for the systematic exploration of the physicochemical properties of porphyrinoids as a function of the number, relative orientation, and degree of macrocycle-π-conjugation of β-substituents, in general, and for dilactone-based porphyrinic chromophores, in particular.

Bacterio- and Isobacteriodilactones by Stepwise or Direct Oxidations of meso-Tetrakis(pentafluorophenyl)porphyrin

Porpholactones are porphyrinoids in which one or more β,β'-bonds of the parent chromophore were replaced by lactone moieties. Accessible to varying degrees by direct and nonselective oxidations of porphyrins, the rational syntheses of all five dilactone isomers along stepwise, controlled, and high-yielding routes via porphyrin → tetrahydroxyisobacteriochlorin metal complexes → isobacteriochlorindilactone metal complexes or porphyrin → tetrahydroxybacteriochlorin → bacteriochlorindilactone (and related) pathways, respectively, are described. A major benefit of these complementary routes over established methods is the simplicity of the isolation of the dilactones because of the reduced number of side products formed. In an alternative approach we report the direct and selective conversion of free base meso-tetrakis(pentafluorophenyl)porphyrin to all isomers of free base isobacteriodilactones using the oxidant cetyltrimethylN⁺MnO₄^-. The solid-state structures of some of the isomers and their precursors are reported, providing data on the conformational modulation induced by the derivatizations. We also rationalize computationally their differing thermodynamic stability and electronic properties. In making new efficient routes toward these dilactone isomers available, we enable the further study of this diverse class of porphyrinoids.

Atomically precise Au144(SR)60 nanoclusters (R = Et, Pr) are capped by 12 distinct ligand types of 5-fold equivalence and display gigantic diastereotopic effects

For two decades, Au144(SR)60 has been one of the most studied and used thiolate (SR) protected gold nanoclusters. In many ways, however, it proved to be a challenging and elusive case, also because of the difficulties in solving its structure by single-crystal X-ray crystallography. We used very short thiols and could prepare Au144(SC2H5)60 and Au144(SC3H7)60 in a very pure form, which was confirmed by UV-vis absorption spectroscopy and very regular electrochemistry patterns. Inductively coupled plasma and electrospray ionization mass spectrometries gave definite proof of the Au144(SR)60 stoichiometry. High-resolution 1D and 2D NMR spectroscopy in the solution phase provided the result of assessing the presence of 12 ligand types in exactly the same amount (5-fold equivalence). Equally important, we found that the two protons belonging to each methylene group along the thiolate chain are diastereotopic. For the α-CH2 protons, the diastereotopic effect can be indeed gigantic, as it reaches chemical-shift differences of 2.9 ppm. DFT calculations provided insights into the relationship between structure and NMR results. In particular, the 12 ligand types and corresponding diastereotopic effects may be explained by considering the presence of C-H···S hydrogen bonds. These results thus provide fundamental insights into the structure of the thiolate layer capping this long-studied gold nanocluster.

Carotenoid-Chlorophyll Interactions in a Photosynthetic Antenna Protein: A Supramolecular QM/MM Approach

Multichromophoric interactions control the initial events of energy capture and transfer in the light harvesting peridinin-chlorophyll a protein (PCP) from marine algae dinoflagellates. Due to the van der Waals association of the carotenoid peridinin (Per) with chlorophyll a in a unique 4:1 stoichiometric ratio, supramolecular quantum mechanical/molecular mechanical (QM/MM) calculations are essential to accurately describe structure, spectroscopy, and electronic coupling. We show that, by enabling inter-chromophore electronic coupling, substantial effects arise in the nature of the transition dipole moment and the absorption spectrum. We further hypothesize that inter-protein domain Per-Per interactions are not negligible, and are needed to explain the experimental reconstruction features of the spectrum in wild-type PCP.

The limited extent of the electronic modulation of chlorins and bacteriochlorins through chromene-annulation

Optical data (UV-vis absorption and fluorescence emission spectra, including fluorescence yields and lifetimes) and electrochemical measurements are used to quantify the modulation of the electronic properties of meso-tetrakis(pentafluorophenyl)-chlorin diol and -bacteriochlorin tetraols upon intramolecular chromene-annulation, including the investigation of regio- and stereoisomers. The small modulations of the frontier orbitals of the porphyrinoids are rationalized using DFT computations and can be traced to small electronic effects due to the co-planarized meso-aryl groups in combination with conformational effects.

Dual Modifications of α-Galactosylceramide Synergize to Promote Activation of Human Invariant Natural Killer T Cells and Stimulate Anti-tumor Immunity

Glycosylceramides that activate CD1d-restricted invariant natural killer T (iNKT) cells have potential therapeutic applications for augmenting immune responses against cancer and infections. Previous studies using mouse models identified sphinganine variants of α-galactosylceramide as promising iNKT cell activators that stimulate cytokine responses with a strongly proinflammatory bias. However, the activities of sphinganine variants in mice have generally not translated well to studies of human iNKT cell responses. Here, we show that strongly proinflammatory and anti-tumor iNKT cell responses were achieved in mice by a variant of α-galactosylceramide that combines a sphinganine base with a hydrocinnamoyl ester on C6″ of the sugar. Importantly, the activities observed with this variant were largely preserved for human iNKT cell responses. Structural and in silico modeling studies provided a mechanistic basis for these findings and suggested basic principles for capturing useful properties of sphinganine analogs of synthetic iNKT cell activators in the design of immunotherapeutic agents.

Light Harvesting by Equally Contributing Mechanisms in a Photosynthetic Antenna Protein

We report supramolecular quantum mechanics/molecular mechanics simulations on the peridinin-chlorophyll a protein (PCP) complex from the causative algal species of red tides. These calculations reproduce for the first time quantitatively the distinct peridinin absorptions, identify multichromophoric molecular excitations, and elucidate the mechanisms regulating the strongly allowed S₀ (1¹A_g^-) → S₂ (1¹B_u⁺) absorptions of the bound peridinins that span a 58 nm spectral range in the region of maximal solar irradiance. We discovered that protein binding site-imposed conformations, local electrostatics, and electronic coupling contribute equally to the spectral inhomogeneity. Electronic coupling causes coherent excitations among the densely packed pigments. Complementary pairing of tuning mechanisms is the result of a competition between pigment-pigment and pigment-environment interactions. We found that the aqueous solvent works in concert with the charge distribution of PCP to produce a strong correlation between peridinin spectral bathochromism and the local dielectric environment.

A magnetic look into the protecting layer of Au25 clusters

The field of molecular metal clusters protected by organothiolates is experiencing a very rapid growth. So far, however, a clear understanding of the fine interactions between the cluster core and the capping monolayer has remained elusive, despite the importance of the latter in interfacing the former to the surrounding medium. Here, we describe a very sensitive methodology that enables comprehensive assessment of these interactions. Pulse electron nuclear double resonance (ENDOR) was employed to study the interaction of the unpaired electron with the protons of the alkanethiolate ligands in four structurally related paramagnetic Au25(SR)018 clusters (R = ethyl, propyl, butyl, 2-methylpropyl). Whereas some of these structures were known, we present the first structural description of the highly symmetric Au25(SPr)018 cluster. Through knowledge of the structural data, the ENDOR signals could be successfully related to the types of ligand and the distance of the relevant protons from the central gold core. We found that orbital distribution affects atoms that can be as far as 6 Å from the icosahedral core. Simulations of the spectra provided the values of the hyperfine coupling constants. The resulting information was compared with that provided by 1H NMR spectroscopy, and molecular dynamics calculations provided useful hints to understanding differences between the ENDOR and NMR results. It is shown that the unpaired electron can be used as a very precise probe of the main structural features of the interface between the metal core and the capping ligands.

The MOD-QM/MM Method: Applications to Studies of Photosystem II and DNA G-Quadruplexes

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and catalytic sites embedded in macrobiomolecules (eg, proteins and nucleic acids). QM/MM methodologies are highly efficient since they implement quantum chemistry methods for modeling only the portion of the system involving bond-breaking/forming processes (QM layer), as influenced by the surrounding molecular environment described in terms of molecular mechanics force fields (MM layer). Some of the limitations of QM/MM methods when polarization effects are not explicitly considered include the approximate treatment of electrostatic interactions between QM and MM layers. Here, we review recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in biomacromolecules and structural refinement, beyond the common limitations of QM/MM hybrid methods. We focus on photosystem II (PSII) with emphasis on the description of the oxygen-evolving complex (OEC) and its high-resolution extended X-ray absorption fine structure spectra (EXAFS) in conjunction with Monte Carlo structural refinement. Furthermore, we review QM/MM structural refinement studies of DNA G4 quadruplexes with embedded monovalent cations and direct comparisons to NMR data.

QM/MM Models of the O2-Evolving Complex of Photosystem II

This paper introduces structural models of the oxygen-evolving complex of photosystem II (PSII) in the dark-stable S1 state, as well as in the reduced S0 and oxidized S2 states, with complete ligation of the metal-oxo cluster by amino acid residues, water, hydroxide, and chloride. The models are developed according to state-of-the-art quantum mechanics/molecular mechanics (QM/MM) hybrid methods, applied in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus, recently reported at 3.5 Å resolution. Manganese and calcium ions are ligated consistently with standard coordination chemistry assumptions, supported by biochemical and spectroscopic data. Furthermore, the calcium-bound chloride ligand is found to be bound in a position consistent with pulsed electron paramagnetic resonance data obtained from acetate-substituted PSII. The ligation of protein ligands includes monodentate coordination of D1-D342, CP43-E354, and D1-D170 to Mn(1), Mn(3), and Mn(4), respectively; η(2) coordination of D1-E333 to both Mn(3) and Mn(2); and ligation of D1-E189 and D1-H332 to Mn(2). The resulting QM/MM structural models are consistent with available mechanistic data and also are compatible with X-ray diffraction models and extended X-ray absorption fine structure measurements of PSII. It is, therefore, conjectured that the proposed QM/MM models are particularly relevant to the development and validation of catalytic water-oxidation intermediates.

Au₂₅(SEt)₁₈, a nearly naked thiolate-protected Au₂₅ cluster: structural analysis by single crystal X-ray crystallography and electron nuclear double resonance

X-ray crystallography has been fundamental in discovering fine structural features of ultrasmall gold clusters capped by thiolated ligands. For still unknown structures, however, new tools capable of providing relevant structural information are sought. We prepared a 25-gold atom nanocluster protected by the smallest ligand ever used, ethanethiol. This cluster displays the electrochemistry, mass spectrometry, and UV-vis absorption spectroscopy features of similar Au25 clusters protected by 18 thiolated ligands. The anionic and the neutral form of Au25(SEt)18 were fully characterized by (1)H and (13)C NMR spectroscopy, which confirmed the monolayer's properties and the paramagnetism of neutral Au25(SEt)18(0). X-ray crystallography analysis of the latter provided the first known structure of a gold cluster protected by a simple, linear alkanethiolate. Here, we also report the direct observation by electron nuclear double resonance (ENDOR) of hyperfine interactions between a surface-delocalized unpaired electron and the gold atoms of a nanocluster. The advantages of knowing the exact molecular structure and having used such a small ligand allowed us to compare the experimental values of hyperfine couplings with DFT calculations unaffected by structure's approximations or omissions.

Stability of N-glycosidic bond of (5'S)-8,5'-cyclo-2'-deoxyguanosine

8,5'-Cyclopurine deoxynucleosides are unique tandem lesions containing an additional covalent bond between the base and the sugar. These mutagenic and genotoxic lesions are repaired only by nucleotide excision repair. The N-glycosidic (or C1'-N9) bond of 2'-deoxyguanosine (dG) derivatives is usually susceptible to acid hydrolysis, but even after cleavage of this bond of the cyclopurine lesions, the base would remain attached to the sugar. Here, the stability of the N-glycosidic bond and the products formed by formic acid hydrolysis of (5'S)-8,5'-cyclo-2'-deoxyguanosine (S-cdG) were investigated. For comparison, the stability of the N-glycosidic bond of 8,5'-cyclo-2',5'-dideoxyguanosine (ddcdG), 8-methyl-2'-deoxyguanosine (8-Me-dG), 7,8-dihydro-8-oxo-2'-deoxyguanosine (8-Oxo-dG), and dG was also studied. In various acid conditions, S-cdG and ddcdG exhibited similar stability to hydrolysis. Likewise, 8-Me-dG and dG showed comparable stability, but the half-lives of the cyclic dG lesions were at least 5-fold higher than those of dG or 8-Me-dG. NMR studies were carried out to investigate the products formed after the cleavage of the C1'-N9 bond. 2-Deoxyribose generated α and β anomers of deoxyribopyranose and deoxyribopyranose oligomers following acid treatment. S-cdG gave α- and β-deoxyribopyranose linked guanine as the major products, but α and β anomers of deoxyribofuranose linked guanine and other products were also detected. The N-glycosidic bond of 8-Oxo-dG was found exceptionally stable in acid. Computational studies determined that both the protonation of the N7 atom and the rate constant in the bond breaking step control the overall kinetics of hydrolysis, but both varied for the molecules studied indicating a delicate balance between the two steps. Nevertheless, the computational approach successfully predicted the trend observed experimentally. For 8-Oxo-dG, the low pK(a) of O(8) and N3 prevented appreciable protonation, making the free energy for N-glycosidic bond cleavage in the subsequent step very high.

Recognition and reactivity in the binding between Raf kinase inhibitor protein and its small-molecule inhibitor locostatin

The present work is aimed to provide detail on the binding process between Raf kinase inhibitor protein (RKIP) and locostatin, the only exogenous compound known to alter the function of RKIP. Understanding the basis of RKIP inhibition for use in pharmacological applications is of considerable interest, as dysregulated RKIP expression has the potential to contribute to pathophysiological processes. Herein, we report a series of atomistic models to describe the protein-ligand recognition step and the subsequent reactivity steps. Modeling approaches include ligand docking, molecular dynamics, and quantum mechanics/molecular mechanics calculations. We expect that such a computational assay will serve to study similar complexes in which potency is associated with recognition and reactivity. Although previous data suggested a single amino acid residue (His86) to be involved in the binding of locostatin, the actual ligand conformation and the steps involved in the reactivity process remain elusive from a detailed atomistic description. We show that the first reaction step, consisting of a nucleophilic attack of the nitrogen (Nε) of His86 at the sp(2)-hybridized carbon (C2) of locostatin, presents a late transition state (almost identical to the product). The reaction is followed by a hydrogen abstraction and hydrolysis. The theoretically predicted overall rate constant (6 M(-1) s(-1)) is in a very good agreement with the experimentally determined rate constant (13 M(-1) s(-1)).

Handedness enantioselection of carbon nanotubes using helical assemblies of flavin mononucleotide

In order to truly unlock advanced applications of single-walled carbon nanotubes (SWNTs), one needs to separate them according to both chirality and handedness. Here we show that the chiral D-ribityl phosphate chain of flavin mononucleotide (FMN) induces a right-handed helix that enriches the left-handed SWNTs for all suspended (n,m) species. Such enantioselectivity stems from the sp(3) hybridization of the N atom anchoring the sugar moiety to the flavin ring. This produces two FMN conformations (syn and anti) analogous to DNA. Electrostatic interactions between the neighboring uracil moiety and the 2'-OH group of the side chain provide greater stability to the anti-FMN conformation that leads to a right-handed FMN helix. The right-handed twist that the FMN helix imposes to the underlying nanotube, similar to "Indian burn", causes diameter dilation of only the left-handed SWNTs, whose improved intermolecular interactions with the overlaying FMN helix, impart enantioselection.

QM/MM Prediction of the Stark Shift in the Active Site of a Protein

Recent developments in the biophysical characterization of proteins have provided a means of directly measuring electrostatic fields by introducing a probe molecule to the system of interest and interpreting photon absorption in the context of the Stark effect. To fully account for this effect, the development of accurate atomistic models is of paramount importance. However, suitable computational protocols for evaluating Stark shifts in proteins are yet to be established. In this work, we present a comprehensive computational method to predict the change in absorption frequency of a probe functional group as a direct result of a perturbation in its surrounding electrostatic field created by a protein environment, i.e., the Stark shift. We apply the method to human aldose reductase, a key protein enzyme that catalyzes the reduction of monosaccharides. We develop a protocol based on a combination of molecular dynamics and moving-domain QM/MM methods, which achieves quantitative agreement with experiment. We outline the difficulties in predicting localized electrostatic field changes within a protein environment, and by extension the Stark shift, due to a protein site mutation. Furthermore, the combined use of Stark effect spectroscopy and computational modeling is used to predict the protonation state of ionizable residues in the vicinity of the electrostatic probe.

Effect of the charge state (z = -1, 0, +1) on the nuclear magnetic resonance of monodisperse Au25[S(CH2)2Ph]18(z) clusters

Monodisperse Au(25)L(18)(0) (L = S(CH(2))(2)Ph) and [n-Oct(4)N(+)][Au(25)L(18)(-)] clusters were synthesized in tetrahydrofuran. An original strategy was then devised to oxidize them: in the presence of bis(pentafluorobenzoyl) peroxide, the neutral or the negatively charged clusters react as efficient electron donors in a dissociative electron-transfer (ET) process, in the former case yielding [Au(25)L(18)(+)][C(6)F(5)CO(2)(-)]. As opposed to other reported redox methods, this dissociative ET approach is irreversible, easily controllable, and clean, particularly for NMR purposes, as no hydrogen atoms are introduced. By using this approach, the -1, 0, and +1 charge states of Au(25)L(18) could be fully characterized by (1)H and (13)C NMR spectroscopy, using one- and two-dimensional techniques, in various solvents, and as a function of temperature. For all charge states, the NMR results and analysis nicely match recent structural findings about the presence of two different ligand populations in the capping monolayer, each resonance of the two ligand families displaying distinct NMR patterns. The radical nature of Au(25)L(18)(0) is particularly evident in the (1)H and (13)C NMR patterns of the inner ligands. The NMR behavior of radical Au(25)L(18)(0) was also simulated by DFT calculations, and the interplay between theory and experiments revealed a fundamental paramagnetic contribution coming from Fermi contact shifts. Interestingly, the NMR patterns of Au(25)L(18)(-) and Au(25)L(18)(+) were found to be quite similar, pointing to the latter cluster form as a diamagnetic species.

An extrapolation method for computing protein solvation energies based on density fragmentation of a graphical surface tessellation

Modeling chemical events inside proteins often require the incorporation of solvent effects via continuum polarizable models. One of these approaches is based on the assumption that the interface between solute and solvent acts as a conductor. Image charges are added on the molecular surface to satisfy the appropriate conductor boundary conditions in the presence of solute charges. As in the case of other polarizable continuum models that are based on surface tessellation, the simplest implementation of this approach is often limited to several hundred atoms due to a matrix inversion, which scales as the cube of the number or tesserae. For larger systems, approaches that use iterative matrix solvers coupled to fast summation methods must be used. In the present work, we develop a self-consistent approach to obtain conductor-like screening charges suitable for applications in proteins. The approach is based on a density fragmentation of a graphical surface tessellation. This method, although approximate, provides a straightforward scheme of parallelization, which can in principle be added to existing linear scaling implementations of conductor-like models. We implement this method in conjunction with a fixed charge model for the protein, as well as with a moving domain QM/MM description of the protein. In the latter case, the overall result leads to a charge distribution within the protein determined by self-polarization and polarization due to solvent.

Locostatin Disrupts Association of Raf Kinase Inhibitor Protein With Binding Proteins by Modifying a Conserved Histidine Residue in the Ligand-Binding Pocket

Raf kinase inhibitor protein (RKIP) interacts with a number of different proteins and regulates multiple signaling pathways. Here, we show that locostatin, a small molecule that covalently binds RKIP, not only disrupts interactions of RKIP with Raf-1 kinase, but also with G protein-coupled receptor kinase 2. In contrast, we found that locostatin does not disrupt binding of RKIP to two other proteins: inhibitor of κB kinase α and transforming growth factor β-activated kinase 1. These results thus imply that different proteins interact with different regions of RKIP. Locostatin's mechanism of action involves modification of a nucleophilic residue on RKIP. We observed that after binding RKIP, part of locostatin is slowly hydrolyzed, leaving a smaller RKIP-butyrate adduct. We identified the residue alkylated by locostatin as His86, a highly conserved residue in RKIP's ligand-binding pocket. Computational modeling of the binding of locostatin to RKIP suggested that the recognition interaction between small molecule and protein ensures that locostatin's electrophilic site is poised to react with His86. Furthermore, binding of locostatin would sterically hinder binding of other ligands in the pocket. These data provide a basis for understanding how locostatin disrupts particular interactions of RKIP with RKIP-binding proteins and demonstrate its utility as a probe of specific RKIP interactions and functions.

Computational and experimental investigations of mono-septanoside binding by Concanavalin A: correlation of ligand stereochemistry to enthalpies of binding

Structure-energy relationships for a small group of pyranose and septanose mono-saccharide ligands are developed for binding to Concanavalin A (ConA). The affinity of ConA for methyl "manno"β-septanoside 7 was found to be higher than any of the previously reported mono-septanoside ligands. Isothermal titration calorimetry (ITC) in conjunction with docking simulations and quantum mechanics/molecular mechanics (QM/MM) modeling established the specific role of binding enthalpy in the structure-energy relations of ConA bound to natural mono-saccharides and unnatural mono-septanosides. An important aspect in the differential binding among ligands is the deformation energy required to reorganize internal hydroxyl groups upon binding of the ligand to ConA.

QM/MM approaches in medicinal chemistry research

One of the goals of medicinal chemistry concerns the ability to compute protein-ligand interactions based on the structural knowledge of the receptor. To this end, the majority of current approaches incorporate classical force field potentials to describe receptor-ligand interactions. One of the most critical problems of standard molecular mechanics (MM) force fields is their fixed-charge treatment of electrostatic interactions. Two problems are derived from this approximation, polarization and charge transfer. As an immediate step in computational complexity, it seems natural to incorporate Quantum Mechanics (QM) within a hybrid QM/MM approach, which has shown to be a useful tool to describe structural and mechanistic aspects of chromophores and prosthetic residues in proteins. In this review, we describe specifically the role of QM/MM methods and their various applications to computational drug design and medicinal chemistry research in general.

The MoD-QM/MM methodology for structural refinement of photosystem II and other biological macromolecules

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and structural refinement of macrobiomolecules (e.g., proteins and nucleic acids). These methods are highly efficient, since they implement quantum chemistry techniques for modeling only the small part of the system (QM layer) that undergoes chemical modifications, charge transfer, etc., under the influence of the surrounding environment. The rest of the system (MM layer) is described in terms of molecular mechanics force fields, assuming that its influence on the QM layer can be roughly decomposed in terms of electrostatic interactions and steric hindrance. Common limitations of QM/MM methods include inaccuracies in the MM force fields, when polarization effects are not explicitly considered, and the approximate treatment of electrostatic interactions at the boundaries between QM and MM layers. This article reviews recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in extended systems beyond the common limitations of QM/MM hybrid methods. We focus on the moving-domain QM/MM (MoD-QM/MM) methodology that partitions the system into many molecular domains and obtains the electrostatic and structural properties of the whole system from an iterative self-consistent treatment of the constituent molecular fragments. We illustrate the MoD-QM/MM method as applied to the description of photosystem II as well as in conjunction with the application of spectroscopically constrained QM/MM optimization methods, based on high-resolution spectroscopic data (extended X-ray absorption fine structure spectra, and exchange coupling constants).

Computational insights into the O2-evolving complex of photosystem II

Mechanistic investigations of the water-splitting reaction of the oxygen-evolving complex (OEC) of photosystem II (PSII) are fundamentally informed by structural studies. Many physical techniques have provided important insights into the OEC structure and function, including X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) spectroscopy as well as mass spectrometry (MS), electron paramagnetic resonance (EPR) spectroscopy, and Fourier transform infrared spectroscopy applied in conjunction with mutagenesis studies. However, experimental studies have yet to yield consensus as to the exact configuration of the catalytic metal cluster and its ligation scheme. Computational modeling studies, including density functional (DFT) theory combined with quantum mechanics/molecular mechanics (QM/MM) hybrid methods for explicitly including the influence of the surrounding protein, have proposed chemically satisfactory models of the fully ligated OEC within PSII that are maximally consistent with experimental results. The inorganic core of these models is similar to the crystallographic model upon which they were based, but comprises important modifications due to structural refinement, hydration, and proteinaceous ligation which improve agreement with a wide range of experimental data. The computational models are useful for rationalizing spectroscopic and crystallographic results and for building a complete structure-based mechanism of water-splitting in PSII as described by the intermediate oxidation states of the OEC. This review summarizes these recent advances in QM/MM modeling of PSII within the context of recent experimental studies.

QM/MM computational studies of substrate water binding to the oxygen-evolving centre of photosystem II

This paper reports computational studies of substrate water binding to the oxygen-evolving centre (OEC) of photosystem II (PSII), completely ligated by amino acid residues, water, hydroxide and chloride. The calculations are based on quantum mechanics/molecular mechanics hybrid models of the OEC of PSII, recently developed in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus. The model OEC involves a cuboidal Mn3CaO4Mn metal cluster with three closely associated manganese ions linked to a single mu4-oxo-ligated Mn ion, often called the 'dangling manganese'. Two water molecules bound to calcium and the dangling manganese are postulated to be substrate molecules, responsible for dioxygen formation. It is found that the energy barriers for the Mn(4)-bound water agree nicely with those of model complexes. However, the barriers for Ca-bound waters are substantially larger. Water binding is not simply correlated to the formal oxidation states of the metal centres but rather to their corresponding electrostatic potential atomic charges as modulated by charge-transfer interactions. The calculations of structural rearrangements during water exchange provide support for the experimental finding that the exchange rates with bulk 18 O-labelled water should be smaller for water molecules coordinated to calcium than for water molecules attached to the dangling manganese. The models also predict that the S1-->S2 transition should produce opposite effects on the two water-exchange rates.

A model of the oxygen-evolving center of photosystem II predicted by structural refinement based on EXAFS simulations

A refined computational structural model of the oxygen-evolving complex (OEC) of photosystem II (PSII) is introduced. The model shows that the cuboidal core Mn3CaO4 with a "dangler" Mn ligated to a corner mu4-oxide ion is maximally consistent with the positioning of the amino acids around the metal cluster as characterized by XRD models and high-resolution spectroscopic data, including polarized EXAFS of oriented single crystals and isotropic EXAFS. It is, therefore, natural to expect that the proposed structural model should be particularly useful to establish the structure of the OEC, consistently with high-resolution spectroscopic data, and for elucidating the mechanism of water-splitting in PSII as described by the intermediate oxidation states of the EC along the catalytic cycle.

Quantum mechanics/molecular mechanics structural models of the oxygen-evolving complex of photosystem II

The annual production of 260 Gtonnes of oxygen, during the process of photosynthesis, sustains life on earth. Oxygen is produced in the thylakoid membranes of green-plant chloroplasts and the internal membranes of cyanobacteria by photocatalytic water oxidation at the oxygen-evolving complex (OEC) of photosystem II (PSII). Recent breakthroughs in X-ray crystallography and advances in quantum mechanics/molecular mechanics (QM/MM) hybrid methods have enabled the construction of chemically sensible models of the OEC of PSII. The resulting computational structural models suggest the complete ligation of the catalytic center by amino acid residues, water, hydroxide and chloride, as determined from the intrinsic electronic properties of the oxomanganese core and the perturbational influence of the surrounding protein environment. These structures are found to be consistent with available mechanistic data, and are also compatible with X-ray diffraction models and extended X-ray absorption fine structure measurements. It is therefore conjectured that these OEC models are particularly relevant for the elucidation of the catalytic mechanism of water oxidation.