B1 Domain Of Protein Research Articles

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Sidechain chi1 distribution in B3 domain of protein G from extensive sets of residual dipolar couplings

Site-specific protonation kinetics of acidic side chains in proteins determined by pH-dependent carboxyl (13)C NMR relaxation.

Proton-transfer dynamics plays a critical role in many biochemical processes, such as proton pumping across membranes and enzyme catalysis. The large majority of enzymes utilize acid-base catalysis and proton-transfer mechanisms, where the rates of proton transfer can be rate limiting for the overall reaction. However, measurement of proton-exchange kinetics for individual side-chain carboxyl groups in proteins has been achieved in only a handful of cases, which typically have involved comparative analysis of mutant proteins in the context of reaction network modeling. Here we describe an approach to determine site-specific protonation and deprotonation rate constants (kon and koff, respectively) of carboxyl side chains, based on (13)C NMR relaxation measurements as a function of pH. We validated the method using an extensively studied model system, the B1 domain of protein G, for which we measured rate constants koff in the range (0.1-3) × 10(6) s(-1) and kon in the range (0.6-300) × 10(9) M(-1) s(-1), which correspond to acid-base equilibrium dissociation constants (Ka) in excellent agreement with previous results determined by chemical shift titrations. Our results further reveal a linear free-energy relationship between log kon and pKa, which provides information on the free-energy landscape of the protonation reaction, showing that the variability among residues in these parameters arises primarily from the extent of charge stabilization of the deprotonated state by the protein environment. We find that side-chain carboxyls with extreme values of koff or kon are involved in hydrogen bonding, thus providing a mechanistic explanation for the observed stabilization of the protonated or deprotonated state.

Quinary structure modulates protein stability in cells

Protein quinary interactions organize the cellular interior and its metabolism. Although the interactions stabilizing secondary, tertiary, and quaternary protein structure are well defined, details about the protein-matrix contacts that comprise quinary structure remain elusive. This gap exists because proteins function in the crowded cellular environment, but are traditionally studied in simple buffered solutions. We use NMR-detected H/D exchange to quantify quinary interactions between the B1 domain of protein G and the cytosol of Escherichia coli. We demonstrate that a surface mutation in this protein is 10-fold more destabilizing in cells than in buffer, a surprising result that firmly establishes the significance of quinary interactions. Remarkably, the energy involved in these interactions can be as large as the energies that stabilize specific protein complexes. These results will drive the critical task of implementing quinary structure into models for understanding the proteome.

Complex Pathways in Folding of Protein G Explored by Simulation and Experiment

The B1 domain of protein G has been a classic model system of folding for decades, the subject of numerous experimental and computational studies. Most of the experimental work has focused on whether the protein folds via an intermediate, but the evidence is mostly limited to relatively slow kinetic observations with a few structural probes. In this work we observe folding on the submillisecond timescale with microfluidic mixers using a variety of probes including tryptophan fluorescence, circular dichroism, and photochemical oxidation. We find that each probe yields different kinetics and compare these observations with a Markov State Model constructed from large-scale molecular dynamics simulations and find a complex network of states that yield different kinetics for different observables. We conclude that there are many folding pathways before the final folding step and that these paths do not have large free energy barriers.

Tunable, mixed-resolution modeling using library-based Monte Carlo and graphics processing units.

Building on our recently introduced library-based Monte Carlo (LBMC) approach, we describe a flexible protocol for mixed coarse-grained (CG)/all-atom (AA) simulation of proteins and ligands. In the present implementation of LBMC, protein side chain configurations are pre-calculated and stored in libraries, while bonded interactions along the backbone are treated explicitly. Because the AA side chain coordinates are maintained at minimal run-time cost, arbitrary sites and interaction terms can be turned on to create mixed-resolution models. For example, an AA region of interest such as a binding site can be coupled to a CG model for the rest of the protein. We have additionally developed a hybrid implementation of the generalized Born/surface area (GBSA) implicit solvent model suitable for mixed-resolution models, which in turn was ported to a graphics processing unit (GPU) for faster calculation. The new software was applied to study two systems: (i) the behavior of spin labels on the B1 domain of protein G (GB1) and (ii) docking of randomly initialized estradiol configurations to the ligand binding domain of the estrogen receptor (ERα). The performance of the GPU version of the code was also benchmarked in a number of additional systems.

