Conformational characterization of disulfide bonds: A tool for protein classification

Background

Throughout evolution, mutations in particular regions of some protein structures have resulted in extra covalent bonds that increase the overall robustness of the fold: disulfide bonds. The two strategically placed cysteines can also have a more direct role in protein function, either by assisting thiol or disulfide exchange, or through allosteric effects. In this work, we verified how the structural similarities between disulfides can reflect functional and evolutionary relationships between different proteins. We analyzed the conformational patterns of the disulfide bonds in a set of disulfide-rich proteins that included twelve SCOP superfamilies: thioredoxin-like and eleven superfamilies containing small disulfide-rich proteins (SDP).

Results

The twenty conformations considered in the present study were characterized by both structural and energetic parameters. The corresponding frequencies present diverse patterns for the different superfamilies. The least-strained conformations are more abundant for the SDP superfamilies, while the “catalytic” +/−RHook is dominant for the thioredoxin-like superfamily. The “allosteric” -RHSaple is moderately abundant for BBI, Crisp and Thioredoxin-like superfamilies and less frequent for the remaining superfamilies. Using a hierarchical clustering analysis we found that the twelve superfamilies were grouped in biologically significant clusters.

Conclusions

In this work, we carried out an extensive statistical analysis of the conformational motifs for the disulfide bonds present in a set of disulfide-rich proteins. We show that the conformational patterns observed in disulfide bonds are sufficient to group proteins that share both functional and structural patterns and can therefore be used as a criterion for protein classification.     

Two reaction pathways for transformation of high potential cytochrome b559 of PS II into the intermediate potential form

This study describes an analysis of different treatments that influence the relative content and the midpoint potential of HP Cyt b559 in PS II membrane fragments from higher plants. Two basically different types of irreversible modification effects are distinguished: the HP form of Cyt b559 is either predominantly affected when the heme group is oxidized (“O-type” effects) or when it is reduced (“R-type” effects). Transformation of HP Cyt b559 to lower potential redox forms (IP and LP forms) by the “O-type” mechanism is induced by high pH and detergent treatments. In this case the effects consist of a gradual decrease in the relative content of HP Cyt b559 while its midpoint potential remains unaffected. Transformation of HP Cyt b559 via an “R-type” mechanism is caused by a number of exogenous compounds denoted L: herbicides, ADRY reagents and tetraphenylboron. These compounds are postulated to bind to the PS II complex at a quinone binding site designated as Q_C which interacts with Cyt b559 and is clearly not the Q_B site. Binding of compounds L to the Q_C site when HP Cyt b559 is oxidized gives rise to a gradual decrease in the E_m of HP Cyt b559 with increasing concentration of L (up to 10 K_ox(L) values) while the relative content of HP Cyt b559 is unaffected. Higher concentrations of compounds L required for their binding to Q_C site when HP Cyt b559 is reduced (described by K_red(L)) induce a conversion of HP Cyt b559 to lower potential redox forms (“R-type” transformation). Two reaction pathways for transitions of Cyt b559 between the different protein conformations that are responsible for the HP and IP/LP redox forms are proposed and new insights into the functional regulation of Cyt b559 via the Q_C site are discussed.     

Structural diversity in a human antibody germline library

To support antibody therapeutic development, the crystal structures of a set of 16 germline variants composed of 4 different kappa light chains paired with 4 different heavy chains have been determined. All four heavy chains of the antigen-binding fragments (Fabs) have the same complementarity-determining region (CDR) H3 that was reported in an earlier Fab structure. The structure analyses include comparisons of the overall structures, canonical structures of the CDRs and the VH:VL packing interactions. The CDR conformations for the most part are tightly clustered, especially for the ones with shorter lengths. The longer CDRs with tandem glycines or serines have more conformational diversity than the others. CDR H3, despite having the same amino acid sequence, exhibits the largest conformational diversity. About half of the structures have CDR H3 conformations similar to that of the parent; the others diverge significantly. One conclusion is that the CDR H3 conformations are influenced by both their amino acid sequence and their structural environment determined by the heavy and light chain pairing. The stem regions of 14 of the variant pairs are in the ‘kinked’ conformation, and only 2 are in the extended conformation. The packing of the VH and VL domains is consistent with our knowledge of antibody structure, and the tilt angles between these domains cover a range of 11 degrees. Two of 16 structures showed particularly large variations in the tilt angles when compared with the other pairings. The structures and their analyses provide a rich foundation for future antibody modeling and engineering efforts.     

