Search PubMed⌕ Search

Biomedical subjects

William F DeGrado

Publications and source records attributed to William F DeGrado.

At least 19 recordsLinked to original sources
PubMed2006

Amide vibrations are delocalized across the hydrophobic interface of a transmembrane helix dimer.

The tertiary interactions between amide-I vibrators on the separate helices of transmembrane helix dimers were probed by ultrafast 2D vibrational photon echo spectroscopy. The 2D IR approach proves to be a useful structural method for the study of membrane-bound structures. The 27-residue human erythrocyte protein Glycophorin A transmembrane peptide sequence: KKITLIIFG(79)VMAGVIGTILLISWG(94)IKK was labeled at G(79) and G(94) with (13)C=(16)O or (13)C=(18)O. The isotopomers and their 50:50 mixtures formed helical dimers in SDS micelles whose 2D IR spectra showed components from homodimers when both helices had either (13)C=(16)O or (13)C=(18)O substitution and a heterodimer when one had (13)C=(16)O substitution and the other had (13)C=(18)O substitution. The cross-peaks in the pure heterodimer 2D IR difference spectrum and the splitting of the homodimer peaks in the linear IR spectrum show that the amide-I mode is delocalized across a pair of helices. The excitation exchange coupling in the range 4.3-6.3 cm(-1) arises from through-space interactions between amide units on different helices. The angle between the two Gly(79) amide-I transition dipoles, estimated at 103 degrees from linear IR spectroscopy and 110 degrees from 2D IR spectroscopy, combined with the coupling led to a structural picture of the hydrophobic interface that is remarkably consistent with results from NMR on helix dimers. The helix crossing angle in SDS is estimated at 45 degrees. Two-dimensional IR spectroscopy also sets limits on the range of geometrical parameters for the helix dimers from an analysis of the coupling constant distribution.

Amides↗
PubMed2006

Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations.

Almost half of the proteome of living organisms is constituted of metalloproteins. Unfortunately, the ability of the current generation of molecular dynamics pairwise-additive forcefields to properly describe metal pockets is severely lacking due to the intrinsic difficulty of handling polarization and charge transfer contributions. In order to improve the description of metalloproteins, a simple reparameterization strategy is proposed herein that does not involve artificial constraints. Specifically, a non-bonded quantum mechanical-based model is used to capture the mean polarization and charge transfer contributions to the interatomic forces within the metal site. The present approach is demonstrated to provide enough accuracy to maintain the integrity of the metal pocket for a variety of metalloproteins during extended (multi-nanosecond) molecular dynamics simulations. The method enables the sampling of small conformational changes and the relaxation of local frustrations in NMR structures.

Binding Sites↗
PubMed2006

Truncation of a cross-linked GCN4-p1 coiled coil leads to ultrafast folding.

Structural perturbation has been extensively used in protein folding studies because it yields valuable conformational information regarding the folding process. Here we have used N-terminal truncation on a cross-linked variant of the GCN4-p1 leucine zipper, aiming to develop a better understanding of the folding mechanism of the coiled-coil motif. Our results indicate that removing the first heptad repeat in this cross-linked GCN4-p1 coiled coil significantly decreases the folding free energy barrier and results in a maximum folding rate of (2.0 +/- 0.3 micros)(-1), which is approximately 50 times faster than that of the full-length protein. Therefore, these results suggest that a set of native or nativelike tertiary interactions, distributed throughout the entire sequence, collectively stabilize the folding transition state of the GCN4-p1 coiled coil. While stable subdomains or triggering sequences have been shown to be critical to the stability of GCN4 coiled coils, our results suggest that the folding of such a subdomain does not seem to dictate the overall folding kinetics.

Amino Acid Sequence↗
PubMed2006

Arylamide derivatives as allosteric inhibitors of the integrin alpha2beta1/type I collagen interaction.

We herein report a group of allosteric inhibitors of integrin alpha(2)beta(1) based on an arylamide scaffold. Compound 4 showed an IC(50) of 4.80 microM in disrupting integrin I-domain/collagen binding in an ELISA. These arylamide compounds are able to block collagen binding to integrin alpha(2)beta(1) on the platelet surface. Further we find that compound 4 recognizes a hydrophobic cleft on the side of the alpha(2) I-domain, suggesting an alternative targeting site for drug development.

