Search

Molecular excitons play a major role within dye aggregates and hold significant potential for (opto)electronic and photovoltaic applications. Numerous studies have documented alterations in the spectral properties of dye homoaggregates, but only limited work has been reported for heteroaggregates. In this article, dimeric dye stacks were constructed from azobenzene-like dyes with identical or distinct structures, and their excitonic features were computationally investigated. Our results show that strong exciton coupling is not limited to identical chromophores, as often assumed, based on a recently made available Frenkel Exciton Hamiltonian and multiconfigurational plus second-order perturbation theory energetics methodology. Heteroaggregate stacks were found to exhibit different absorption features from the corresponding interacting monomers, indicating considerable coupling interactions between units. We analyzed how such coupling may vary according to various aspects, such as the relative positions of the interacting monomers or the differences in their energetics. Such qualitative and semiquantitative analyses allow the evaluation of the excitonic behavior of these dye aggregates to encourage further efforts toward a deeper understanding of the excitonic properties of tailored dye heteroaggregate systems.

Oxygenic photosynthesis begins in the reaction center (RC) of the protein complex photosystem II (PSII). PSII has an intriguing, nearly symmetrical arrangement of cofactors within its RC. Despite this symmetry, evolution has favored only one of the two branches of PSII for efficient electron transfer. Current spectroscopic experiments explore the electronic dynamics during the picoseconds after energy has entered the RC and until the electron transfers to the pheophytin of the first branch. We present state-of-the-art multiconfigurational multireference calculations of the excitation energies or site energies of the four chlorophyll pigments of the RC without protein environment considerations. We see a significant variation that breaks the apparent symmetry of the RC. The inner chlorophyll of the productive RC branch possessed the lowest excitation energy of the four central chlorophylls. Our computational method used here is expensive; thus, geometry optimization of the crystal structure is currently not possible. In future work, charge and energy dynamics within the RC will be included as well as a dynamic description of the protein environment and its coupling to the RC. Other state-of-the-art studies of the RC, at lower levels of electronic structure, include a static treatment of the protein environment. These almost unanimously report that the outer chlorophyll of the active branch had the lowest excitation energy. Future work is needed to reconcile this discrepancy., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.

We present a computational protocol for the fast and automated screening of excited-state hybrid quantum mechanics/molecular mechanics (QM/MM) models of rhodopsins to be used as fluorescent probes based on the automatic rhodopsin modeling protocol ( a -ARM). Such " a -ARM fluorescence screening protocol" is implemented through a general Python-based driver, PyARM, that is also proposed here. The implementation and performance of the protocol are benchmarked using different sets of rhodopsin variants whose absorption and, more relevantly, emission spectra have been experimentally measured. We show that, despite important limitations that make unsafe to use it as a black-box tool, the protocol reproduces the observed trends in fluorescence and it is capable of selecting novel potentially fluorescent rhodopsins. We also show that the protocol can be used in mechanistic investigations to discern fluorescence enhancement effects associated with a near degeneracy of the S 1 /S 2 states or, alternatively, with a barrier generated via coupling of the S 0 /S 1 wave functions.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter., (© 2022. The Author(s).)

In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties., (© 2022. The Author(s).)

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments., (© 2021 American Society for Photobiology.)

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

This article introduces Web-ARM, a specialized tool, online available, designed to build quantum mechanical/molecular mechanical models of rhodopsins, a widely spread family of light-responsive proteins. Web-ARM allows the rapidly building of models of rhodopsins with a documented quality and the prediction of trends in UV-vis absorption maximum wavelengths, based on their excitation energies computed at the CASPT2//CASSCF/Amber level of theory. Web-ARM builds upon the recently reported, python-based a -ARM protocol [ J. Chem. Theory Comput. , 2019 , 15, 3134-3152] and, as such, necessitates only a crystallographic structure or a comparative model in PDB format and a very basic knowledge of the studied rhodopsin system. The user-friendly web interface uses such input to generate congruous, gas-phase models of rhodopsins and, if requested, their mutants. We present two possible applications of Web-ARM, which showcase how the interface can be employed to assist both research and educational activities in fields at the interface between chemistry and biology. The first application shows how, through Web-ARM, research projects (e.g., rhodopsin and rhodopsin mutant screening) can be carried out in significantly less time with respect to using the required computational photochemistry tools via a command line. The second application documents the use of Web-ARM in a real-life educational/training activity, through a hands-on experience illustrating the concepts of rhodopsin color tuning.

