Peter J. Winn

Comparison of the dynamics of substrate access channels in three cytochrome P450s reveals different opening mechanisms and a novel functional role for a buried arginine.

Understanding the mechanism and specificity of substrate binding in the cytochrome P450 (P450) superfamily is an important step toward explaining its key role in drug metabolism, toxicity, xenobiotic degradation, and several biosynthetic pathways. Here we investigate the ligand exit pathways and mechanisms of P450cam (CYP101), P450BM-3 (CYP102), and P450eryF (CYP107A1) by using random expulsion molecular dynamics and classical molecular dynamics simulations. Although several different pathways are found for each protein, one pathway is common to all three. The mechanism of ligand exit along this pathway is, however, quite different in the three different proteins. For P450cam, small backbone conformational changes, in combination with aromatic side chain rotation, allow for the passage of the rather rigid, compact, and hydrophobic substrate, camphor. In P450BM-3, larger transient backbone changes are observed on ligand exit. R47, situated at the entrance to the channel, appears important in guiding negatively charged fatty acid substrates in and out of the active site. In P450eryF, an isolated buried arginine, R185, stabilized by four hydrogen bonds to backbone carbonyl oxygen atoms, is located in the exit channel and is identified as having a particularly unusual functionality, dynamically gating channel opening. The results for these three P450s suggest that the channel opening mechanisms are adjusted to the physico-chemical properties of the substrate and can kinetically modulate protein-substrate specificity.

Do mammalian cytochrome P450s show multiple ligand access pathways and ligand channelling

Understanding substrate binding and product release in cytochrome P450 (CYP) enzymes is important for explaining their key role in drug metabolism, toxicity, xenobiotic degradation and biosynthesis. Here, molecular simulations of substrate and product exit from the buried active site of a mammalian P450, the microsomal CYP2C5, identified a dominant exit channel, termed pathway (pw) 2c. Previous simulations with soluble bacterial P450s showed a different dominant egress channel, pw2a. Combining these, we propose two mechanisms in CYP2C5: (i) a one‐way route by which lipophilic substrates access the enzyme from the membrane by pw2a and hydroxylated products egress along pw2c; and (ii) a two‐way route for access and egress, along pw2c, for soluble compounds. The proposed differences in substrate access and product egress routes between membrane‐bound mammalian P450s and soluble bacterial P450s highlight the adaptability of the P450 fold to the requirements of differing cellular locations and substrate specificity profiles.

ProMetCS: An Atomistic Force Field for Modeling Protein-Metal Surface Interactions in a Continuum Aqueous Solvent.

In order to study protein-inorganic surface association processes, we have developed a physics-based energy model, the ProMetCS model, which describes protein-surface interactions at the atomistic level while treating the solvent as a continuum. Here, we present an approach to modeling the interaction of a protein with an atomically flat Au(111) surface in an aqueous solvent. Protein-gold interactions are modeled as the sum of van der Waals, weak chemisorption, and electrostatic interactions, as well as the change in free energy due to partial desolvation of the protein and the metal surface upon association. This desolvation energy includes the effects of water-protein, water-surface, and water-water interactions and has been parametrized using molecular dynamics (MD) simulations of water molecules and a test atom at a gold-water interface. The proposed procedure for computing the energy terms is mostly grid-based and is therefore efficient for application to long-time simulations of protein binding processes. The approach was tested for capped amino acid residues whose potentials of mean force for binding to a gold surface were computed and compared with those obtained previously in MD simulations with water treated explicitly. Calculations show good quantitative agreement with the results from MD simulations for all but one amino acid (Trp), as well as correspondence with available experimental data on the adhesion properties of amino acids.

A conserved motif flags acyl carrier proteins for β-branching in polyketide synthesis.