Allosteric Regulation↗
PubMed2006

Activation of individual alphaIIbbeta3 integrin molecules by disruption of transmembrane domain interactions in the absence of clustering.

We used laser tweezers-based force spectroscopy to measure the binding strength between fibrinogen molecules covalently bound to latex beads and either wild-type alphaIIbbeta3 molecules or alphaIIbbeta3 molecules containing the transmembrane domain mutations beta3 G708N or alphaIIb G972N expressed on Chinese hamster ovary cells. As we demonstrated previously for alphaIIbbeta3 on agonist-stimulated platelets and for purified alphaIIbbeta3 molecules incubated with Mn(2+), two regimes of rupture forces were present when wild-type alphaIIbbeta3 was activated by the monoclonal antibody PT25-2: rupture forces of 20-60 pN with an exponentially decreasing probability of detection and rupture forces in the range of 60-150 pN with a maximum at approximately 70-80 pN. Both rupture force regimes were specific for fibrinogen binding to the activated conformation of alphaIIbbeta3 because they were inhibited by alphaIIbbeta3-specific antagonists. Identical rupture force regimes were present constitutively when cells expressing the alphaIIb and beta3 transmembrane domain mutants were studied, confirming that these mutations induced an active alphaIIbbeta3 conformation. Moreover, there were no significant differences in the yield strength of the low-to-moderate and strong force regimes when alphaIIbbeta3 was activated by PT25-2 or the transmembrane domain mutations, implying that there was no fundamental difference in the way these forms of activated alphaIIbbeta3 interacted with fibrinogen. Thus, the two-step pathway of the interaction of alphaIIbbeta3 with fibrinogen we have identified appears to be a fundamental property of the high-affinity state of alphaIIbbeta3 and is identical regardless of whether this affinity state is achieved by intracellular, extracellular, or membrane-associated events.

Animals↗
PubMed2006

Controlling the shape and flexibility of arylamides: a combined ab initio, ab initio molecular dynamics, and classical molecular dynamics study.

Using quantum chemistry plus ab initio molecular dynamics and classical molecular dynamics methods, we address the relationship between molecular conformation and the biomedical function of arylamide polymers. Specifically, we have developed new torsional parameters for a class of these polymers and applied them in a study of the interaction between a representative arylamide and one of its biomedical targets, the anticoagulant drug heparin. Our main finding is that the torsional barrier of a C(aromatic)-C(carbonyl) bond increases significantly upon addition of an o-OCH2CH2NH3+ substituent on the benzene ring. Our molecular dynamics studies that are based on the original general AMBER force field (GAFF) and GAFF modified to include our newly developed torsional parameters show that the binding mechanism between the arylamide and heparin is very sensitive to the choice of torsional potentials. Ab initio molecular dynamics simulation of the arylamide independently confirms the degree of flexibility we obtain by classical molecular dynamics when newly developed torsional potentials are used.

Acetanilides↗
PubMed2006

Characterization of nonbiological antimicrobial polymers in aqueous solution and at water-lipid interfaces from all-atom molecular dynamics.

We have applied molecular dynamics to investigate the structural properties and activity of recently synthesized amphiphilic polymethacrylate derivatives, designed to mimic the antimicrobial activity of natural peptides. The composition, molecular weight, and hydrophobicity (ratio of hydrophobic and cationic units) of these short copolymers can be modulated to achieve structural diversity, which is crucial in controlling the antimicrobial activity. We have carried out all-atom molecular dynamics to systematically investigate the conformations adopted by these copolymers in water and at the water-lipid interface as a function of sequence and the chemical nature of the monomers. For two sequences, we observe partial insertion into the bilayer. Formation of strong interactions between the lipid headgroups and the amine groups of the polymers assists in the initial association with the lipids. However, the primary driving force for the observed partial insertion appears to be the hydrophobic effect. Our results indicate sensitive dependence of the overall shape on the sequence, suggesting that experimentally observed changes in activity can be correlated with particular sequences, providing an avenue for rational design.

