K. Anton Feenstra

Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems

A systematic analysis is performed on the effectiveness of removing degrees of freedom from hydrogen atoms and/or increasing hydrogen masses to increase the efficiency of molecular dynamics simulations of hydrogen‐rich systems such as proteins in water. In proteins, high‐frequency bond‐angle vibrations involving hydrogen atoms limit the time step to 3 fs, which is already a factor of 1.5 beyond the commonly used time step of 2 fs. Removing these degrees of freedom from the system by constructing hydrogen atoms as dummy atoms, allows the time step to be increased to 7 fs, a factor of 3.5 compared with 2 fs. Additionally, a gain in simulation stability can be achieved by increasing the masses of hydrogen atoms with remaining degrees of freedom from 1 to 4 u. Increasing hydrogen mass without removing the high‐frequency degrees of freedom allows the time step to be increased only to 4 fs, a factor of two, compared with 2 fs. The net gain in efficiency of sampling configurational space may be up to 15% lower than expected from the increase in time step due to the increase in viscosity and decrease in diffusion constant. In principle, introducing dummy atoms and increasing hydrogen mass do not influence thermodynamical properties of the system and dynamical properties are shown to be influenced only to a moderate degree. Comparing the maximum time step attainable with these methods (7 fs) to the time step of 2 fs that is routinely used in simulation, and taking into account the increase in viscosity and decrease in diffusion constant, we can say that a net gain in simulation efficiency of a factor of 3 to 3.5 can be achieved. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 786–798, 1999

PRALINEtm: a strategy for improved multiple alignment of transmembrane proteins.

MOTIVATION Membrane-bound proteins are a special class of proteins. The regions that insert into the cell-membrane have a profoundly different hydrophobicity pattern compared with soluble proteins. Multiple alignment techniques use scoring schemes tailored for sequences of soluble proteins and are therefore in principle not optimal to align membrane-bound proteins. RESULTS Transmembrane (TM) regions in protein sequences can be reliably recognized using state-of-the-art sequence prediction techniques. Furthermore, membrane-specific scoring matrices are available. We have developed a new alignment method, called PRALINETM, which integrates these two features to enhance multiple sequence alignment. We tested our algorithm on the TM alignment benchmark set by Bahr et al. (2001), and showed that the quality of TM alignments can be significantly improved compared with the quality produced by a standard multiple alignment technique. The results clearly indicate that the incorporation of these new elements into current state-of-the-art alignment methods is crucial for optimizing the alignment of TM proteins. AVAILABILITY A webserver is available at http://www.ibi.vu.nl/programs/pralinewww.

Sequence comparison by sequence harmony identifies subtype-specific functional sites

Multiple sequence alignments are often used to reveal functionally important residues within a protein family. They can be particularly useful for the identification of key residues that determine functional differences between protein subfamilies. We present a new entropy-based method, Sequence Harmony (SH) that accurately detects subfamily-specific positions from a multiple sequence alignment. The SH algorithm implements a novel formula, able to score compositional differences between subfamilies, without imposing conservation, in a simple manner on an intuitive scale. We compare our method with the most important published methods, i.e. AMAS, TreeDet and SDP-pred, using three well-studied protein families: the receptor-binding domain (MH2) of the Smad family of transcription factors, the Ras-superfamily of small GTPases and the MIP-family of integral membrane transporters. We demonstrate that SH accurately selects known functional sites with higher coverage than the other methods for these test-cases. This shows that compositional differences between protein subfamilies provide sufficient basis for identification of functional sites. In addition, SH selects a number of sites of unknown function that could be interesting candidates for further experimental investigation.

Metabolic Regio- and Stereoselectivity of Cytochrome P450 2D6 towards 3,4-Methylenedioxy-N-alkylamphetamines: in Silico Predictions and Experimental Validation