Light-harvesting systems 2 and 3 (LH2 and LH3) act as antennas for the initial light capture by photosynthetic purple bacteria, thus initiating the conversion of solar energy into chemical energy. The main absorbers are carotenoids and bacteriochlorophylls (BChls), which harvest different parts of the solar spectrum. The first two optical transitions in BChl produce the Q y and Q x absorption bands. The large size of BChl molecules has prevented accurate computational determination of the electronic structures for the relevant states, until we recently succeeded in obtaining the excitation energies and transition dipole moments of the first (Q y ) transition for all BChls in LH2 and LH3 using multi-state multiconfigurational second-order perturbation theory calculations. In this work, we go one step further, compute the corresponding values for the Q x transition, in line with previous work [ J. Am. Chem. Soc . 2017 , 139 , 7558 - 7567 ], and compare and assess our data against excitation energies obtained through time-dependent density functional theory methods. Interestingly, we find that the two transitions respond differently to BChls' geometrical factors, such as the macrocycle ring curvature and the dihedral torsion of the acetyl moiety. These findings will aid the unraveling of structure-function relationships for absorption and energy transfer processes in purple bacteria, and once again this demonstrates the viability of multireference quantum chemical methods as computational tools for the photophysics of biomolecules.

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol -1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol -1 (0.04 eV).

A methodology for the automatic production of quantum mechanical/molecular mechanical (QM/MM) models of retinal-binding rhodopsin proteins and subsequent prediction of their spectroscopic properties has been proposed recently by some of the authors. The technology employed for the evaluation of the excitation energies is called Automatic Rhodopsin Modeling (ARM), and it involves the use of the complete active space self-consistent field (CASSCF) method followed by a multiconfiguration second-order perturbation theory (in particular, CASPT2) calculation of external correlation energies. Although it was shown that ARM is capable of successfully reproducing and predicting spectroscopic property trends in chromophore-embedding protein sets, practical applications of such technology are limited by the high computational costs of the multiconfiguration perturbation theory calculations. In the present work we benchmark the more affordable multiconfiguration pair-density functional theory (MC-PDFT) method whose accuracy has been recently validated for retinal chromophores in the gas phase, indicating that MC-PDFT could potentially be used to analyze large (e.g., few hundreds) sets of rhodopsin proteins. Here, we test this theory for a set of rhodopsin QM/MM models whose experimental absorption maxima (λ a max ) have been measured. The results indicate that MC-PDFT may be employed to calculate λ a max values for this important class of photoresponsive proteins.

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

Natural light-harvesting is performed by pigment-protein complexes, which collect and funnel the solar energy at the start of photosynthesis. The identity and arrangement of pigments largely define the absorption spectrum of the antenna complex, which is further regulated by a palette of structural factors. Small alterations are induced by pigment-protein interactions. In light-harvesting systems 2 and 3 from Rhodoblastus acidophilus , the pigments are arranged identically, yet the former has an absorption peak at 850 nm that is blue-shifted to 820 nm in the latter. While the shift has previously been attributed to the removal of hydrogen bonds, which brings changes in the acetyl moiety of the bacteriochlorophyll, recent work has shown that other mechanisms are also present. Using computational and modeling tools on the corresponding crystal structures, we reach a different conclusion: The most critical factor for the shift is the curvature of the macrocycle ring. The bending of the planar part of the pigment is identified as the second-most important design principle for the function of pigment-protein complexes-a finding that can inspire the design of novel artificial systems., Competing Interests: The authors declare no conflict of interest.

Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.

