Impact of Peripheral Hydrogen Bond on Electronic Properties of the Primary Acceptor Chlorophyll in the Reaction Center of Photosystem I.

Photosystem I (PS I) is a photosynthetic pigment-protein complex that absorbs light and uses the absorbed energy to initiate electron transfer. Electron transfer has been shown to occur concurrently along two (A- and B-) branches of reaction center (RC) cofactors. The electron transfer chain originates from a special pair of chlorophyll a molecules (P700), followed by two chlorophylls and one phylloquinone in each branch (denoted as A-1, A0, A1, respectively), converging in a single iron-sulfur complex Fx. While there is a consensus that the ultimate electron donor-acceptor pair is P700+A0-, the involvement of A-1 in electron transfer, as well as the mechanism of the very first step in the charge separation sequence, has been under debate. To resolve this question, multiple groups have targeted electron transfer cofactors by site-directed mutations. In this work, the peripheral hydrogen bonds to keto groups of A0 chlorophylls have been disrupted by mutagenesis. Four mutants were generated: PsaA-Y692F; PsaB-Y667F; PsaB-Y667A; and a double mutant PsaA-Y692F/PsaB-Y667F. Contrary to expectations, but in agreement with density functional theory modeling, the removal of the hydrogen bond by Tyr → Phe substitution was found to have a negligible effect on redox potentials and optical absorption spectra of respective chlorophylls. In contrast, Tyr → Ala substitution was shown to have a fatal effect on the PS I function. It is thus inferred that PsaA-Y692 and PsaB-Y667 residues have primarily structural significance, and their ability to coordinate respective chlorophylls in electron transfer via hydrogen bond plays a minor role.

Capturing CO2 in Quadrupolar Binding Pockets: Broadband Microwave Spectroscopy of Pyrimidine-(CO2)n, n = 1,2.

Pyrimidine has two in-plane CH(δ+)/N̈(δ-)/CH(δ+) binding sites that are complementary to the (δ-/2δ+/δ-) quadrupole moment of CO2. We recorded broadband microwave spectra over the 7.5-17.5 GHz range for pyrimidine-(CO2)n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO2 complex is planar, with CO2 occupying one of the quadrupolar binding sites, forming a structure in which the CO2 is stabilized in the plane by interactions with the C-H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO2 complex studied previously by (Doran, J. L. J. Mol. Struct. 2012, 1019, 191-195). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C2v symmetry in which the second CO2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO2 is -13.7 kJ mol-1, including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO2)2 is almost exactly double that size, indicating little interaction between the two CO2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

Detangling Solvatochromic Effects by the Effective Fragment Potential Method.

Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.

Quantifying the Nearly Random Microheterogeneity of Aqueous tert-Butyl Alcohol Solutions Using Vibrational Spectroscopy.

The microheterogeneous structure of aqueous tert-butyl alcohol (TBA) solutions is quantified by combining experimental, simulations, and theoretical results. Experimental Raman multivariate curve resolution (Raman-MCR) C-H frequency shift measurements are compared with predictions obtained using combined quantum mechanical and effective fragment potential (QM/EFP) calculations, as well as with molecular dynamics (MD), random mixture (RM), and finite lattice (FL) predictions. The results indicate that the microheterogeneous aggregation in aqueous TBA solutions is slightly less than that predicted by MD simulations performed using either CHARMM generalized force field (CGenFF) or optimized parameters for liquid simulations all atom (OPLS-AA) force fields but slightly more than that in a self-avoiding RM of TBA-like molecules. The results imply that the onset of microheterogeneity in aqueous solutions occurs when solute contact free energies are about an order of magnitude smaller than thermal fluctuations, thus suggesting a fundamental bound of relevance to biological self-assembly.

Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example.

Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.

Ion-Size Dependent Adsorption Crossover on the Surface of a Water Droplet.

The adsorption of ionic and neutral spherical solutes on the surface of a liquid water droplet are investigated using molecular dynamics simulations and theoretical analyses. The results reveal a crossover in the sign of the adsorption free energy as a function of ion size, with ions larger than iodide predicted to be increasingly surface active. Adsorption free energies are decomposed into competing energetic and entropic contributions arising from direct solute-water interaction energy and its fluctuations. The entropically driven surface activity of large ions is predicted to increase with ion size, while small ions are typically driven away from the interface by a more delicate balance of energetic and entropic contributions, with a nonmonotonic ion size dependence linked to the ion's hydration-shell structure and stability. The physical interpretation of the results is illuminated by comparisons with dielectric linear response and cavity formation predictions and implications to interfacial acidity and enhanced chemical reactivity are discussed.

