Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5

Disulfide bridge formation to increase thermostability of DFPase enzyme: A computational study

Organophosphate compounds bioremediation by use of organophosphorus degradation enzymes such as DFPase is a developing interest in industry and medicine. The most important problem with the bio-catalytic enzymes is their instability on high temperatures. This work carried out to find suitable locations for introducing disulfide bridges in DFPase enzyme. We employed some computational approaches to design the disulfide bridges and evaluate their roles in the enzyme structural thermostability. According to the in silico results, mutant 6 (V24C, C76) increased the enzyme thermostability relative to wild-type.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

New Dynamic Rotamer Libraries: Data-Driven Analysis of Side-Chain Conformational Propensities

Most rotamer libraries are generated from subsets of the PDB and do not fully represent the conformational scope of protein side chains. Previous attempts to rectify this sparse coverage of conformational space have involved application of weighting and smoothing functions. We resolve these limitations by using physics-based molecular dynamics simulations to determine more accurate frequencies of rotameric states. This work forms part of our Dynameomics initiative and uses a set of 807 proteins selected to represent 97% of known autonomous protein folds, thereby eliminating the bias toward common topologies found within the PDB. Our Dynameomics derived rotamer libraries encompass 4.8 × 10(9) rotamers, sampled from at least 51,000 occurrences of each of 93,642 residues. Here, we provide a backbone-dependent rotamer library, based on secondary structure ϕ/ψ regions, and an update to our 2011 backbone-independent library that addresses the doubling of our dataset since its original publication.

Biotechnology of cold-active proteases

The bulk of Earth's biosphere is cold (<5 °C) and inhabited by psychrophiles. Biocatalysts from psychrophilic organisms (psychrozymes) have attracted attention because of their application in the ongoing efforts to decrease energy consumption. Proteinases as a class represent the largest category of industrial enzymes. There has been an emphasis on employing cold-active proteases in detergents because this allows laundry operations at ambient temperatures. Proteases have been used in environmental bioremediation, food industry and molecular biology. In view of the present limited understanding and availability of cold-active proteases with diverse characteristics, it is essential to explore Earth's surface more in search of an ideal cold-active protease. The understanding of molecular and mechanistic details of these proteases will open up new avenues to tailor proteases with the desired properties. A detailed account of the developments in the production and applications of cold-active proteases is presented in this review.

Protein folds and protein folding

The classification of protein folds is necessarily based on the structural elements that distinguish domains. Classification of protein domains consists of two problems: the partition of structures into domains and the classification of domains into sets of similar structures (or folds). Although similar topologies may arise by convergent evolution, the similarity of their respective folding pathways is unknown. The discovery and the characterization of the majority of protein folds will be followed by a similar enumeration of available protein folding pathways. Consequently, understanding the intricacies of structural domains is necessary to understanding their collective folding pathways. We review the current state of the art in the field of protein domain classification and discuss methods for the systematic and comprehensive study of protein folding across protein fold space via atomistic molecular dynamics simulation. Finally, we discuss our large-scale Dynameomics project, which includes simulations of representatives of all autonomous protein folds.

Dynameomics: protein dynamics and unfolding across fold space

All currently known structures of proteins together define ‘protein fold space’. To increase the general understanding of protein dynamics and protein folding, we selected a set of 807 proteins and protein domains that represent 95% of the currently known autonomous folded domains present in globular proteins. Native state and unfolding simulations of these representatives are now complete and accessible via a novel database containing over 11 000 simulations. Because protein folding is a microscopically reversible process, these simulations effectively sample protein folding across all of protein fold space. Here, we give an overview of how the representative proteins were selected and how the simulations were performed and validated. We then provide examples of different types of analyses that can be performed across our large set of simulations, made possible by the database approach. We further show how the unfolding simulations can be used to compare unfolding of structural elements in isolation and in different structural contexts, using as an example a short, triple stranded β-sheet that forms the WW domain and is present in several larger unrelated proteins.

Conformational changes below the Tm: molecular dynamics studies of the thermal pretransition of ribonuclease A

