Biomedical subjects

Jeffery G Saven

Publications and source records attributed to Jeffery G Saven.

17 recordsLinked to original sources

PubMed2007

Combinatorial protein design strategies using computational methods.

Computational methods continue to facilitate efforts in protein design. Most of this work has focused on searching sequence space to identify one or a few sequences compatible with a given structure and functionality. Probabilistic computational methods provide information regarding the range of amino acid variability permitted by desired functional and structural constraints. Such methods may be used to guide the construction of both individual sequences and combinatorial libraries of proteins.

Amino Acid Sequence↗

PubMed2006

Design of functional ferritin-like proteins with hydrophobic cavities.

Ferritin four-helix bundle subunits self-assemble to create a stable multimer with a large central hydrophilic cavity where metal ions bind. To explore the versatility of this reaction vessel, computational design was used to generate cavities with increasingly apolar surface areas inside a dodecameric ferritin-like protein, Dps. Cavity mutants, in which as many as 120 surface accessible hydrophilic residues were replaced with hydrophobic amino acids, were shown to still assemble properly using size-exclusion chromatography and dynamic light scattering measurements. Wild-type Dps exhibited highly cooperative subunit folding and assembly, which was monitored by changes in Trp fluorescence and UV circular dichroism. The hydrophobic cavity mutants showed distinctly less cooperative unfolding behavior, with one mutant forming a partially assembled intermediate upon guanidine denaturation. Although the stability of Dps to such denaturation decreased with increasing apolar surface area, all proteins exhibited high melting temperatures, T(m) = 74-90 degrees C. Despite the large number of mutations, near-native ability to mineralize iron was maintained. This work illustrates the versatility of the ferritin scaffold for engineering large protein cavities with novel properties.

Bacterial Proteins↗

PubMed2006

Limitations of yeast surface display in engineering proteins of high thermostability.

Engineering proteins that can fold to unique structures remains a challenge. Protein stability has previously been engineered via the observed correlation between thermal stability and eukaryotic secretion level. To explore the limits of an expression-based approach, variants of the highly thermostable three-helix bundle protein alpha3D were studied using yeast surface display. A library of alpha3D mutants was created to explore the possible correlation of protein stability and fold with expression level. Five efficiently expressed mutants were then purified and further studied biochemically. Despite their differences in stability, most mutants expressed at levels comparable with that of wild-type alpha3D. Two other related sequences (alpha3A and alpha3B) that form collapsed, stable molten globules but lack a uniquely folded structure were similarly expressed at high levels by yeast display. Together these observations suggest that the quality control system in yeast is unable to discriminate between well-folded proteins of high stability and molten globules. The present study, therefore, suggests that an optimization of the surface display efficiency on yeast may yield proteins that are thermally and chemically stable yet are poorly folded.

Amino Acid Sequence↗

PubMed2006

Ultrafast folding of a computationally designed Trp-cage mutant: Trp2-cage.

Miniproteins provide useful model systems for understanding the principles of protein folding and design. These proteins also serve as useful test cases for theories of protein folding, and their small size and ultrafast folding kinetics put them in a regime of size and time scales that is now becoming accessible to molecular dynamics simulations. Previous estimates have suggested the "speed limit" for folding is on the order of 1 mus. Here a computationally designed mutant of the 20-residue Trp-cage miniprotein, Trp2-cage, is presented. The Trp2-cage has greater stability than the parent and folds on the ultrafast time scale of 1 mICROs at room temperature, as determined from infrared temperature-jump experiments.

Computer Simulation↗

PubMed2005

Simulation of pH-dependent edge strand rearrangement in human beta-2 microglobulin.