Molecular and Structural Characterization of Isomeric Compounds in Atmospheric Organic Aerosol Using Ion Mobility-Mass Spectrometry.

Secondary organic aerosol (SOA) formed through multiphase atmospheric chemistry makes up a large fraction of airborne particles. The chemical composition and molecular structures of SOA constituents vary between different emission sources and aging processes in the atmosphere, which complicates their identification. In this work, we employ drift tube ion mobility spectrometry with quadrupole time-of-flight mass spectrometry (IM-MS) detection for rapid gas-phase separation and multidimensional characterization of isomers in two biogenic SOAs produced from ozonolysis of isomeric monoterpenes, d-limonene (LSOA) and α-pinene (PSOA). SOA samples were ionized using electrospray ionization (ESI) and characterized using IM-MS in both positive and negative ionization modes. The IM-derived collision cross sections in nitrogen gas (DTCCSN2 ) for individual SOA components were obtained using multifield and single-field measurements. A novel application of IM multiplexing/high-resolution demultiplexing methodology was employed to increase sensitivity, improve peak shapes, and augment mobility baseline resolution, which revealed several isomeric structures for the measured ions. For LSOA and PSOA samples, we report significant structural differences of the isomer structures. Molecular structural calculations using density functional theory combined with the theoretical modeling of CCS values provide insights into the structural differences between LSOA and PSOA constituents. The average DTCCSN2 values for monomeric SOA components observed as [M + Na]+ ions are 3-6% higher than those of their [M - H]- counterparts. Meanwhile, dimeric and trimeric isomer components in both samples showed an inverse trend with the relevant values of [M - H]- ions being 3-7% higher than their [M + Na]+ counterparts, respectively. The results indicate that the structures of Na+-coordinated oligomeric ions are more compact than those of the corresponding deprotonated species. The coordination with Na+ occurs on the oxygen atoms of the carbonyl groups leading to a compact configuration. Meanwhile, deprotonated molecules have higher DTCCSN2 values due to their elongated structures in the gas phase. Therefore, DTCCSN2 values of isomers in SOA mixtures depend strongly on the mode of ionization in ESI. Additionally, PSOA monomers and dimers exhibit larger DTCCSN2 values (1-4%) than their LSOA counterparts owing to more rigid structures. A cyclobutane ring is present with functional groups pointing in opposite directions in PSOA compounds, as compared to noncyclic flexible LSOA structures, forming more compact ions in the gas phase. Lastly, we investigated the effects of direct photolysis on the chemical transformations of selected individual PSOA components. We use IM-MS to reveal structural changes associated with aerosol aging by photolysis. This study illustrates the detailed molecular and structural descriptors for the detection and annotation of structural isomers in complex SOA mixtures.

Copper(I)-Pyrazolate Complexes as Solid-State Phosphors: Deep-Blue Emission through a Remote Steric Effect.

We describe a novel manifestation of rigidochromic behavior in a series of tetranuclear Cu(I)-pyrazolate (Cu4pz4) macrocycles, with implications for solid-state luminescence at deep-blue wavelengths (<460 nm). The Cu4pz4 emissions are remarkably sensitive to structural effects far from the luminescent core: when 3,5-di-tert-butylpyrazoles are used as bridging ligands, adding a C4 substituent can induce a blue shift of more than 100 nm. X-ray crystal and computational analyses reveal that C4 units influence the conformational behavior of adjacent tert-butyl groups, with a subsequent impact on the global conformation of the Cu4pz4 complex. Emissions are mediated primarily through a cluster-centered triplet (3CC) state; compression of the Cu4 cluster into a nearly close-packed geometry prevents the reorganization of its excited-state structure and preserves the 3CC energy at a high level. The remote steric effect may thus offer alternative strategies toward the design of phosphors with rigid excited-state geometries.

Expulsion of Hydroxide Ions from Methyl Hydration Shells.

