DoOPSearch: a web-based tool for finding and analysing common conserved motifs in the promoter regions of different chordate and plant genes.

BACKGROUND: The comparative genomic analysis of a large number of orthologous promoter regions of the chordate and plant genes from the DoOP databases shows thousands of conserved motifs. Most of these motifs differ from any known transcription factor binding site (TFBS). To identify common conserved motifs, we need a specific tool to be able to search amongst them. Since conserved motifs from the DoOP databases are linked to genes, the result of such a search can give a list of genes that are potentially regulated by the same transcription factor(s). RESULTS: We have developed a new tool called DoOPSearch http://doopsearch.abc.hu for the analysis of the conserved motifs in the promoter regions of chordate or plant genes. We used the orthologous promoters of the DoOP database to extract thousands of conserved motifs from different taxonomic groups. The advantage of this approach is that different sets of conserved motifs might be found depending on how broad the taxonomic coverage of the underlying orthologous promoter sequence collection is (consider e.g. primates vs. mammals or Brassicaceae vs. Viridiplantae). The DoOPSearch tool allows the users to search these motif collections or the promoter regions of DoOP with user supplied query sequences or any of the conserved motifs from the DoOP database. To find overrepresented gene ontologies, the gene lists obtained can be analysed further using a modified version of the GeneMerge program. CONCLUSION: We present here a comparative genomics based promoter analysis tool. Our system is based on a unique collection of conserved promoter motifs characteristic of different taxonomic groups. We offer both a command line and a web-based tool for searching in these motif collections using user specified queries. These can be either short promoter sequences or consensus sequences of known transcription factor binding sites. The GeneMerge analysis of the search results allows the user to identify statistically overrepresented Gene Ontology terms that might provide a clue on the function of the motifs and genes.

Proton-driven amide bond-cleavage pathways of gas-phase peptide ions lacking mobile protons.

The mobile proton model (Dongre, A. R., Jones, J. L., Somogyi, A. and Wysocki, V. H. J. Am. Chem. Soc. 1996, 118 , 8365-8374) of peptide fragmentation states that the ionizing protons play a critical role in the gas-phase fragmentation of protonated peptides upon collision-induced dissociation (CID). The model distinguishes two classes of peptide ions, those with or without easily mobilizable protons. For the former class mild excitation leads to proton transfer reactions which populate amide nitrogen protonation sites. This enables facile amide bond cleavage and thus the formation of b and y sequence ions. In contrast, the latter class of peptide ions contains strongly basic functionalities which sequester the ionizing protons, thereby often hindering formation of sequence ions. Here we describe the proton-driven amide bond cleavages necessary to produce b and y ions from peptide ions lacking easily mobilizable protons. We show that this important class of peptide ions fragments by different means from those with easily mobilizable protons. We present three new amide bond cleavage mechanisms which involve salt-bridge, anhydride, and imine enol intermediates, respectively. All three new mechanisms are less energetically demanding than the classical oxazolone b(n)-y(m) pathway. These mechanisms offer an explanation for the formation of b and y ions from peptide ions with sequestered ionizing protons which are routinely fragmented in large-scale proteomics experiments.

Long-distance proton transfer with a break in the bacteriorhodopsin active site.

Bacteriorhodopsin is a proton-pumping membrane protein found in the plasma membrane of the archaeon Halobacterium salinarium. Light-induced isomerization of the retinal chromophore from all-trans to 13-cis leads to a sequence of five conformation-coupled proton transfer steps and the net transport of one proton from the cytoplasmic to the extracellular side of the membrane. The mechanism of the long-distance proton transfer from the primary acceptor Asp85 to the extracellular proton release group during the O --> bR is poorly understood. Experiments suggest that this long-distance transfer could involve a transient state [O] in which the proton resides on the intermediate carrier Asp212. To assess whether the transient protonation of Asp212 participates in the deprotonation of Asp85, we performed hybrid Quantum Mechanics/Molecular Mechanics proton transfer calculations using different protein structures and with different retinal geometries and active site water molecules. The structural models were assessed by computing UV-vis excitation energies and C=O vibrational frequencies. The results indicate that a transient [O] conformer with protonated Asp212 could indeed be sampled during the long-distance proton transfer to the proton release group. Our calculations suggest that, in the starting proton transfer state O, the retinal is strongly twisted and at least three water molecules are present in the active site.

