Formation of transient protein complexes

Solar water splitting Pt-nanoparticle photosystem I thylakoid systems: Catalyst identification, location and oligomeric structure

Photosynthetic conversion of light energy into chemical energy occurs in sheet-like membrane-bound compartments called thylakoids and is mediated by large integral membrane protein-pigment complexes called reaction centers (RCs). Oxygenic photosynthesis of higher plants, cyanobacteria and algae requires the symbiotic linking of two RCs, photosystem II (PSII) and photosystem I (PSI), to split water and assimilate carbon dioxide. Worldwide there is a large research investment in developing RC-based hybrids that utilize the highly evolved solar energy conversion capabilities of RCs to power catalytic reactions for solar fuel generation. Of particular interest is the solar-powered production of H₂, a clean and renewable energy source that can replace carbon-based fossil fuels and help provide for ever-increasing global energy demands. Recently, we developed thylakoid membrane hybrids with abiotic catalysts and demonstrated that photosynthetic Z-scheme electron flow from the light-driven water oxidation at PSII can drive H₂ production from PSI. One of these hybrid systems was created by self-assembling Pt-nanoparticles (PtNPs) with the stromal subunits of PSI that extend beyond the membrane plane in both spinach and cyanobacterial thylakoids. Using PtNPs as site-specific probe molecules, we report the electron microscopic (EM) imaging of oligomeric structure, location and organization of PSI in thylakoid membranes and provide the first direct visualization of photosynthetic Z-scheme solar water-splitting biohybrids for clean H₂ production.

Dynamics and mechanistic interpretations of nonribosomal peptide synthetase cyclization domains

Ox-/thiazoline groups in nonribosomal peptides are formed by a variant of peptide-forming condensation domains called heterocyclization (Cy) domains and appear in a range of pharmaceutically important natural products and virulence factors. Recent cryo-EM, crystallographic, and NMR studies of Cy domains make it opportune to revisit outstanding questions regarding their molecular mechanisms. This review covers structural and dynamical findings about Cy domains that will inform future bioengineering efforts and our understanding of natural product synthesis.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

Cofactor Dynamics Couples the Protein Surface to the Heme in Cytochrome c, Facilitating Electron Transfer

Electron transport through biomolecules and in biological transport networks is of great importance to bioenergetics and biocatalysis. More generally, it is of crucial importance to understand how the pathways that connect buried metallocofactors to other cofactors, and to protein surfaces, affect the biological chemistry of metalloproteins. In terms of electron transfer (ET), the strongest coupling pathways usually comprise covalent and hydrogen bonded networks, with a limited number of through-space contacts. Herein, we set out to determine the relative roles of hydrogen bonds involved in ET via an established heme-to-surface tunneling pathway in cytochrome (cyt) c (i.e., heme-W59-D60-E61-N62). A series of cyt c variants were produced where a ruthenium tris(diimine) photooxidant was placed at position 62 via covalent modification of the N62C residue. Surprisingly, variants where the H-bonding residues W59 and D60 were replaced (i.e., W59F and D60A) showed no change in ET rate from the ferrous heme to Ru(III). In contrast, changing the composition of an alternative tunneling pathway (i.e., heme-M64-N63-C62) with the M64L substitution shows a factor of 2 decrease in the rate of heme-to-Ru ET. This pathway involves a through-space tunneling step between the heme and M64 residue, and such steps are usually disfavored. To rationalize why the heme-M64-N63-C62 is preferred, molecular dynamics (MD) simulations and Pathways analysis were employed. These simulations show that the change in heme-Ru ET rates is attributed to different conformations with compressed donor-acceptor distances, by ∼2 Å in pathway distance, in the M64-containing protein as compared to the M64L protein. The change in distance is correlated with changes in the electronic coupling that are in accord with the experimentally observed heme-Ru ET rates. Remarkably, the M64L variation at the core of the protein translates to changes in cofactor dynamics at the protein surface. The surface changes identified by MD simulations include dynamic anion-π and dipole-dipole interactions. These interactions influence the strength of tunneling pathways and ET rates by facilitating decreases in through-space tunneling distances in key coupling pathways.

Are transient protein-protein interactions more dispensable?

Protein-protein interactions (PPIs) are key drivers of cell function and evolution. While it is widely assumed that most permanent PPIs are important for cellular function, it remains unclear whether transient PPIs are equally important. Here, we estimate and compare dispensable content among transient PPIs and permanent PPIs in human. Starting with a human reference interactome mapped by experiments, we construct a human structural interactome by building three-dimensional structural models for PPIs, and then distinguish transient PPIs from permanent PPIs using several structural and biophysical properties. We map common mutations from healthy individuals and disease-causing mutations onto the structural interactome, and perform structure-based calculations of the probabilities for common mutations (assumed to be neutral) and disease mutations (assumed to be mildly deleterious) to disrupt transient PPIs and permanent PPIs. Using Bayes' theorem we estimate that a similarly small fraction (<~20%) of both transient and permanent PPIs are completely dispensable, i.e., effectively neutral upon disruption. Hence, transient and permanent interactions are subject to similarly strong selective constraints in the human interactome.

Impact of distant peptide substrate residues on enzymatic activity of SlyD