The affinity of hydroxide ions for methyl hydration shells is assessed using a combined experimental and theoretical analysis of tert-butyl alcohol (TBA) dissolved in pure water and aqueous NaOH and NaI. The experimental results are obtained using Raman multivariate curve resolution (Raman-MCR) and a new three-component total least squares (Raman-TLS) spectral decomposition strategy used to highlight vibrational perturbations resulting from interactions between TBA and aqueous ions. The experiments are interpreted and extended with the aid of effective fragment potential molecular dynamics (EFP-MD) simulations, as well as Kirkwood-Buff calculations and octanol/water partition measurements, to relate TBA-ion distribution functions to TBA solubility changes. The combined experimental and simulation results reveal that methyl group hydration shells more strongly expel hydroxide than iodide anions, whose populations near the methyl groups of TBA are predicted to be correlated with sodium counterion localization near the TBA hydroxyl group.

Hydration and Seamless Integration of Hydrogen Peroxide in Water.

Raman multivariate curve resolution is used to decompose the vibrational spectra of aqueous hydrogen peroxide (H2O2) into pure water, dilute H2O2, and concentrated H2O2 spectral components. The dilute spectra reveal four sub-bands in the OH stretch region, assigned to the OH stretch and Fermi resonant bend overtone of H2O2, and two nonequivalent OH groups on water molecules that donate a hydrogen bond to H2O2. At high concentrations, a spectral component resembling pure H2O2 emerges. Our results further demonstrate that H2O2 perturbs the structure of water significantly less than either methanol or sodium chloride of the same concentration, as evidenced by comparing the hydration-shell spectra of tert-butyl alcohol dissolved in the three aqueous solutions.

Interaction of Polymers with Enzalutamide Nanodroplets-Impact on Droplet Properties and Induction Times.

Amorphous solid dispersions (ASDs), which consist of a drug dispersed in a polymeric matrix, are increasingly being applied to improve the in vivo performance of poorly water-soluble drugs delivered orally. The polymer is a critical component, playing several roles including facilitating drug release from the ASD, as well as delaying crystallization from the supersaturated solution generated upon dissolution. Certain ASD formulations dissolve to produce amorphous drug-rich nanodroplets. The interaction of the polymer with these nanodroplets is poorly understood but is thought to be important for inhibiting crystallization in these systems. In this study, the impact of ionic polymers on the crystallization kinetics of enzalutamide from supersaturated solutions containing different amounts of amorphous nanodroplets was evaluated by determination of nucleation induction times. The amount of the polymer associated with the drug nanodroplets was also determined. When comparing two polymers, hydroxypropylmethyl cellulose acetate succinate (HPMCAS) and Eudragit E PO, it was found that the crystallization tendency and physical properties of the drug nanodroplets varied in the presence of these two polymers. Both polymers distributed between the aqueous phase and the drug-rich nanodroplets. A greater amount of Eudragit E PO was associated with the drug-rich nanodroplets. Despite this, Eudragit E PO was a less-effective crystallization inhibitor than HPMCAS in systems containing nanodroplets. In conclusion, in supersaturated solutions containing amorphous nanodroplets, the extent of association of a polymer with the drug nanodroplet does not solely predict crystallization inhibition.

Effective Fragment Potentials for Flexible Molecules: Transferability of Parameters and Amino Acid Database.

An accurate but efficient description of noncovalent interactions is a key to predictive modeling of biological and materials systems. The effective fragment potential (EFP) is an ab initio-based force field that provides a physically meaningful decomposition of noncovalent interactions of a molecular system into Coulomb, polarization, dispersion, and exchange-repulsion components. An EFP simulation protocol consists of two steps, preparing parameters for molecular fragments by a series of ab initio calculations on each individual fragment, and calculation of interaction energy and properties of a total molecular system based on the prepared parameters. As the fragment parameters (distributed multipoles, polarizabilities, localized wave function, etc.) depend on a fragment geometry, straightforward application of the EFP method requires recomputing parameters of each fragment if its geometry changes, for example, during thermal fluctuations of a molecular system. Thus, recomputing fragment parameters can easily become both computational and human bottlenecks and lead to a loss of efficiency of a simulation protocol. An alternative approach, in which fragment parameters are adjusted to different fragment geometries, referred to as "flexible EFP", is explored here. The parameter adjustment is based on translations and rotations of local coordinate frames associated with fragment atoms. The protocol is validated on extensive benchmark of amino acid dimers extracted from molecular dynamics snapshots of a cryptochrome protein. A parameter database for standard amino acids is developed to automate flexible EFP simulations in proteins. To demonstrate applicability of flexible EFP in large-scale protein simulations, binding energies and vertical electron ionization and electron attachment energies of a lumiflavin chromophore of the cryptochrome protein are computed. The results obtained with flexible EFP are in a close agreement with the standard EFP procedure but provide a significant reduction in computational cost.