Amide temperature coefficients in the protein G B1 domain

Temperature coefficients have been measured for backbone amide (1)H and (15)N nuclei in the B1 domain of protein G (GB1), using temperatures in the range 283-313K, and pH values from 2.0 to 9.0. Many nuclei display pH-dependent coefficients, which were fitted to one or two pK(a) values. (1)H coefficients showed the expected behaviour, in that hydrogen-bonded amides have less negative values, but for those amides involved in strong hydrogen bonds in regular secondary structure there is a negative correlation between strength of hydrogen bond and size of temperature coefficient. The best correlation to temperature coefficient is with secondary shift, indicative of a very approximately uniform thermal expansion. The largest pH-dependent changes in coefficient are for amides in loops adjacent to sidechain hydrogen bonds rather than the amides involved directly in hydrogen bonds, indicating that the biggest determinant of the temperature coefficient is temperature-dependent loss of structure, not hydrogen bonding. Amide (15)N coefficients have no clear relationship with structure.

A Bayesian Approach for Determining Protein Side-Chain Rotamer Conformations Using Unassigned NOE Data

A major bottleneck in protein structure determination via nuclear magnetic resonance (NMR) is the lengthy and laborious process of assigning resonances and nuclear Overhauser effect (NOE) cross peaks. Recent studies have shown that accurate backbone folds can be determined using sparse NMR data, such as residual dipolar couplings (RDCs) or backbone chemical shifts. This opens a question of whether we can also determine the accurate protein side-chain conformations using sparse or unassigned NMR data. We attack this question by using unassigned nuclear Overhauser effect spectroscopy (NOESY) data, which records the through-space dipolar interactions between protons nearby in three-dimensional (3D) space. We propose a Bayesian approach with a Markov random field (MRF) model to integrate the likelihood function derived from observed experimental data, with prior information (i.e., empirical molecular mechanics energies) about the protein structures. We unify the side-chain structure prediction problem with the side-chain structure determination problem using unassigned NMR data, and apply the deterministic dead-end elimination (DEE) and A* search algorithms to provably find the global optimum solution that maximizes the posterior probability. We employ a Hausdorff-based measure to derive the likelihood of a rotamer or a pairwise rotamer interaction from unassigned NOESY data. In addition, we apply a systematic and rigorous approach to estimate the experimental noise in NMR data, which also determines the weighting factor of the data term in the scoring function derived from the Bayesian framework. We tested our approach on real NMR data of three proteins: the FF Domain 2 of human transcription elongation factor CA150 (FF2), the B1 domain of Protein G (GB1), and human ubiquitin. The promising results indicate that our algorithm can be applied in high-resolution protein structure determination. Since our approach does not require any NOE assignment, it can accelerate the NMR structure determination process.

Simulation of Chaperonin Effect on Protein Folding: A Shift from Nucleation–Condensation to Framework Mechanism

The iterative annealing mechanism (IAM) of chaperonin-assisted protein folding is explored in a framework of a well-established coarse-grained protein modeling tool, which enables the study of protein dynamics in a time-scale well beyond classical all-atom molecular mechanics. The chaperonin mechanism of action is simulated for two paradigm systems of protein folding, B domain of protein A (BdpA) and B1 domain of protein G (GB1), and compared to chaperonin-free simulations presented here for BdpA and recently published for GB1. The prediction of the BdpA transition state ensemble (TSE) is in perfect agreement with experimental findings. It is shown that periodic distortion of the polypeptide chains by hydrophobic chaperonin interactions can promote rapid folding and leads to a decrease in folding temperature. It is also demonstrated how chaperonin action prevents kinetically trapped conformations and modulates the observed folding mechanisms from nucleation–condensation to a more framework-like.