Amyloid fibrils formed from unrelated proteins often share morphological similarities, suggesting common biophysical mechanisms for amyloidogenesis. Biochemical studies of human beta-2 microglobulin (beta2M) have shown that its transition from a water-soluble protein to insoluble aggregates can be triggered by low pH. Additionally, biophysical measurements of beta2M using NMR have identified residues of the protein that participate in the formation of amyloid fibrils. The crystal structure of monomeric human beta2M determined at pH 5.7 shows that one of its edge beta-strands (strand D) adopts a conformation that differs from other structures of the same protein obtained at higher pH. This alternate beta-strand arrangement lacks a beta-bulge, which may facilitate protein aggregation through intermolecular beta-sheet association. To explore whether the pH change may yield the observed conformational difference, molecular dynamics simulations of beta2M were performed. The effects of pH were modeled by specifying the protonation states of Asp, Glu, and His, as well as the C terminus of the main chain. The bulged conformation of strand D is preferred at medium pH (pH 5-7), whereas at low pH (pH < 4) the straight conformation is observed. Therefore, low pH may stabilize the straight conformation of edge strand D and thus increase the amyloidogenicity of beta2M.

Amyloid↗

PubMed2005

Molecular Dynamics Simulation of WSK-3, a Computationally Designed, Water-Soluble Variant of the Integral Membrane Protein KcsA.

Poor solubility and low expression levels often make membrane proteins difficult to study. An alternative to the use of detergents to solubilize these aggregation-prone proteins is the partial redesign of the sequence so as to confer water solubility. Recently, computationally assisted membrane protein solubilization (CAMPS) has been reported, where exposed hydrophobic residues on a protein's surface are computationally redesigned. Herein, the structure and fluctuations of a designed, water-soluble variant of KcsA (WSK-3) were studied using molecular dynamics simulations. The root mean square deviation of the protein from its starting structure, where the backbone coordinates are those of KcsA, was 1.8 angstroms. The structure of salt bridges involved in structural specificity and solubility were examined. The preferred configuration of ions and water in the selectivity filter of WSK-3 was consistent with the reported preferences for KcsA. The structure of the selectivity filter was maintained, which is consistent with WSK-3 having an affinity for agitoxin2 comparable to that of wild-type KcsA. In contrast to KcsA, the central cavity's side chains were observed to reorient, allowing water diffusion through the side of the cavity wall. These simulations provide an atomistic analysis of the CAMPS strategy and its implications for further investigations of membrane proteins.

Bacterial Proteins↗

PubMed2005

Statistical theory for protein ensembles with designed energy landscapes.

Combinatorial protein libraries provide a promising route to investigate the determinants and features of protein folding and to identify novel folding amino acid sequences. A library of sequences based on a pool of different monomer types are screened for folding molecules, consistent with a particular foldability criterion. The number of sequences grows exponentially with the length of the polymer, making both experimental and computational tabulations of sequences infeasible. Herein a statistical theory is extended to specify the properties of sequences having particular values of global energetic quantities that specify their energy landscape. The theory yields the site-specific monomer probabilities. A foldability criterion is derived that characterizes the properties of sequences by quantifying the energetic separation of the target state from low-energy states in the unfolded ensemble and the fluctuations of the energies in the unfolded state ensemble. For a simple lattice model of proteins, excellent agreement is observed between the theory and the results of exact enumeration. The theory may be used to provide a quantitative framework for the design and interpretation of combinatorial experiments.

Models, Statistical↗

PubMed2005

Statistical and molecular dynamics studies of buried waters in globular proteins.

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.

Amino Acids↗

PubMed2005

De novo design of a redox-active minimal rubredoxin mimic.

Metal-binding sites in metalloproteins frequently occur at the interfaces of elements of secondary structure, which has enabled the retrostructural analysis of natural proteins and the de novo design of helical bundles that bind metal ion cofactors. However, the design of metalloproteins containing beta-structure is less well developed, despite the frequent occurrence of beta-conformations in natural metalloproteins. Here, we describe the design and construction of a beta-protein, RM1, that forms a stable, redox-active 4-Cys thiolate Fe(II/III) site analogous to the active site of rubredoxin. The protein folds into a beta-structure in the presence and absence of metal ions and binds Fe(II/III) to form a redox-active site that is stable to repeated cycles of oxidation and reduction, even in an aerobic environment.

Binding Sites↗

PubMed2005

Modulating the DNA affinity of Elk-1 with computationally selected mutations.