Structural refinement of membrane proteins by restrained molecular dynamics and solvent accessibility data

We present an approach for incorporating solvent accessibility data from electron paramagnetic resonance experiments in the structural refinement of membrane proteins through restrained molecular dynamics simulations. The restraints have been parameterized from oxygen (ΠO₂) and nickel-ethylenediaminediacetic acid (ΠNiEdda) collision frequencies, as indicators of lipid or aqueous exposed spin-label sites. These are enforced through interactions between a pseudoatom representation of the covalently attached Nitroxide spin-label and virtual “solvent” particles corresponding to O₂ and NiEdda in the surrounding environment. Interactions were computed using an empirical potential function, where the parameters have been optimized to account for the different accessibilities of the spin-label pseudoatoms to the surrounding environment. This approach, “pseudoatom-driven solvent accessibility refinement”, was validated by refolding distorted conformations of the Streptomyces lividans potassium channel (KcsA), corresponding to a range of 2-30 Å root mean-square deviations away from the native structure. Molecular dynamics simulations based on up to 58 electron paramagnetic resonance restraints derived from spin-label mutants were able to converge toward the native structure within 1-3 Å root mean-square deviations with minimal computational cost. The use of energy-based ranking and structure similarity clustering as selection criteria helped in the convergence and identification of correctly folded structures from a large number of simulations. This approach can be applied to a variety of integral membrane protein systems, regardless of oligomeric state, and should be particularly useful in calculating conformational changes from a known reference crystal structure.     

Predicting antibody hypervariable loop conformations. II: Minimization and molecular dynamics studies of MCPC603 from many randomly generated loop conformations

We describe a method for predicting the conformations of loops in proteins and its application to four of the complementarity determining regions [CDRs] in the crystallographically determined structure of MCPC603. The method is based on the generation of a large number of randomly generated conformations for the backbone of the loop being studied, followed by either minimization or molecular dynamics followed by minimization starting from these random structures. The details of the algorithm for the generation of the loops are presented in the first paper in this series (Shenkin et al. [submitted]). The results of minimization and molecular dynamics applied to these loops is presented here. For the two shortest CDRs studied (H1 and L2, which are five and seven amino acids long), minimizations and dynamics simulations which ignore interactions of the loop amino acids beyond the carbon beta replicate the conformation of the crystal structure closely. This suggests that these loops fold independently of sequence variation. For the third CDR (L3, which is nine amino acids), those portions of the CDR near its base which are hydrogen bonded to framework are well replicated by our procedures, but the top of the loop shows significant conformational variability. This variability persists when side chain interactions for the MCPC603 sequence are included. For a fourth CDR (H3, which is 11 amino acids long), new low-energy backbone conformations are found; however, only those which are close to the crystal are compatible with the sequence when side chain interactions are taken into account. Results from minimization and dynamics on single CDRs with all other CDRs removed are presented. These allow us to explore the extent to which individual CDR conformations are determined by interactions with framework only.     

Backbone Flexibility of CDR3 and Immune Recognition of Antigens

