The 1.59Å resolution structure of the minor pseudopilin EpsH of Vibrio cholerae reveals a long flexible loop.

The type II secretion complex exports folded proteins from the periplasm to the extracellular milieu. It is used by the pathogenic bacterium Vibrio cholerae to export several proteins, including its major virulence factor, cholera toxin. The pseudopilus is an essential component of the type II secretion system and likely acts as a piston to push the folded proteins across the outer membrane through the secretin pore. The pseudopilus is composed of the major pseudopilin, EpsG, and four minor pseudopilins, EpsH, EpsI, EpsJ and EpsK. We determined the x-ray crystal structure of the head domain of EpsH at 1.59Å resolution using molecular replacement with the previously reported EpsH structure, 2qv8, as the template. Three additional N-terminal amino acids present in our construct prevent an artifactual conformation of residues 160-166, present in one of the two monomers of the 2qv8 structure. Additional crystal contacts stabilize a long flexible loop comprised of residues 104-135 that is more disordered in the 2qv8 structure but is partially observed in our structure in very different positions for the two EpsH monomers in the asymmetric unit. In one of the conformations the loop is highly extended. Modeling suggests the highly charged loop is capable of contacting EpsG and possibly secreted protein substrates, suggesting a role in specificity of pseudopilus assembly or secretion function.

Microsecond unfolding kinetics of sheep prion protein reveals an intermediate that correlates with susceptibility to classical scrapie.

The microsecond folding and unfolding kinetics of ovine prion proteins (ovPrP) were measured under various solution conditions. A fragment comprising residues 94-233 of the full-length ovPrP was studied for four variants with differing susceptibilities to classical scrapie in sheep. The observed biexponential unfolding kinetics of ovPrP provides evidence for an intermediate species. However, in contrast to previous results for human PrP, there is no evidence for an intermediate under refolding conditions. Global analysis of the kinetic data, based on a sequential three-state mechanism, quantitatively accounts for all folding and unfolding data as a function of denaturant concentration. The simulations predict that an intermediate accumulates under both folding and unfolding conditions, but is observable only in unfolding experiments because the intermediate is optically indistinguishable from the native state. The relative population of intermediates in two ovPrP variants, both transiently and under destabilizing equilibrium conditions, correlates with their propensities for classical scrapie. The variant susceptible to classical scrapie has a larger population of the intermediate state than the resistant variant. Thus, the susceptible variant should be favored to undergo the PrP(C) to PrP(Sc) conversion and oligomerization.

A general polymer model of unfolded proteins under folding conditions.

There is increasing evidence that a polypeptide chain in solvent conditions that favor folding may have transient structure and is significantly more compact than a fully denatured chain. However, there is no sequence-dependent model to capture such interactions. In this work, we present a simple and computationally inexpensive model based on a wormlike chain with excluded volume. The probability distribution of millions of such chains is reweighted to bias compact conformations in which residues of similar hydrophobicity are located near each other. This model, which has one adjustable parameter, is fit to measured values of intramolecular contact formation, which has been shown to be extremely sensitive to various models of intrachain distances. We show that under various denaturant concentrations, there is good convergence of the model for several different sequences with a wide range of dynamics. We also show that this model quantitatively predicts paramagnetic resonance enhancement (PRE) measurements with no adjustable parameters. Therefore a simple, probabilistic model that accounts for sequence-specific interactions may give a more realistic starting point for predictions of protein folding.

Homology-based modeling of the Erwinia amylovora type III secretion chaperone DspF used to identify amino acids required for virulence and interaction with the effector DspE.

The structure of DspF, a type III secretion system (T3SS) chaperone required for virulence of the fruit tree pathogen Erwinia amylovora, was modeled based on predicted structural homology to characterized T3SS chaperones. This model guided the selection of 11 amino acid residues that were individually mutated to alanine via site-directed mutagenesis. Each mutant was assessed for its effect on virulence complementation, dimerization and interaction with the N-terminal chaperone-binding site of DspE. Four amino acid residues were identified that did not complement the virulence defect of a dspF knockout mutant, and three of these residues were required for interaction with the N-terminus of DspE. This study supports the significance of the predicted beta-sheet helix-binding groove in DspF chaperone function.