Over the last decade, we have investigated and exploited the photophysical properties of triangulenium dyes. Azadioxatriangulenium (ADOTA) and diazaoxatriangulenium (DAOTA), in particular, have features that make them useful in various fluorescence-based technologies (e.g., bioimaging). Through our work with ADOTA and DAOTA, we became aware that the reported fluorescence quantum yields (ϕ fl ) for these dyes are lower than their actual values. We thus set out to further investigate the fundamental structure-property relationships in these unique conjugated cationic systems. The nonradiative processes in the systems were explored using transient absorption spectroscopy and time-resolved emission spectroscopy in combination with computational chemistry. The influence of molecular oxygen on the fluorescence properties was explored, and the singlet oxygen sensitization efficiencies of ADOTA and DAOTA were determined. We conclude that, for these dyes, the amount of nonradiative deactivation of the first excited singlet state (S 1 ) of the azaoxa-triangulenium fluorophores is low, that the rate of such deactivation is slower than what is observed in common cationic dyes, that there are no observable radiative transitions occurring from the first excited triplet state (T 1 ) of these dyes, and that the efficiency of sensitized singlet oxygen production is low (ϕ Δ ≤ 10%). These photophysical results provide a solid base upon which technological applications of these fluorescent dyes can be built., Competing Interests: The authors declare no competing financial interest.

The modeling of vibrations in optical spectra relies heavily on the simplifications brought about by using harmonic oscillators. However, realistic molecular systems can deviate substantially from this description. We develop two methods which show that the extension to arbitrarily shaped potential energy surfaces is not only straightforward, but also efficient. These methods are applied to an electronic two-level system with potential energy surfaces of polynomial form and used to study anharmonic features such as the zero-phonon line shape and mirror-symmetry breaking between absorption and fluorescence spectra. The first method, which constructs vibrational wave functions as linear combinations of the harmonic oscillator wave functions, is shown to be extremely robust and can handle large anharmonicities. The second method uses the cumulant expansion, which is readily solved, even at high orders, thanks to an ideally suited matrix theorem.

Benign prostatic hyperplasia is a common disease in men aged over 50 years old, with an incidence increasing to more than 80% over the age of 70, that is increasingly going to attract pharmaceutical interest. Within conventional therapies, such as α -adrenoreceptor antagonists and 5 α -reductase inhibitor, there is a large requirement for treatments with less adverse events on, e.g., blood pressure and sexual function: phytotherapy may be the right way to fill this need. Serenoa repens standardized extract has been widely studied and its ability to reduce lower urinary tract symptoms related to benign prostatic hyperplasia is comprehensively described in literature. An innovative investigation on the mechanism of inhibition of 5 α -reductase by Serenoa repens extract active principles is proposed in this work through computational methods, performing molecular docking simulations on the crystal structure of human liver 5 β -reductase. The results confirm that both sterols and fatty acids can play a role in the inhibition of the enzyme, thus, suggesting a competitive mechanism of inhibition. This work proposes a further confirmation for the rational use of herbal products in the management of benign prostatic hyperplasia, and suggests computational methods as an innovative, low cost, and non-invasive process for the study of phytocomplex activity toward proteic targets., Competing Interests: The authors declare there are no competing interests.

In this study a number of chromophores based on boron subphthalocyanines are investigated for use in the future design of organic photovoltaic devices based on molecular triad systems. The computational study is performed at the TD-DFT CAM-B3LYP/6-311G(d) level of theory. The absorption spectra of these chromophores are simulated using TD-DFT and compared to experimental results. All investigated chromophores absorb light in the visible range and thus are suitable for absorption of sunlight in solar cell applications. On the basis of energy-level alignments, suitable combinations of moieties for a molecular triad system are proposed. The molecular triads will be used in future work as the functional part of organic photovoltaic devices, where the chromophore will be used both to absorb the incoming solar radiation and to increase the distance between the separated charges on donor and acceptor units to increase the lifetime of the charge-separated state.