Conformational entropy is an important component of protein–protein interactions; however, there is no reliable method for computing this parameter. We have developed a statistical measure of residual backbone entropy in folded proteins by using the ?–ψ distributions of the 20 amino acids in common secondary structures. The backbone entropy patterns of amino acids within helix, sheet or coil form clusters that recapitulate the branching and hydrogen bonding properties of the side chains in the secondary structure type. The same types of residues in coil and sheet have identical backbone entropies, while helix residues have much smaller conformational entropies. We estimated the backbone entropy change for immunoglobulin complementarity-determining regions (CDRs) from the crystal structures of 34 low-affinity T-cell receptors and 40 high-affinity Fabs as a result of the formation of protein complexes. Surprisingly, we discovered that the computed backbone entropy loss of only the CDR3, but not all CDRs, correlated significantly with the kinetic and affinity constants of the 74 selected complexes. Consequently, we propose a simple algorithm to introduce proline mutations that restrict the conformational flexibility of CDRs and enhance the kinetics and affinity of immunoglobulin interactions. Combining the proline mutations with rationally designed mutants from a previous study led to 2400-fold increase in the affinity of the A6 T-cell receptor for Tax-HLAA2. However, this mutational scheme failed to induce significant binding changes in the already-high-affinity C225–Fab/huEGFR interface. Our results will serve as a roadmap to formulate more effective target functions to design immune complexes with improved biological functions.     

Length-independent structural similarities enrich the antibody CDR canonical class model

Complementarity-determining regions (CDRs) are antibody loops that make up the antigen binding site. Here, we show that all CDR types have structurally similar loops of different lengths. Based on these findings, we created length-independent canonical classes for the non-H3 CDRs. Our length variable structural clusters show strong sequence patterns suggesting either that they evolved from the same original structure or result from some form of convergence. We find that our length-independent method not only clusters a larger number of CDRs, but also predicts canonical class from sequence better than the standard length-dependent approach.

To demonstrate the usefulness of our findings, we predicted cluster membership of CDR-L3 sequences from 3 next-generation sequencing datasets of the antibody repertoire (over 1,000,000 sequences). Using the length-independent clusters, we can structurally classify an additional 135,000 sequences, which represents a ～20% improvement over the standard approach. This suggests that our length-independent canonical classes might be a highly prevalent feature of antibody space, and could substantially improve our ability to accurately predict the structure of novel CDRs identified by next-generation sequencing.     

Extensive conformational transitions are required to turn on ATP hydrolysis in myosin

Conventional myosin is representative of biomolecular motors in which the hydrolysis of adenosine triphosphate (ATP) is coupled to large-scale structural transitions both in and remote from the active site. The mechanism that underlies such “mechanochemical coupling,” especially the causal relationship between hydrolysis and allosteric structural changes, has remained elusive despite extensive experimental and computational analyses. In this study, using combined quantum mechanical and molecular mechanical simulations and different conformations of the myosin motor domain, we provide evidence to support that regulation of ATP hydrolysis activity is not limited to residues in the immediate environment of the phosphate. Specifically, we illustrate that efficient hydrolysis of ATP depends not only on the proper orientation of the lytic water but also on the structural stability of several nearby residues, especially the Arg238-Glu459 salt bridge (the numbering of residues follows myosin II in Dictyostelium discoideum) and the water molecule that spans this salt bridge and the lytic water. More importantly, by comparing the hydrolysis activities in two motor conformations with very similar active-site (i.e., Switches I and II) configurations, which distinguished this work from our previous study, the results clearly indicate that the ability of these residues to perform crucial electrostatic stabilization relies on the configuration of residues in the nearby N-terminus of the relay helix and the “wedge loop.” Without the structural support from those motifs, residues in a closed active site in the post-rigor motor domain undergo subtle structural variations that lead to consistently higher calculated ATP hydrolysis barriers than in the pre-powerstroke state. In other words, starting from the post-rigor state, turning on the ATPase activity requires not only displacement of Switch II to close the active site but also structural transitions in the N-terminus of the relay helix and the “wedge loop,” which have been proposed previously to be ultimately coupled to the rotation of the converter subdomain 40 Å away.     

The crystal structure analysis of group B Streptococcus sortase C1: a model for the "lid" movement upon substrate binding