Ab Initio Calculations of Deuterium Isotope Effects on Chemical Shifts of Salt-Bridged Lysines

Deuterium isotope effects measure the change in chemical shift on substitution of a proton by deuterium. They have been calculated by direct treatment of the H/D nuclear quantum effect using a multicomponent ab initio molecular orbital method based on a non-Born-Oppenheimer approximation. This method enables the determination of both the electronic and the protonic (deuteronic) wave functions simultaneously and can directly calculate the geometrical difference induced by H/D isotope effects. The calculations show that the one-bond deuterium isotope effects on (15)N nuclear shielding, (1)Δ(15)N(D), in ammonium and amines decrease as a counterion or water molecule moves closer to the nitrogen. (1)Δ(15)N(D) and (2)Δ(1)H(D) of the NH(3)(+) groups of lysine residues in the B1 domain of protein G have been calculated using truncated side chains and also determined experimentally by NMR. Comparisons show that the structures in solution are different from those in the crystal and that solvation plays an important role in weakening the hydrogen bonds.

Molecular Dynamics Simulations of a β‐Hairpin Fragment of Protein G by Means of Atom‐Bond Electronegativity Equalization Method Fused into Molecular Mechanics (ABEEMδπ/MM)

AbstractThere are some controversial opinions about the origin of folding β‐hairpin stability in aqueous solution. In this study, the structural and dynamic behavior of a 16‐residue β‐hairpin from B1 domain of protein G has been investigated at 280, 300, 350 and 450 K using molecular dynamics (MD) simulations by means of Atom‐Bond Electronegativity Equalization Method Fused into Molecular Mechanics i.e., ABEEMδπ/MM and the explicit ABEEM‐7P water solvent model. In addition, a 300 K simulation of one mutant having the aromatic residues substituted with alanines has been performed. The hydrophobic surface area, hydrophilic surface area and some structural properties have been used to measure the role of the hydrophobic interactions. It is found that the aromatic residues substituted with alanines have shown an evident destabilization of the structure and unfolding started after 1.5 ns. It is also found that the number of the main chain hydrogen bonds have different distributions through three different simulations. All above demonstrate that the hydrophobic interactions and the main chain hydrogen bonds play an important role in the stability of the folding structure of β‐hairpin in solution. Furthermore, through the structural analyses of the β‐hairpin structures from four temperature simulations and the comparison with other MD simulations of β‐hairpin peptides, the new ABEEMδπ force field can reproduce the structural data in good agreement with the experimental data.

The Impact of Hydrogen Bonding on Amide 1H Chemical Shift Anisotropy Studied by Cross-Correlated Relaxation and Liquid Crystal NMR Spectroscopy

Site-specific 1H chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide moieties in the B3 domain of protein G (GB3). Experimental input data include residual chemical shift anisotropy (RCSA), measured in six mutants that align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension, and cross-correlated relaxation rates between the 1HN CSA tensor and either the 1H−15N, the 1H−13C′, or the 1H−13Cα dipolar interactions. Analyses with the assumption that the 1HN CSA tensor is symmetric with respect to the peptide plane (three-parameter fit) or without this premise (five-parameter fit) yield very similar results, confirming the robustness of the experimental input data, and that, to a good approximation, one of the principal components orients orthogonal to the peptide plane. 1HN CSA tensors are found to deviate strongly from axial symmetry, with the most shielded tensor component roughly parallel to the N−H vector, and the least shielded component orthogonal to the peptide plane. DFT calculations on pairs of N-methyl acetamide and acetamide in H-bonded geometries taken from the GB3 X-ray structure correlate with experimental data and indicate that H-bonding effects dominate variations in the 1HN CSA. Using experimentally derived 1HN CSA tensors, the optimal relaxation interference effect needed for narrowest 1HN TROSY line widths is found at ∼1200 MHz.