A series of 3,4-methylenedioxy-N-alkylamphetamines (MDAAs) were automatically docked and subjected to molecular dynamics (MD) simulations in a cytochrome P450 2D6 (CYP2D6) protein model. The predicted substrate binding orientations, sites of oxidation, and relative reactivities were compared to the experimental data of wild-type and Phe120Ala mutant CYP2D6. Automated docking results were not sufficient to accurately rationalize experimental binding orientations of 3,4-methylenedioxy-N-methylamphetamine (MDMA) in the two enzymes as measured with spin lattice relaxation NMR. Nevertheless, the docking results could be used as starting structures for MD simulations. Predicted binding orientations of MDMA and sites of oxidation of the MDAAs derived from MD simulations matched well with the experimental data. It appeared the experimental results were best described in MD simulations considering the nitrogen atoms of the MDAAs in neutral form. Differences in regioselectivity and stereoselectivity in the oxidative metabolism of the MDAAs by the Phe120Ala mutant CYP2D6 were correctly predicted, and the effects of the Phe120Ala mutation could be rationalized as well.

Interaction of 14-3-3 proteins with the Estrogen Receptor Alpha F domain provides a drug target interface

Estrogen receptor alpha (ERα) is involved in numerous physiological and pathological processes, including breast cancer. Breast cancer therapy is therefore currently directed at inhibiting the transcriptional potency of ERα, either by blocking estrogen production through aromatase inhibitors or antiestrogens that compete for hormone binding. Due to resistance, new treatment modalities are needed and as ERα dimerization is essential for its activity, interference with receptor dimerization offers a new opportunity to exploit in drug design. Here we describe a unique mechanism of how ERα dimerization is negatively controlled by interaction with 14-3-3 proteins at the extreme C terminus of the receptor. Moreover, the small-molecule fusicoccin (FC) stabilizes this ERα/14-3-3 interaction. Cocrystallization of the trimeric ERα/14-3-3/FC complex provides the structural basis for this stabilization and shows the importance of phosphorylation of the penultimate Threonine (ERα-T594) for high-affinity interaction. We confirm that T594 is a distinct ERα phosphorylation site in the breast cancer cell line MCF-7 using a phospho-T594–specific antibody and by mass spectrometry. In line with its ERα/14-3-3 interaction stabilizing effect, fusicoccin reduces the estradiol-stimulated ERα dimerization, inhibits ERα/chromatin interactions and downstream gene expression, resulting in decreased cell proliferation. Herewith, a unique functional phosphosite and an alternative regulation mechanism of ERα are provided, together with a small molecule that selectively targets this ERα/14-3-3 interface.

Executing Multicellular Differentiation: Quantitative Predictive Modelling of C. elegans Vulval Development

MOTIVATION Understanding the processes involved in multi-cellular pattern formation is a central problem of developmental biology, hopefully leading to many new insights, e.g. in the treatment of various diseases. Defining suitable computational techniques for development modelling, able to perform in silico simulation experiments, is an open and challenging problem. RESULTS Previously, we proposed a coarse-grained, quantitative approach based on the basic Petri net formalism, to mimic the behaviour of the biological processes during multicellular differentiation. Here, we apply our modelling approach to the well-studied process of Caenorhabditis elegans vulval development. We show that our model correctly reproduces a large set of in vivo experiments with statistical accuracy. It also generates gene expression time series in accordance with recent biological evidence. Finally, we modelled the role of microRNA mir-61 during vulval development and predict its contribution in stabilizing cell pattern formation.

Hard-wired heterogeneity in blood stem cells revealed using a dynamic regulatory network model