Unfolded-state dynamics and structure of protein L characterized by simulation and experiment.

While several experimental techniques now exist for characterizing protein unfolded states, all-atom simulation of unfolded states has been challenging due to the long time scales and conformational sampling required. We address this problem by using a combination of accelerated calculations on graphics processor units and distributed computing to simulate tens of thousands of molecular dynamics trajectories each up to approximately 10 mus (for a total aggregate simulation time of 127 ms). We used this approach in conjunction with Trp-Cys contact quenching experiments to characterize the unfolded structure and dynamics of protein L. We employed a polymer theory method to make quantitative comparisons between high-temperature simulated and chemically denatured experimental ensembles and find that reaction-limited quenching rates calculated from simulation agree remarkably well with experiment. In both experiment and simulation, we find that unfolded-state intramolecular diffusion rates are very slow compared to highly denatured chains and that a single-residue mutation can significantly alter unfolded-state dynamics and structure. This work suggests a view of the unfolded state in which surprisingly low diffusion rates could limit folding and opens the door for all-atom molecular simulation to be a useful predictive tool for characterizing protein unfolded states along with experiments that directly measure intramolecular diffusion.

Expression, purification, crystallization and preliminary X-ray studies of Vibrio cholerae pseudopilin EpsH.

EpsH is a minor pseudopilin protein of the Vibrio cholerae type II secretion system. A truncated form of EpsH with a C-terminal noncleavable His tag was constructed and expressed in Escherichia coli, purified and crystallized by sitting-drop vapor diffusion. A complete data set was collected to 1.71 A resolution. The crystals belonged to space group P2(1)2(1)2(1), with unit-cell parameters a = 53.39, b = 71.11, c = 84.64 A. There were two protein molecules in the asymmetric unit, which gave a Matthews coefficient V(M) of 2.1 A(3) Da(-1), corresponding to 41.5% solvent content.

Ruggedness in the folding landscape of protein L.

By exploring the folding pathways of the B1 domain of protein L with a series of equilibrium and rapid kinetic experiments, we have found its unfolded state to be more complex than suggested by two-state folding models. Using an ultrarapid mixer to initiate protein folding within approximately 2-4 microseconds, we observe folding kinetics by intrinsic tryptophan fluorescence and fluorescence resonance energy transfer. We detect at least two processes faster than 100 mus that would be hidden within the burst phase of a stopped-flow instrument measuring tryptophan fluorescence. Previously reported measurements of slow intramolecular diffusion are commensurate with the slower of the two observed fast phases. These results suggest that a multidimensional energy landscape is necessary to describe the folding of protein L, and that the dynamics of the unfolded state is dominated by multiple small energy barriers.

Dynamic similarity of the unfolded states of proteins L and G.

The formation of specific intramolecular contacts has been studied under a range of denaturing conditions in single domains of the immunoglobulin-binding proteins L and G. Although they share no significant sequence similarity and have dissimilar folding pathways, the two domains have a similar native fold. Our measurements show that the rates of forming corresponding contacts in the unfolded states of both proteins are remarkably similar and even exhibit similar dependence on denaturant concentration. The unfolded proteins were modeled using Szabo, Schulten, and Schulten (SSS) theory as wormlike chains with excluded volume; when combined with our experimental data, the SSS analysis suggests that the unfolded state becomes uniformly more compact and less diffusive (i.e., rearranges more slowly) with decreasing denaturant concentrations.

Ab initio modeling of the herpesvirus VP26 core domain assessed by CryoEM density.