Automatic detection of exonic splicing enhancers (ESEs) using SVMs.

BACKGROUND: Exonic splicing enhancers (ESEs) activate nearby splice sites and promote the inclusion (vs. exclusion) of exons in which they reside, while being a binding site for SR proteins. To study the impact of ESEs on alternative splicing it would be useful to have a possibility to detect them in exons. Identifying SR protein-binding sites in human DNA sequences by machine learning techniques is a formidable task, since the exon sequences are also constrained by their functional role in coding for proteins. RESULTS: The choice of training examples needed for machine learning approaches is difficult since there are only few exact locations of human ESEs described in the literature which could be considered as positive examples. Additionally, it is unclear which sequences are suitable as negative examples. Therefore, we developed a motif-oriented data-extraction method that extracts exon sequences around experimentally or theoretically determined ESE patterns. Positive examples are restricted by heuristics based on known properties of ESEs, e.g. location in the vicinity of a splice site, whereas negative examples are taken in the same way from the middle of long exons. We show that a suitably chosen SVM using optimized sequence kernels (e.g., combined oligo kernel) can extract meaningful properties from these training examples. Once the classifier is trained, every potential ESE sequence can be passed to the SVM for verification. Using SVMs with the combined oligo kernel yields a high accuracy of about 90 percent and well interpretable parameters. CONCLUSION: The motif-oriented data-extraction method seems to produce consistent training and test data leading to good classification rates and thus allows verification of potential ESE motifs. The best results were obtained using an SVM with the combined oligo kernel, while oligo kernels with oligomers of a certain length could be used to extract relevant features.

Sequence-scrambling fragmentation pathways of protonated peptides.

The gas-phase structures and fragmentation pathways of the N-terminal b and a fragments of YAGFL-NH(2), AGLFY-NH(2), GFLYA-NH(2), FLYAG-NH(2), and LYAGF-NH(2) were investigated using collision-induced dissociation (CID) and detailed molecular mechanics and density functional theory (DFT) calculations. Our combined experimental and theoretical approach allows probing of the scrambling and rearrangement reactions that take place in CID of b and a ions. It is shown that low-energy CID of the b(5) fragments of the above peptides produces nearly the same dissociation patterns. Furthermore, CID of protonated cyclo-(YAGFL) generates the same fragments with nearly identical ion abundances when similar experimental conditions are applied. This suggests that rapid cyclization of the primarily linear b(5) ions takes place and that the CID spectrum is indeed determined by the fragmentation behavior of the cyclic isomer. This can open up at various amide bonds, and its fragmentation behavior can be understood only by assuming a multitude of fragmenting linear structures. Our computational results fully support this cyclization-reopening mechanism by showing that protonated cyclo-(YAGFL) is energetically favored over the linear b(5) isomers. Furthermore, the cyclization-reopening transition structures are energetically less demanding than those of conventional bond-breaking reactions, allowing fast interconversion among the cyclic and linear isomers. This chemistry can lead in principle to complete loss of sequence information upon CID, as documented for the b(5) ion of FLYAG-NH(2). CID of the a(5) ions of the above peptides produces fragment ion distributions that can be explained by assuming b-type scrambling of their parent population and a --> a*-type rearrangement pathways ( Vachet , R. W. , Bishop , B. M. , Erickson , B. W. , and Glish , G. L. J. Am. Chem. Soc. 1997, 119, 5481 ). While a ions easily undergo cyclization, the resulting macrocycle predominantly reopens to regenerate the original linear structure. Computational data indicate that the a --> a*-type rearrangement pathways of the linear a isomers involve post-cleavage proton-bound dimer intermediates in which the fragments reassociate and the originally C-terminal fragment is transferred to the N-terminus.

Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study.

A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.

Key role of active-site water molecules in bacteriorhodopsin proton-transfer reactions.

