Elucidation of the ribonuclease A aggregation process mediated by 3D domain swapping: a computational approach reveals possible new multimeric structures

By lyophilization from 40% acetic acid solutions, bovine pancreatic ribonuclease A forms several three-dimensional (3D) domain-swapped oligomers: dimers, trimers, tetramers, pentamers, hexamers, and traces of high-order oligomers, purifiable by cation-exchange chromatography. Each oligomeric species consists of at least two conformers displaying different basicity density, and/or exposure of positive charges. The structures of the two dimers and one trimer have been solved. Plausible models have been proposed for a second RNase A trimer and four tetramers, but not all the models are certainly assignable to the tetramers purified. Further studies have also been made on the pentameric and hexameric species, again without reaching structurally clear-cut results. This work is focused on the detailed modeling of the tetrameric RNase A species, using four different approaches to possibly clarify unknown structural aspects. The results obtained do not confirm the validity of one tetrameric model previously proposed, but allow the proposal of a novel tetrameric structure displaying new interfaces that are absent in the other known conformers. New details concerning other tetrameric structures are also described. RNase A multimers larger than tetramers, i.e., pentamers, hexamers, octamers, nonamers, up to dodecamers, are also modeled, with the proposal of novel domain-swapped structures, and the confirmation of what had previously been inferred. Finally, the propensity of RNase A to possibly form high-order supramolecular multimers is analyzed starting from the large number of domain-swapped RNase A conformers modeled.     

Propensity for C-terminal domain swapping correlates with increased regional flexibility in the C-terminus of RNase A

Domain swapping is a type of oligomerization in which monomeric proteins exchange a structural element, resulting in oligomers whose subunits recapitulate the native, monomeric fold. It has been implicated as a potential mechanism for protein aggregation, which provides a strong impetus to understand the structural determinants and folding mechanisms that trigger domain swapping. Bovine pancreatic ribonuclease A (RNase A) is a well-studied protein known to domain swap under extreme conditions, such as lyophilization from acetic acid. The major domain-swapped dimer form of RNase A exchanges a β-strand at its C-terminus to form a C-terminal domain-swapped dimer. To study the mechanism by which C-terminal swapping occurs, we used a variant of RNase A containing a P114G mutation that readily domain swaps under physiological conditions. Using NMR and hydrogen-deuterium exchange, we find that the P114G variant has decreased protection from hydrogen exchange compared to the wild-type protein near the C-terminal hinge region. Our results suggest that domain swapping occurs via a local high-energy fluctuation at the C-terminus.     

Substitution of asparagine 76 by a tyrosine residue induces domain swapping in Helicobacter pylori phosphopantetheine adenylyltransferase

Phosphopantetheine adenylyltransferase (PPAT) catalyses the penultimate step in coenzyme A biosynthesis in bacteria and is therefore a candidate target for antibacterial drug development. We randomly mutated the residues in the Helicobacter pylori PPAT sequence to identify those that govern protein folding and ligand binding, and we describe the crystal structure of one of these mutants (I4V/N76Y) that contains the mutations I4?→?V and N76?→?Y. Unlike other PPATs, which are homohexamers, I4V/N76Y is a domain-swapped homotetramer. The protomer structure of this mutant is an open conformation in which the 65 C-terminal residues are intertwined with those of a neighbouring protomer. Despite structural differences between wild-type PPAT and IV4/N76Y, they had similar ligand-binding properties. ATP binding to these two proteins was enthalpically driven, whereas that for Escherichia coli PPAT is entropically driven. The structural packing of the subunits may affect the thermal denaturation of wild-type PPAT and I4V/N76Y. Mutations in hinge regions often induce domain swapping, i.e. the spatial exchange of portions of adjacent protomers, but residues 4 and 76 of H. pylori PPAT are not located in or near to the hinge region. However, one or both of these residues is responsible for the large conformational change in the C-terminal region of each protomer. To identify the residue(s) responsible, we constructed the single-site mutant, N76Y, and found a large displacement of α-helix 4, which indicated that its flexibility allowed the domain swap to occur.     