In order to regulate gene expression, transcription factors must first bind their target DNA sequences. The affinity of this binding is determined by both the network of interactions at the interface and the entropy change associated with the complex formation. To study the role of structural fluctuation in fine-tuning DNA affinity, we performed molecular dynamics simulations of two highly homologous proteins, Elk-1 and SAP-1, that exhibit different sequence specificity. Simulation studies show that several residues in Elk have significantly higher main-chain root-mean-square deviations than their counterparts in SAP. In particular, a single residue, D69, may contribute to Elk's lower DNA affinity for P(c-fos) by structurally destabilizing the carboxy terminus of the recognition helix. While D69 does not contact DNA directly, the increased mobility in the region may contribute to its weaker binding. We measured the ability of single point mutants of Elk to bind P(c-fos) in a reporter assay, in which D69 of wild-type Elk has been mutated to other residues with higher helix propensity in order to stabilize the local conformation. The gains in transcriptional activity and the free energy of binding suggested from these measurements correlate well with stability gains computed from helix propensity and charge-macrodipole interactions. The study suggests that residues that are distal to the binding interface may indirectly modulate the binding affinity by stabilizing the protein scaffold required for efficient DNA interaction.

Amino Acid Sequence↗

PubMed2005

Computational de novo design and characterization of a four-helix bundle protein that selectively binds a nonbiological cofactor.

We report the complete de novo design of a four-helix bundle protein that selectively binds the nonbiological DPP-Fe(III) metalloporphyrin cofactor (DPP-Fe(III) = 5, 15-Di[(4-carboxymethyleneoxy)phenyl]porphinato iron(III)). A tetrameric, D2-symmetric backbone scaffold was constructed to encapsulate two DPP-Fe(III) units through bis(His) coordination. The complete sequence was determined with the aid of the statistical computational design algorithm SCADS. The 34-residue peptide was chemically synthesized. UV-vis and CD spectroscopy, size-exclusion chromatography, and analytical ultracentrifugation indicated the peptide undergoes a transition from a predominantly random coil monomer to an alpha-helical tetramer upon binding DPP-Fe(III). EPR spectroscopy studies indicated the axial imidazole ligands were oriented in a perpendicular fashion, as defined by second-shell interactions that were included in the design. The 1-D 1H NMR spectrum of the assembled protein displayed features of a well-packed interior. The assembled protein possessed functional redox properties different from those of structurally similar systems containing the heme cofactor. The designed peptide demonstrated remarkable cofactor selectivity with a significantly weaker binding affinity for the natural heme cofactor. These findings open a path for the selective incorporation of more elaborate cofactors into designed scaffolds for constructing molecularly well-defined nanoscale materials.

Amino Acid Sequence↗

PubMed2004

Computational design of water-soluble analogues of the potassium channel KcsA.

Although the interiors of membrane and water-soluble proteins are similar in their physicochemical properties, membrane proteins differ in having larger fractions of hydrophobic residues on their exteriors. Thus, it should be possible to water-solubilize membrane proteins by mutating their lipid-contacting side chains to more polar groups. Here, a computational approach was used to generate water-soluble variants of the potassium channel KcsA. As a probe of the correctness of the fold, the proteins contain an agitoxin2 binding site from a mammalian homologue of the channel. The resulting proteins express in high yield in Escherichia coli and share the intended functional and structural properties with KcsA, including secondary structure, tetrameric quaternary structure, and tight specific binding to both agitoxin2 and a small molecule channel blocker.

Amino Acid Sequence↗

PubMed2004

Advances in computational protein design.

Computational protein design continues to experience a variety of methodological advances. Several improvements have been suggested for the objective functions used to quantify sequence/structure compatibility. Disparate design strategies based upon dead-end elimination, simulated annealing and statistical design have each recently yielded striking successes involving de novo designed proteins with sizes on the order of 100 residues or greater. Such methods may be used to design new proteins, as well as to redesign natural proteins to facilitate structural and biophysical studies.

Computational Biology↗

PubMed2003