The functional mechanism of the light-driven proton pump protein bacteriorhodopsin depends on the location of water molecules in the active site at various stages of the photocycle and on their roles in the proton-transfer steps. Here, free energy computations indicate that electrostatic interactions favor the presence of a cytoplasmic-side water molecule hydrogen bonding to the retinal Schiff base in the state preceding proton transfer from the retinal Schiff base to Asp85. However, the nonequilibrium nature of the pumping process means that the probability of occupancy of a water molecule in a given site depends both on the free energies of insertion of the water molecule in this and other sites during the preceding photocycle steps and on the kinetic accessibility of these sites on the time scale of the reaction steps. The presence of the cytoplasmic-side water molecule has a dramatic effect on the mechanism of proton transfer: the proton is channeled on the Thr89 side of the retinal, whereas the transfer on the Asp212 side is hindered. Reaction-path simulations and molecular dynamics simulations indicate that the presence of the cytoplasmic-side water molecule permits a low-energy bacteriorhodopsin conformer in which the water molecule bridges the twisted retinal Schiff base and the proton acceptor Asp85. From this low-energy conformer, proton transfer occurs via a concerted mechanism in which the water molecule participates as an intermediate proton carrier.

Vibrational spectroscopy and conformational structure of protonated polyalanine peptides isolated in the gas phase.

The conformational structures of protonated polyalanine peptides, Ala(n)H(+), have been investigated in the gas phase for n = 3, 4, 5, and 7 using a combination of resonant-infrared multiphoton dissociation (R-IRMPD) spectroscopy in the NH and OH stretch regions and quantum chemical calculations. Agreement between theoretical IR and experimental R-IRMPD spectral features has enabled the assignment of specific hydrogen-bonded conformational motifs in the short protonated peptides and revealed their conformational evolution under elevated-temperature conditions, as a function of increasing chain length. The shortest peptide, Ala(3)H(+), adopts a mixture of extended and cyclic chain conformations, protonated, respectively, at a backbone carbonyl or the N-terminus. The longer peptides adopt folded, cyclic, and globular charge-solvated conformations protonated at the N-terminus, consistent with previous ion-mobility studies. The longest peptide, Ala(7)H(+), adopts a globular conformation with the N-terminus completely charge-solvated, demonstrating the emergence of "physiologically relevant" intramolecular interactions in the peptide backbone. The computed conformational relative free energies highlight the importance of entropic contributions in these peptides.

SERpredict: detection of tissue- or tumor-specific isoforms generated through exonization of transposable elements.

BACKGROUND: Transposed elements (TEs) are known to affect transcriptomes, because either new exons are generated from intronic transposed elements (this is called exonization), or the element inserts into the exon, leading to a new transcript. Several examples in the literature show that isoforms generated by an exonization are specific to a certain tissue (for example the heart muscle) or inflict a disease. Thus, exonizations can have negative effects for the transcriptome of an organism. RESULTS: As we aimed at detecting other tissue- or tumor-specific isoforms in human and mouse genomes which were generated through exonization of a transposed element, we designed the automated analysis pipeline SERpredict (SER = Specific Exonized Retroelement) making use of Bayesian Statistics. With this pipeline, we found several genes in which a transposed element formed a tissue- or tumor-specific isoform. CONCLUSION: Our results show that SERpredict produces relevant results, demonstrating the importance of transposed elements in shaping both the human and the mouse transcriptomes. The effect of transposed elements on the human transcriptome is several times higher than the effect on the mouse transcriptome, due to the contribution of the primate-specific Alu elements.

Epigenetic reactivation of tumor suppressor genes by a novel small-molecule inhibitor of human DNA methyltransferases.

DNA methylation regulates gene expression in normal and malignant cells. The possibility to reactivate epigenetically silenced genes has generated considerable interest in the development of DNA methyltransferase inhibitors. Here, we provide a detailed characterization of RG108, a novel small molecule that effectively blocked DNA methyltransferases in vitro and did not cause covalent enzyme trapping in human cell lines. Incubation of cells with low micromolar concentrations of the compound resulted in significant demethylation of genomic DNA without any detectable toxicity. Intriguingly, RG108 caused demethylation and reactivation of tumor suppressor genes, but it did not affect the methylation of centromeric satellite sequences. These results establish RG108 as a DNA methyltransferase inhibitor with fundamentally novel characteristics that will be particularly useful for the experimental modulation of epigenetic gene regulation.