A unique feature of the class-C-type sortases, enzymes essential for Gram-positive pilus biogenesis, is the presence of a flexible “lid” anchored in the active site. However, the mechanistic details of the “lid” displacement, suggested to be a critical prelude for enzyme catalysis, are not yet known. This is partly due to the absence of enzyme-substrate and enzyme-inhibitor complex crystal structures. We have recently described the crystal structures of the Streptococcus agalactiae SAG2603 V/R sortase SrtC1 in two space groups (type II and type III) and that of its “lid” mutant and proposed a role of the “lid” as a protector of the active-site hydrophobic environment. Here, we report the crystal structures of SAG2603 V/R sortase C1 in a different space group (type I) and that of its complex with a small-molecule cysteine protease inhibitor. We observe that the catalytic Cys residue is covalently linked to the small-molecule inhibitor without lid displacement. However, the type I structure provides a view of the sortase SrtC1 lid displacement while having structural elements similar to a substrate sorting motif suitably positioned in the active site. We propose that these major conformational changes seen in the presence of a substrate mimic in the active site may represent universal features of class C sortase substrate recognition and enzyme activation.     

Thermal and chemical unfolding and refolding of a eukaryotic sodium channel

Voltage-gated sodium channels are dynamic membrane proteins essential for signaling in nervous and muscular systems. They undergo substantial conformational changes associated with the closed, open and inactivated states. However, little information is available regarding their conformational stability. In this study circular dichroism spectroscopy was used to investigate the changes in secondary structure accompanying chemical and thermal denaturation of detergent-solubilised sodium channels isolated from Electrophorus electricus electroplax. The proteins appear to be remarkably resistant to either type of treatment, with “denatured” channels, retaining significant helical secondary structure even at 77 °C or in 10% SDS. Further retention of helical secondary structure at high temperature was observed in the presence of the channel-blocking tetrodotoxin. It was possible to refold the thermally-denatured (but not chemically-denatured) channels in vitro. The correctly refolded channels were capable of undergoing the toxin-induced conformational change indicative of ligand binding. In addition, flux measurements in liposomes showed that the thermally-denatured (but not chemically-denatured) proteins were able to re-adopt native, active conformations. These studies suggest that whilst sodium channels must be sufficiently flexible to undergo major conformational changes during their functional cycle, the proteins are highly resistant to unfolding, a feature that is important for maintaining structural integrity during dynamic processes.     

Zipping and Unzipping of Adenylate Kinase: Atomistic Insights into the Ensemble of Open ↔ Closed Transitions

Adenylate kinase (AdK), a phosphotransferase enzyme, plays an important role in cellular energy homeostasis. It undergoes a large conformational change between an open and a closed state, even in the absence of substrate. We investigate the apo-AdK transition at the atomic level both with free-energy calculations and with our new dynamic importance sampling (DIMS) molecular dynamics method. DIMS is shown to sample biologically relevant conformations as verified by comparing an ensemble of hundreds of DIMS transitions to AdK crystal structure intermediates. The simulations reveal in atomic detail how hinge regions partially and intermittently unfold during the transition. Conserved salt bridges are seen to have important structural and dynamic roles; in particular, four ionic bonds that open in a sequential, zipper-like fashion and, thus, dominate the free-energy landscape of the transition are identified. Transitions between the closed and open conformations only have to overcome moderate free-energy barriers. Unexpectedly, the closed state and the open state encompass broad free-energy basins that contain conformations differing in domain hinge motions by up to 40°. The significance of these extended states is discussed in relation to recent experimental Förster resonance energy transfer measurements. Taken together, these results demonstrate how a small number of cooperative key interactions can shape the overall dynamics of an enzyme and suggest an “all-or-nothing” mechanism for the opening and closing of AdK. Our efficient DIMS molecular dynamics computer simulation approach can provide a detailed picture of a functionally important macromolecular transition and thus help to interpret and suggest experiments to probe the conformational landscape of dynamic proteins such as AdK.     

The diversity of H3 loops determines the antigen‐binding tendencies of antibody CDR loops