Efforts in structural biology have targeted the systematic determination of all protein structures through experimental determination or modeling. In recent years, 3-D electron cryomicroscopy (cryoEM) has assumed an increasingly important role in determining the structures of these large macromolecular assemblies to intermediate resolutions (6-10 A). While these structures provide a snapshot of the assembly and its components in well-defined functional states, the resolution limits the ability to build accurate structural models. In contrast, sequence-based modeling techniques are capable of producing relatively robust structural models for isolated proteins or domains. In this work, we developed and applied a hybrid modeling approach, utilizing cryoEM density and ab initio modeling to produce a structural model for the core domain of a herpesvirus structural protein, VP26. Specifically, this method, first tested on simulated data, utilizes the cryoEM density map as a geometrical constraint in identifying the most native-like models from a gallery of models generated by ab initio modeling. The resulting model for the core domain of VP26, based on the 8.5-A resolution herpes simplex virus type 1 (HSV-1) capsid cryoEM structure and mutational data, exhibited a novel fold. Additionally, the core domain of VP26 appeared to have a complementary interface to the known upper-domain structure of VP5, its cognate binding partner. While this new model provides for a better understanding of the assembly and interactions of VP26 in HSV-1, the approach itself may have broader applications in modeling the components of large macromolecular assemblies.

An approximately 140-kb deletion associated with feline spinal muscular atrophy implies an essential LIX1 function for motor neuron survival.

The leading genetic cause of infant mortality is spinal muscular atrophy (SMA), a clinically and genetically heterogeneous group of disorders. Previously we described a domestic cat model of autosomal recessive, juvenile-onset SMA similar to human SMA type III. Here we report results of a whole-genome scan for linkage in the feline SMA pedigree using recently developed species-specific and comparative mapping resources. We identified a novel SMA gene candidate, LIX1, in an approximately140-kb deletion on feline chromosome A1q in a region of conserved synteny to human chromosome 5q15. Though LIX1 function is unknown, the predicted secondary structure is compatible with a role in RNA metabolism. LIX1 expression is largely restricted to the central nervous system, primarily in spinal motor neurons, thus offering explanation of the tissue restriction of pathology in feline SMA. An exon sequence screen of 25 human SMA cases, not otherwise explicable by mutations at the SMN1 locus, failed to identify comparable LIX1 mutations. Nonetheless, a LIX1-associated etiology in feline SMA implicates a previously undetected mechanism of motor neuron maintenance and mandates consideration of LIX1 as a candidate gene in human SMA when SMN1 mutations are not found.

LRRK2 in Parkinson's disease: protein domains and functional insights.

Parkinson's disease (PD) is the most common motor neurodegenerative disease. Mutations in the gene encoding leucine-rich repeat kinase 2 (LRRK2) have been linked recently with autosomal-dominant parkinsonism that is clinically indistinguishable from typical, idiopathic, late-onset PD. Thus, the protein LRRK2 has emerged as a promising therapeutic target for treatment of PD. LRRK2 is extraordinarily large and complex, with multiple enzymatic and protein-interaction domains, each of which is targeted by pathogenic mutations in familial PD. This review places the PD-associated mutations of LRRK2 in a structural and functional framework, with the ultimate aim of deciphering the molecular basis of LRRK2-associated pathogenesis. This, in turn, should advance our understanding and treatment of familial and idiopathic PD.

The AidB component of the Escherichia coli adaptive response to alkylating agents is a flavin-containing, DNA-binding protein.

Upon exposure to alkylating agents, Escherichia coli increases expression of aidB along with three genes (ada, alkA, and alkB) that encode DNA repair proteins. In order to begin to identify the role of AidB in the cell, the protein was purified to homogeneity, shown to possess stoichiometric amounts of flavin adenine dinucleotide (FAD), and confirmed to have low levels of isovaleryl-coenzyme A (CoA) dehydrogenase activity. A homology model of an AidB homodimer was constructed based on the structure of a four-domain acyl-CoA oxidase. The predicted structure revealed a positively charged groove connecting the two active sites and a second canyon of positive charges in the C-terminal domain, both of which could potentially bind DNA. Three approaches were used to confirm that AidB binds to double-stranded DNA. On the basis of its ability to bind DNA and its possession of a redox-active flavin, AidB is predicted to catalyze the direct repair of alkylated DNA.