GOPET: a tool for automated predictions of Gene Ontology terms.

BACKGROUND: Vast progress in sequencing projects has called for annotation on a large scale. A Number of methods have been developed to address this challenging task. These methods, however, either apply to specific subsets, or their predictions are not formalised, or they do not provide precise confidence values for their predictions. DESCRIPTION: We recently established a learning system for automated annotation, trained with a broad variety of different organisms to predict the standardised annotation terms from Gene Ontology (GO). Now, this method has been made available to the public via our web-service GOPET (Gene Ontology term Prediction and Evaluation Tool). It supplies annotation for sequences of any organism. For each predicted term an appropriate confidence value is provided. The basic method had been developed for predicting molecular function GO-terms. It is now expanded to predict biological process terms. This web service is available via http://genius.embnet.dkfz-heidelberg.de/menu/biounit/open-husar CONCLUSION: Our web service gives experimental researchers as well as the bioinformatics community a valuable sequence annotation device. Additionally, GOPET also provides less significant annotation data which may serve as an extended discovery platform for the user.

CAFTAN: a tool for fast mapping, and quality assessment of cDNAs.

BACKGROUND: The German cDNA Consortium has been cloning full length cDNAs and continued with their exploitation in protein localization experiments and cellular assays. However, the efficient use of large cDNA resources requires the development of strategies that are capable of a speedy selection of truly useful cDNAs from biological and experimental noise. To this end we have developed a new high-throughput analysis tool, CAFTAN, which simplifies these efforts and thus fills the gap between large-scale cDNA collections and their systematic annotation and application in functional genomics. RESULTS: CAFTAN is built around the mapping of cDNAs to the genome assembly, and the subsequent analysis of their genomic context. It uses sequence features like the presence and type of PolyA signals, inner and flanking repeats, the GC-content, splice site types, etc. All these features are evaluated in individual tests and classify cDNAs according to their sequence quality and likelihood to have been generated from fully processed mRNAs. Additionally, CAFTAN compares the coordinates of mapped cDNAs with the genomic coordinates of reference sets from public available resources (e.g., VEGA, ENSEMBL). This provides detailed information about overlapping exons and the structural classification of cDNAs with respect to the reference set of splice variants. The evaluation of CAFTAN showed that is able to correctly classify more than 85% of 5950 selected "known protein-coding" VEGA cDNAs as high quality multi- or single-exon. It identified as good 80.6 % of the single exon cDNAs and 85 % of the multiple exon cDNAs. The program is written in Perl and in a modular way, allowing the adoption of this strategy to other tasks like EST-annotation, or to extend it by adding new classification rules and new organism databases as they become available. We think that it is a very useful program for the annotation and research of unfinished genomes. CONCLUSION: CAFTAN is a high-throughput sequence analysis tool, which performs a fast and reliable quality prediction of cDNAs. Several thousands of cDNAs can be analyzed in a short time, giving the curator/scientist a first quick overview about the quality and the already existing annotation of a set of cDNAs. It supports the rejection of low quality cDNAs and helps in the selection of likely novel splice variants, and/or completely novel transcripts for new experiments.

Genomic analysis of Xenopus organizer function.