Structures of the two 3D domain-swapped RNase A trimers

When concentrated in mildly acidic solutions, bovine pancreatic ribonuclease (RNase A) forms long-lived oligomers including two types of dimer, two types of trimer, and higher oligomers. In previous crystallographic work, we found that the major dimeric component forms by a swapping of the C-terminal beta-strands between the monomers, and that the minor dimeric component forms by swapping the N-terminal alpha-helices of the monomers. On the basis of these structures, we proposed that a linear RNase A trimer can form from a central molecule that simultaneously swaps its N-terminal helix with a second RNase A molecule and its C-terminal strand with a third molecule. Studies by dissociation are consistent with this model for the major trimeric component: the major trimer dissociates into both the major and the minor dimers, as well as monomers. In contrast, the minor trimer component dissociates into the monomer and the major dimer. This suggests that the minor trimer is cyclic, formed from three monomers that swap their C-terminal beta-strands into identical molecules. These conclusions are supported by cross-linking of lysyl residues, showing that the major trimer swaps its N-terminal helix, and the minor trimer does not. We verified by X-ray crystallography the proposed cyclic structure for the minor trimer, with swapping of the C-terminal beta-strands. This study thus expands the variety of domain-swapped oligomers by revealing the first example of a protein that can form both a linear and a cyclic domain-swapped oligomer. These structures permit interpretation of the enzymatic activities of the RNase A oligomers on double-stranded RNA.     

Three-dimensional domain swapping as a mechanism to lock the active conformation in a super-active octamer of SARS-CoV main protease

Proteolytic processing of viral polyproteins is indispensible for the lifecycle of coronaviruses. The main protease (M^pro) of SARS-CoV is an attractive target for anti-SARS drug development as it is essential for the polyprotein processing. M^pro is initially produced as part of viral polyproteins and it is matured by autocleavage. Here, we report that, with the addition of an N-terminal extension peptide, M^pro can form a domain-swapped dimer. After complete removal of the extension peptide from the dimer, the mature M^pro self-assembles into a novel super-active octamer (AO-M^pro). The crystal structure of AO-M^pro adopts a novel fold with four domain-swapped dimers packing into four active units with nearly identical conformation to that of the previously reported M^pro active dimer, and 3D domain swapping serves as a mechanism to lock the active conformation due to entanglement of polypeptide chains. Compared with the previously well characterized form of M^pro, in equilibrium between inactive monomer and active dimer, the stable AO-M^pro exhibits much higher proteolytic activity at low concentration. As all eight active sites are bound with inhibitors, the polyvalent nature of the interaction between AO-M^pro and its polyprotein substrates with multiple cleavage sites, would make AO-M^pro functionally much more superior than the M^pro active dimer for polyprotein processing. Thus, during the initial period of SARS-CoV infection, this novel active form AO-M^pro should play a major role in cleaving polyproteins as the protein level is extremely low. The discovery of AO-M^pro provides new insights about the functional mechanism of M^pro and its maturation process.     

Crowding agents and osmolytes provide insight into the formation and dissociation of RNase A oligomers

RNase A forms 3D domain-swapped oligomers with novel enzymatic and biological activities. We study how crowding agents and osmolytes affect the formation and dissociation of RNase A oligomers. The crowding agents Ficoll and dextran were found to enhance oligomer formation, whereas the stabilizers sodium sulfate, glycine and trimethylammonium oxide (TMAO) do not. In contrast, TMAO significantly slows RNase A dimer dissociation, while the effect of Ficoll is small. These results lead us to propose that the mechanisms of oligomer formation and dissociation are different. In the RNase A “C-dimer”, the C-terminal β-strand is swapped between two subunits. The loop preceding this β-strand adopts a β-sheet which has been proposed to resemble amyloid structurally. Hydrogen/deuterium (H/D) exchange of the RNase A C-dimer reveal that the H-bonds formed between the swapped C-terminal β-strand and the other subunit are strong. Their rupture may be crucial for C-dimer dissociation. In contrast, H-bonds formed by Asn 113 in the novel β-sheet adopted by the hinge loop in the C-dimer are not strongly protected. Besides the fundamental insights obtained, the results represent a technical advance for obtaining increased oligomer yields and storage lifetimes.     

