Gerrit Groenhof

GROMACS: Fast, flexible, and free

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum‐chemical packages (MOPAC, GAMES‐UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

Femtosecond structural dynamics drives the trans/cis isomerization in photoactive yellow protein.

Visualizing a response to light Many biological processes depend on detecting and responding to light. The response is often mediated by a structural change in a protein that begins when absorption of a photon causes isomerization of a chromophore bound to the protein. Pande et al. used x-ray pulses emitted by a free electron laser source to conduct time-resolved serial femtosecond crystallography in the time range of 100 fs to 3 ms. This allowed for the real-time tracking of the trans-cis isomerization of the chromophore in photoactive yellow protein and the associated structural changes in the protein. Science, this issue p. 725 The trans-to-cis isomerization of a key chromophore is characterized on ultrafast time scales. A variety of organisms have evolved mechanisms to detect and respond to light, in which the response is mediated by protein structural changes after photon absorption. The initial step is often the photoisomerization of a conjugated chromophore. Isomerization occurs on ultrafast time scales and is substantially influenced by the chromophore environment. Here we identify structural changes associated with the earliest steps in the trans-to-cis isomerization of the chromophore in photoactive yellow protein. Femtosecond hard x-ray pulses emitted by the Linac Coherent Light Source were used to conduct time-resolved serial femtosecond crystallography on photoactive yellow protein microcrystals over a time range from 100 femtoseconds to 3 picoseconds to determine the structural dynamics of the photoisomerization reaction.

Visualizing a protein quake with time-resolved X-ray scattering at a free-electron laser

We describe a method to measure ultrafast protein structural changes using time-resolved wide-angle X-ray scattering at an X-ray free-electron laser. We demonstrated this approach using multiphoton excitation of the Blastochloris viridis photosynthetic reaction center, observing an ultrafast global conformational change that arises within picoseconds and precedes the propagation of heat through the protein. This provides direct structural evidence for a protein quake: the hypothesis that proteins rapidly dissipate energy through quake-like structural motions.

Signal transduction in the photoactive yellow protein. I. Photon absorption and the isomerization of the chromophore

Molecular dynamics simulation techniques together with time‐dependent density functional theory calculations have been used to investigate the effect of photon absorption by a 4‐hydroxy‐cinnamic acid chromophore on the structural properties of the photoactive yellow protein (PYP) from Ectothiorodospira halophila. The calculations suggest that the protein not only modifies the absorption spectrum of the chromophore but also regulates the subsequent isomerization of the chromophore by stabilizing the isomerization transition state. Although signaling from PYP is thought to involve partial unfolding of the protein, the mechanical effects accompanying isomerization do not appear to directly destabilize the protein. Proteins 2002;48:202–211.

Signal transduction in the photoactive yellow protein. II. Proton transfer initiates conformational changes

Molecular dynamics simulation techniques, together with semiempirical PM3 calculations, have been used to investigate the effect of photoisomerization of the 4‐hydroxy‐cinnamic acid chromophore on the structural properties of the photoactive yellow protein (PYP) from Ectothiorodospira halophila. In this bacteria, exposure to blue light leads to a negative photoactic response. The calculations suggest that the isomerization does not directly destabilize the protein. However, because of the isomerization, a proton transfer from a glutamic acid residue (Glu46) to the phenolate oxygen atom of the chromophore becomes energetically favorable. The proton transfer initiates conformational changes within the protein, which are in turn believed to lead to signaling. Proteins 2002;48:212–219.

Combined QM/MM Investigation on the Light-Driven Electron-Induced Repair of the (6-4) Thymine Dimer Catalyzed by DNA Photolyase

The (6-4) photolyases are blue-light-activated enzymes that selectively bind to DNA and initiate splitting of mutagenic thymine (6-4) thymine photoproducts (T(6-4)T-PP) via photoinduced electron transfer from flavin adenine dinucleotide anion (FADH(-)) to the lesion triggering repair. In the present work, the repair mechanism after the initial electron transfer and the effect of the protein/DNA environment are investigated theoretically by means of hybrid quantum mechanical/molecular mechanical (QM/MM) simulations using X-ray structure of the enzyme-DNA complex. By comparison of three previously proposed repair mechanisms, we found that the lowest activation free energy is required for the pathway in which the key step governing the repair photocycle is electron transfer coupled with the proton transfer from the protonated histidine, His365, to the N3 nitrogen of the pyrimidone thymine. The transfer simultaneously occurs with concerted intramolecular OH transfer without formation of an oxetane or isolated water molecule intermediate. In contrast to previously suggested mechanisms, this newly identified pathway requires neither a subsequent two-photon process nor electronic excitation of the photolesion.

New QM/MM implementation of the DFTB3 method in the gromacs package

The approximate density‐functional tight‐binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle–mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free‐energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron‐induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.

Thermodynamics of hydronium and hydroxide surface solvation

The concentration of hydronium and hydroxide at the water–air interface has been debated for a long time. Recent evidence from a range of experiments and theoretical calculations strongly suggests the water surface to be somewhat acidic. Using novel polarizable models we have performed potential of mean force calculations of a hydronium ion, a hydroxide ion and a water molecule in a water droplet and a water slab and we were able to rationalize that hydronium, but not hydroxide, is slightly enriched at the surface for two reasons. First, because the hydrogen bond acceptance capacity of hydronium is weaker than water and it is more favorable to have the hydronium oxygen on the surface. Second, hydroxide ions are expelled from the surface of the droplets, due to the entropy being lower when a hydroxide ion is hydrated on the surface. As a result, the water dissociation constant pKw increases slightly near the surface. The results are corroborated by calculations of surface tension of NaOH solutions that are in reasonable agreement with the experiment. The structural and thermodynamic interpretation of hydronium and hydroxide hydration provided by these calculations opens the route to a better understanding of atmospheric and surface chemistry.

The room temperature crystal structure of a bacterial phytochrome determined by serial femtosecond crystallography

Phytochromes are a family of photoreceptors that control light responses of plants, fungi and bacteria. A sequence of structural changes, which is not yet fully understood, leads to activation of an output domain. Time-resolved serial femtosecond crystallography (SFX) can potentially shine light on these conformational changes. Here we report the room temperature crystal structure of the chromophore-binding domains of the Deinococcus radiodurans phytochrome at 2.1u2009Å resolution. The structure was obtained by serial femtosecond X-ray crystallography from microcrystals at an X-ray free electron laser. We find overall good agreement compared to a crystal structure at 1.35u2009Å resolution derived from conventional crystallography at cryogenic temperatures, which we also report here. The thioether linkage between chromophore and protein is subject to positional ambiguity at the synchrotron, but is fully resolved with SFX. The study paves the way for time-resolved structural investigations of the phytochrome photocycle with time-resolved SFX.

Hydrogen Bond Fluctuations Control Photochromism in a Reversibly Photo-Switchable Fluorescent Protein.

Reversibly switchable fluorescent proteins (RSFPs) are essential for high-resolution microscopy of biological samples, but the reason why these proteins are photochromic is still poorly understood. To address this problem, we performed molecular dynamics simulations of the fast switching Met159Thr mutant of the RSFP Dronpa. Our simulations revealed a ground state structural heterogeneity in the chromophore pocket that consists of three populations with one, two, or three hydrogen bonds to the phenolate moiety of the chromophore. By means of non-adiabatic quantum mechanics/molecular dynamics simulations, we demonstrated that the subpopulation with a single hydrogen bond is responsible for off-switching through photo-isomerization of the chromophore, whereas two or more hydrogen bonds inhibit the isomerization and promote fluorescence instead. While rational design of new RSFPs has so far focused on structure alone, our results suggest that structural heterogeneity must be considered as well.