Polarization Effects Stabilize Bacteriorhodopsin’s Chromophore Binding Pocket: A Molecular Dynamics Study

Earliest photic zone niches probed by ancestral microbial rhodopsins

For billions of years, life has continuously adapted to dynamic physical conditions near the Earth’s surface. Fossils and other preserved biosignatures in the paleontological record are the most direct evidence for reconstructing the broad historical contours of this adaptive interplay. However, biosignatures dating to Earth’s earliest history are exceedingly rare. Here, we combine phylogenetic inference of primordial rhodopsin proteins with modeled spectral features of the Precambrian Earth environment to reconstruct the paleobiological history of this essential family of photoactive transmembrane proteins. Our results suggest that ancestral microbial rhodopsins likely acted as light-driven proton pumps and were spectrally tuned toward the absorption of green light, which would have enabled their hosts to occupy depths in a water column or biofilm where UV wavelengths were attenuated. Subsequent diversification of rhodopsin functions and peak absorption frequencies was enabled by the expansion of surface ecological niches induced by the accumulation of atmospheric oxygen. Inferred ancestors retain distinct associations between extant functions and peak absorption frequencies. Our findings suggest that novel information encoded by biomolecules can be used as “paleosensors” for conditions of ancient, inhabited niches of host organisms not represented elsewhere in the paleontological record. The coupling of functional diversification and spectral tuning of this taxonomically diverse protein family underscores the utility of rhodopsins as universal testbeds for inferring remotely detectable biosignatures on inhabited planetary bodies.

Understanding Colour Tuning Rules and Predicting Absorption Wavelengths of Microbial Rhodopsins by Data-Driven Machine-Learning Approach

The light-dependent ion-transport function of microbial rhodopsin has been widely used in optogenetics for optical control of neural activity. In order to increase the variety of rhodopsin proteins having a wide range of absorption wavelengths, the light absorption properties of various wild-type rhodopsins and their artificially mutated variants were investigated in the literature. Here, we demonstrate that a machine-learning-based (ML-based) data-driven approach is useful for understanding and predicting the light-absorption properties of microbial rhodopsin proteins. We constructed a database of 796 proteins consisting of microbial rhodopsin wildtypes and their variants. We then proposed an ML method that produces a statistical model describing the relationship between amino-acid sequences and absorption wavelengths and demonstrated that the fitted statistical model is useful for understanding colour tuning rules and predicting absorption wavelengths. By applying the ML method to the database, two residues that were not considered in previous studies are newly identified to be important to colour shift.

X-ray structure analysis of bacteriorhodopsin at 1.3 Å resolution

Bacteriorhodopsin (bR) of Halobacterium salinarum is a membrane protein that acts as a light-driven proton pump. bR and its homologues have recently been utilized in optogenetics and other applications. Although the structures of those have been reported so far, the resolutions are not sufficient for elucidation of the intrinsic structural features critical to the color tuning and ion pumping properties. Here we report the accurate crystallographic analysis of bR in the ground state. The influence of X-rays was suppressed by collecting the data under a low irradiation dose at 15 K. Consequently, individual atoms could be separately observed in the electron density map at better than 1.3 Å resolution. Residues from Thr5 to Ala233 were continuously constructed in the model. The twist of the retinal polyene was determined to be different from those in the previous models. Two conformations were observed for the proton release region. We discuss the meaning of these fine structural features.

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

Density functional tight binding: values of semi-empirical methods in an ab initio era

Semi-empirical (SE) methods are derived from Hartree-Fock (HF) or Density Functional Theory (DFT) by neglect and approximation of electronic integrals. Thereby, parameters are introduced which have to be determined from reference calculations and/or by fitting to available experimental data. This leads to computational methods that are about 2-3 orders of magnitude faster than the standard HF/DFT methods using medium sized basis sets while being about 3 orders of magnitude slower than empirical force field methods (Molecular Mechanics: MM). Therefore, SE methods are most appropriate for a specific range of applications. These include the study of systems that contain a large number of atoms and therefore being too large for ab initio or DFT methods and also problems where dynamic or entropic effects are particularly important. In the latter case, the errors made by considering a very limited number of molecular structures or neglecting entropic contributions can be much larger than the accuracy lost due to the use of SE methods. Another area where SE methods are attractive concerns the analysis of systems for which reliable MM models are not readily available. Therefore, even in an era when rapid progress is being made in ab initio methods, there is considerable interest in further developing SE methods. We illustrate this point by focusing on the discussion of recent development and application of the Density Functional Tight Binding method.