Of the complementarity‐determining regions (CDRs) of antibodies, H3 loops, with varying amino acid sequences and loop lengths, adopt particularly diverse loop conformations. The diversity of H3 conformations produces an array of antigen recognition patterns involving all the CDRs, in which the residue positions actually in contact with the antigen vary considerably. Therefore, for a deeper understanding of antigen recognition, it is necessary to relate the sequence and structural properties of each residue position in each CDR loop to its ability to bind antigens. In this study, we proposed a new method for characterizing the structural features of the CDR loops and obtained the antigen‐binding ability of each residue position in each CDR loop. This analysis led to a simple set of rules for identifying probable antigen‐binding residues. We also found that the diversity of H3 loop lengths and conformations affects the antigen‐binding tendencies of all the CDR loops.     

Combined effects of collagen type I alpha1 (COL1A1) Sp1 polymorphism and osteoporosis risk factors on bone mineral density in Turkish postmenopausal women

Identification of risk factors for osteoporosis has been essential for understanding the development of osteoporosis. The collagen type I alpha1 (COL1A1) gene is suggested to be implicated in reduced bone mineral density (BMD) in osteoporosis. In the present study, the investigation of the effects of Sp1 polymorphic variants of COL1A1 gene on BMD values, and the determination of the association between COL1A1 Sp1 gene variants and osteoporosis risk factors in the context of gene–environment interaction in Turkish postmenopausal women were aimed. For the detection of COL1A1 Sp1 polymorphism, PCR-RFLP techniques have been used. BMD for lumbar spine (L1–L4) and hip (femoral neck and total hip) was measured by DXA. This study was carried out using a sample of 254 postmenopausal women. We observed a trend decrease in BMD values in the subjects with “ss” genotype having lower BMD of lumbar spine, femoral neck and total hip than those with “SS” and “Ss” genotype, however the differences did not reach statistical significance (P > 0.05). We also found that the frequencies of the BMD under mean values at the femoral neck (57.5%) and total hip (76.2%) increased considerably in the subjects carrying “Ss/ss” genotypes in combination of having family history of osteoporosis (61.5% for femoral neck) and smoking history (90.0% for total hip). This population-based study indicates that COL1A1 Sp1 polymorphism may contribute to the development of osteoporosis in combination of osteoporosis risk factors in Turkish postmenopausal women.     

On stabilization of a neutral aromatic ligand by π-cation interactions in monoclonal antibodies

It has been shown that anti-PAH mAb can bind a particular cross-reactant by adopting two distinct “red” and “blue” conformations of its binding sites [N.M. Grubor et al. PNAS 102, 2005, 7453-7458]. In the case of red conformation of pyrene (Py)/anti-PAH mAb (with a broad fluorescence (0,0)-band with fwhm ~ 140 cm⁻¹), the central role in complex formation was played by π-π interactions. The nature of the blue-shifted conformation with very narrow fluorescence (0,0)-band (fwhm ~ 75 cm⁻¹) was left unclear due to the lack of suitable data for comparison. In this work, we suggest spectroscopic and modeling results obtained for the blue conformation of Py in several mAb (including 4D5 mAb) are consistent with π-cation interactions, underscoring the importance of π-cation interaction in ligand binding and stabilization in agreement with earlier modeling studies [J-L. Pellequer, et al. J. Mol. Biol. 302, 2000, 691-699]. We propose considerable narrowing of the fluorescence origin band of ligand in the protein environment could be regarded as a simple indicator of π-cation interactions. Since 4D5 mAb forms only the blue-shifted conformation, while anti-PAH and 8E11 mAbs form both blue- and red-shifted conformations, we suggest mAb interactions, with Py molecules lacking H-bonding functionality, may induce distinct conformations of mAb binding sites that allow binding by π-π and/or π-cation interactions.     

Analysis of binding site hot spots on the surface of Ras GTPase