Type I PKSs often utilise programmed β-branching, via enzymes of an “HMG-CoA synthase (HCS) cassette”, to incorporate various side chains at the second carbon from the terminal carboxylic acid of growing polyketide backbones. We identified a strong sequence motif in Acyl Carrier Proteins (ACPs) where β-branching is known. Substituting ACPs confirmed a correlation of ACP type with β-branching specificity. While these ACPs often occur in tandem, NMR analysis of tandem β-branching ACPs indicated no ACP-ACP synergistic effects and revealed that the conserved sequence motif forms an internal core rather than an exposed patch. Modelling and mutagenesis identified ACP Helix III as a probable anchor point of the ACP-HCS complex whose position is determined by the core. Mutating the core affects ACP functionality while ACP-HCS interface substitutions modulate system specificity. Our method for predicting β-carbon branching expands the potential for engineering novel polyketides and lays a basis for determining specificity rules.

Protein–protein docking by simulating the process of association subject to biochemical constraints

We present a computational procedure for modeling protein–protein association and predicting the structures of protein–protein complexes. The initial sampling stage is based on an efficient Brownian dynamics algorithm that mimics the physical process of diffusional association. Relevant biochemical data can be directly incorporated as distance constraints at this stage. The docked configurations are then grouped with a hierarchical clustering algorithm into ensembles that represent potential protein–protein encounter complexes. Flexible refinement of selected representative structures is done by molecular dynamics simulation. The protein–protein docking procedure was thoroughly tested on 10 structurally and functionally diverse protein–protein complexes. Starting from X‐ray crystal structures of the unbound proteins, in 9 out of 10 cases it yields structures of protein–protein complexes close to those determined experimentally with the percentage of correct contacts >30% and interface backbone RMSD <4 Å. Detailed examination of all the docking cases gives insights into important determinants of the performance of the computational approach in modeling protein–protein association and predicting of protein–protein complex structures. Proteins 2008.

Multiple, ligand-dependent routes from the active site of cytochrome P450 2C9.

The active site of liver-specific, drug-metabolizing cytochrome P450 (CYP) monooxygenases is deeply buried in the protein and is connected to the protein surface through multiple tunnels, many of which were found open in different CYP crystal structures. It has been shown that different tunnels could serve as ligand passage routes in different CYPs. However, it is not understood whether one CYP uses multiple routes for substrate access and product release and whether these routes depend on ligand properties. From 300 ns of molecular dynamics simulations of CYP2C9, the second most abundant CYP in the human liver we found four main ligand exit routes, the occurrence of each depending on the ligand type and the conformation of the F-G loop, which is likely to be affected by the CYP-membrane interaction. A non-helical F-G loop favored exit towards the putative membrane-embedded region. Important protein features that direct ligand exit include aromatic residues that divide the active site and whose motions control access to two pathways. The ligands interacted with positively charged residues on the protein surface through hydrogen bonds that appear to select for acidic substrates. The observation of multiple, ligand-dependent routes in a CYP aids understanding of how CYP mutations affect drug metabolism and provides new possibilities for CYP inhibition.

Characterization of mutations in the PAS domain of the EvgS sensor kinase selected by laboratory evolution for acid resistance in Escherichia coli

Laboratory‐based evolution and whole‐genome sequencing can link genotype and phenotype. We used evolution of acid resistance in exponential phase Escherichia coli to study resistance to a lethal stress. Iterative selection at pH 2.5 generated five populations that were resistant to low pH in early exponential phase. Genome sequencing revealed multiple mutations, but the only gene mutated in all strains was evgS, part of a two‐component system that has already been implicated in acid resistance. All these mutations were in the cytoplasmic PAS domain of EvgS, and were shown to be solely responsible for the resistant phenotype, causing strong upregulation at neutral pH of genes normally induced by low pH. Resistance to pH 2.5 in these strains did not require the transporter GadC, or the sigma factor RpoS. We found that EvgS‐dependent constitutive acid resistance to pH 2.5 was retained in the absence of the regulators GadE or YdeO, but was lost if the oxidoreductase YdeP was also absent. A deletion in the periplasmic domain of EvgS abolished the response to low pH, but not the activity of the constitutive mutants. On the basis of these results we propose a model for how EvgS may become activated by low pH.