Multiagent Consensus Equilibrium in Molecular Structure Determination.

Multiagent consensus equilibrium (MACE) is demonstrated for the integration of experimental observables as constraints in molecular structure determination and for the systematic merging of multiple computational architectures. MACE is founded on simultaneously determining the equilibrium point between multiple experimental and/or computational agents; the returned state description (e.g., atomic coordinates for molecular structure) represents the intersection of each manifold and is not equivalent to the average optimum state for each agent. The moment of inertia, determined directly from microwave spectroscopy measurements, serves to illustrate the mechanism through which MACE evaluations merge experimental and quantum chemical modeling. MACE results reported combine gradient descent optimization of each ab initio agent with an agent that predicts the chemical structure based on root-mean-square deviation of the predicted inertia tensor with experimentally measured moments of inertia. Successful model fusion for several small molecules was achieved as well as the larger molecule solketal. Fusing a model of moment of inertia, an underdetermined predictor of structure, with low cost computational methods yielded structure determination performance comparable to standard computational methods such as MP2/cc-pVTZ and greater agreement with experimental observables.

Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding.

Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.

Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex.

High efficiency of light harvesting in photosynthetic pigment-protein complexes is governed by evolutionary-perfected protein-assisted tuning of individual pigment properties and interpigment interactions. Due to the large number of spectrally overlapping pigments in a typical photosynthetic complex, experimental methods often fail to unambiguously identify individual chromophore properties. Here, we report a first-principles-based modeling protocol capable of predicting properties of pigments in protein environment to a high precision. The technique was applied to successfully uncover electronic properties of the Fenna-Matthews-Olson (FMO) pigment-protein complex. Each of the three subunits of the FMO complex contains eight strongly coupled bacteriochlorophyll a (BChl a) pigments. The excitonic structure of FMO can be described by an electronic Hamiltonian containing excitation (site) energies of BChl a pigments and electronic couplings between them. Several such Hamiltonians have been developed in the past based on the information from various spectroscopic measurements of FMO; however, fine details of the excitonic structure and energy transfer in FMO, especially assignments of short-lived high-energy sites, remain elusive. Utilizing polarizable embedding quantum mechanics/molecular mechanics with the effective fragment potentials, we computed the electronic Hamiltonian of FMO that is in general agreement with previously reported empirical Hamiltonians and quantitatively reproduces experimental absorption and circular dichroism spectra of the FMO protein. The developed computational protocol is sufficiently simple and can be utilized for predictive modeling of other wild-type and mutated photosynthetic pigment-protein complexes.

FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method.

In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.

The effects of site asymmetry on near-degenerate state-to-state vibronic mixing in flexible bichromophores.

Laser-induced fluorescence excitation and dispersed fluorescence spectra of a model flexible bichromophore, 1,1-diphenylethane (DPE), have been recorded under jet-cooled conditions in the gas phase in the region near the first pair of near-degenerate excited states (S1 and S2). The S1 and S2 origin transitions have been identified at 37 397 and 37 510 cm-1, a splitting of 113 cm-1. This splitting is four times smaller than the excitonic splitting calculated by ab initio methods at the EOM-CCSD/cc-pVDZ level of theory (410 cm-1), which necessarily relies on the Born-Oppenheimer approximation. Dispersed fluorescence spectra provide a state-to-state picture of the vibronic coupling. These results are compared with the results of a multimode vibronic coupling model capable of treating chromophores in asymmetric environments. This model was used to predict the splitting between S1 and S2 origins close to the experiment, reduced from its pure excitonic value by Franck-Condon quenching. Quantitative accuracy is achieved by the model, lending insight into the state-to-state mixing that occurs between individual S1 and S2 vibronic levels. The S2 origin is determined to be mixed with S1(v) levels by two mechanisms common to internal conversion in almost any setting; namely, (i) mixing involving near-degenerate levels with large vibrational quantum number changes that are not governed by Δv = 1 Herzberg-Teller (HT) selection rules, and (ii) mixing with levels with larger energy gaps that do follow these selection rules. In DPE, the asymmetric ring flapping vibrational mode R¯ dominates the HT coupling.

Polarizable embedding for simulating redox potentials of biomolecules.

Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana. The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (-0.15 V).