BACKGROUND: Studies of the Xenopus organizer have laid the foundation for our understanding of the conserved signaling pathways that pattern vertebrate embryos during gastrulation. The two primary activities of the organizer, BMP and Wnt inhibition, can regulate a spectrum of genes that pattern essentially all aspects of the embryo during gastrulation. As our knowledge of organizer signaling grows, it is imperative that we begin knitting together our gene-level knowledge into genome-level signaling models. The goal of this paper was to identify complete lists of genes regulated by different aspects of organizer signaling, thereby providing a deeper understanding of the genomic mechanisms that underlie these complex and fundamental signaling events. RESULTS: To this end, we ectopically overexpress Noggin and Dkk-1, inhibitors of the BMP and Wnt pathways, respectively, within ventral tissues. After isolating embryonic ventral halves at early and late gastrulation, we analyze the transcriptional response to these molecules within the generated ectopic organizers using oligonucleotide microarrays. An efficient statistical analysis scheme, combined with a new Gene Ontology biological process annotation of the Xenopus genome, allows reliable and faithful clustering of molecules based upon their roles during gastrulation. From this data, we identify new organizer-related expression patterns for 19 genes. Moreover, our data sub-divides organizer genes into separate head and trunk organizing groups, which each show distinct responses to Noggin and Dkk-1 activity during gastrulation. CONCLUSION: Our data provides a genomic view of the cohorts of genes that respond to Noggin and Dkk-1 activity, allowing us to separate the role of each in organizer function. These patterns demonstrate a model where BMP inhibition plays a largely inductive role during early developmental stages, thereby initiating the suites of genes needed to pattern dorsal tissues. Meanwhile, Wnt inhibition acts later during gastrulation, and is essential for maintenance of organizer gene expression throughout gastrulation, a role which may depend on its ability to block the expression of a host of ventral, posterior, and lateral fate-specifying factors.

Discovery of two novel, small-molecule inhibitors of DNA methylation.

DNA methyltransferases are promising targets for cancer therapy. In many cancer cells promoters of tumor suppressor genes are hypermethylated, which results in gene inactivation. It has been shown that DNA methyltransferase inhibitors can suppress tumor growth and have significant therapeutic value. However, the established inhibitors are limited in their application due to their substantial cytotoxicity. To discover novel compounds for the inhibition of human DNA methyltransferases, we have screened a set of small molecules available from the NCI database. Using a 3-dimensional model of the human DNA methyltransferase 1 and a modified docking and scoring procedure, we have identified a small list of molecules with high affinities for the active site of the enzyme. The two highest scoring structures were found to inhibit DNA methyltransferase activity in vitro and in vivo. The newly discovered inhibitors validate our screening procedure and also provide a useful basis for further rational drug development.

Revising the proton affinity scale of the naturally occurring alpha-amino acids.

The proton affinities (PA) of the 20 naturally occurring alpha-amino acids (AA) have been determined computationally by means of density functional theory (DFT) and high-level G2(MP2) calculations. These theoretical PAs, together with data that have appeared since 1997 in the literature, are used to validate the most reasonable currently available PA scale for AAs (Harrison, A. G. Mass Spectrom. Rev. 1997, 16, 201-217.). Significant scatter is observed for the PAs of Ser, Asp, Phe, Asn, Met, Pro, Gln, Glu, Trp, His, Lys, and Arg, many of which have a basic side-chain functionality. Critical review of the available data leads to new consensus PAs for Asn, Gln, Met, and Arg of 222.4, 230.5, 223.7, and 250.2 kcal/mol, respectively.

Isotope labeling and theoretical study of the formation of a3* ions from protonated tetraglycine.

Extensive 13C, 15N, and 2H labeling of tetraglycine was used to investigate the b3+ --> a3* reaction during low-energy collision-induced dissociation (CID) in a quadrupole ion-trap mass spectrometer. The patterns observed with respect to the retention or elimination of the isotope labels demonstrate that the reaction pathway involves elimination of CO and NH3. The ammonia molecule includes 2 H atoms from amide or amino positions, and one from an alpha-carbon position. The loss of NH3 does not involve elimination of the N-terminal amino group but, instead, the N atom of the presumed oxazolone ring in the b3+ ion. The CO molecule eliminated is the carbonyl group of the same oxazolone ring, and the alpha-carbon H atom is transferred from the amino acid adjacent to the oxazolone ring. Quantum chemical calculations indicate a multistep reaction cascade involving CO loss on the b3 --> a3 pathway and loss of NH=CH2 from the a3 ion to form b2. In the postreaction complex of b2 and NH=CH2, the latter can be attacked by the N-terminal amino group of the former. The product of this attack, an isomerized a3 ion, can eliminate NH3 from its N-terminus to form a3*. Calculations suggest that the ammonia and a3* species can form various ion-molecule complexes, and NH3 can initiate relay-type mobilization of the oxazolone H atoms from alpha-carbon positions to form a new oxazolone isomer. This multiple-step reaction scheme clearly explains the isotope labeling results, including unexpected scrambling of H atoms from alpha-carbon positions.