Peptidyl-prolyl isomerases (PPIases) catalyze intrinsically slow and often rate-limiting isomerization of prolyl-peptide bonds in unfolded or partially folded proteins, thereby speeding up the folding process and preventing misfolding. They often possess binding and chaperone domains in addition to the domain carrying the isomerization activity. Although generally, their substrates display no identity in their amino acid sequence upstream and downstream of the proline with 20 possibilities for each residue, PPIases are efficient enzymes. SlyD is a highly efficient PPIase consisting of an isomerase domain and an additional chaperone domain. The binding of peptide substrates to SlyD and its enzymatic activity depend to some extend on the proline-proximal residues, however, the impact of proline-distant residues has not been investigated so far. Here, we introduce a label-free NMR-based method to measure SlyD activity on different peptide substrates and analysed the data in the context of obtained binding affinities and several co-crystal structures. We show that especially charged and aromatic residues up to eight positions downstream and three positions upstream of the proline and outside the canonical region of similar conformations affect the activity and binding, although they rarely display distinct conformations in our crystal structures. We hypothesize that these positions primarily influence the association reaction. In the absence of the chaperone domain the isomerase activity strongly correlates with substrate affinity, whereas additional factors play a role in its presence. The mutual orientation of isomerase and chaperone domains depends on the presence of substrates in both binding sites, implying allosteric regulation of enzymatic activity.

Altered coordination in a blue copper protein upon association with redox partner revealed by carbon-deuterium vibrational probes

Proteins tune the reactivity of metal sites; less understood is the impact of association with a redox partner. We demonstrate the utility of carbon-deuterium labels for selective analysis of delicate...

New Insights into the Evolution of the Electron Transfer from Cytochrome f to Photosystem I in the Green and Red Branches of Photosynthetic Eukaryotes

In cyanobacteria and most green algae of the eukaryotic green lineage, the copper-protein plastocyanin (Pc) alternatively replaces the heme-protein cytochrome c6 (Cc6) as the soluble electron carrier from cytochrome f (Cf) to photosystem I (PSI). The functional and structural equivalence of 'green' Pc and Cc6 has been well established, representing an example of convergent evolution of two unrelated proteins. However, plants only produce Pc, despite having evolved from green algae. On the other hand, Cc6 is the only soluble donor available in most species of the red lineage of photosynthetic organisms, which includes, among others, red algae and diatoms. Interestingly, Pc genes have been identified in oceanic diatoms, probably acquired by horizontal gene transfer from green algae. However, the mechanisms that regulate the expression of a functional Pc in diatoms are still unclear. In the green eukaryotic lineage, the transfer of electrons from Cf to PSI has been characterized in depth. The conclusion is that in the green lineage, this process involves strong electrostatic interactions between partners, which ensure a high affinity and an efficient electron transfer (ET) at the cost of limiting the turnover of the process. In the red lineage, recent kinetic and structural modeling data suggest a different strategy, based on weaker electrostatic interactions between partners, with lower affinity and less efficient ET, but favoring instead the protein exchange and the turnover of the process. Finally, in diatoms the interaction of the acquired green-type Pc with both Cf and PSI may not yet be optimized.

Enhancing the population of the encounter complex affects protein complex formation efficiency

Optimal charge distribution is considered to be important for efficient formation of protein complexes. Electrostatic interactions guide encounter complex formation that precedes the formation of an active protein complex. However, disturbing the optimized distribution by introduction of extra charged patches on cytochrome c peroxidase does not lead to a reduction in productive encounters with its partner cytochrome c. To test whether a complex with a high population of encounter complex is more easily affected by suboptimal charge distribution, the interactions of cytochrome c mutant R13A with wild‐type cytochrome c peroxidase and a variant with an additional negative patch were studied. The complex of the peroxidase and cytochrome c R13A was reported to have an encounter state population of 80%, compared to 30% for the wild‐type cytochrome c. NMR analysis confirms the dynamic nature of the interaction and demonstrates that the mutant cytochrome c samples the introduced negative patch. Kinetic experiments show that productive complex formation is fivefold to sevenfold slower at moderate and high ionic strength values for cytochrome c R13A but the association rate is not affected by the additional negative patch on cytochrome c peroxidase, showing that the total charge on the protein surface can compensate for less optimal charge distribution. At low ionic strength (44 mm), the association with the mutant cytochrome c reaches the same high rates as found for wild‐type cytochrome c, approaching the diffusion limit.

Concerning P450 Evolution: Structural Analyses Support Bacterial Origin of Sterol 14α-Demethylases

Sterol biosynthesis, primarily associated with eukaryotic kingdoms of life, occurs as an abbreviated pathway in the bacterium Methylococcus capsulatus. Sterol 14α-demethylation is an essential step in this pathway and is catalyzed by cytochrome P450 51 (CYP51). In M. capsulatus, the enzyme consists of the P450 domain naturally fused to a ferredoxin domain at the C-terminus (CYP51fx). The structure of M. capsulatus CYP51fx was solved to 2.7 Å resolution and is the first structure of a bacterial sterol biosynthetic enzyme. The structure contained one P450 molecule per asymmetric unit with no electron density seen for ferredoxin. We connect this with the requirement of P450 substrate binding in order to activate productive ferredoxin binding. Further, the structure of the P450 domain with bound detergent (which replaced the substrate upon crystallization) was solved to 2.4 Å resolution. Comparison of these two structures to the CYP51s from human, fungi, and protozoa reveals strict conservation of the overall protein architecture. However, the structure of an "orphan" P450 from nonsterol-producing Mycobacterium tuberculosis that also has CYP51 activity reveals marked differences, suggesting that loss of function in vivo might have led to alterations in the structural constraints. Our results are consistent with the idea that eukaryotic and bacterial CYP51s evolved from a common cenancestor and that early eukaryotes may have recruited CYP51 from a bacterial source. The idea is supported by bioinformatic analysis, revealing the presence of CYP51 genes in >1,000 bacteria from nine different phyla, >50 of them being natural CYP51fx fusion proteins.

Molecular Basis of Multiple Mitochondrial Dysfunctions Syndrome 2 Caused by CYS59TYR BOLA3 Mutation