Light-harvesting system 2 (LH2) of purple bacteria is one of the most popular antenna complexes used to study Nature's way of collecting and channeling solar energy. The dynamics of the absorbed energy is probed by ultrafast spectroscopy. Simulation of these experiments relies on fitting a range of parameters to reproduce the spectra. Here, we present a method that can determine key parameters to chemical accuracy. These will eliminate free variables in the modeling, thus reducing the problem. Using MS-RASPT2/RASSCF calculations, we compute excitation energies and transition dipole moments of all bacteriochlorophylls in LH2. We find that the excitation energies vary among the bacteriochlorophyll monomers and that they are regulated by the curvature of the macrocycle ring and the dihedral angle of an acetyl moiety. Increasing the curvature lifts the ground state energy, which causes a red shift of the excitation energy. Increasing the torsion of the acetyl moiety raises the excited state energy, resulting in a blue shift of the excitation energy. The obtained results mark a giant leap for multiconfigurational multireference quantum chemical methods in the photochemistry of biological systems, which can prove instrumental in exposing the underlying physics of photosynthetic light-harvesting.

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization., (© 2015 Wiley Periodicals, Inc.)

HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636-1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.

Our previously presented method for high throughput computational screening of mutant activity (Hediger et al., 2012) is benchmarked against experimentally measured amidase activity for 22 mutants of Candida antarctica lipase B (CalB). Using an appropriate cutoff criterion for the computed barriers, the qualitative activity of 15 out of 22 mutants is correctly predicted. The method identifies four of the six most active mutants with ≥3-fold wild type activity and seven out of the eight least active mutants with ≤0.5-fold wild type activity. The method is further used to screen all sterically possible (386) double-, triple- and quadruple-mutants constructed from the most active single mutants. Based on the benchmark test at least 20 new promising mutants are identified.

We present a semi-empirical (PM6-based) computational method for systematically estimating the effect of all possible single mutants, within a certain radius of the active site, on the barrier height of an enzymatic reaction. The intent of this method is not a quantitative prediction of the barrier heights, but rather to identify promising mutants for further computational or experimental study. The method is applied to identify promising single and double mutants of Bacillus circulans xylanase (BCX) with increased hydrolytic activity for the artificial substrate ortho-nitrophenyl β-xylobioside (ONPX2). The estimated reaction barrier for wild-type (WT) BCX is 18.5 kcal/mol, which is in good agreement with the experimental activation free energy value of 17.0 kcal/mol extracted from the observed k cat using transition state theory (Joshi et al., 2001). The PM6 reaction profiles for eight single point mutations are recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase in barrier height for all eight mutants while FMO predicts an increase for six of the eight mutants. Both methods predict that the largest change in barrier occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase, respectively. We thus conclude that PM6 is sufficiently accurate to identify promising mutants for further study. We prepared a set of all theoretically possible (342) single mutants in which every amino acid of the active site (except for the catalytically active residues E78 and E172) was mutated to every other amino acid. Based on results from the single mutants we construct a set of 111 double mutants consisting of all possible pairs of single mutants with the lowest barrier for a particular position and compute their reaction profile. None of the mutants have, to our knowledge, been prepared experimentally and therefore present experimentally testable predictions.

Nanowire-based field-effect transistors (FETs) can be used as ultra-sensitive and label-free biosensors for detecting protein-protein interactions. A way to increase the performance of such sensors is to dilute the sensing buffer drastically. However, we show here that this can have an important effect on the function of the proteins. Moreover, it is demonstrated that this dilution significantly affects the pH stability of the sensing buffer, which consequently impacts the charge of the protein and thus the response and signal-to-noise ratio in the sensing experiments. Three model systems are investigated experimentally to illustrate the impact on ligand-protein and protein-protein interactions. Simulations are performed to illustrate the effect on the performance of the sensors. Combining various parameters, the current study provides a means for evaluating and selecting the most appropriate buffer composition for bioFET measurements.