ESTAnnotator: A tool for high throughput EST annotation.

In high throughput sequence analysis, it is often necessary to combine the results of contemporary bioinformatics tools, because no individual tool alone computes all the requested information. ESTAnnotator is a tool for the high throughput annotation of expressed sequence tags (ESTs) by automatically running a collection of bioinformatics applications. In the first step, a quality check is performed and repeats, vector parts and low quality sequences are masked. Then successive steps of database searching and EST clustering are performed. Already known transcripts present within mRNA and genomic DNA reference databases are identified. Subsequently, tools for the clustering of anonymous ESTs, and for further database searches at the protein level, are applied. Finally, the outputs of each individual tool are gathered and the relevant results presented in a descriptive summary. ESTAnnotator was already successfully applied for the systematic identification and characterisation of novel human genes involved in cartilage/bone formation, growth, differentiation and homeostasis. ESTAnnotator is available at http://genome.dkfz-heidelberg.de, contact: genome@dkfz.de.

High-throughput protein analysis integrating bioinformatics and experimental assays.

The wealth of transcript information that has been made publicly available in recent years requires the development of high-throughput functional genomics and proteomics approaches for its analysis. Such approaches need suitable data integration procedures and a high level of automation in order to gain maximum benefit from the results generated. We have designed an automatic pipeline to analyse annotated open reading frames (ORFs) stemming from full-length cDNAs produced mainly by the German cDNA Consortium. The ORFs are cloned into expression vectors for use in large-scale assays such as the determination of subcellular protein localization or kinase reaction specificity. Additionally, all identified ORFs undergo exhaustive bioinformatic analysis such as similarity searches, protein domain architecture determination and prediction of physicochemical characteristics and secondary structure, using a wide variety of bioinformatic methods in combination with the most up-to-date public databases (e.g. PRINTS, BLOCKS, INTERPRO, PROSITE SWISSPROT). Data from experimental results and from the bioinformatic analysis are integrated and stored in a relational database (MS SQL-Server), which makes it possible for researchers to find answers to biological questions easily, thereby speeding up the selection of targets for further analysis. The designed pipeline constitutes a new automatic approach to obtaining and administrating relevant biological data from high-throughput investigations of cDNAs in order to systematically identify and characterize novel genes, as well as to comprehensively describe the function of the encoded proteins.

Mechanism of primary proton transfer in bacteriorhodopsin.

Recent structures of putative intermediates in the bacteriorhodopsin photocycle have provided valuable snapshots of the mechanism by which protons are pumped across the membrane. However, key steps remain highly controversial, particularly the proton transfer occurring immediately after retinal trans-->cis photoisomerization. The gradual release of stored energy is inherently nonequilibrium: which photocycle intermediates are populated depends not only on their energy but also on their interconversion rates. To understand why the photocycle follows a productive (i.e., pumping), rather than some unproductive, relaxation pathway, it is necessary to know the relative energy barriers of individual steps. To discriminate between the many proposed scenarios of this process, we computed all its possible minimum-energy paths. This reveals that not one, but three very different pathways have energy barriers consistent with experiment. This result reconciles the conflicting views held on the mechanism and suggests a strategy by which the protein renders this essential step resilient.

Tuning of retinal twisting in bacteriorhodopsin controls the directionality of the early photocycle steps.

Productive proton pumping by bacteriorhodopsin requires that, after the all-trans to 13-cis photoisomerization of the retinal chromophore, the photocycle proceeds with proton transfer and not with thermal back-isomerization. The question of how the protein controls these events in the active site is addressed here using quantum mechanical/molecular mechanical reaction-path calculations. The results indicate that, while retinal twisting significantly contributes to lowering the barrier for the thermal cis-trans back-isomerization, the rate-limiting barrier for this isomerization is still 5-6 kcal/mol larger than that for the first proton-transfer step. In this way, the retinal twisting is finely tuned so as to store energy to drive the subsequent photocycle while preventing wasteful back-isomerization.