Site-specific backbone amide (15)N chemical shift anisotropy tensors in a small protein from liquid crystal and cross-correlated relaxation measurements.

Site-specific (15)N chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide (15)N nuclei in the B3 domain of protein G (GB3) from residual chemical shift anisotropy (RCSA) measured in six different mutants that retain the native structure but align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension. This information is complemented by measurement of cross-correlated relaxation rates between the (15)N CSA tensor and either the (15)N-(1)H or (15)N-(13)C' dipolar interaction. In agreement with recent solid state NMR measurements, the (15)N CSA tensors exhibit only a moderate degree of variation from averaged values, but have larger magnitudes in alpha-helical (-173 +/- 7 ppm) than in beta-sheet (-162 +/- 6 ppm) residues, a finding also confirmed by quantum computations. The orientations of the least shielded tensor component cluster tightly around an in-peptide-plane vector that makes an angle of 19.6 +/- 2.5 degrees with the N-H bond, with the asymmetry of the (15)N CSA tensor being slightly smaller in alpha-helix (eta = 0.23 +/- 0.17) than in beta-sheet (eta = 0.31 +/- 0.11). The residue-specific (15)N CSA values are validated by improved agreement between computed and experimental (15)N R(1rho) relaxation rates measured for (15)N-{(2)H} sites in GB3, which are dominated by the CSA mechanism. Use of residue-specific (15)N CSA values also results in more uniform generalized order parameters, S(2), and predicts considerable residue-by-residue variations in the magnetic field strengths where TROSY line narrowing is most effective.

Recombinant S-Layer Proteins of Lactobacillus brevis Mediating Antibody Adhesion to Calf Intestine Alleviated Neonatal Diarrhea Syndrome

A chimeric gene encoding enhanced green fluorescent protein (EGFP) and a S-layer protein from Lactobacillus brevis KCTC3102, and/or two copies of the Fc-binding Z-domain, a synthetic analog of the B-domain of protein A, was constructed and expressed in Escherichia coli BL21(DE3). The S-layer fusion proteins produced in a 500-l fermentor were likely to be stable in the range of pH 5 to 8 and 0 degree to 40 degrees . Their adhesive property enabled an easy and rapid immobilization of enzymes or antibodies on solid materials such as plastics, glass, sol-gel films, and intestinal epithelial cells. Owing to their affinity towards intestinal cells and immunoglobulin G, the Slayer fusion proteins enabled the adhesion of antibodies to human epithelial cells. In addition, feeding a mixture of the S-layer fusion proteins and antibodies against neonatal calf diarrhea (coronavirus, rotavirus, Escherichia coli, and Salmonella typhimurium) to Hanwoo calves resulted in 100% prevention of neonatal calf diarrhea syndrome (p<0.01),whereas feeding antibodies only resulted in 56% prevention.

A test of implicit solvent models on the folding simulation of the GB1 peptide

To quantify the effects of implicit solvent model on the folding of protein, thermodynamics calculation was performed for the folding of C-terminal hairpin from B1 domain of protein G using AMBER ff96 force field and different implicit solvent models (GB(OBC) and GBn). In the multiple AMBER force fields, ff96 strongly favors the hairpin conformations whereas the other force fields (e.g., ff94 and ff99) favor the helical conformations. In recent studies ff96 was found to be better balanced for various secondary structures than other force fields. The two implicit solvent models mentioned above were shown to be the most appropriate for beta-hairpin folding [M. S. Shell et al., J. Phys. Chem. B 112, 6878 (2008)]. The free energy landscape as a function of several different reaction coordinates was studied for both solvent models. The comparison of the results from these implicit models shows that the combination of AMBER ff96 with the GB(OBC) model yields stable structures that are in better agreement with experiments and results calculated using the explicit solvent model.

Characterization of Salt Bridges to Lysines in the Protein G B1 Domain