Multiple mitochondrial dysfunctions syndrome (MMDS) is a rare neurodegenerative disorder associated with mutations in genes with a vital role in the biogenesis of mitochondrial [4Fe-4S] proteins. Mutations in one of these genes encoding for BOLA3 protein lead to MMDS type 2 (MMDS2). Recently, a novel phenotype for MMDS2 with complete clinical recovery was observed in a patient containing a novel variant (c.176G > A, p.Cys59Tyr) in compound heterozygosity. In this work, we aimed to rationalize this unique phenotype observed in MMDS2. To do so, we first investigated the structural impact of the Cys59Tyr mutation on BOLA3 by NMR, and then we analyzed how the mutation affects both the formation of a hetero-complex between BOLA3 and its protein partner GLRX5 and the iron-sulfur cluster-binding properties of the hetero-complex by various spectroscopic techniques and by experimentally driven molecular docking. We show that (1) the mutation structurally perturbed the iron-sulfur cluster-binding region of BOLA3, but without abolishing [2Fe-2S]²⁺ cluster-binding on the hetero-complex; (2) tyrosine 59 did not replace cysteine 59 as iron-sulfur cluster ligand; and (3) the mutation promoted the formation of an aberrant apo C59Y BOLA3-GLRX5 complex. All these aspects allowed us to rationalize the unique phenotype observed in MMDS2 caused by Cys59Tyr mutation.

Mimicking the Nucleosomal Context in Peptide-Based Binders of a H3K36me Reader Increases Binding Affinity While Altering the Binding Mode

Targeting of proteins in the histone modification machinery has emerged as a promising new direction to fight disease. The search for compounds that inhibit proteins that readout histone modification has led to several new epigenetic drugs, mostly for proteins involved in recognition of acetylated lysines. However, this approach proved to be a challenging task for methyllysine readers, which typically feature shallow binding pockets. Moreover, reader proteins of trimethyllysine K36 on the histone H3 (H3K36me3) not only bind the methyllysine but also the nucleosomal DNA. Here, we sought to find peptide-based binders of H3K36me3 reader PSIP1, which relies on DNA interactions to tightly bind H3K36me3 modified nucleosomes. We designed several peptides that mimic the nucleosomal context of H3K36me3 recognition by including negatively charged Glu-rich regions. Using a detailed NMR analysis, we find that addition of negative charges boosts binding affinity up to 50-fold while decreasing binding to the trimethyllysine binding pocket. Since screening and selection of compounds for reader domains is typically based solely on affinity measurements due to their lack of enzymatic activity, our case highlights the need to carefully control for the binding mode, in particular for the challenging case of H3K36me3 readers.

Efficient Encounter Complex Formation and Electron Transfer to Cytochrome c Peroxidase with an Additional, Distant Electrostatic Binding Site

Electrostatic interactions can strongly increase the efficiency of protein complex formation. The charge distribution in redox proteins is often optimized to steer a redox partner to the electron transfer active binding site. To test whether the optimized distribution is more important than the strength of the electrostatic interactions, an additional negative patch was introduced on the surface of cytochrome c peroxidase, away from the stereospecific binding site, and its effect on the encounter complex as well as the rate of complex formation was determined. Monte Carlo simulations and paramagnetic relaxation enhancement NMR experiments indicate that the partner, cytochrome c, interacts with the new patch. Unexpectedly, the rate of the active complex formation was not reduced, but rather slightly increased. The findings support the idea that for efficient protein complex formation the strength of the electrostatic interaction is more critical than an optimized charge distribution.

An overview of the factors playing a role in cytochrome P450 monooxygenase and ferredoxin interactions

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins that are ubiquitously present in organisms, including non-living entities such as viruses. With the exception of self-sufficient P450s, all other P450 enzymes need electrons to perform their enzymatic activity and these electrons are supplied by P450 redox proteins. Different types of P450 redox proteins can be found in organisms and are classified into different classes. Bacterial P450s (class I) receive electrons from ferredoxins which are iron-sulfur cluster proteins. The presence of more than one copy and different types of ferredoxins within a bacterial species poses fundamental questions about the selectivity of P450s and ferredoxins in relation to each other. Apart from transferring electrons, ferredoxins have also been found to modulate P450 functions. Achieving an understanding of the interaction between ferredoxins and P450s is required to harness their biotechnological potential for designing a universal electron transfer protein. A brief overview of factors playing a role in ferredoxin and P450 interactions is presented in this review article.

Protein-Protein Binding as a Two-Step Mechanism: Preselection of Encounter Poses during the Binding of BPTI and Trypsin

Biomolecular recognition between proteins follows complex mechanisms, the understanding of which can substantially advance drug discovery efforts. Here, we track each step of the binding process in atomistic detail with molecular dynamics simulations using trypsin and its inhibitor bovine pancreatic trypsin inhibitor (BPTI) as a model system. We use umbrella sampling to cover a range of unbinding pathways. Starting from these simulations, we subsequently seed classical simulations at different stages of the process and combine them to a Markov state model. We clearly identify three kinetically separated states (an unbound state, an encounter state, and the final complex) and describe the mechanisms that dominate the binding process. From our model, we propose the following sequence of events. The initial formation of the encounter complex is driven by long-range interactions because opposite charges in trypsin and BPTI draw them together. The encounter complex features the prealigned binding partners with binding sites still partially surrounded by solvation shells. Further approaching leads to desolvation and increases the importance of van der Waals interactions. The native binding pose is adopted by maximizing short-range interactions. Thereby side-chain rearrangements ensure optimal shape complementarity. In particular, BPTI’s P1 residue adapts to the S1 pocket and prime site residues reorient to optimize interactions. After the paradigm of conformation selection, binding-competent conformations of BPTI and trypsin are already present in the apo ensembles and their probabilities increase during this proposed two-step association process. This detailed characterization of the molecular forces driving the binding process includes numerous aspects that have been discussed as central to the binding of trypsin and BPTI and protein complex formation in general. In this study, we combine all these aspects into one comprehensive model of protein recognition. We thereby contribute to enhance our general understanding of this fundamental mechanism, which is particularly critical as the development of biopharmaceuticals continuously gains significance.