We have recently discovered an allosteric switch in Ras, bringing an additional level of complexity to this GTPase whose mutants are involved in nearly 30% of cancers. Upon activation of the allosteric switch, there is a shift in helix 3/loop 7 associated with a disorder to order transition in the active site. Here, we use a combination of multiple solvent crystal structures and computational solvent mapping (FTMap) to determine binding site hot spots in the “off” and “on” allosteric states of the GTP-bound form of H-Ras. Thirteen sites are revealed, expanding possible target sites for ligand binding well beyond the active site. Comparison of FTMaps for the H and K isoforms reveals essentially identical hot spots. Furthermore, using NMR measurements of spin relaxation, we determined that K-Ras exhibits global conformational dynamics very similar to those we previously reported for H-Ras. We thus hypothesize that the global conformational rearrangement serves as a mechanism for allosteric coupling between the effector interface and remote hot spots in all Ras isoforms. At least with respect to the binding sites involving the G domain, H-Ras is an excellent model for K-Ras and probably N-Ras as well. Ras has so far been elusive as a target for drug design. The present work identifies various unexplored hot spots throughout the entire surface of Ras, extending the focus from the disordered active site to well-ordered locations that should be easier to target.     

Electrochemical detection of Pseudomonas aeruginosa 16S rRNA using a biosensor based on immobilized stem-loop structured probe

We have designed an electrochemical DNA biosensor based on stem-loop structured probes for enzymatic detection of Pseudomonas aeruginosa 16S ribosomal RNA (rRNA) in composting degradation. The probe modified with a thiol at its 5′ end and a biotin at its 3′ end was immobilized on a gold electrode through self-assembly. The stem-loop structured probes were “closed” when target was absent, then the hybridization of the target induced the conformational changes to “open”, along with the biotin at its 3′ end binding with streptavidin-horseradish peroxidase (HRP), and subsequent quanti?cation of the target was detected via electrochemical detecting the enzymatic product in the presence of substrate. Under the optimum experiment conditions, the amperometric current response to HRP-catalyzed reaction was linearly related to the logarithm of the target nucleic acid concentration, ranging from 0.3 and 600 pg/μL, with the detection limit of 0.012 pg/μL. A correlation coefficient of 0.9960 was identified. The 16S rRNA extracted from P. aeruginosa was analyzed by this proposed sensor. The results were in agreement with the reference values deduced from UV spectrometric data. The biosensor was indicative of good precision, stability, sensitivity, and selectivity.     

Does conformational free energy distinguish loop conformations in proteins?

Limitations in protein homology modeling often arise from the inability to adequately model loops. In this paper we focus on the selection of loop conformations. We present a complete computational treatment that allows the screening of loop conformations to identify those that best fit a molecular model. The stability of a loop in a protein is evaluated via computations of conformational free energies in solution, i.e., the free energy difference between the reference structure and the modeled one. A thermodynamic cycle is used for calculation of the conformational free energy, in which the total free energy of the reference state (i.e., gas phase) is the CHARMm potential energy. The electrostatic contribution of the solvation free energy is obtained from solving the finite-difference Poisson-Boltzmann equation. The nonpolar contribution is based on a surface area-based expression. We applied this computational scheme to a simple but well-characterized system, the antibody hypervariable loop (complementarity-determining region, CDR). Instead of creating loop conformations, we generated a database of loops extracted from high-resolution crystal structures of proteins, which display geometrical similarities with antibody CDRs. We inserted loops from our database into a framework of an antibody; then we calculated the conformational free energies of each loop. Results show that we successfully identified loops with a "reference-like" CDR geometry, with the lowest conformational free energy in gas phase only. Surprisingly, the solvation energy term plays a confusing role, sometimes discriminating "reference-like" CDR geometry and many times allowing "non-reference-like" conformations to have the lowest conformational free energies (for short loops). Most "reference-like" loop conformations are separated from others by a gap in the gas phase conformational free energy scale. Naturally, loops from antibody molecules are found to be the best models for long CDRs (> or = 6 residues), mainly because of a better packing of backbone atoms into the framework of the antibody model.