Anti-Infective Agents↗
PubMed2006

Computational design of heterochiral peptides against a helical target.

Polypeptides incorporating D-amino acids occasionally occur in nature and are an important class of pharmaceutical molecules. With the use of heterochiral Monte Carlo (HCMC), a method inspired by the de novo design of proteins, we develop peptide scaffolds for interacting with a molecular target, a left-handed alpha-helix. The HCMC approach concurrently seeks to optimize a peptide sequence, its internal conformation, and its docked conformation with a target surface. Several major classes of interactions are observed: (1) homochiral interactions between two alphaL helices, (2) heterochiral interactions between an alphaL and an alphaR helix, and (3) heterochiral interactions between the alphaL target and novel nonhelical structures. We explore the application of HCMC to simulating the preferential enantioselectivity of heterochiral complexes. Implications for biomimetic design in molecular recognition are discussed.

Alanine↗
PubMed2005

Response of a designed metalloprotein to changes in metal ion coordination, exogenous ligands, and active site volume determined by X-ray crystallography.

The de novo protein DF1 is a minimal model for diiron and dimanganese metalloproteins, such as soluble methane monooxygenase. DF1 is a homodimeric four-helix bundle whose dinuclear center is formed by two bridging Glu side chains, two chelating Glu side chains, and two monodentate His ligands. Here, we report the di-Mn(II) and di-Co(II) derivatives of variants of this protein. Together with previously solved structures, 23 crystallographically independent four-helix bundle structures of DF1 variants have been determined, which differ in the bound metal ions and size of the active site cavity. For the di-Mn(II) derivatives, as the size of the cavity increases, the number and polarity of exogenous ligands increases. This collection of structures was analyzed to determine the relationship between protein conformation and the geometry of the active site. The primary mode of backbone movement involves a coordinated tilting and sliding of the first helix in the helix-loop-helix motif. Sliding depends on crystal-packing forces, the steric bulk of a critical residue that determines the dimensions of the active site access cavity, and the intermetal distance. Additionally, a torsional motion of the bridging carboxylates modulates the intermetal distance. This analysis provides a critical evaluation of how conformation, flexibility, and active site accessibility affect the geometry and ligand-binding properties of a metal center. The geometric parameters defining the DF structures were compared to natural diiron proteins; DF proteins have a restricted active site cavity, which may have implications for substrate recognition and chemical stability.

Binding Sites↗
PubMed2005

Spectroscopic and computational studies of the de novo designed protein DF2t: correlation to the biferrous active site of ribonucleotide reductase and factors that affect O2 reactivity.

DF2t, a de novo designed protein that mimics the active-site structure of many non-heme biferrous enzymes, has been studied using a combination of circular dichroism (CD), magnetic circular dichroism (MCD), and variable-temperature variable-field (VTVH) MCD. The active site of DF2t is found to have one five-coordinate iron and one four-coordinate iron, which are weakly antiferromagnetically coupled through a mu-1,3 carboxylate bridge. These results bear a strong resemblance to the spectra of Escherichia coli ribonucleotide reductase (R2), and density functional theory calculations were conducted on the W48F/D84E R2 mutant in order to determine the energetics of formation of a monodentate end-on-bound O2 to one iron in the binuclear site. The mu-1,3 carboxylate bridges found in O2-activating enzymes lack efficient superexchange pathways for the second electron transfer (i.e., the OH/oxo bridge in hemerythrin), and simulations of the binding of O2 in a monodentate end-on manner revealed that the bridging carboxylate ligands do not appear capable of transferring an electron to O2 from the remote Fe. Comparison of the results from previous studies of the mu-1,2 biferric-peroxo structure, which bridges both irons, finds that the end-on superoxide mixed-valent species is considerably higher in energy than the bridging peroxo-diferric species. Thus, one of the differences between O2-activating and O2-binding proteins appears to be the ability of O2 to bridge both Fe centers to generate a peroxo intermediate capable of further reactivity.

Binding Sites↗
PubMed2005

Alpha-alpha linking motifs and interhelical orientations.