The Transient Complex of Cytochrome c and Cytochrome c Peroxidase: Insights into the Encounter Complex from Multifrequency EPR and NMR Spectroscopy

We present a novel approach to study transient protein‐protein complexes with standard, 9 GHz, and high‐field, 95 GHz, electron paramagnetic resonance (EPR) and paramagnetic NMR at ambient temperatures and in solution. We apply it to the complex of yeast mitochondrial iso‐1‐cytochrome c (Cc) with cytochrome c peroxidase (CcP) with the spin label [1‐oxyl‐2,2,5,5‐tetramethyl‐Δ3‐pyrroline‐3‐methyl)‐methanethiosulfonate] attached at position 81 of Cc (SL−Cc). A dissociation constant K_D of 20±4×10⁻⁶ M (EPR and NMR) and an equal amount of stereo‐specific and encounter complex (NMR) are found. The EPR spectrum of the fully bound complex reveals that the encounter complex has a significant population (60 %) that shares important features, such as the Cc‐interaction surface, with the stereo‐specific complex.

Unlocking latent kinetic information from label-free binding

Transient affinity binding interactions are central to life, composing the fundamental elements of biological networks including cell signaling, cell metabolism and gene regulation. Assigning a defined reaction mechanism to affinity binding interactions is critical to our understanding of the associated structure-function relationship, a cornerstone of biophysical characterization. Transient kinetics are currently measured using low throughput methods such as nuclear magnetic resonance, or stop-flow spectrometry-based techniques, which are not practical in many settings. In contrast, label-free biosensors measure reaction kinetics through direct binding, and with higher throughout, impacting life sciences with thousands of publications each year. Here we have developed a methodology enabling label-free biosensors to measure transient kinetic interactions towards providing a higher throughput approach suitable for mechanistic understanding of these processes. The methodology relies on hydrodynamic dispersion modeling of a smooth analyte gradient under conditions that maintain the quasi-steady-state boundary layer assumption. A transient peptide-protein interaction of relevance to drug discovery was analyzed thermodynamically using transition state theory and numerical simulations validated the approach over a wide range of operating conditions. The data establishes the technical feasibility of this approach to transient kinetic analyses supporting further development towards higher throughput applications in life science.

Deamidation disrupts native and transient contacts to weaken the interaction between UBC13 and RING-finger E3 ligases

The deamidase OspI from enteric bacteria Shigella flexneri deamidates a glutamine residue in the host ubiquitin-conjugating enzyme UBC13 and converts it to glutamate (Q100E). Consequently, its polyubiquitination activity in complex with the RING-finger ubiquitin ligase TRAF6 and the downstream NF-κB inflammatory response is silenced. The precise role of deamidation in silencing the UBC13/TRAF6 complex is unknown. We report that deamidation inhibits the interaction between UBC13 and TRAF6 RING-domain (TRAF6^RING) by perturbing both the native and transient interactions. Deamidation creates a new intramolecular salt-bridge in UBC13 that competes with a critical intermolecular salt-bridge at the native UBC13/TRAF6^RING interface. Moreover, the salt-bridge competition prevents transient interactions necessary to form a typical UBC13/RING complex. Repulsion between E100 and the negatively charged surface of RING also prevents transient interactions in the UBC13/RING complex. Our findings highlight a mechanism wherein a post-translational modification perturbs the conformation and stability of transient complexes to inhibit protein-protein association.

Shigella is a highly infectious group of bacteria that attack the human digestive tract, causing severe and often deadly diarrhoea, especially in children. There is currently no vaccine to protect against the disease, and some strains are also now resistant to antibiotics. People get infected by eating or drinking contaminated foods and water. After passing through the stomach, Shigella invades and then multiplies in the lining of the intestine, eventually causing tissue damage and irritation.

During this process, Shigella ‘hides’ from its host’s immune system by blocking how intestinal cells respond to infection. Normally, infected cells send out chemical signals that act like a call for help, attracting specialised immune cells to clear the infection. In intestinal cells, two proteins called UBC13 and TRAF6 work together to switch on this response. Specifically, TRAF6 needs to bind to UBC13 for the switch to turn on.

Like many proteins, UBC13 is formed of thousands of atoms; some of these are organized in ‘functional groups’, a collection of atoms joined in a specific manner and with special chemical properties. During Shigella infection, the bacteria produce an enzyme that changes a single functional group (an amino group) at a specific location within UBC13 for a different one (an hydroxyl group).

Previous research showed that this could stop the immune response in intestinal cells, but the mechanism remained unknown. Mohanty et al. therefore set out to determine exactly how a change of so few atoms could have such a dramatic effect.

Biochemical studies using purified proteins revealed that Shigella’s alteration to UBC13 did not change its overall structure. However, the altered protein could no longer bind to its partner TRAF6. Theoretical analysis and computer simulations revealed that the normal binding process relies on a positively charged amino acid (one of the protein’s building blocks) in UBC13 and a negatively charged one in TRAF6 being attracted to each other. Shigella’s substitution, however, introduces a second negatively charged amino acid in UBC13. This ‘steals’ the positively charged amino acid that would normally interact with TRAF6: the electrical attraction between the two proteins is disrupted, and this stops them from binding.

The work by Mohanty et al. reveals the exact mechanism Shigella uses to dampen its host’s immune response during infection. In the future, this knowledge could be used to develop more effective drugs that would help control outbreaks of diarrhoea.

