Stephen W. Michnick

All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins †

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

Atypical Membrane Topology and Heteromeric Function of Drosophila Odorant Receptors In Vivo

Drosophila olfactory sensory neurons (OSNs) each express two odorant receptors (ORs): a divergent member of the OR family and the highly conserved, broadly expressed receptor OR83b. OR83b is essential for olfaction in vivo and enhances OR function in vitro, but the molecular mechanism by which it acts is unknown. Here we demonstrate that OR83b heterodimerizes with conventional ORs early in the endomembrane system in OSNs, couples these complexes to the conserved ciliary trafficking pathway, and is essential to maintain the OR/OR83b complex within the sensory cilia, where odor signal transduction occurs. The OR/OR83b complex is necessary and sufficient to promote functional reconstitution of odor-evoked signaling in sensory neurons that normally respond only to carbon dioxide. Unexpectedly, unlike all known vertebrate and nematode chemosensory receptors, we find that Drosophila ORs and OR83b adopt a novel membrane topology with their N-termini and the most conserved loops in the cytoplasm. These loops mediate direct association of ORs with OR83b. Our results reveal that OR83b is a universal and integral part of the functional OR in Drosophila. This atypical heteromeric and topological design appears to be an insect-specific solution for odor recognition, making the OR/OR83b complex an attractive target for the development of highly selective insect repellents to disrupt olfactory-mediated host-seeking behaviors of insect disease vectors.

An in Vivo Map of the Yeast Protein Interactome

Protein interactions regulate the systems-level behavior of cells; thus, deciphering the structure and dynamics of protein interaction networks in their cellular context is a central goal in biology. We have performed a genome-wide in vivo screen for protein-protein interactions in Saccharomyces cerevisiae by means of a protein-fragment complementation assay (PCA). We identified 2770 interactions among 1124 endogenously expressed proteins. Comparison with previous studies confirmed known interactions, but most were not known, revealing a previously unexplored subspace of the yeast protein interactome. The PCA detected structural and topological relationships between proteins, providing an 8-nanometer–resolution map of dynamically interacting complexes in vivo and extended networks that provide insights into fundamental cellular processes, including cell polarization and autophagy, pathways that are evolutionarily conserved and central to both development and human health.

PKB/Akt modulates TGF-beta signalling through a direct interaction with Smad3

Transforming growth factor β (TGF-β) has a major role in cell proliferation, differentiation and apoptosis in many cell types. Integration of the TGF-β pathway with other signalling cascades that control the same cellular processes may modulate TGF-β responses. Here we report the discovery of a new functional link between TGF-β and growth factor signalling pathways, mediated by a physical interaction between the serine-threonine kinase PKB (protein kinase B)/Akt and the transcriptional activator Smad3. Formation of the complex is induced by insulin, but inhibited by TGF-β stimulation, placing PKB–Smad3 at a point of convergence between these two pathways. PKB inhibits Smad3 by preventing its phosphorylation, binding to Smad4 and nuclear translocation. In contrast, Smad3 does not inhibit PKB. Inhibition of Smad3 by PKB occurs through a kinase-activity-independent mechanism, resulting in a decrease in Smad3-mediated transcription and protection of cells against TGF-β-induced apoptosis. Consistently, knockdown of the endogenous PKB gene with small-interfering RNA (siRNA) has the opposite effect. Our results suggest a very simple mechanism for the integration of signals arising from growth-factor- and TGF-β-mediated pathways.

β-Lactamase protein fragment complementation assays as in vivo and in vitro sensors of protein-protein interactions

We have previously described a strategy for detecting protein–protein interactions based on protein interaction–assisted folding of rationally designed fragments of enzymes. We call this strategy the protein fragment complementation assay (PCA). Here we describe PCAs based on the enzyme TEM-1 β-lactamase (EC: 3.5.2.6), which include simple colorimetric in vitro assays using the cephalosporin nitrocefin and assays in intact cells using the fluorescent substrate CCF2/AM (ref. 6). Constitutive protein–protein interactions of the GCN4 leucine zippers and of apoptotic proteins Bcl2 and Bad, and the homodimerization of Smad3, were tested in an in vitro assay using cell lysates. With the same in vitro assay, we also demonstrate interactions of protein kinase PKB with substrate Bad. The in vitro assay is facile and amenable to high-throughput modes of screening with signal-to-background ratios in the range of 10:1 to 250:1, which is superior to other PCAs developed to date. Furthermore, we show that the in vitro assay can be used for quantitative analysis of a small molecule–induced protein interaction, the rapamycin-induced interaction of FKBP and yeast FRB (the FKBP-rapamycin binding domain of TOR (target of rapamycin)). The assay reproduces the known dissociation constant and number of sites for this interaction. The combination of in vitro colorimetric and in vivo fluorescence assays of β-lactamase in mammalian cells suggests a wide variety of sensitive and high-throughput large-scale applications, including in vitro protein array analysis of protein–protein or enzyme–protein interactions and in vivo applications such as clonal selection for cells expressing interacting protein partners.

