Luca De Vico

Software news and update MOLCAS 7 : The Next Generation

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.

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

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.

Theoretical study of the chemiluminescent decomposition of dioxetanone.

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.

New general tools for constrained geometry optimizations

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.

Quantifying signal changes in nano-wire based biosensors

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.

Effects of buffer composition and dilution on nanowire field-effect biosensors

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.

A Computational Methodology to Screen Activities of Enzyme Variants

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.

Multireference Excitation Energies for Bacteriochlorophylls A within Light Harvesting System 2.

Light-harvesting system 2 (LH2) of purple bacteria is one of the most popular antenna complexes used to study Natures 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 Silico Prediction of Mutant HIV-1 Proteases Cleaving a Target Sequence

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.

Intermolecular Modes between LH2 Bacteriochlorophylls and Protein Residues: The Effect on the Excitation Energies

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.