NMR investigations have been carried out on the B1 domain of protein G. This protein has six lysine residues, of which three are consistently found to form surface-exposed salt bridges in crystal structures, while the other three are not. The Nzeta and Hzeta chemical shifts of all six lysines are similar and are not affected significantly by pH titration of the carboxylate groups in the protein, except for a relatively small titration of K39 Nzeta. Deuterium isotope effects on nitrogen and proton are of the size expected for a simple hydrated amine (a result supported by density functional theory calculations), and also do not titrate with the carboxylates. The line shapes of the J-coupled (15)N signals suggest rapid internal reorientation of all NH(3)(+) groups. pK(a) values have been measured for all charged side chains except Glu50 and do not show the perturbations expected for salt bridge formation, except that E35 has a Hill coefficient of 0.84. The main differential effect seen is that the lysines that are involved in salt bridges in the crystal display faster exchange of the amine protons with the solvent, an effect attributed to general base catalysis by the carboxylates. This explanation is supported by varying buffer composition, which demonstrates reduced electrostatic shielding at low concentration. In conclusion, the study demonstrates that the six surface-exposed lysines in protein G are not involved in significant salt bridge interactions, even though such interactions are found consistently in crystal structures. However, the intrahelical E35-K39 (i,i+4) interaction is partially present.

A Second-Harmonic-Active Unnatural Amino Acid as a Structural Probe of Biomolecules on Surfaces

Second-harmonic generation (SHG) is highly sensitive to the net, average orientation of SH-active molecules on surfaces and has recently emerged as a technique for detecting biomolecules and their conformational changes. As most biomolecules are not intrinsically SH-active, they must be labeled with probes to render them detectable. To date, exogenous probes have been used to do this, but second-harmonic-active unnatural amino acids offer important advantages for the long-range goal of precisely and directly determining structural changes in real time and may be used for both buried and surface sites. Results of the first known SH-active unnatural amino acid, Aladan, are presented here. Aladan is found to be SH-active by detecting it at an interface, both alone and incorporated into the B1 domain of protein G (GB1), a globular immunoglobulin-binding protein, at both buried and exposed sites. The tilt angle of Aladan alone adsorbed on a mica surface is determined by polarization experiments, and its nonlinear polarizability alpha((2)) is found to be ca. 10(-30) esu. Aladan GB1 mutants are detectable by SHG, either when coupled covalently to a derivatized glass surface or bound to IgG immobilized via protein A. Addition of an Fc domain to this GB1 complex causes a small but defined change in the SH signal when Aladan is incorporated at site Ala(24), but not at Leu(7), consistent with a local conformational change of GB1. This structural change is not apparent in either X-ray crystallography or fluorescence studies, demonstrating that SHG can detect subtle orientational changes, including protein-protein interactions in which no significant rearrangements occur.

Evolutional selection of a combinatorial phage library displaying randomly-rearranged various single domains of immunoglobulin (Ig)-binding proteins (IBPs) with four kinds of Ig molecules

BackgroundProtein A, protein G and protein L are three well-defined immunoglobulin (Ig)-binding proteins (IBPs), which show affinity for specific sites on Ig of mammalian hosts. Although the precise functions of these molecules are not fully understood, it is thought that they play an important role in pathogenicity of bacteria. The single domains of protein A, protein G and protein L were all demonstrated to have function to bind to Ig. Whether combinations of Ig-binding domains of various IBPs could exhibit useful novel binding is interesting.ResultsWe used a combinatorial phage library which displayed randomly-rearranged various-peptide-linked molecules of D and A domains of protein A, designated PA(D) and PA(A) respectively, B2 domain of protein G (PG) and B3 domain of protein L (PL) for affinity selection with human IgG (hIgG), human IgM (hIgM), human IgA (hIgA) and recombinant hIgG1-Fc as bait respectively. Two kinds of novel combinatorial molecules with characteristic structure of PA(A)-PG and PA(A)-PL were obtained in hIgG (hIgG1-Fc) and hIgM (hIgA) post-selection populations respectively. In addition, the linking peptides among all PA(A)-PG and PA(A)-PL structures was strongly selected, and showed interestingly divergent and convergent distribution. The phage binding assays and competitive inhibition experiments demonstrated that PA(A)-PG and PA(A)-PL combinations possess comparable binding advantages with hIgG/hIgG1-Fc and hIgM/hIgA respectively.ConclusionIn this work, a combinatorial phage library displaying Ig-binding domains of protein A, protein G, or protein L joined by various random linking peptides was used to conducted evolutional selection in vitro with four kinds of Ig molecules. Two kinds of novel combinations of Ig-binding domains, PA(A)-PG and PA(A)-PL, were obtained, and demonstrate the novel Ig binding properties.