We present a fast computational method to efficiently screen enzyme activity. In the presented method, the effect of mutations on the barrier height of an enzyme-catalysed reaction can be computed within 24 hours on roughly 10 processors. The methodology is based on the PM6 and MOZYME methods as implemented in MOPAC2009, and is tested on the first step of the amide hydrolysis reaction catalyzed by the Candida Antarctica lipase B (CalB) enzyme. The barrier heights are estimated using adiabatic mapping and shown to give barrier heights to within 3 kcal/mol of B3LYP/6-31G(d)//RHF/3-21G results for a small model system. Relatively strict convergence criteria (0.5 kcal/(molÅ)), long NDDO cutoff distances within the MOZYME method (15 Å) and single point evaluations using conventional PM6 are needed for reliable results. The generation of mutant structures and subsequent setup of the semiempirical calculations are automated so that the effect on barrier heights can be estimated for hundreds of mutants in a matter of weeks using high performance computing.

We present a web interface which allows us to conveniently set up calculations based on the BioFET-SIM model. With the interface, the signal of a BioFET sensor can be calculated depending on its parameters, as well as the signal dependence on pH. As an illustration, two case studies are presented. In the first case, a generic peptide with opposite charges on both ends is inverted in orientation on a semiconducting nanowire surface leading to a corresponding change in sign of the computed sensitivity of the device. In the second case, the binding of an antibody/antigen complex on the nanowire surface is studied in terms of orientation and analyte/nanowire surface distance. We demonstrate how the BioFET-SIM web interface can aid in the understanding of experimental data and postulate alternative ways of antibody/antigen orientation on the nanowire surface.

A single charge screening model of surface charge sensors in liquids (De Vico et al., Nanoscale, 2011, 3, 706-717) is extended to multiple charges to model the effect of the charge distributions of analyte proteins on FET sensor response. With this model we show that counter-intuitive signal changes (e.g. a positive signal change due to a net positive protein binding to a p-type conductor) can occur for certain combinations of charge distributions and Debye lengths. The new method is applied to interpret published experimental data on Streptavidin (Ishikawa et al., ACS Nano, 2009, 3, 3969-3976) and Nucleocapsid protein (Ishikawa et al., ACS Nano, 2009, 3, 1219-1224).

In this work, we present a computational methodology for predicting the change in signal (conductance sensitivity) of a nano-BioFET sensor (a sensor based on a biomolecule binding another biomolecule attached to a nano-wire field effect transistor) upon binding its target molecule. The methodology is a combination of the screening model of surface charge sensors in liquids developed by Brandbyge and co-workers [Sørensen et al., Appl. Phys. Lett., 2007, 91, 102105], with the PROPKA method for predicting the pH-dependent charge of proteins and protein-ligand complexes, developed by Jensen and co-workers [Li et al., Proteins: Struct., Funct., Bioinf., 2005, 61, 704-721, Bas et al., Proteins: Struct., Funct., Bioinf., 2008, 73, 765-783]. The predicted change in conductance sensitivity based on this methodology is compared to previously published data on nano-BioFET sensors obtained by other groups. In addition, the conductance sensitivity dependence from various parameters is explored for a standard wire, representative of a typical experimental setup. In general, the experimental data can be reproduced with sufficient accuracy to help interpret them. The method has the potential for even more quantitative predictions when key experimental parameters (such as the charge carrier density of the nano-wire or receptor density on the device surface) can be determined (and reported) more accurately.

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies., (Copyright 2009 Wiley Periodicals, Inc.)