Model of a Kinetically Driven Crosstalk between Paralogous Protein Encounter Complexes

Proteins interact with one another across a broad spectrum of affinities. Our understanding of the low end of this spectrum, as characterized by millimolar dissociation constants, relies on a handful of cases in which weak encounters have experimentally been identified. These weak interactions away from the specific target binding site can lead toward a higher-affinity complex. Recently, we detected weak encounters between two paralogous phosphotransferase pathways of Escherichia coli, which regulate various metabolic processes and stress responses. In addition to encounters that are known to occur between cognate proteins, i.e., those that can exchange phosphate groups with each other, surprisingly, encounters involving noncognates were also observed. It is not clear whether these "futile" encounters have a cooperative or competitive role. Using agent-based simulations, we find that the encounter complexes can be cooperative or competitive so as to increase or lower the effective binding affinity of the specific complex under different circumstances. This finding invites further questions into how organisms might exploit such low affinities to connect their signaling components.

Heterogeneous and Highly Dynamic Interface in Plastocyanin-Cytochrome f Complex Revealed by Site-Specific 2D-IR Spectroscopy

Transient protein complexes are crucial for sustaining dynamic cellular processes. The complexes of electron-transfer proteins are a notable example, such as those formed by plastocyanin (Pc) and cytochrome f (cyt f) in the photosynthetic apparatus. The dynamic and heterogeneous nature of these complexes, however, makes their study challenging. To better elucidate the complex of Nostoc Pc and cyt f, 2D-IR spectroscopy coupled to site-specific labeling with cyanophenylalanine infrared (IR) probes was employed to characterize how the local environments at sites along the surface of Pc were impacted by cyt f binding. The results indicate that Pc most substantially engages with cyt f via the hydrophobic patch around the copper redox site. Complexation with cyt f led to an increase in inhomogeneous broadening of the probe absorptions, reflective of increased heterogeneity of interactions with their environment. Notably, most of the underlying states interconverted very rapidly (1 to 2 ps), suggesting a complex with a highly mobile interface. The data support a model of the complex consisting of a large population of an encounter complex. Additionally, the study demonstrates the application of 2D-IR spectroscopy with site-specifically introduced probes to reveal new quantitative insight about dynamic biochemical systems.

Assessing protein conformational landscapes: integration of DEER data in Maximum Occurrence analysis

The properties of the conformational landscape of a biomolecule are of capital importance to understand its function. It is widely accepted that a statistical ensemble is far more representative than a single structure, especially for proteins with disordered regions. While experimental data provide the most important handle on the conformational variability that the system is experiencing, they usually report on either time or ensemble averages. Since the available conformations largely outnumber the (independent) available experimental data, the latter can be equally well reproduced by a variety of ensembles. We have proposed the Maximum Occurrence (MaxOcc) approach to provide an upper bound of the statistical weight of each conformation. This method is expected to converge towards the true statistical weights by increasing the number of independent experimental datasets. In this paper we explore the ability of DEER (Double Electron Electron Resonance) data, which report on the distance distribution between two spin labels attached to a biomolecule, to restrain the MaxOcc values and its complementarity to previously introduced experimental techniques such as NMR and Small-Angle X-ray Scattering. We here present the case of Ca2+ bound calmodulin (CaM) as a test case and show that DEER data impose a sizeable reduction of the conformational space described by high MaxOcc conformations.

Solution structure of the cytochrome P450 reductase-cytochrome c complex determined by neutron scattering

Electron transfer in all living organisms critically relies on formation of complexes between the proteins involved. The function of these complexes requires specificity of the interaction to allow for selective electron transfer but also a fast turnover of the complex, and they are therefore often transient in nature, making them challenging to study. Here, using small-angle neutron scattering with contrast matching with deuterated protein, we report the solution structure of the electron transfer complex between cytochrome P450 reductase (CPR) and its electron transfer partner cytochrome c This is the first reported solution structure of a complex between CPR and an electron transfer partner. The structure shows that the interprotein interface includes residues from both the FMN- and FAD-binding domains of CPR. In addition, the FMN is close to the heme of cytochrome c but distant from the FAD, indicating that domain movement is required between the electron transfer steps in the catalytic cycle of CPR. In summary, our results reveal key details of the CPR catalytic mechanism, including interactions of two domains of the reductase with cytochrome c and motions of these domains relative to one another. These findings shed light on interprotein electron transfer in this system and illustrate a powerful approach for studying solution structures of protein-protein complexes.

A Structural Model of a P450-Ferredoxin Complex from Orientation-Selective Double Electron-Electron Resonance Spectroscopy

Cytochrome P450 (CYP) monooxygenases catalyze the oxidation of chemically inert carbon-hydrogen bonds in diverse endogenous and exogenous organic compounds by atmospheric oxygen. This C-H bond oxy-functionalization activity has huge potential in biotechnological applications. Class I CYPs receive the two electrons required for oxygen activation from NAD(P)H via a ferredoxin reductase and ferredoxin. The interaction of Class I CYPs with their cognate ferredoxin is specific. In order to reconstitute the activity of diverse CYPs, structural characterization of CYP-ferredoxin complexes is necessary, but little structural information is available. Here we report a structural model of such a complex (CYP199A2-HaPux) in frozen solution derived from distance and orientation restraints gathered by the EPR technique of orientation-selective double electron-electron resonance (os-DEER). The long-lived oscillations in the os-DEER spectra were well modeled by a single orientation of the CYP199A2-HaPux complex. The structure is different from the two known Class I CYP-Fdx structures: CYP11A1-Adx and CYP101A1-Pdx. At the protein interface, HaPux residues in the [Fe₂S₂] cluster-binding loop and the α3 helix and the C-terminus residue interact with CYP199A2 residues in the proximal loop and the C helix. These residue contacts are consistent with biochemical data on CYP199A2-ferredoxin binding and electron transfer. Electron-tunneling calculations indicate an efficient electron-transfer pathway from the [Fe₂S₂] cluster to the heme. This new structural model of a CYP-Fdx complex provides the basis for tailoring CYP enzymes for which the cognate ferredoxin is not known, to accept electrons from HaPux and display monooxygenase activity.