While the geometry and sequence preferences of turns that link two beta-strands have been exhaustively explored, the corresponding preferences for sequences that link helical structures have been less well studied. Here we examine the interhelical geometry of two connected helices as a function of their link's length. The interhelical geometry of a helical pair appears to be significantly influenced by the number of linking residues. Furthermore, for relatively short link lengths, a very limited number of predominant conformations are observed, which can be categorized by their phi/psi angles. No more than two predominant linking backbone conformations are observed for a given link length, and some linking backbone conformations correlate strongly with distinctive interhelical geometric parameters. In this study, sequence and hydrogen-bonding patterns were defined for predominant interhelical link motifs. These results should assist in both protein structure prediction and de novo protein design.

Amino Acid Motifs↗
PubMed2005

Artificial di-iron proteins: solution characterization of four helix bundles containing two distinct types of inter-helical loops.

Peptide-based models have an enormous impact for the development of metalloprotein models, as they seem appropriate candidates to mimic both the structural characteristics and reactivity of the natural systems. Through the de novo design of four-helix bundles, we developed the DF (Due Ferri) family of artificial proteins, as models of di-iron and di-manganese metalloproteins. The goal of our research is to elucidate how the electrostatic environment, polarity and solvent accessibility of the metal-binding site, influence the functional properties of di-iron proteins. The first two subsets of the DF protein family, DF1 and DF2, consist of two non-covalently associated helix-loop-helix motifs, which bind the di-metal cofactor near the center of the structure. The DF2 subset was designed to improve the properties of DF1: DF2 and DF2t have several changes in their sequences to improve solubility and metal ion access, as well as a change in the loop connecting the two helices. In order to evaluate how these changes affect the overall structure of the model proteins, we solved the NMR structures of the di-Zn(II) complexes of DF2 and DF2t, and compared these structures with those recently obtained from X-ray crystallography. Further, we examined the thermodynamic consequences associated with the mutations, by measuring the stability of DF2t in the presence of different metal ions, and comparing the results with the data already obtained for DF2. Taken together, analysis of all the data showed the importance of the turn conformation in the design and stability of four-helix bundle.

Amino Acid Sequence↗
PubMed2005

T-jump infrared study of the folding mechanism of coiled-coil GCN4-p1.

Partially folded intermediates have been frequently observed in equilibrium and kinetic protein folding studies. However, folding intermediates that exist at the native side of the rate-limiting step are rather difficult to study because they often evade detection by conventional folding kinetic methods. Here, we demonstrated that a laser-induced temperature-jump method can potentially be used to identify the existence of such post-transition or hidden intermediates. Specifically, we studied two cross-linked variants of GCN4-p1 coiled-coil. The GCN4 leucine zipper has been studied extensively and most of these studies have regarded it as a two-state folder. Our static circular dichroism and infrared data also indicate that the thermal unfolding of these two monomeric coiled-coils can be adequately described by an apparent two-state model. However, their temperature-jump-induced relaxation kinetics exhibit non-monoexponential behavior, dependent upon sequence and temperature. Taken together, our results support a folding mechanism wherein at least one folding intermediate populates behind the main rate-limiting step.

Basic-Leucine Zipper Transcription Factors↗
PubMed2005

Structural organization and interactions of transmembrane domains in tetraspanin proteins.