The unimolecular chemiluminescent decomposition of unsubstituted dioxetanone was studied at the complete active space self-consistent field level of theory combined with the multistate second-order multiconfigurational perturbation theory energy correction. The calculations revealed interesting features. Two transition states, two conical intersections, and one intermediate stable biradical structure along the lowest energy reaction path were identified. It was noted that the conical intersections are found at or in very close proximity to the transition states. The first and second transition states correspond to O-O and C-C cleavages, respectively. In particular, a planar structure is supported by the (1)(sigma,sigma*) state during the O-O dissociation up to the first transition state and conical intersection. At this point the (1)(sigma,sigma*) state dissociation path bifurcates, corresponding to a torsion of the O-C-C-O angle. Simultaneously, the (1)(n,sigma*) state becomes lower in energy while still favoring a planar structure. As the lowest-energy reaction path proceeds toward the second transition state and conical intersection, the (1)(n,sigma*), (3)(n,sigma*), and (1)(sigma,sigma*) states are close in energy. This work suggests that the vibrational distribution at the first conical intersection and the interactions among the states as the reaction proceeds between the two transition states are the origin of the population of the chemiluminescent (n,sigma*) states.

Photochemical reactions rationalization is a key aspect for the understanding and setup of novel experiment and novel photoinitiated pathways. In this respect, the relationship between minimum energy paths over an excited-state and the intersection to lower potential energy surfaces is fundamental. In order to help the understanding of this relationship, in this study we present a novel kind of constraint for geometry optimizations, namely, an "orthogonality" constraint. Its possible applications are described. A complete example on how to retrieve the direct relationship between a minimum energy path over an excited-state potential energy surface and a conical intersection seam is given for C2v symmetry constrained formaldehyde. The advantages of using the novel constraint when rationalizing a (photo)chemical reaction are presented.

The UV photodissociation (<5 eV) of diiodomethane (CH(2)I(2)) is investigated by spin-orbit ab initio calculations. The experimentally observed photodissociation channels in the gas and condensed phases are clearly assigned by multi-state second-order multiconfigurational perturbation theory in conjunction with spin-orbit interaction through complete active space-state interaction potential energy curves. The calculated results indicate that the fast dissociations of the first two singlet states of CH(2)I(2) and CH(2)I--I lead to geminate-radical products, CH(2)I (.)+I((2)P(3/2)) or CH(2)I (.)+ I*((2)P(1/2)). The recombination process from CH(2)I--I to CH(2)I(2) is explained by an isomerization process and a secondary photodissociation reaction of CH(2)I--I. Finally, the study reveals that spin-orbits effects are significant in the quantitative analysis of the electronic spectrum of the CH(2)I--I species.

A modification of the constrained geometry optimization method by Anglada and Bofill (Anglada, J. M.; Bofill, J. M. J. Comput. Chem. 1997, 18, 992-1003) is designed and implemented. The changes include the choice of projection, quasi-line-search, and the use of a Rational Function optimization approach rather than a reduced-restricted-quasi-Newton-Raphson method in the optimization step. Furthermore, we show how geometrical constrains can be implemented in an approach based on nonredundant curvilinear coordinates avoiding the inclusion of the constraints in the set of redundant coordinates used to define the internal coordinates. The behavior of the new implementation is demonstrated in geometry optimizations featuring single or multiple geometrical constraints (bond lengths, angles, etc.), optimizations on hyperspherical cross sections (as in the computation of steepest descent paths), and location of energy minima on the intersection subspace of two potential energy surfaces (i.e. minimum energy crossing points). In addition, a novel scheme to determine the crossing point geometrically nearest to a given molecular structure is proposed.