A single mutation induces amyloid aggregation in the alpha-spectrin SH3 domain: analysis of the early stages of fibril formation

The Src-homology region 3 domain of chicken alpha-spectrin (Spc-SH3) is a small two-state folding protein, which has never been described to form amyloid fibrils under any condition investigated so far. We show here that the mutation of asparagine 47 to alanine at the distal loop, which destabilises similarly the native and folding transition states of the domain, induces the formation of amyloid fibrils under mild acid conditions. Amyloid aggregation of the mutant is enhanced by the increase in temperature, protein concentration and NaCl concentration. The early stages of amyloid formation have been monitored as a function of time and temperature using a variety of biophysical methods. Differential scanning calorimetry experiments under conditions of amyloid formation have allowed the identification of different thermal transitions corresponding to conformational and aggregation processes as well as to the high-temperature disaggregation and unfolding of the amyloid fibrils. Aggregation is preceded by a rapid conformational change in the monomeric domain involving about 40% of the global unfolding enthalpy, considerable change in secondary structure, large loss of tertiary structure and exposure of hydrophobic patches to the solvent. The conformational change is followed by formation of a majority of oligomeric species with apparent hydrodynamic radius between 2.5 nm and 10nm, depending on temperature, together with the appearance and progressive growth of protofibrillar aggregates. After these early aggregation stages, long and curved fibrils of up to several micrometers start to develop by elongation of the protofibrils. The calorimetric data indicate that the specific enthalpy of fibril disaggregation and unfolding is relatively low, suggesting a low density of interactions within the fibril structure as compared to the native protein and a main entropy contribution to the stability of the amyloid fibrils.     

Loss of CYP3A7 gene induction by 1,25-dihydroxyvitamin D3 is caused by less binding of VDR to the proximal ER6 in CYP3A7 gene

Cytochrome P450 3A4 and 3A7 (CYP3A4 and CYP3A7, respectively) are predominant forms in the human adult and fetal liver, respectively. 1,25-Dihydroxyvitamin D(3) (1,25(OH)(2)D(3)) is known to be a potent inducer of CYP3A4 in human colon carcinoma Caco-2 via vitamin D receptor (VDR). However, whether CYP3A7 is inducible by 1,25(OH)(2)D(3) has not yet been elucidated. In the present study, we examined the effect of 1,25(OH)(2)D(3) on CYP3A7 gene expression in Caco-2 cells, which express CYP3A4 and CYP3A7 mRNAs. 1,25(OH)(2)D(3) hardly induced the expression of CYP3A7 mRNA in contrast to the marked induction of CYP3A4 mRNA. Reporter assay using 5'-franking region CYP3A4 and CYP3A7 genes also revealed that 1,25(OH)(2)D(3) activates CYP3A4 promoter, but not CYP3A7 promoter, which has two mutations in the proximal ER6 site compared with CYP3A4 promoter. In addition, we found that the binding of VDR to the proximal ER6 in CYP3A7 gene was markedly less than that to the proximal ER6 in CYP3A4 gene using gel shift assay. Taken together, the decrease of VDR binding to the proximal ER6 caused by the mutation results in the loss of CYP3A7 gene activation by 1,25(OH)(2)D(3).     

The design of a hyperstable mutant of the Abp1p SH3 domain by sequence alignment analysis

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 degrees C and a deltaGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 degrees C with concomitant increases in deltaGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 degrees C and a deltaGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.     

A phosphoinositide-binding sequence is shared by PH domain target molecules—a model for the binding of PH domains to proteins