Interaction and electron transfer between ferredoxin-NADP+ oxidoreductase and its partners: structural, functional, and physiological implications

Ferredoxin-NADP⁺ reductase (FNR) catalyzes the last step of linear electron transfer in photosynthetic light reactions. The FAD cofactor of FNR accepts two electrons from two independent reduced ferredoxin molecules (Fd) in two sequential steps, first producing neutral semiquinone and then the fully anionic reduced, or hydroquinone, form of the enzyme (FNR_hq). FNR_hq transfers then both electrons in a single hydride transfer step to NADP⁺. We are presenting the recent progress in studies focusing on Fd:FNR interaction and subsequent electron transfer processes as well as on interaction of FNR with NADP⁺/H followed by hydride transfer, both from the structural and functional point of views. We also present the current knowledge about the physiological role(s) of various FNR isoforms present in the chloroplasts of higher plants and the functional impact of subchloroplastic location of FNR. Moreover, open questions and current challenges about the structure, function, and physiology of FNR are discussed.

Insights into the Dynamics and Dissociation Mechanism of a Protein Redox Complex Using Molecular Dynamics

Leishmania major peroxidase (LmP) is structurally and functionally similar to the well-studied yeast Cytochrome c peroxidase (CCP). A recent Brownian dynamics study showed that L. major Cytochrome c (LmCytc) associates with LmP by forming an initial complex with the N-terminal helix A of LmP, followed by a movement toward the electron transfer (ET) site observed in the LmP-LmCytc crystal structure. Critical to forming the active electron transfer complex is an intermolecular Arg-Asp ion pair at the center of the interface. If the dissociation reaction is effectively the reverse of the association reaction, then rupture of the Asp-Arg ion pair should be followed by movement of LmCytc back toward LmP helix A. To test this possibility, we have performed multiple molecular dynamics (MD) simulations of the LmP-LmCytc complex. In five separate simulations, LmCytc is observed to indeed move toward helix A, and in two of the simulations, the Asp-Arg ion pair breaks, which frees LmCytc to fully associate with the LmP helix A secondary binding site. These results support the "bind and crawl" or "velcro" mechanism of association, wherein LmCytc forms a nonspecific electrostatic complex with LmP helix A, followed by a "crawl" toward the ET-active site, where the Asp-Arg ion pair holds the LmCytc in position for rapid ET. These simulations also point to Tyr134_LmP as being important in the association/dissociation reactions. Experimentally mutating Tyr134 to Phe was found to decrease K_m by 3.6-fold, which is consistent with its predicted role in complex formation by MD simulations.

Identification of productive and futile encounters in an electron transfer protein complex

Well-defined, stereospecific states in protein complexes are often in exchange with an ensemble of more dynamic orientations: the encounter states. The structure of the stereospecific complex between cytochrome P450cam and putidaredoxin was solved recently by X-ray diffraction as well as paramagnetic NMR spectroscopy. Other than the stereospecific complex, the NMR data clearly show the presence of additional states in the complex in solution. In these encounter states, populated for a small percentage of the time, putidaredoxin assumes multiple orientations and samples a large part of the surface of cytochrome P450cam. To characterize the nature of the encounter states, an extensive paramagnetic NMR dataset has been analyzed using the Maximum Occurrence of Regions methodology. The analysis reveals the location and maximal spatial extent of the additional states needed to fully explain the NMR data. Under the assumption of sparsity of the size of the conformational ensemble, several minor states can be located quite precisely. The distribution of these minor states correlates with the electrostatic potential map around cytochrome P450cam. Whereas some minor states are on isolated positively charged patches, others are connected to the stereospecific site via positively charged paths. The existence of electrostatically favorable pathways between the stereospecific interaction site and the different minor states or lack thereof suggests a means to discriminate between productive and futile encounter states.

Stability of a Transient Protein Complex in a Charged Aqueous Droplet with Variable pH

Electrospray ionization mass spectrometry (ESI-MS) has the potential to become a high-throughput robust experimental method for the detection of protein-protein equilibrium constants. Poorly understood processes that affect the stability of weak noncovalent protein complexes in the intervening droplet environment are a significant factor that precludes the advancement of the method. We use molecular dynamics to study the stability of a ubiquitin and ubiquitin-associated domain complex (RCSB PDB code 2MRO ) in an aqueous droplet with changing size and charge concentration. We present evidence that a weak protein complex changes conformation and may dissociate in shrinking droplets. Then, the droplets containing these dissociated proteins divide. Our findings suggest that in some cases ESI-MS does not measure the correct association constants. The study intends to stimulate research for systematic development of experimental protocols that stabilize weakly bound protein interfaces in droplets.

Mapping Lipid Bilayer Recognition Sites of Metalloproteinases and Other Prospective Peripheral Membrane Proteins

Peripheral binding of proteins to lipid bilayers is critical not only in intracellular signaling but also in metalloproteinase shedding of signaling proteins from cell surfaces. Assessment of how proteins recognize fluid bilayers peripherally using crystallography or structure-based predictions has been important but incomplete. Assay of dynamic protein-bilayer interactions in solution has become feasible and reliable using paramagnetic NMR and site-directed fluor labeling. Details of preparations and assay protocols for these spectroscopic measurements of bilayer proximity or contact, respectively, are described.