Computational design and characterization of a monomeric helical dinuclear metalloprotein.

The de novo design of di-iron proteins is an important step towards understanding the diversity of function among this complex family of metalloenzymes. Previous designs of due ferro (DF) proteins have resulted in tetrameric and dimeric four-helix bundles having crystallographically well-defined structures and active-site geometries. Here, the design and characterization of DFsc, a 114 residue monomeric four-helix bundle, is presented. The backbone was modeled using previous oligomeric structures and appropriate inter-helical turns. The identities of 26 residues were predetermined, including the primary and secondary ligands in the active site, residues involved in active site accessibility, and the gamma beta gamma beta turn between helices 2 and 3. The remaining 88 amino acid residues were determined using statistical computer aided design, which is based upon a recent statistical theory of protein sequences. Rather than sampling sequences, the theory directly provides the site-specific amino acid probabilities, which are then used to guide sequence design. The resulting sequence (DFsc) expresses well in Escherichia coli and is highly soluble. Sedimentation studies confirm that the protein is monomeric in solution. Circular dichroism spectra are consistent with the helical content of the target structure. The protein is structured in both the apo and the holo forms, with the metal-bound form exhibiting increased stability. DFsc stoichiometrically binds a variety of divalent metal ions, including Zn(II), Co(II), Fe(II), and Mn(II), with micromolar affinities. 15N HSQC NMR spectra of both the apo and Zn(II) proteins reveal excellent dispersion with evidence of a significant structural change upon metal binding. DFsc is then a realization of complete de novo design, where backbone structure, activity, and sequence are specified in the design process.

Amino Acid Sequence↗

PubMed2003

Probabilistic approach to the design of symmetric protein quaternary structures.

Probabilistic methods have been developed that estimate the site-specific probabilities of the amino acids in sequences likely to fold to a particular target structure, and such information can be used to guide the de novo design of proteins and to probe sequence variability. An extension of these methods for the design of symmetric homo-oligomeric quaternary structures is presented. The theory is in excellent agreement with the results of studies on exactly solvable lattice models. Application to an atomically detailed representation of proteins verifies the utility of a symmetry assumption, which greatly simplifies and accelerates the calculations. The method may be applied to a wide variety of symmetric and periodic protein structures.

Amino Acids↗

PubMed2002

Designing gene libraries from protein profiles for combinatorial protein experiments.

Protein combinatorial libraries provide new ways to probe the determinants of folding and to discover novel proteins. Such libraries are often constructed by expressing an ensemble of partially random gene sequences. Given the intractably large number of possible sequences, some limitation on diversity must be imposed. A non-uniform distribution of nucleotides can be used to reduce the number of possible sequences and encode peptide sequences having a predetermined set of amino acid probabilities at each residue position, i.e., the amino acid sequence profile. Such profiles can be determined by inspection, multiple sequence alignment or physically-based computational methods. Here we present a computational method that takes as input a desired sequence profile and calculates the individual nucleotide probabilities among partially random genes. The calculated gene library can be readily used in the context of standard DNA synthesis to generate a protein library with essentially the desired profile. The fidelity between the desired profile and the calculated one coded by these partially random genes is quantitatively evaluated using the linear correlation coefficient and a relative entropy, each of which provides a measure of profile agreement at each position of the sequence. On average, this method of identifying such codon frequencies performs as well or better than other methods with regard to fidelity to the original profile. Importantly, the method presented here provides much better yields of complete sequences that do not contain stop codons, a feature that is particularly important when all or large fractions of a gene are subject to combinatorial mutation.

Amino Acids↗

PubMed2002

Combinatorial protein design.

Combinatorial protein libraries permit the examination of a wide range of sequences. Such methods are being used for denovo design and to investigate the determinants of protein folding. The exponentially large number of possible sequences, however, necessitates restrictions on the diversity of sequences in a combinatorial library. Recently, progress has been made in developing theoretical tools to bias and characterize the ensemble of sequences that fold into a given structure - tools that can be applied to the design and interpretation of combinatorial experiments.

Amino Acids↗