A molecular model of chromatin organisation and transcription: how a multi‐RNA polymerase II machine transcribes and remodels the β‐globin locus during development

We present a molecular model of eukaryotic gene transcription. For the β‐globin locus, we hypothesise that a transcription machine composed of multiple RNA polymerase II (PolII) assembles using the locus control region as a foundation. Transcription and locus remodelling can be achieved by pulling DNA through this multi‐PolII ‘reading head’. Once a transcription complex is formed, it may engage an active gene in several rounds of transcription. Observed intergenic sense and antisense transcripts may be the result of PolII pulling the DNA through the reading head whilst searching for the promoter of a gene. Support for this hypothesis is provided using various data from the literature. In the model, DNA is packed in a 30‐nm chromatin fibre, thus gene regulatory regions separated by kilobases are close in space. This, and the need to store transcription‐induced supercoiling, may explain why functionally interacting regions are often separated by many kilobases.

The Radical Ion Chemistry of S-Nitrosylated Peptides

The radical ion chemistry of a suite of S-nitrosopeptides has been investigated. Doubly and triply-protonated ions of peptides NYCGLPGEYWLGNDK, NYCGLPGEYWLGNDR, NYCGLPGERWLGNDR, NACGAPGEKWAGNDK, NYCGLPGEKYLGNDK, NYGLPGCEKWYGNDK and NYGLPGEKWYGCNDK were subjected to electron capture dissociation (ECD), and collision-induced dissociation (CID). The peptide sequences were selected such that the effect of the site of S-nitrosylation, the nature and position of the basic amino acid residues, and the nature of the other amino acid side chains, could be interrogated. The ECD mass spectra were dominated by a peak corresponding to loss of •NO from the charge-reduced precursor, which can be explained by a modified Utah-Washington mechanism. Some backbone fragmentation in which the nitrosyl modification was preserved was also observed in the ECD of some peptides. Molecular dynamics simulations of peptide ion structure suggest that the ECD behavior was dependent on the surface accessibility of the protonated residue. CID of the S-nitrosylated peptides resulted in homolysis of the S–N bond to form a long-lived radical with loss of •NO. The radical peptide ions were isolated and subjected to ECD and CID. ECD of the radical peptide ions provided an interesting comparison to ECD of the unmodified peptides. The dominant process was electron capture without further dissociation (ECnoD). CID of the radical peptide ions resulted in cysteine, leucine, and asparagine side chain losses, and radical-induced backbone fragmentation at tryptophan, tyrosine, and asparagine residues, in addition to charge-directed backbone fragmentation.

Structural and electrostatic properties of ubiquitination and related pathways.

Post-translational modification by ubiquitin and ubiquitin-like (UBL) proteins is a key mechanism for cellular control. The specificity of the enzymes of ubiquitination and their close paralogs is dependent on their molecular electrostatic potentials. For example, analysis of molecular electrostatic potentials and electrostatically key residues can account for the selectivity of different E1s (activating enzymes) and of different SUMO proteases. The molecular interactions of the ubiquitin conjugating enzymes, the ubiquitin family proteins (UFP) and UBL domains are discussed in detail. An interesting observation is that the Non Canonical Ubiquitin Conjugating Enzymes (NCUBEs) have electrostatic potentials that are more similar to the UBC9 orthologs, the SUMO conjugating enzymes, than they are to other ubiquitin conjugating enzymes. It had previously been suggested that UBC9 may select for SUMO based on its difference in electrostatic potential as compared to other E2s but the NCUBE exception suggests that this may not be the case. The web site http://www.ubiquitin-resource.org/ allows users to find the E2s most electrostatically similar to a query E2. Where possible, models have been made for all E2 domains in the SMART database (http://smart.embl-heidelberg.de/). A brief overview of molecular electrostatic potentials and their application to understanding protein function is also given.