Calmodulin Gates Aquaporin 0 Permeability through a Positively Charged Cytoplasmic Loop

Aquaporin 0 (AQP0), the major intrinsic protein of the eye lens, plays a vital role in maintaining lens clarity by facilitating the transport of water across lens fiber cell membranes. AQP0 reduces its osmotic water permeability constant (P_f) in response to increases in the external calcium concentration, an effect that is mediated by an interaction with the calcium-binding messenger protein, calmodulin (CaM), and phosphorylation of the CaM-binding site abolishes calcium sensitivity. Despite recent structural characterization of the AQP0-CaM complex, the mechanism by which CaM modulates AQP0 remains poorly understood. By combining atomistic molecular dynamics simulations and oocyte permeability assays, we conclude that serine phosphorylation of AQP0 does not inhibit CaM binding to the whole AQP0 protein. Instead, AQP0 phosphorylation alters calcium sensitivity by modifying the AQP0-CaM interaction interface, particularly at an arginine-rich loop that connects the fourth and fifth transmembrane helices. This previously unexplored loop, which sits outside of the canonical CaM-binding site on the AQP0 cytosolic face, mechanically couples CaM to the pore-gating residues of the second constriction site. We show that this allosteric loop is vital for CaM regulation of the channels, facilitating cooperativity between adjacent subunits and regulating factors such as serine phosphorylation. Similar allosteric interactions may also mediate CaM modulation of the properties of other CaM-regulated proteins.

A Cytochrome P450 3A4 Biosensor Based on Generation 4.0 PAMAM Dendrimers for the Detection of Caffeine

Cytochromes P450 (CYP, P450) are a large family of heme-active-site proteins involved in many catalytic processes, including steroidogenesis. In humans, four primary enzymes are involved in the metabolism of almost all xenobiotics. Among these enzymes, CYP3A4 is responsible for the inactivation of the majority of used drugs which makes this enzyme an interesting target for many fields of research, especially pharmaceutical research. Since the late 1970s, attempts have been made to construct and develop electrochemical sensors for the determination of substrates. This paper is concerned with the establishment of such a CYP3A4-containing biosensor. The sensor was constructed by adsorption of alternating layers of sub-nanometer gold particle-modified PAMAM (poly-amido-amine) dendrimers of generation 4.0, along with the enzyme by a layer-by-layer assembly technique. Atomic force microscopy (AFM), quartz crystal microbalance (QCM), and Fourier-transformed infrared spectroscopy (FTIR) were employed to elucidate the sensor assembly. Additionally, the biosensor was tested by cyclic voltammetry using caffeine as a substrate.

Weak self-association of cytochrome c peroxidase molecules observed by paramagnetic NMR

There is growing experimental evidence that many proteins exhibit a tendency for (ultra)weak homo- or hetero- oligomerization interactions. With the development of paramagnetic relaxation enhancement NMR spectroscopy it has become possible to characterize weak complexes experimentally and even detect complexes with affinities in the 1-25 mM range. We present evidence for a weak complex between cytochrome c peroxidase (CcP) molecules. In a previous study, we attached nitroxide based spin labels at three positions on CcP with the intent of observing intramolecular PRE effects. However, several intermolecular PRE effects were also observed suggesting a weak self-association between CcP molecules. The CcP-CcP complex was characterized using paramagnetic NMR and protein docking. The interaction occurs between the surface that is also part of the stereo-specific binding site for its physiological partner, cytochrome c (Cc), and several small, positively charged patches on the "back" of CcP. The CcP-CcP complex is not a stereo-specific complex. It is a dynamic ensemble of orientations, characteristic of an encounter state. The contact areas resemble those observed for CcP molecules in crystals. The CcP-CcP complex formation competes with that of the CcP-Cc complex. However, the affinity for Cc is much larger and thus it is expected that, under physiological conditions, auto-inhibition will be limited. A weak self-association between cytochrome c peroxidase molecules was characterized using paramagnetic NMR.

Path to Collagenolysis: COLLAGEN V TRIPLE-HELIX MODEL BOUND PRODUCTIVELY AND IN ENCOUNTERS BY MATRIX METALLOPROTEINASE-12

Collagenolysis is essential in extracellular matrix homeostasis, but its structural basis has long been shrouded in mystery. We have developed a novel docking strategy guided by paramagnetic NMR that positions a triple-helical collagen V mimic (synthesized with nitroxide spin labels) in the active site of the catalytic domain of matrix metalloproteinase-12 (MMP-12 or macrophage metalloelastase) primed for catalysis. The collagenolytically productive complex forms by utilizing seven distinct subsites that traverse the entire length of the active site. These subsites bury ∼1,080 Å(2)of surface area, over half of which is contributed by the trailing strand of the synthetic collagen V mimic, which also appears to ligate the catalytic zinc through the glycine carbonyl oxygen of its scissile G∼VV triplet. Notably, the middle strand also occupies the full length of the active site where it contributes extensive interfacial contacts with five subsites. This work identifies, for the first time, the productive and specific interactions of a collagen triple helix with an MMP catalytic site. The results uniquely demonstrate that the active site of the MMPs is wide enough to accommodate two strands from collagen triple helices. Paramagnetic relaxation enhancements also reveal an extensive array of encounter complexes that form over a large part of the catalytic domain. These transient complexes could possibly facilitate the formation of collagenolytically active complexes via directional Brownian tumbling.

Visualizing the Ensemble Structures of Protein Complexes Using Chemical Cross-Linking Coupled with Mass Spectrometry

Graphical Abstract

Abstract