Coupling density functional theory to polarizable force fields for efficient and accurate Hamiltonian molecular dynamics simulations

Hybrid molecular dynamics (MD) simulations, in which the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules, pose a challenge. A corresponding computational approach should guarantee energy conservation, exclude artificial distortions of the electron density at the interface between the DFT and PMM fragments, and should treat the long-range electrostatic interactions within the hybrid simulation system in a linearly scaling fashion. Here we describe a corresponding Hamiltonian DFT/(P)MM implementation, which accounts for inducible atomic dipoles of a PMM environment in a joint DFT/PMM self-consistency iteration. The long-range parts of the electrostatics are treated by hierarchically nested fast multipole expansions up to a maximum distance dictated by the minimum image convention of toroidal boundary conditions and, beyond that distance, by a reaction field approach such that the computation scales linearly with the number of PMM atoms. Short-range over-polarization artifacts are excluded by using Gaussian inducible dipoles throughout the system and Gaussian partial charges in the PMM region close to the DFT fragment. The Hamiltonian character, the stability, and efficiency of the implementation are investigated by hybrid DFT/PMM-MD simulations treating one molecule of the water dimer and of bulk water by DFT and the respective remainder by PMM.

Polarizable water models from mixed computational and empirical optimization

Here we suggest a mixed computational and empirical approach serving to optimize the parameters of complex and polarizable molecular mechanics (PMM) models for complicated liquids. The computational part of the parameter optimization relies on hybrid calculations combining density functional theory (DFT) for a solute molecule with a PMM treatment of its solvent environment at well-defined thermodynamic conditions. As an application we have developed PMM models for water featuring ν = 3, 4, and 5 points of force action, a Gaussian inducible dipole and a Buckingham potential at the oxygen, the experimental liquid phase geometry, the experimental gas phase polarizability α(exp)(g) = 1.47 ų, and, for ν = 4 and 5, the gas phase value μ(exp)(g) = 1.855 D for the static dipole moment. The widths of the Gaussian dipoles and, for ν = 4 and 5, also the electrostatic geometries of these so-called TLνP models are derived from self-consistent DFT/PMM calculations, and the parameters of the Buckingham potentials (and the static TL3P dipole moment) are estimated from molecular dynamics (MD) simulations. The high quality of the resulting models is demonstrated for the observables targeted during optimization (potential energy per molecule, pressure, radial distribution functions) and a series of predicted properties (quadrupole moments, density at constant pressure, dielectric constant, diffusivity, viscosity, compressibility, heat capacity) at certain standard conditions. Remaining deficiencies and possible ways for their removal are discussed.

QM and QM/MM simulations of proteins

Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.

Building a model of the blue cone pigment based on the wild type rhodopsin structure with QM/MM methods

Understanding the mechanism of color tuning of the retinal chromophore by its host protein became one of the key issues in the research on rhodopsins. While early mutation studies addressed its genetic origin, recent studies advanced to investigate its structural origin, based on X-ray crystallographic structures. For the human cone pigments, no crystal structures have been produced, and homology models were employed to elucidate the origin of its blue-shifted absorption. In this theoretical study, we take a different route to establish a structural model for human blue. Starting from the well-resolved structure of bovine rhodopsin, we derive multiple mutant models by stepwise mutation and equilibration using molecular dynamics simulations in a hybrid quantum mechanics/molecular mechanics framework. Our 30fold mutant reproduces the experimental UV-vis absorption shift of 0.45 eV and provides new insights about both structural and genetic factors that affect the excitation energy. Electrostatic effects of individual amino acids and collaborative structural effects are analyzed using semiempirical (OM2/MRCI) and ab initio (SORCI) multireference approaches.

DFT/MM description of flavin IR spectra in BLUF domains

