Simulated annealing with restrained molecular dynamics using CONGEN: energy refinement of the NMR solution structures of epidermal and type-alpha transforming growth factors.

The new functionality of the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168; Bassolino-Klimas D et al., 1996, Protein Sci 5:593-603) has been applied for energy refinement of two previously determined solution NMR structures, murine epidermal growth factor (mEGF) and human type-alpha transforming growth factor (hTGF alpha). A summary of considerations used in converting experimental NMR data into distance constraints for CONGEN is presented. A general protocol for simulated annealing with restrained molecular dynamics is applied to generate NMR solution structures using CONGEN together with real experimental NMR data. A total of 730 NMR-derived constraints for mEGF and 424 NMR-derived constraints for hTGF alpha were used in these energy-refinement calculations. Different weighting schemes and starting conformations were studied to check and/or improve the sampling of the low-energy conformational space that is consistent with all constraints. The results demonstrate that loosened (i.e., "relaxed") sets of the EGF and hTGF alpha internuclear distance constraints allow molecules to overcome local minima in the search for a global minimum with respect to both distance restraints and conformational energy. The resulting energy-refined structures of mEGF and hTGF alpha are compared with structures determined previously and with structures of homologous proteins determined by NMR and X-ray crystallography.     

Combined use of 13C chemical shift and 1Hα−13Cα heteronuclear NOE data in monitoring a protein NMR structure refinement

Summary A large portion of the ¹³C resonance assignments for murine epidermal growth factor (mEGF) at pH 3.1 and 28°C has been determined at natural isotope abundance. Sequence-specific ¹³C assignments are reported for 100% of the assignable C, 96% of the C, 86% of the aromatic and 70% of the remaining peripheral aliphatic resonances of mEGF. A good correlation was observed between experimental and back-calculated C chemical shifts for regions of regular -sheet structure. These assignments also provide the basis for interpreting ¹H–¹³C heteronuclear NOE (HNOE) values in mEGF at natural isotope abundance. Some of the backbone polypeptide segments with high internal mobility, indicated by these ¹H–¹³C HNOE measurements, correlate with locations of residues involved in the putative mEGF-receptor binding site. Using four families of mEGF structures obtained over the last few years, we demonstrate that standard deviations between experimental and back-calculated C values can be used to monitor the refinement of this protein's structure, particularly for -sheet regions. Improved agreement between calculated and observed values of C is correlated with other measures of structure quality, including lowered values of residual constraint violations and more negative values of conformational energy. These results support the view that experimental conformation-dependent chemical shifts, C, can provide a reliable source of information for monitoring the process of protein structure refinement and are potentially useful restraints for driving the refinement.Abbreviations HSQC heteronuclear single-quantum coherence spectroscopy - PFG pulsed-field gradient - TOCSY ¹H-¹H total correlation spectroscopy - EGF epidermal growth factor - mEGF murine EGF - hEGF human EGF - hTGF human type- transforming growth factor - DIPSI spm-locking pulse sequence - NOE nuclear Overhauser effect - HNOE heteronuclear Overhauser effect     

Human epidermal growth factor. Distinct roles of tyrosine 37 and arginine 41 in receptor binding as determined by site-directed mutagenesis and nuclear magnetic resonance spectroscopy

Site-directed mutagenesis was employed to examine the function of two highly conserved residues, Tyr-37 and Arg-41, of human EGF (hEGF) in receptor binding. Both a conservative change to phenylalanine and a semi-conservative change to histidine at position 37 yield proteins with receptor affinity similar to wild-type hEGF. A non-conservative change to alanine results in a molecule with about 40% of the receptor affinity, indicating that an aromatic residue is not essential at this position. Both conservative (to lysine) and non-conservative (to alanine) substitutions at position 41 drastically reduced receptor binding to less than 0.5% of the wild-type activity. 1D-NMR data indicate that the replacement of Arg-41 by lysine does not significantly alter the native protein conformation. Thus, Arg-41 may be directly involved in ligand receptor interaction, whereas the side chain of Tyr-37, although possibly important structurally, is not essential for receptor binding.     

The New York Consortium on Membrane Protein Structure (NYCOMPS): a high-throughput platform for structural genomics of integral membrane proteins

The New York Consortium on Membrane Protein Structure (NYCOMPS) was formed to accelerate the acquisition of structural information on membrane proteins by applying a structural genomics approach. NYCOMPS comprises a bioinformatics group, a centralized facility operating a high-throughput cloning and screening pipeline, a set of associated wet labs that perform high-level protein production and structure determination by x-ray crystallography and NMR, and a set of investigators focused on methods development. In the first three years of operation, the NYCOMPS pipeline has so far produced and screened 7,250 expression constructs for 8,045 target proteins. Approximately 600 of these verified targets were scaled up to levels required for structural studies, so far yielding 24 membrane protein crystals. Here we describe the overall structure of NYCOMPS and provide details on the high-throughput pipeline.     