BACKGROUND: Proteins of the tetraspanin family contain four transmembrane domains (TM1-4) linked by two extracellular loops and a short intracellular loop, and have short intracellular N- and C-termini. While structure and function analysis of the larger extracellular loop has been performed, the organization and role of transmembrane domains have not been systematically assessed. RESULTS: Among 28 human tetraspanin proteins, the TM1-3 sequences display a distinct heptad repeat motif (abcdefg)n. In TM1, position a is occupied by structurally conserved bulky residues and position d contains highly conserved Asn and Gly residues. In TM2, position a is occupied by conserved small residues (Gly/Ala/Thr), and position d has a conserved Gly and two bulky aliphatic residues. In TM3, three a positions of the heptad repeat are filled by two leucines and a glutamate/glutamine residue, and two d positions are occupied by either Phe/Tyr or Val/Ile/Leu residues. No heptad motif is apparent in TM4 sequences. Mutations of conserved glycines in human CD9 (Gly25 and Gly32 in TM1; Gly67 and Gly74 in TM2) caused aggregation of mutant proteins inside the cell. Modeling of the TM1-TM2 interface in CD9, using a novel algorithm, predicts tight packing of conserved bulky residues against conserved Gly residues along the two helices. The homodimeric interface of CD9 was mapped, by disulfide cross-linking of single-cysteine mutants, to the vicinity of residues Leu14 and Phe17 in TM1 (positions g and c) and Gly77, Gly80 and Ala81 in TM2 (positions d, g and a, respectively). Mutations of a and d residues in both TM1 and TM2 (Gly25, Gly32, Gly67 and Gly74), involved in intramolecular TM1-TM2 interaction, also strongly diminished intermolecular interaction, as assessed by cross-linking of Cys80. CONCLUSION: Our results suggest that tetraspanin intra- and intermolecular interactions are mediated by conserved residues in adjacent, but distinct regions of TM1 and TM2. A key structural element that defines TM1-TM2 interaction in tetraspanins is the specific packing of bulky residues against small residues.

Algorithms↗
PubMed2005

Automated use of mutagenesis data in structure prediction.

In the absence of experimental structural determination, numerous methods are available to indirectly predict or probe the structure of a target molecule. Genetic modification of a protein sequence is a powerful tool for identifying key residues involved in binding reactions or protein stability. Mutagenesis data is usually incorporated into the modeling process either through manual inspection of model compatibility with empirical data, or through the generation of geometric constraints linking sensitive residues to a binding interface. We present an approach derived from statistical studies of lattice models for introducing mutation information directly into the fitness score. The approach takes into account the phenotype of mutation (neutral or disruptive) and calculates the energy for a given structure over an ensemble of sequences. The structure prediction procedure searches for the optimal conformation where neutral sequences either have no impact or improve stability and disruptive sequences reduce stability relative to wild type. We examine three types of sequence ensembles: information from saturation mutagenesis, scanning mutagenesis, and homologous proteins. Incorporating multiple sequences into a statistical ensemble serves to energetically separate the native state and misfolded structures. As a result, the prediction of structure with a poor force field is sufficiently enhanced by mutational information to improve accuracy. Furthermore, by separating misfolded conformations from the target score, the ensemble energy serves to speed up conformational search algorithms such as Monte Carlo-based methods.

Algorithms↗
PubMed2005

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins.

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.

Amino Acids↗
PubMed2005

Synergistic interactions between aqueous and membrane domains of a designed protein determine its fold and stability.

Membrane-spanning proteins contain both aqueous and membrane-spanning regions, both of which contribute to folding and stability. To explore the interplay between these two domains we have designed and studied the assembly of coiled-coil peptides that span from the membrane into the aqueous phase. The membrane-spanning segment is based on MS1, a transmembrane coiled coil that contains a single Asn at a buried a position of a central heptad in its sequence. This Asn has been shown to drive assembly of the monomeric peptide in a membrane environment to a mixture of dimers and trimers. The coiled coil has now been extended into the aqueous phase by addition of water-soluble helical extensions. Although too short to fold in isolation, these helical extensions were expected to interact synergistically with the transmembrane domain and modulate its stability as well as its conformational specificity for forming dimers versus trimers. One design contains Asn at a position of the aqueous helical extension, which was expected to specify a dimeric state; a second peptide, which contains Val at this position, was expected to form trimers. The thermodynamics of assembly of the hybrid peptides were studied in micelles by sedimentation equilibrium ultracentrifugation. The aqueous helical extensions indeed conferred additional stability and conformational specificity to MS1 in the expected manner. These studies highlight the delicate interplay between membrane-spanning and water-soluble regions of proteins, and demonstrate how these different environments define the thermodynamics of a given specific interaction. In this case, an Asn in the transmembrane domain provided a strong driving force for folding but failed to specify a unique oligomerization state, while an Asn in the water-soluble domain was able to define specificity for a specific aggregation state as well as modulate stability.

Amino Acid Sequence↗