Recent work suggests that some native conformations of proteins can vary with temperature. To obtain an atomic-level description of this structural and conformational variation, we have performed all-atom, explicit-solvent molecular dynamics simulations of bovine pancreatic ribonuclease A (RNase A) up to its melting temperature (Tm approximately 337 K). RNase A has a thermal pretransition near 320 K [Stelea, S. D., Pancoska, P., Benight, A. S., and Keiderling, T. A. (2001) Protein Sci. 10, 970-978]. Our simulations identify a conformational change that coincides with this pretransition. Between 310 and 320 K, there is a small but significant decrease in the number of native contacts, beta-sheet hydrogen bonding, and deviation of backbone conformation from the starting structure, and an increase in the number of nonnative contacts. Native contacts are lost in beta-sheet regions and in alpha1, partially due to movement of alpha1 away from the beta-sheet core. At 330 and 340 K, a nonnative helical segment of residues 15-20 forms, corresponding to a helix observed in the N-terminal domain-swapped dimer [Liu, Y. S., Hart, P. J., Schulnegger, M. P., and Eisenberg, D. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 3437-3432]. The conformations observed at the higher temperatures possess nativelike topology and overall conformation, with many native contacts, but they have a disrupted active site. We propose that these conformations may represent the native state at elevated temperature, or the N' state. These simulations show that subtle, functionally important changes in protein conformation can occur below the Tm.

Multi-template approach to modeling engineered disulfide bonds

The key issue for disulfide bond engineering is to select the most appropriate location in the protein. By surveying the structure of experimentally engineered disulfide bonds, we found about half of them that have geometry incompatible with any native disulfide bond geometry. To improve the current prediction methods that tend to apply either ideal geometrical or energetical criteria to single three-dimensional structures, we have combined a novel computational protocol with the usage of multiple protein structures to take into account protein backbone flexibility. The multiple structures can be selected from either independently determined crystal structures for identical proteins, models of nuclear magnetic resonance experiments, or crystal structures of homology-related proteins. We have validated our approach by comparing the predictions with known disulfide bonds. The accuracy of prediction for native disulfide bonds reaches 99.6%. In a more stringent test on the reported engineered disulfide bonds, we have obtained a success rate of 93%. Our protocol also determines the oxido-reduction state of a predicted disulfide bond and the corresponding mutational cost. From the energy ranking, the user can easily choose top predicted sites for mutagenesis experiments. Our method provides information about local stability of the engineered disulfide bond surroundings.

Enhancing the stability of microsomal cytochrome b5: a rational approach informed by comparative studies with the outer mitochondrial membrane isoform

The outer mitochondrial membrane isoform of mammalian cytochrome b5 (OM b5) is much less prone to lose heme than the microsomal isoform (Mc b5), with a conserved difference at position 71 (leucine versus serine) playing a major role. We replaced Ser71 in Mc b5 with Leu, with the prediction that it would retard heme loss by diminishing polypeptide expansion accompanying rupture of the histidine to iron bonds. The strategy was partially successful in that it slowed dissociation of heme from its less stable orientation in bMc b5 (B). Heme dissociation from orientation A was accelerated to a similar extent, however, apparently owing to increased binding pocket dynamic mobility related to steric strain. A second mutation (L32I) guided by results of previous comparative studies of Mc and OM b5s diminished the steric strain, but much greater relief was achieved by replacing heme with iron deuteroporphyrin IX (FeDPIX). Indeed, the stability of the Mc(S71L) b5 FeDPIX complex is similar to that of the FeDPIX complex of OM b5. The results suggest that maximizing heme binding pocket compactness in the apo state is a useful general strategy for increasing the stability of engineered or designed proteins.

Insertion of the cytochrome b5 heme-binding loop into an SH3 domain. Effects on structure and stability, and clues about the cytochrome's architecture

Under native conditions, apocytochrome b(5) exhibits a stable core and a disordered heme-binding region that refolds upon association with the cofactor. The termini of this flexible region are in close proximity, suggesting that loop closure may contribute to the thermodynamic properties of the apocytochrome. A chimeric protein containing 43 residues encompassing the cytochrome loop was constructed using the cyanobacterial photosystem I accessory protein E (PsaE) from Synechococcus sp. PCC 7002 as a structured scaffold. PsaE has the topology of an SH3 domain, and the insertion was engineered to replace its 14-residue CD loop. NMR and optical spectroscopies showed that the hybrid protein (named EbE1) was folded under native conditions and that it retained the characteristics of an SH3 domain. NMR spectroscopy revealed that structural and dynamic differences were confined near the site of loop insertion. Variable-temperature 1D NMR spectra of EbE1 confirmed the presence of a kinetic unfolding barrier. Thermal and chemical denaturations of PsaE and EbE1 demonstrated cooperative, two-state transitions; the stability of the PsaE scaffold was found only moderately compromised by the insertion, with a DeltaT(m) of 8.3 degrees C, a DeltaC(m) of 1.5 M urea, and a DeltaDeltaG degrees of 4.2 kJ/mole. The data implied that the penalty for constraining the ends of the inserted region was lower than the approximately 6.4 kJ/mole calculated for a self-avoiding chain. Extrapolation of these results to cytochrome b(5) suggested that the intrinsic stability of the folded portion of the apoprotein reflected only a small detrimental contribution from the large heme-binding domain.