CASPT2//CASSCF/6-31G photochemical reaction path computations for two 4-cis-nona-2,4,6,8-tetraeniminium cation derivatives, with the 4-cis double bond embedded in a seven- and eight-member ring, are carried out to model the reactivity of the corresponding ring-locked retinal chromophores. The comparison of the excited state branches of the two reaction paths with that of the native chromophore, is used to unveil the factors responsible for the remarkably short (60 fs) excited state (S(1)) lifetime observed when an artificial rhodopsin containing an eight member ring-locked retinal is photoexcited. Indeed, it is shown that the strain imposed by the eight-member ring on the chromophore backbone leads to a dramatic change in the shape of the S(1) energy surface. Our models are also used to investigate the nature of the primary photoproducts observed in different artificial rhodopsins. It is seen that only the eight member ring-locked retinal model can access a shallow energy minimum on the ground state. This result implies that the primary, photorhodopsin-like, transient observed in artificial rhodopsins could correspond to a shallow excited state minimum. Similarly, the second, bathorhodopsin-like, transient species could be assigned to a ground state structure displaying a nearly all-trans conformation.

Thio amino acids can be integrated into the backbone of peptides without significantly perturbing their structure. In this contribution we use ultrafast infrared and visible spectroscopy as well as state-of-the-art ab initio computations to investigate the photoisomerization of the trans form of N-methylthioacetamide (NMTAA) as a model conformational photoswitch. Following the S2 excitation of trans-NMTAA in water, the return of the molecule into the trans ground state and the formation of the cis isomer is observed on a dual time scale, with a fast component of 8-9 ps and a slow time constant of approximately 250 ps. On both time scales the probability of isomerization to the cis form is found to be 30-40%, independently of excitation wavelength. Ab initio CASPT2//CASSCF photochemical reaction path calculations indicate that, in vacuo, the trans-->cis isomerization event takes place on the S1 and/or T1 triplet potential energy surfaces and is controlled by very small energy barriers, in agreement with the experimentally observed picosecond time scale. Furthermore, the calculations identify one S2/S1 and four nearly isoenergetic S1/S0 conical intersection decay channels. In line with the observed isomerization probability, only one of the S1/S0 conical intersections yields the cis conformation upon S1-->S0 decay. A substantially equivalent excited-state relaxation results from four T1/S0 intersystem crossing points.

The results of new and recently reported CASSCF/6-31G* photoisomerization path computations of a series of models of the 11-cis retinal chromophore of the visual pigment rhodopsin are discussed. The results indicate that, with respect to the chromophore in vacuo, certain structural, intramolecular and environmental factors are capable of speeding up the excited-state decay associated with the cis --> trans isomerization motion. Using suitable protonated Schiff-base models, it is shown that three structural factors can potentially speed up the isomerization: (i) reducing the length of the conjugated chain, (ii) twisting of the hydrocarbon end of the conjugated chain with respect to the protonated Schiff-base end and (iii) ring locking of the conjugated chain with an eight-membered ring. All these factors operate through increasing the slope of the excited-state energy surface and enhancing the coupling between stretching and torsional modes. We argue that the protein catalysis seen in rhodopsin may, at least partly, exploit the same principles.

Multiconfigurational second-order perturbation theory computations and reaction path mapping for the retinal protonated Schiff base models all-trans-nona-2,4,6,8-tetraeniminium and 2-cis-nona-2,4,6,8-tetraeniminium cation demonstrate that, in isolated conditions, retinal chromophores exhibit at least three competing excited-state double bond isomerization paths. These paths are associated with the photoisomerization of the double bonds in positions 9, 11, and 13, respectively, and are controlled by barriers that favor the position 11. The computations provide a basis for the understanding of the observed excited-state lifetime in both naturally occurring and synthetic chromophores in solution and, tentatively, in the protein environment. In particular, we provide a rationalization of the excited-state lifetimes observed for a group of locked retinal chromophores which suggests that photoisomerization in bacteriorhodopsin is the result of simultaneous specific "catalysis" (all-trans --> 13-cis path) accompanied by specific "inhibition" (all-trans --> 11-cis path). The nature of the S(1) --> S(0) decay channel associated with the three paths has also been investigated at the CASSCF level of theory. It is shown that the energy surfaces in the vicinity of the conical intersection for the photoisomerization about the central double bond of retinal (position 11) and the two corresponding lateral double bonds (positions 9 and 13) are structurally different.