Computational Design of a Protein-Based Enzyme Inhibitor

While there has been considerable progress in designing protein–protein interactions, the design of proteins that bind polar surfaces is an unmet challenge. We describe the computational design of a protein that binds the acidic active site of hen egg lysozyme and inhibits the enzyme. The design process starts with two polar amino acids that fit deep into the enzyme active site, identifies a protein scaffold that supports these residues and is complementary in shape to the lysozyme active-site region, and finally optimizes the surrounding contact surface for high-affinity binding. Following affinity maturation, a protein designed using this method bound lysozyme with low nanomolar affinity, and a combination of NMR studies, crystallography, and knockout mutagenesis confirmed the designed binding surface and orientation. Saturation mutagenesis with selection and deep sequencing demonstrated that specific designed interactions extending well beyond the centrally grafted polar residues are critical for high-affinity binding.     

Solution NMR structure and folding dynamics of the N terminus of a rat non-muscle alpha-tropomyosin in an engineered chimeric protein

Tropomyosin is an alpha-helical coiled-coil protein that aligns head-to-tail along the length of the actin filament and regulates its function. The solution structure of the functionally important N terminus of a short 247-residue non-muscle tropomyosin was determined in an engineered chimeric protein, GlyTM1bZip, consisting of the first 19 residues of rat short alpha-tropomyosin and the last 18 residues of the GCN4 leucine zipper. A gene encoding GlyTM1bZip was synthesized, cloned and expressed in Escherichia coli. Triple resonance NMR spectra were analyzed with the program AutoAssign to assign its backbone resonances. Multidimensional nuclear Overhauser effect spectra, X-filtered spectra and (3)J(H(N)-H(alpha)) scalar coupling were analyzed using AutoStructure. This is the first application of this new program to determine the three-dimensional structure of a symmetric homodimer and a structure not previously reported. Residues 7-35 in GlyTM1bZip form a coiled coil, but neither end is helical. Heteronuclear (15)N-(1)H nuclear Overhauser effect data showed that the non-helical N-terminal residues are flexible. The (13)C' chemical shifts of the coiled-coil backbone carbonyl groups in GlyTM1bZip showed a previously unreported periodicity, where resonances arising from residues at the coiled-coil interface in a and d positions of the heptad repeat were displaced relatively upfield and those arising from residues in c positions were displaced relatively downfield. Heteronuclear single quantum coherence spectra, collected as a function of temperature, showed that cross-peaks arising from the alpha-helical backbone and side-chains at the coiled-coil interface broadened or shifted with T(M) values approximately 20 degrees C lower than the loss of alpha-helix measured by circular dichroism, suggesting the presence of a folding intermediate. The side-chain of Ile14, a residue essential for binding interactions, exhibited multiple conformations. The conformational flexibility of the N termini of short tropomyosins may be important for their binding specificity.     

Cold-shock induced high-yield protein production in Escherichia coli

Overexpression of proteins in Escherichia coli at low temperature improves their solubility and stability. Here, we apply the unique features of the cspA gene to develop a series of expression vectors, termed pCold vectors, that drive the high expression of cloned genes upon induction by cold-shock. Several proteins were produced with very high yields, including E. coli EnvZ ATP-binding domain (EnvZ-B) and Xenopus laevis calmodulin (CaM). The pCold vector system can also be used to selectively enrich target proteins with isotopes to study their properties in cell lysates using NMR spectroscopy. We have cloned 38 genes from a range of prokaryotic and eukaryotic organisms into both pCold and pET14 (ref. 3) systems, and found that pCold vectors are highly complementary to the widely used pET vectors.     

Primer Prim'r: A web based server for automated primer design

The Northeast Structural Genomics Consortium (NESG) is one of nine NIH-funded pilot projects created to develop technologies needed for structural studies of proteins on a genome-wide scale. One of the most challenging aspects of this emerging field is the production of protein samples amenable to structural determination. To do this efficiently, all steps in the protein production pipeline must be automated. Here we describe the Primer program (linked from http://www-nmr.cabm.rutgers.edu/bioinformatics, www-nmr.cabm.rutgers.edu/bioinformatics, a web-based primer design program freely available to the scientific community, which was created to automate this time consuming and laborious task. This program has the ability to simultaneously calculate plasmid specific primer sets for multiple open reading frame (ORF) targets, including 96-well and greater formats. Primer includes a library of commonly used plasmid systems and possesses the ability to upload user-defined plasmid systems. In addition to calculating gene-specific annealing regions for each target, the program also adds appropriate restriction endonuclease recognition or viral recombination sites while preserving a reading frame with plasmid based fusions. Primer has several useful features such as sorting calculated primer sets by target size, facilitating interpretation of PCR amplifications by agarose gel electrophoresis, as well as supplying the molecular biologist with many important characteristics of each target such as the expected size of the PCR amplified DNA fragment and internal restriction sites. The NESG has cloned over 1500 genes using oligonucleotide primers designed by Primer.