Chemical cross-linking coupled with mass spectrometry (CXMS) identifies protein residues that are close in space, and has been increasingly used for modeling the structures of protein complexes. Here we show that a single structure is usually sufficient to account for the intermolecular cross-links identified for a stable complex with sub-µmol/L binding affinity. In contrast, we show that the distance between two cross-linked residues in the different subunits of a transient or fleeting complex may exceed the maximum length of the cross-linker used, and the cross-links cannot be fully accounted for with a unique complex structure. We further show that the seemingly incompatible cross-links identified with high confidence arise from alternative modes of protein-protein interactions. By converting the intermolecular cross-links to ambiguous distance restraints, we established a rigid-body simulated annealing refinement protocol to seek the minimum set of conformers collectively satisfying the CXMS data. Hence we demonstrate that CXMS allows the depiction of the ensemble structures of protein complexes and elucidates the interaction dynamics for transient and fleeting complexes.

Electronic supplementary material

The online version of this article (doi:10.1007/s41048-015-0015-y) contains supplementary material, which is available to authorized users.

Design and engineering of a man-made diffusive electron-transport protein

Maquettes are man-made cofactor-binding oxidoreductases designed from first principles with minimal reference to natural protein sequences. Here we focus on water-soluble maquettes designed and engineered to perform diffusive electron transport of the kind typically carried out by cytochromes, ferredoxins and flavodoxins and other small proteins in photosynthetic and respiratory energy conversion and oxido-reductive metabolism. Our designs were tested by analysis of electron transfer between heme maquettes and the well-known natural electron transporter, cytochrome c. Electron-transfer kinetics were measured from seconds to milliseconds by stopped-flow, while sub-millisecond resolution was achieved through laser photolysis of the carbon monoxide maquette heme complex. These measurements demonstrate electron transfer from the maquette to cytochrome c, reproducing the timescales and charge complementarity modulation observed in natural systems. The ionic strength dependence of inter-protein electron transfer from 9.7×10(6) M(-1) s(-1) to 1.2×10(9) M(-1) s(-1) follows a simple Debye-Hückel model for attraction between +8 net charged oxidized cytochrome c and -19 net charged heme maquette, with no indication of significant protein dipole moment steering. Successfully recreating essential components of energy conversion and downstream metabolism in man-made proteins holds promise for in vivo clinical intervention and for the production of fuel or other industrial products. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Transient protein-protein interactions visualized by solution NMR

Proteins interact with each other to establish their identities in cell. The affinities for the interactions span more than ten orders of magnitude, and KD values in μM-mM regimen are considered transient and are important in cell signaling. Solution NMR including diamagnetic and paramagnetic techniques has enabled atomic-resolution depictions of transient protein-protein interactions. Diamagnetic NMR allows characterization of protein complexes with KD values up to several mM, whereas ultraweak and fleeting complexes can be modeled with the use of paramagnetic NMR especially paramagnetic relaxation enhancement (PRE). When tackling ever-larger protein complexes, PRE can be particularly useful in providing long-range intermolecular distance restraints. As NMR measurements are averaged over the ensemble of complex structures, structural information for dynamic protein-protein interactions besides the stereospecific one can often be extracted. Herein the protein interaction dynamics are exemplified by encounter complexes, alternative binding modes, and coupled binding/folding of intrinsically disordered proteins. Further integration of NMR with other biophysical techniques should allow better visualization of transient protein-protein interactions. In particular, single-molecule data may facilitate the interpretation of ensemble-averaged NMR data. Though same structures of proteins and protein complexes were found in cell as in diluted solution, we anticipate that the dynamics of transient protein protein-protein interactions be different, which awaits awaits exploration by NMR. This article is part of a Special Issue entitled: Physiological Enzymology and Protein Functions. This article is part of a Special Issue entitled: Physiological Enzymology and Protein Functions.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.

Synergic role of nucleophosmin three-helix bundle and a flanking unstructured tail in the interaction with G-quadruplex DNA

Nucleophosmin (NPM1) is a nucleocytoplasmic shuttling protein, mainly localized at nucleoli, that plays a number of functions in ribosome biogenesis and export, cell cycle control, and response to stress stimuli. NPM1 is the most frequently mutated gene in acute myeloid leukemia; mutations map to the C-terminal domain of the protein and cause its denaturation and aberrant cytoplasmic translocation. NPM1 C-terminal domain binds G-quadruplex regions at ribosomal DNA and at gene promoters, including the well characterized sequence from the nuclease-hypersensitive element III region of the c-MYC promoter. These activities are lost by the leukemic variant. Here we analyze the NPM1/G-quadruplex interaction, focusing on residues belonging to both the NPM1 terminal three-helix bundle and a lysine-rich unstructured tail, which has been shown to be necessary for high affinity recognition. We performed extended site-directed mutagenesis and measured binding rate constants through surface plasmon resonance analysis. These data, supported by molecular dynamics simulations, suggest that the unstructured tail plays a double role in the reaction mechanism. On the one hand, it facilitates the formation of an encounter complex through long range electrostatic interactions; on the other hand, it directly contacts the G-quadruplex scaffold through multiple and transient electrostatic interactions, significantly enlarging the contact surface.

The cytochrome c peroxidase and cytochrome c encounter complex: the other side of the story

Formation of an encounter complex is important for efficient protein complex formation. The encounter state consists of an ensemble of orientations of two proteins in the complex. Experimental description of such ensembles inherently suffers from insufficient data availability. We have measured paramagnetic relaxation enhancements (PRE) on cytochrome c peroxidase (CcP) caused by its partner cytochrome c (Cc) carrying a spin label. The data complement earlier PRE data of spin labelled CcP, identifying several new interactions. This work demonstrates the need of obtaining as many independent data sets as possible to achieve the most accurate description of an encounter complex.