Characterization of peptide folding nuclei by hydrogen/deuterium exchange-mass spectrometry.

Covalently linked pairs of well-chosen peptides can be good model systems for protein folding studies because they can adopt stable secondary, side-chain, and tertiary structure under certain conditions. We demonstrate a method for characterizing the structure in such peptide pairs by hydrogen/deuterium exchange of individual amide groups analyzed by collision-induced dissociation tandem mass spectrometry, in concert with circular dichroism spectroscopy. We apply the method to two peptides (and their three possible pairs) from bovine pancreatic trypsin inhibitor to address specific hypotheses regarding the stabilization of local secondary structure by long-range interactions.

Contributions of amino acid side chains to the kinetics and thermodynamics of the bivalent binding of protein L to Ig kappa light chain.

Protein L is a bacterial surface protein with 4-5 immunoglobulin (Ig)-binding domains (B1-B5), each of which appears to have two binding sites for Ig, corresponding to the two edges of its beta-sheet. To verify these sites biochemically and to probe their relative contributions to the protein L-Ig kappa light chain (kappa) interaction, we compared the binding of PLW (the Y47W mutant of the B1 domain) to that of mutants designed to disrupt binding to sites 1 and 2, using gel filtration, BIAcore surface plasmon resonance, fluorescence titration, and solid-phase radioimmunoassays. Gel filtration experiments show that PLW binds kappa both in 1:1 complexes and multivalently, consistent with two binding sites. Covalent dimers of the A20C and V51C mutants of PLW were prepared to eliminate site 1 and site 2 binding, respectively; both the A20C and V51C dimers bind kappa in 1:1 complexes and multivalently, indicating that neither site 1 nor site 2 is solely responsible for kappa binding. The A20R mutant was designed computationally to eliminate site 1 binding while preserving site 2 binding; consistent with this design, the A20R mutant binds kappa in 1:1 complexes but not multivalently. To probe the contributions of amino acid side chains to binding, we prepared 75 point mutants spanning nearly every residue of PLW; BIAcore studies of these mutants revealed two binding-energy "hot spots" consistent with sites 1 and 2. These data indicate that PLW binds kappa at both sites with similar affinities (high nanomolar), with the strongest contributions to the binding energy from Tyr34 (site 2) and Tyr36 (site 1). Compared to other protein-protein complexes, the binding is insensitive to amino acid substitutions at these sites, consistent with the large number of main chain interactions relative to side chain interactions. The strong binding of protein L to Ig kappa light chains of various species may result from the ambidextrous binding of the B1-B5 domains and the unimportance of specific side chain interactions.

Efficient minimization of angle-dependent potentials for polypeptides in internal coordinates.

Angular potentials play an important role in the refinement of protein structures through angle-dependent restraints (e.g., those determined by cross-correlated relaxations, residual dipolar couplings, and hydrogen bonds). Analytic derivatives of such angular potentials with respect to the dihedral angles of proteins would be useful for optimizing such restraints and other types of angular potentials (i.e., such as we are now introducing into protein structure prediction) but have not been described. In this article, analytic derivatives are calculated for four types of angular potentials and integrated with the efficient recursive derivative calculation methods of Go and coworkers. The formulas are implemented in publicly available software and illustrated by refining a low-resolution protein structure with idealized vector-angle, dipolar-coupling, and hydrogen-bond restraints. The method is now being used routinely to optimize hydrogen-bonding potentials in ROSETTA.

Rosetta predictions in CASP5: successes, failures, and prospects for complete automation.

We describe predictions of the structures of CASP5 targets using Rosetta. The Rosetta fragment insertion protocol was used to generate models for entire target domains without detectable sequence similarity to a protein of known structure and to build long loop insertions (and N-and C-terminal extensions) in cases where a structural template was available. Encouraging results were obtained both for the de novo predictions and for the long loop insertions; we describe here the successes as well as the failures in the context of current efforts to improve the Rosetta method. In particular, de novo predictions failed for large proteins that were incorrectly parsed into domains and for topologically complex (high contact order) proteins with swapping of segments between domains. However, for the remaining targets, at least one of the five submitted models had a long fragment with significant similarity to the native structure. A fully automated version of the CASP5 protocol produced results that were comparable to the human-assisted predictions for most of the targets, suggesting that automated genomic-scale, de novo protein structure prediction may soon be worthwhile. For the three targets where the human-assisted predictions were significantly closer to the native structure, we identify the steps that remain to be automated.