Quantifying the Structural Requirements of the Folding Transition State of Protein A and Other Systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of ψ analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational ϕ analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger α/β proteins, that the TSEs of most small proteins achieve ∼ 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

Mechanism and Energy Landscape of Domain Swapping in the B1 Domain of Protein G

Three-dimensional domain swapping has emerged as a ubiquitous process for homo-oligomer formation in many unrelated proteins, but the molecular mechanism of this process is still poorly understood. Here we present a mechanism for the swapping reaction in the B1 domain of the immunoglobulin G binding protein from group G of Streptococcus (GB1). This is a particularly attractive system for investigating the swapping process, as the swapped dimer formed by the quadruple mutant (L5V/F30V/Y33F/A34F) of GB1 was recently shown to exist in equilibrium with a monomer-like conformation over time scales of minutes. According to our mechanism, swapping in GB1 starts from the C-terminus of the polypeptide chain and progresses by exchanging an increasing portion of the chains until a stable conformational state is reached. This exchange process does not involve unfolding. Rather, the conformational changes of individual monomers and their association are tightly coupled to minimize solvent exposure and maximize the total number of native contacts at all times, thereby closely approximating the minimum energy path of the reaction. Using detailed atomic descriptions, we compute the complete free-energy profiles of the exchange reaction for the GB1 quadruple mutant that forms swapped dimers and for the wild-type protein, which is monomeric. In both GB1 forms, intermediates sample a surprisingly wide range of nearly isoenergetic association modes and hinge conformations, indicating that the exchange reaction is a non-specific process akin to encounter complex formation where the amino acid sequence plays a marginal role. The main role of the mutations in the swapping process is to destabilize the GB1 monomer state, while stabilizing the swapped dimer conformation, with non-native intersubunit interactions, fostered by mutant side chains, contributing significantly to this stabilization. Our findings are rationalized in terms of a generic swapping mechanism that involves the association of activated molecular species, and it is argued that a similar mechanism may apply to swapping in other protein systems.

A systematic analysis of backbone amide assignments achieved via combinatorial selective labelling of amino acids

With the advent of high-yield cell-free expressions systems, many researchers are exploiting selective isotope labelling of amino acids to increase the efficiency and accuracy of the NMR assignment process. We developed recently a combinatorial selective labelling (CSL) method capable of yielding large numbers of residue-type and sequence-specific backbone amide assignments, which involves comparing cross-peak intensities in 1H- 15N HSQC and 2D 1H- 15N HNCO spectra collected for five samples containing different combinations of 13C- and 15N-labelled amino acids [Parker MJ, Aulton-Jones M, Hounslow A, Craven C J (2004) J Am Chem Soc 126:5020-5021]. In this paper we develop a robust method for establishing the reliability of these assignments. We have performed a detailed statistical analysis of the CSL data collected for a model system (the B1 domain of protein G from Streptococcus), developing a scoring method which allows the confidence in assignments to be assessed, and which enables the effects of overlap on assignment fidelity to be predicted. To further test the scoring method and also to assess the performance of CSL in relation to sample quality, we have applied the method to the CSL data collected for GFP in our previous study.