Pleckstrin homology (PH) domains have been proven to bind phosphoinositides (PI) and inositolphosphates (IP). On the other hand, a binding of PH domains to proteins is still a matter of debate. The goal of this work was to identify potential PH domain protein target sites and to build a model for PH domain–protein binding. A candidate sequence, called HIKE, was identified by sequence homology analysis of the proteins that are considered the strongest PH binding candidates, i.e., G_β, PKC, and Akt. HIKE contains a PI binding sequence and fulfills several criteria for a potential PH-binding site, i.e., it is present in other PH-binding candidates, lies in regulatory regions independently predicted to bind PH domains, and is conserved in 3-D structure among different molecules. These findings and the similarities with the mode of binding of PTB and PDZ domains suggest a β strand–β strand coordination model for PH–protein binding. The HIKE model predicts that membrane anchoring of PH domains and their targets could be a critical step in their interaction, which would consistently explain why PH–protein binding has only been detected in the presence of PI. Proteins 31:1–9, 1998. © 1998 Wiley-Liss, Inc.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

A sphingosine-dependent protein kinase that specifically phosphorylates 14-3-3 (SDK1) is identified as the kinase domain of PKCdelta: a preliminary note

A specific protein kinase that phosphorylates Ser60, Ser59, or Ser58 of 14-3-3beta, eta, or zeta, respectively, only in the presence of sphingosine (Sph) or N,N-dimethyl-Sph (DMS), was termed "sphingosine-dependent protein kinase-1" (SDK1) [J. Biol. Chem. 273(34) (1998) 21834]. We have now identified SDK1 as a protein having the same amino acid sequence as in the C-terminal-half kinase domain of PKCdelta, with approximately 40 kDa molecular mass, based on large-scale purification of a protein from rat liver, and partial sequence using three different combinations of LC-MS or LC-MS/MS with respective search engine. PKCdelta did not display any SDK1 activity and PKCdelta activity was inhibited by Sph and DMS. However, strong SDK1 activity, only in the presence of Sph or DMS, became detectable when PKCdelta was incubated with caspase-3, which releases the approximately 40 kDa kinase domain.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Electron transfer in the complex of membrane-bound human cytochrome P450 3A4 with the flavin domain of P450BM-3: The effect of oligomerization of the heme protein and intermittent modulation of the spin equilibrium

We studied the kinetics of NADPH-dependent reduction of human CYP3A4 incorporated into Nanodiscs (CYP3A4-ND) and proteoliposomes in order to probe the effect of P450 oligomerization on its reduction. The flavin domain of cytochrome P450-BM3 (BMR) was used as a model electron donor partner. Unlike CYP3A4 oligomers, where only 50% of the enzyme was shown to be reducible by BMR, CYP3A4-ND could be reduced almost completely. High reducibility was also observed in proteoliposomes with a high lipid-to-protein ratio (L/P = 910), where the oligomerization equilibrium is displaced towards monomers. In contrast, the reducibililty in proteoliposomes with L/P = 76 did not exceed 55 ± 6%. The effect of the surface density of CYP3A4 in proteoliposomes on the oligomerization equilibrium was confirmed with a FRET-based assay employing a cysteine-depleted mutant labeled on Cys-468 with BODIPY iodoacetamide. These results confirm a pivotal role of CYP3A4 oligomerization in its functional heterogeneity. Furthermore, the investigation of the initial phase of the kinetics of CYP3A4 reduction showed that the addition of NADPH causes a rapid low-to-high-spin transition in the CYP3A4-BMR complex, which is followed by a partial slower reversal. This observation reveals a mechanism whereby the CYP3A4 spin equilibrium is modulated by the redox state of the bound flavoprotein.     

Kinetics of electron transfer in the complex of cytochrome P450 3A4 with the flavin domain of cytochrome P450BM-3 as evidence of functional heterogeneity of the heme protein

We used a rapid scanning stop-flow technique to study the kinetics of reduction of cytochrome P450 3A4 (CYP3A4) by the flavin domain of cytochrome P450-BM3 (BMR), which was shown to form a stoichiometric complex (K_D = 0.48 μM) with CYP3A4. In the absence of substrates only about 50% of CYP3A4 was able to accept electrons from BMR. Whereas the high-spin fraction was completely reducible, the reducibility of the low-spin fraction did not exceed 42%. Among four substrates tested (testosterone, 1-pyrenebutanol, bromocriptine, or α-naphthoflavone (ANF)) only ANF is capable of increasing the reducibility of the low-spin fraction to 75%. Our results demonstrate that the pool of CYP3A4 is heterogeneous, and not all P450 is competent for electron transfer in the complex with reductase. The increase in the reducibility of the enzyme in the presence of ANF may represent an important element of the mechanism of action of this activator.     