A highly sensitive protein-protein interaction assay based on Gaussia luciferase

Protein-fragment complementation assays (PCAs) provide a general strategy to study the dynamics of protein-protein interactions in vivo and in vitro. The full potential of PCA requires assays that are fully reversible and sensitive at subendogenous protein expression levels. We describe a new assay that meets these criteria, based on the Gaussia princeps luciferase enzyme, demonstrating chemical reversal, and induction and inhibition of a key interaction linking insulin and TGFβ signaling.

Universal strategies in research and drug discovery based on protein-fragment complementation assays.

Changes in the interactions among proteins that participate in a biochemical pathway can reflect the immediate regulatory responses to intrinsic or extrinsic perturbations of the pathway. Thus, methods that allow for the direct detection of the dynamics of protein–protein interactions can be used to probe the effects of any perturbation on any pathway of interest. Here we describe experimental strategies — based on protein-fragment complementation assays (PCAs) — that can achieve this. PCA-based strategies can be used with or instead of traditional target-based drug discovery strategies to identify novel pathway-component proteins of therapeutic interest, to increase the quantity and quality of information about the actions of potential drugs, and to gain insight into the intricate networks that make up the molecular machinery of living cells.

Solution structure of FKBP, a rotamase enzyme and receptor for FK506 and rapamycin.

Immunophilins, when complexed to immunosuppressive ligands, appear to inhibit signal transduction pathways that result in exocytosis and transcription. The solution structure of one of these, the human FK506 and rapamycin binding protein (FKBP), has been determined by nuclear magnetic resonance (NMR). FKBP has a previously unobserved antiparallel beta-sheet folding topology that results in a novel loop crossing and produces a large cavity lined by a conserved array of aromatic residues; this cavity serves as the rotamase active site and drug-binding pocket. There are other significant structural features (such as a protruding positively charged loop and an apparently flexible loop) that may be involved in the biological activity of FKBP.

Chromatin- and Transcription-Related Factors Repress Transcription from within Coding Regions throughout the Saccharomyces cerevisiae Genome

Previous studies in Saccharomyces cerevisiae have demonstrated that cryptic promoters within coding regions activate transcription in particular mutants. We have performed a comprehensive analysis of cryptic transcription in order to identify factors that normally repress cryptic promoters, to determine the amount of cryptic transcription genome-wide, and to study the potential for expression of genetic information by cryptic transcription. Our results show that a large number of factors that control chromatin structure and transcription are required to repress cryptic transcription from at least 1,000 locations across the S. cerevisiae genome. Two results suggest that some cryptic transcripts are translated. First, as expected, many cryptic transcripts contain an ATG and an open reading frame of at least 100 codons. Second, several cryptic transcripts are translated into proteins. Furthermore, a subset of cryptic transcripts tested is transiently induced in wild-type cells following a nutritional shift, suggesting a possible physiological role in response to a change in growth conditions. Taken together, our results demonstrate that, during normal growth, the global integrity of gene expression is maintained by a wide range of factors and suggest that, under altered genetic or physiological conditions, the expression of alternative genetic information may occur.

An in vivo library-versus-library selection of optimized protein-protein interactions.

We describe a rapid and efficient in vivo library-versus-library screening strategy for identifying optimally interacting pairs of heterodimerizing polypeptides. Two leucine zipper libraries, semi-randomized at the positions adjacent to the hydrophobic core, were genetically fused to either one of two designed fragments of the enzyme murine dihydrofolate reductase (mDHFR), and cotransformed into Escherichia coli. Interaction between the library polypeptides reconstituted enzymatic activity of mDHFR, allowing bacterial growth. Analysis of the resulting colonies revealed important biases in the zipper sequences relative to the original libraries, which are consistent with selection for stable, heterodimerizing pairs. Using more weakly associating mDHFR fragments, we increased the stringency of selection. We enriched the best-performing leucine zipper pairs by multiple passaging of the pooled, selected colonies in liquid culture, as the best pairs allowed for better bacterial propagation. This competitive growth allowed small differences among the pairs to be amplified, and different sequence positions were enriched at different rates. We applied these selection processes to a library-versus-library sample of 2.0 × 106 combinations and selected a novel leucine zipper pair that may be appropriate for use in further in vivo heterodimerization strategies.