Motivation: Combinatorial interactions of transcription factors with cis-regulatory elements control the dynamic progression through successive cellular states and thus underpin all metazoan development. The construction of network models of cis-regulatory elements, therefore, has the potential to generate fundamental insights into cellular fate and differentiation. Haematopoiesis has long served as a model system to study mammalian differentiation, yet modelling based on experimentally informed cis-regulatory interactions has so far been restricted to pairs of interacting factors. Here, we have generated a Boolean network model based on detailed cis-regulatory functional data connecting 11 haematopoietic stem/progenitor cell (HSPC) regulator genes. Results: Despite its apparent simplicity, the model exhibits surprisingly complex behaviour that we charted using strongly connected components and shortest-path analysis in its Boolean state space. This analysis of our model predicts that HSPCs display heterogeneous expression patterns and possess many intermediate states that can act as ‘stepping stones’ for the HSPC to achieve a final differentiated state. Importantly, an external perturbation or ‘trigger’ is required to exit the stem cell state, with distinct triggers characterizing maturation into the various different lineages. By focusing on intermediate states occurring during erythrocyte differentiation, from our model we predicted a novel negative regulation of Fli1 by Gata1, which we confirmed experimentally thus validating our model. In conclusion, we demonstrate that an advanced mammalian regulatory network model based on experimentally validated cis-regulatory interactions has allowed us to make novel, experimentally testable hypotheses about transcriptional mechanisms that control differentiation of mammalian stem cells. Contact: j.heringa@vu.nl or ioannis.xenarios@isb-sib.ch or bg200@cam.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Multi-Harmony: detecting functional specificity from sequence alignment

Many protein families contain sub-families with functional specialization, such as binding different ligands or being involved in different protein–protein interactions. A small number of amino acids generally determine functional specificity. The identification of these residues can aid the understanding of protein function and help finding targets for experimental analysis. Here, we present multi-Harmony, an interactive web sever for detecting sub-type-specific sites in proteins starting from a multiple sequence alignment. Combining our Sequence Harmony (SH) and multi-Relief (mR) methods in one web server allows simultaneous analysis and comparison of specificity residues; furthermore, both methods have been significantly improved and extended. SH has been extended to cope with more than two sub-groups. mR has been changed from a sampling implementation to a deterministic one, making it more consistent and user friendly. For both methods Z-scores are reported. The multi-Harmony web server produces a dynamic output page, which includes interactive connections to the Jalview and Jmol applets, thereby allowing interactive analysis of the results. Multi-Harmony is available at http://www.ibi.vu.nl/ programs/shmrwww.

In vitro evolution of styrene monooxygenase from Pseudomonas putida CA-3 for improved epoxide synthesis

The styAB genes from Pseudomonas putida CA-3, which encode styrene monooxygenase, were subjected to three rounds of in vitro evolution using error-prone polymerase chain reaction with a view to improving the rate of styrene oxide and indene oxide formation. Improvements in styrene monooxygenase activity were monitored using an indole to indigo conversion assay. Each round of random mutagenesis generated variants improved in indigo formation with third round variants improved nine- to 12-fold over the wild type enzyme. Each round of in vitro evolution resulted in two to three amino acid substitutions in styrene monooxygenase. While the majority of mutations occurred in styA (oxygenase), mutations were also observed in styB (reductase). A mutation resulting in the substitution of valine with isoleucine at amino acid residue 303 occurred near the styrene and flavin adenine dinucleotide binding site of styrene monooxygenase. One mutation caused a shift in the reading frame in styA and resulted in a StyA variant that is 19 amino acids longer than the wild-type protein. Whole cells expressing the best styrene monooxygenase variants (round 3) exhibited eight- and 12-fold improvements in styrene and indene oxidation rates compared to the wild-type enzyme. In all cases, a single enantiomer, (S)-styrene oxide, was formed from styrene while (1S,2R)-indene oxide was the predominant enantiomer (e.e. 97%) formed from indene. The average yield of styrene oxide and indene oxide from their respective alkene substrates was 65% and 90%, respectively.

Sequence harmony: detecting functional specificity from alignments

Multiple sequence alignments are often used for the identification of key specificity-determining residues within protein families. We present a web server implementation of the Sequence Harmony (SH) method previously introduced. SH accurately detects subfamily specific positions from a multiple alignment by scoring compositional differences between subfamilies, without imposing conservation. The SH web server allows a quick selection of subtype specific sites from a multiple alignment given a subfamily grouping. In addition, it allows the predicted sites to be directly mapped onto a protein structure and displayed. We demonstrate the use of the SH server using the family of plant mitochondrial alternative oxidases (AOX). In addition, we illustrate the usefulness of combining sequence and structural information by showing that the predicted sites are clustered into a few distinct regions in an AOX homology model. The SH web server can be accessed at www.ibi.vu.nl/programs/seqharmwww.