Galpha(olf) is necessary for coupling D1 and A2a receptors to adenylyl cyclase in the striatum

In the brain, dopamine and adenosine stimulate cyclic AMP (cAMP) production through D1 and A2a receptors, respectively. Using mutant mice deficient in the olfactory isoform of the stimulatory GTP-binding protein alpha subunit, Galpha(olf), we demonstrate here the obligatory role of this protein in the adenylyl cyclase responses to dopamine and adenosine in the caudate putamen. Responses to dopamine were also dramatically decreased in the nucleus accumbens but remained unaffected in the prefrontal cortex. Moreover, in the caudate putamen of mice heterozygous for the mutation, the amounts of Galpha(olf) were half of the normal levels, and the efficacy of dopamine- and CGS 21680 A(2) agonist-stimulated cAMP production was decreased. Together, these results identify Galpha(olf) as a critical parameter in the responses to dopamine and adenosine in the basal ganglia.     

Hepatitis C virus NS5A is able to competitively displace c‐Myc from the Bin1 SH3 domain in vitro

We studied the interaction of the SH3 domain of Bin1 with a 15‐mer peptide of HCV NS5A and show its potency to competitively displace a 15‐mer human c‐Myc fragment, which is a physiological ligand of Bin1, using NMR spectroscopy. Fluorescence spectroscopy and ITC were employed to determine the affinity of Bin1 SH3 to NS5A(347–361), yielding a submicromolar affinity to NS5A. Our study compares the binding dynamics and affinities of the relevant regions for binding of c‐Myc and NS5A to Bin1 SH3. The result gives further insights into the potential role of NS5A in Bin1‐mediated apoptosis. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

Bioengineering anabolic vitamin D-25-hydroxylase activity into the human vitamin D catabolic enzyme, cytochrome P450 CYP24A1, by a V391L mutation

CYP24A1 is a mitochondrial cytochrome P450 (CYP) that catabolizes 1α,25-dihydroxyvitamin D(3) (1α,25-(OH)(2)D(3)) to different products: calcitroic acid or 1α,25-(OH)(2)D(3)-26,23-lactone via multistep pathways commencing with C24 and C23 hydroxylation, respectively. Despite the ability of CYP24A1 to catabolize a wide range of 25-hydroxylated analogs including 25-hydroxyvitamin D(3), the enzyme is unable to metabolize the synthetic prodrug, 1α-hydroxyvitamin D(3) (1α-OH-D(3)), presumably because it lacks a C25-hydroxyl. In the current study we show that a single V391L amino acid substitution in the β3a-strand of human CYP24A1 converts this enzyme from a catabolic 1α,25-(OH)(2)D(3)-24-hydroxylase into an anabolic 1α-OH-D(3)-25-hydroxylase, thereby forming the hormone, 1α,25-(OH)(2)D(3). Furthermore, because the mutant enzyme retains its basal ability to catabolize 1α,25-(OH)(2)D(3) via C24 hydroxylation, it can also make calcitroic acid. Previous work has shown that an A326G mutation is responsible for the regioselectivity differences observed between human (primarily C24-hydroxylating) and opossum (C23-hydroxylating) CYP24A1. When the V391L and A326G mutations were combined (V391L/A326G), the mutant enzyme continued to form 1α,25-(OH)(2)D(3) from 1α-OH-D(3), but this initial product was diverted via the C23 hydroxylation pathway into the 26,23-lactone. The relative position of Val-391 in the β3a-strand of a homology model and the crystal structure of rat CYP24A1 is consistent with hydrophobic contact of Val-391 and the substrate side chain near C21. We interpret that the substrate specificity of V391L-modified human CYP24A1 toward 1α-OH-D(3) is enabled by an altered contact with the substrate side chain that optimally positions C25 of the 1α-OH-D(3) above the heme for hydroxylation.