Proline isomerization in bovine pancreatic ribonuclease A. 2. Folding conditions.

The kinetics of cis-trans isomerization of individual X-Pro peptide groups is used to study the backbone dynamics of bovine pancreatic ribonuclease A (RNase A). We previously developed and validated a fluorescence method for monitoring the cis-trans isomerization of the Tyr92-Pro93 and Asn113-Pro114 peptide groups of RNase A under unfolding conditions [Juminaga, D., Wedemeyer, W. J., and Scheraga, H. A. (1998) Biochemistry 37, 11614-11620]. The essence of this method is to introduce a fluorescent residue (Tyr or Trp) in a position adjacent to the isomerizing proline (if one is not already present) and to eliminate the fluorescence of other such residues adjacent to prolines by mutating them to phenylalanine. Here, we extend this method to observe the cis-trans isomerization of these peptide groups under folding conditions using two site-directed mutants (Y92F and Y115F) of RNase A. Both isomerizations decelerate with increasing concentrations of GdnHCl, with nearly identical m values (1.11 and 1.19 M(-1), respectively) and extrapolated zero-GdnHCl time constants (42 and 32 s, respectively); by contrast, under unfolding conditions, the cis-trans isomerizations of both Pro93 and Pro114 are independent of GdnHCl concentration. Remarkably, the isomerization rates under folding conditions at GdnHCl concentrations above 1 M are significantly slower than those measured under unfolding conditions. The temperature dependence of the Pro114 isomerization under folding conditions is also unusual; whereas Pro93 exhibits an activation energy typical of proline isomerization (19.4 kcal/mol), Pro114 exhibits a sharply reduced activation energy of 5.7 kcal/mol. A structurally plausible model accounts for these results and, in particular, shows that folding conditions strongly accelerate the cis-trans isomerization of both peptide groups to their native cis conformation, suggesting the presence of flickering local structure in their beta-hairpins.

Proline cis-trans isomerization and protein folding.

Proline cis-trans isomerization plays a key role in the rate-determining steps of protein folding. The energetic origin of this isomerization process is summarized, and the folding and unfolding of disulfide-intact bovine pancreatic ribonuclease A is used as an example to illustrate the kinetics and structural features of conformational changes from the heterogeneous unfolded state (consisting of cis and trans isomers of X-Pro peptide groups) to the native structure in which only one set of proline isomers is present.

Protein structure prediction in 2002.

Central issues concerning protein structure prediction have been highlighted by the recently published summary of the fourth community-wide protein structure prediction experiment (CASP4). Although sequence/structure alignment remains the bottleneck in comparative modeling, there has been substantial progress in fully automated remote homolog detection and in de novo structure prediction. Significant further progress will probably require improvements in high-resolution modeling.

Exact solutions for chemical bond orientations from residual dipolar couplings.

New methods for determining chemical structures from residual dipolar couplings are presented. The fundamental dipolar coupling equation is converted to an elliptical equation in the principal alignment frame. This elliptical equation is then combined with other angular or dipolar coupling constraints to form simple polynomial equations that define discrete solutions for the unit vector(s). The methods are illustrated with residual dipolar coupling data on ubiquitin taken in a single anisotropic medium. The protein backbone is divided into its rigid groups (namely, its peptide planes and Calpha frames), which may be solved for independently. A simple procedure for recombining these independent solutions results in backbone dihedral angles phi and psi that resemble those of the known native structure. Subsequent refinement of these phi-psi angles by the ROSETTA program produces a structure of ubiquitin that agrees with the known native structure to 1.1 A Calpha rmsd.