A class of photoreceptors occurring in various organisms consists of domains that are blue light sensing using flavin (BLUF). The vibrational spectra of the flavin chromophore are spectroscopically well characterized for the dark-adapted resting states and for the light-adapted signaling states of BLUF domains in solution. Here we present a theoretical analysis of such spectra by applying density functional theory (DFT) to the flavin embedded in molecular mechanics (MM) models of its protein and solvent environment. By DFT/MM we calculate flavin spectra for seven different X-ray and NMR structures of BLUF domains occurring in the transcriptional antirepressor AppA and in the blue light receptor B (BlrB) of the purple bacterium Rb. Sphaeroides as well as in the phototaxis photoreceptor Slr1694 of the cyanobacterium Synechocystis. By considering the dynamical stabilities of associated all-atom simulation models and by comparing calculated with observed vibrational spectra, we show that two of the considered structures (both AppA) are obviously erroneous and that specific features of two further crystal structures (BlrB and Slr1694) cannot represent the states of the respective BLUF domains in solution. Thereby, the conformational transitions elicited by solvation are identified. In this context we demonstrate how hydrogen bonds of varying strengths can tune in BLUF domains the C═O stretching frequencies of the flavin chromophore. Furthermore we show that the DFT/MM spectra of the flavin calculated for two different AppA BLUF conformations, which are called Trp(in) and Met(in), fit very well to the spectroscopic data observed for the dark and light states, respectively, if (i) polarized MM force fields, which are calculated by an iterative DFT/MM procedure, are employed for the flavin binding pockets and (ii) the calculated frequencies are properly scaled. Although the associated analysis indicates that the Trp(in) conformation belongs to the dark state, no clear light vs dark distinction emerges for the Met(in) conformation. In this connection, a number of methodological issues relevant for such complex computations are thoroughly discussed showing, in particular, why our current descriptions could not decide the light vs dark question for Met(in).

Structure of the chromophore binding pocket in the Pr state of plant phytochrome phyA

A homology structural model was generated for plant phytochrome phyA utilizing the crystal structure of the sensory module of cyanobacterial phytochrome Cph1 (Cph1Δ2). As chromophores, either the native phytochromobilin cofactor (PΦB) or phycocyanobilin (PCB), the natural cofactor in Cph1, was incorporated. These homology models were further optimized by molecular dynamics (MD) simulations revealing a satisfying overall agreement with the crystal structure of Cph1Δ2. Notable differences in the PΦB adduct of phyA result from a restructuring of the small helical segment α(7) that leads to displacements of a few amino acids away from the cofactor. This repositioning of residues also include aspartate 218 such that, instead of its carbonyl function as in Cph1Δ2, an additional water molecule forms hydrogen bonds with the ring B and C NH groups. To validate the phyA structural model in the chromophore binding pocket, Raman spectra of the cofactor were calculated by means of the quantum mechanics/molecular mechanics (QM/MM) hybrid methodology and compared with the experimental resonance Raman (RR) spectra. The satisfactory overall agreement between calculated and experimental spectra is taken as an indication for the good quality of the structural model. Moreover, the methine bridge stretching modes and the effects of isotopic labeling at selected positions of the chromophore are very well reproduced to allow confirming even details of the methine bridge geometry as predicted by the homology model. Specifically, it is demonstrated that the experimental RR spectra are consistent with a torsional angle of ring D with respect to ring C that is distinctly higher for phyA-PCB (45°) and phyA-PΦB (42°) than for Cph1Δ2 (30°). Raman spectra calculated from different points of the MD trajectory display variations of the mode frequencies and intensities reflecting the structural fluctuations from snapshot to snapshot. The snapshot spectrum of the lowest energy structure and the sum of all snapshot spectra afford an equally good description of the experimental data. Particularly large variations between the snapshots are noted for the N-H in-plane bending mode of the pyrrole rings B and C, which reflect alterations of the hydrogen bond interactions brought about by fluctuations of water molecules in the cofactor cavity. This overestimation of the water molecule mobility is a consequence of the deficiency of the current QM/MM methodology that, due to the lack of appropriate protein force fields, cannot adequately account for the electrostatics in the cofactor pocket.

On easy implementation of a variant of the replica exchange with solute tempering in GROMACS

To reduce the number of replicas required in the conventional replica exchange method for huge systems, recently the replica exchange with solute tempering (REST) method was proposed. Here we showed that a variant of REST realized by rescaling the force-field parameters can be performed with GROMACS 4 without changing the code. We tested the variant REST for alanine dipeptide and an N-terminal peptide from p53 confirming its performance nearly equal to the original REST.