House fly cytochrome b5 exhibits kinetically trapped hemin and selectivity in hemin binding

We report that cytochrome b(5) (cyt b(5)) from Musca domestica (house fly) is more thermally stable than all other microsomal (Mc) cytochromes b(5) that have been examined to date. It also exhibits a much higher barrier to equilibration of the two isomeric forms of the protein, which differ by a 180 degrees rotation about the alpha-gamma-meso axis of hemin (ferric heme). In fact, hemin is kinetically trapped in a nearly statistical 1.2:1 ratio of rotational forms in freshly expressed protein. The equilibrium ratio (5.5:1) is established only upon incubation at temperatures above 37 degrees C. House fly Mc cyt b(5) is only the second b-hemoprotein that has been shown to exhibit kinetically trapped hemin at room temperature or above, the first being cyt b(5) from the outer membrane of rat liver mitochondria (rat OM cyt b(5)). Finally, we show that the small excess of one orientational isomer over the other in freshly expressed protein results from selective binding of hemin by the apoprotein, a phenomenon that has not heretofore been established for any apocyt b(5).

Conformational flexibility in sigma70 region 2 during transcription initiation

Prokaryotic RNA polymerase holoenzyme is composed of core subunits (alpha(2)betabeta'omega) plus a sigma factor that confers promoter specificity allowing for regulation of gene expression. Holoenzyme is known to undergo several conformational changes during the multiple steps of transcription initiation. However, the effects of these changes on the functions of specific regions have not been well characterized. In this work, we addressed the role of possible conformational change in region 2 of Escherichia coli sigma(70) by engineering disulfide bonds that "lock" region 2.1 with region 2.2 and region 2.2 with region 2.3. When these mutant holoenzymes were characterized for gross defects in multiple-round transcription, we found that insertion of either disulfide bond did not result in a fundamental block, indicating that the disulfide-containing holoenzymes are active. However, both disulfide-containing holoenzymes exhibited defects in formation and stability of the open complex. Our results suggest that conformational flexibility within sigma(70) region 2 facilitates open complex formation and transcription initiation.

Using disulfide bond engineering to study conformational changes in the beta'260-309 coiled-coil region of Escherichia coli RNA polymerase during sigma(70) binding

RNA polymerase of Escherichia coli is the sole enzyme responsible for mRNA synthesis in the cell. Upon binding of a sigma factor, the holoenzyme can direct transcription from specific promoter sequences. We have previously defined a region of the beta' subunit (beta'260-309, amino acids 260 to 309) which adopts a coiled-coil conformation shown to interact with sigma(70) both in vitro and in vivo. However, it was not known if the coiled-coil conformation was maintained upon binding to sigma(70). In this work, we engineered a disulfide bond within beta'240-309 that locks the beta' coiled-coil region in the coiled-coil conformation, and we show that this "locked" peptide is able to bind to sigma(70). We also show that the locked coiled-coil is capable of inducing a conformational change within sigma(70) that allows recognition of the -10 nontemplate strand of DNA. This suggests that the coiled-coil does not adopt a new conformation upon binding sigma(70) or upon recognition of the -10 nontemplate strand of DNA.

Structure and characterization of Ectothiorhodospira vacuolata cytochrome b(558), a prokaryotic homologue of cytochrome b(5)

A soluble cytochrome b(558) from the purple phototropic bacterium Ectothiorhodospira vacuolata was completely sequenced by a combination of automated Edman degradation and mass spectrometry. The protein, with a measured mass of 10,094.7 Da, contains 90 residues and binds a single protoheme. Unexpectedly, the sequence shows homology to eukaryotic cytochromes b(5). As no prokaryotic homologue had been reported so far, we developed a protocol for the expression, purification, and crystallization of recombinant cytochrome b(558). The structure was solved by molecular replacement to a resolution of 1.65 A. It shows that cytochrome b(558) is indeed the first bacterial cytochrome b(5) to be characterized and differs from its eukaryotic counterparts by the presence of a disulfide bridge and a four-residue insertion in front of the sixth ligand (histidine). Eukaryotes contain a variety of b(5) homologues, including soluble and membrane-bound multifunctional proteins as well as multidomain enzymes such as sulfite oxidase, fatty-acid desaturase, nitrate reductase, and lactate dehydrogenase. A search of the Mycobacterium tuberculosis genome showed that a previously unidentified gene encodes a fatty-acid desaturase with an N-terminal b(5) domain. Thus, it may provide another example of a bacterial b(5) homologue.