An Ensemble of Rapidly Interconverting Orientations in Electrostatic Protein-Peptide Complexes Characterized by NMR Spectroscopy

Erv1 and Cytochrome c Mediate Rapid Electron Transfer via A Collision-Type Interaction

Being essential for oxidative protein folding in the mitochondrial intermembrane space, the mitochondrial disulfide relay relies on the electron transfer (ET) from the sulfhydryl oxidase Erv1 to cytochrome c (Cc). Using solution NMR spectroscopy, we demonstrate that while the yeast Cc-Erv1 system is functionally active, no observable binding of the protein partners takes place. The transient interaction between Erv1 and Cc can be rationalized by molecular modeling, suggesting that a large surface area of Erv1 can sustain a fast ET to Cc via a collision-type mechanism, without the need for a canonical protein complex formation. We suggest that, by preventing the direct ET to molecular oxygen (O₂), the collision-type Cc-Erv1 interaction plays a role in protecting the organism against reactive oxygen species.

Investigating energy-based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure-based and physics-based atomistic force field with an efficient sampling strategy is adopted to simulate a model di-domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low-energy structures and the minimum-size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small-angle X-ray scattering data. It is illustrated that the regularizations of energy and ensemble-size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high-energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure-ensemble optimizations with a topology-based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.

Kinetic and Structural Characterization of the Effects of Membrane on the Complex of Cytochrome b 5 and Cytochrome c

Cytochrome b₅ (cytb₅) is a membrane protein vital for the regulation of cytochrome P450 (cytP450) metabolism and is capable of electron transfer to many redox partners. Here, using cyt c as a surrogate for cytP450, we report the effect of membrane on the interaction between full-length cytb₅ and cyt c for the first time. As shown through stopped-flow kinetic experiments, electron transfer capable cytb₅ - cyt c complexes were formed in the presence of bicelles and nanodiscs. Experimentally measured NMR parameters were used to map the cytb₅-cyt c binding interface. Our experimental results identify differences in the binding epitope of cytb₅ in the presence and absence of membrane. Notably, in the presence of membrane, cytb₅ only engaged cyt c at its lower and upper clefts while the membrane-free cytb₅ also uses a distal region. Using restraints generated from both cytb₅ and cyt c, a complex structure was generated and a potential electron transfer pathway was identified. These results demonstrate the importance of studying protein-protein complex formation in membrane mimetic systems. Our results also demonstrate the successful preparation of novel peptide-based lipid nanodiscs, which are detergent-free and possesses size flexibility, and their use for NMR structural studies of membrane proteins.

(S)Pinning down protein interactions by NMR

Protein molecules are highly diverse communication platforms and their interaction repertoire stretches from atoms over small molecules such as sugars and lipids to macromolecules. An important route to understanding molecular communication is to quantitatively describe their interactions. These types of analyses determine the amounts and proportions of individual constituents that participate in a reaction as well as their rates of reactions and their thermodynamics. Although many different methods are available, there is currently no single method able to quantitatively capture and describe all types of protein reactions, which can span orders of magnitudes in affinities, reaction rates, and lifetimes of states. As the more versatile technique, solution NMR spectroscopy offers a remarkable catalogue of methods that can be successfully applied to the quantitative as well as qualitative descriptions of protein interactions. In this review we provide an easy-access approach to NMR for the non-NMR specialist and describe how and when solution state NMR spectroscopy is the method of choice for addressing protein ligand interaction. We describe very briefly the theoretical background and illustrate simple protein-ligand interactions as well as typical strategies for measuring binding constants using NMR spectroscopy. Finally, this review provides examples of caveats of the method as well as the options to improve the outcome of an NMR analysis of a protein interaction reaction.

Structural and Dynamic Insights into a Glycine-Mediated Short Analogue of a Designed Peptide in Lipopolysaccharide Micelles: Correlation Between Compact Structure and Anti-Endotoxin Activity

In this study, we report an interaction study of a 13-residue analogue peptide VG13P (VARGWGRKCPLFG), derived from a designed VG16KRKP peptide (VARGWKRKCPLFGKGG), with a Lys6Gly mutation and removal of the last three residues Lys¹⁴-Gly¹⁵-Gly¹⁶, in lipopolysaccharide (LPS), a major component of the outer membrane of Gram-negative bacteria and responsible for sepsis or septic shock. VG13P displays an enhanced anti-endotoxin property as evident from significant reduction in LPS-induced TNF-α gene expression levels in a monocytic cell line, while it retains almost unchanged antimicrobial activity as its parent VG16KRKP against Gram-negative bacterial as well as fungal pathogens. In addition, in vitro LPS binding properties of VG13P in comparison to its parent VG16KRKP also remained unhindered, suggesting that the flexible C-terminal end of VG16KRKP may not play a major role in its observed antibacterial and LPS binding properties. An NMR-resolved solution structure of VG13P in LPS reveals two consecutive β-turns: one at the N-terminus, followed by another at the central region, closely resembling a rocking chair. The crucial Lys6Gly mutation along with C-terminal truncation from VG16KRKP reorients the hydrophobic hub in VG13P in a unique way so as to fold the N-terminal end back on itself, forming a turn and allowing Val1 and Ala2 to interact with Leu11 and Phe12 to bring the hydrophobic residues closer together to form a more compact hub compared to its parent. The hub is further strengthened via CH-π interaction between Gly4 and Phe12. This accounts for its improved anti-endotoxin activity as well as to its uninterrupted antimicrobial activity.

Structural and Biological Interaction of hsc-70 Protein with Phosphatidylserine in Endosomal Microautophagy

hsc-70 (HSPA8) is a cytosolic molecular chaperone, which plays a central role in cellular proteostasis, including quality control during protein refolding and regulation of protein degradation. hsc-70 is pivotal to the process of macroautophagy, chaperone-mediated autophagy, and endosomal microautophagy. The latter requires hsc-70 interaction with negatively charged phosphatidylserine (PS) at the endosomal limiting membrane. Herein, by combining plasmon resonance, NMR spectroscopy, and amino acid mutagenesis, we mapped the C terminus of the hsc-70 LID domain as the structural interface interacting with endosomal PS, and we estimated an hsc-70/PS equilibrium dissociation constant of 4.7 ± 0.1 μm. This interaction is specific and involves a total of 4-5 lysine residues. Plasmon resonance and NMR results were further experimentally validated by hsc-70 endosomal binding experiments and endosomal microautophagy assays. The discovery of this previously unknown contact surface for hsc-70 in this work elucidates the mechanism of hsc-70 PS/membrane interaction for cytosolic cargo internalization into endosomes.

Capsaicin interaction with TRPV1 channels in a lipid bilayer: molecular dynamics simulation

Transient receptor potential vanilloid subtype 1 (TRPV1) is a heat-sensitive ion channel also involved in pain sensation, and is the receptor for capsaicin, the active ingredient of hot chili peppers. The recent structures of TRPV1 revealed putative ligand density within the S1 to S4 voltage-sensor-like domain of the protein. However, questions remain regarding the dynamic role of the lipid bilayer in ligand binding to TRPV1. Molecular dynamics simulations were used to explore behavior of capsaicin in a 1-palmitoyl-2-oleoyl phosphatidylcholine bilayer and with the target S1–S4 transmembrane helices of TRPV1. Equilibrium simulations reveal a preferred interfacial localization for capsaicin. We also observed a capsaicin molecule flipping from the extracellular to the intracellular leaflet, and subsequently able to access the intracellular TRPV1 binding site. Calculation of the potential of mean force (i.e., free energy profile) of capsaicin along the bilayer normal confirms that it prefers an interfacial localization. The free energy profile indicates that there is a nontrivial but surmountable barrier to the flipping of capsaicin between opposing leaflets of the bilayer. Molecular dynamics of the S1–S4 transmembrane helices of the TRPV1 in a lipid bilayer confirm that Y511, known to be crucial to capsaicin binding, has a distribution along the bilayer normal similar to that of the aromatic group of capsaicin. Simulations were conducted of the TRPV1 S1–S4 transmembrane helices in the presence of capsaicin placed in the aqueous phase, in the lipid, or docked to the protein. No stable interaction between ligand and protein was seen for simulations initiated with capsaicin in the bilayer. However, interactions were seen between TRPV1 and capsaicin starting from the cytosolic aqueous phase, and capsaicin remained stable in the majority of simulations from the docked pose. We discuss the significance of capsaicin flipping from the extracellular to the intracellular leaflet and mechanisms of binding site access by capsaicin.

Interactions of the Chemokine CCL5/RANTES with Medium-Sized Chondroitin Sulfate Ligands

Interactions of the chemokine CCL5 (RANTES) with glycosaminoglycans (GAGs) are crucial to the CCL5-mediated inflammation process. However, structural information on interactions between CCL5 and longer GAG fragments is lacking. In this study, the interactions between oligosaccharides derived from chondroitin sulfate and a dimeric variant of CCL5 were investigated using solution nuclear magnetic resonance. The data indicate that, in addition to the BBXB motif in the 40s loop, GAGs also contact residues in the N loop in a manner similar to interactions between chemokine and the receptor N terminus, leading to possible stabilization of the dimer. Using 2,2,6,6-tetramethylpiperidin-1-yl)oxidanyl-tagged hexasaccharides, the binding orientation of the hexasaccharides was shown to be highly dependent on the sulfation pattern of the N-acetyl galactosamine groups. Finally, a model of the CCL5 dimer complexed to chondroitin sulfate hexasaccharides was constructed using paramagnetic relaxation enhancement and intra- and intermolecular nuclear Overhauser effect constraints.

Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos

Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.

Molecular landscape of the interaction between the urease accessory proteins UreE and UreG

Urease, the most efficient enzyme so far discovered, depends on the presence of nickel ions in the catalytic site for its activity. The transformation of inactive apo-urease into active holo-urease requires the insertion of two Ni(II) ions in the substrate binding site, a process that involves the interaction of four accessory proteins named UreD, UreF, UreG and UreE. This study, carried out using calorimetric and NMR-based structural analysis, is focused on the interaction between UreE and UreG from Sporosarcina pasteurii, a highly ureolytic bacterium. Isothermal calorimetric protein-protein titrations revealed the occurrence of a binding event between SpUreE and SpUreG, entailing two independent steps with positive cooperativity (Kd1=42±9μM; Kd2=1.7±0.3μM). This was interpreted as indicating the formation of the (UreE)2(UreG)2 hetero-oligomer upon binding of two UreG monomers onto the pre-formed UreE dimer. The molecular details of this interaction were elucidated using high-resolution NMR spectroscopy. The occurrence of SpUreE chemical shift perturbations upon addition of SpUreG was investigated and analyzed to establish the protein-protein interaction site. The latter appears to involve the Ni(II) binding site as well as mobile portions on the C-terminal and the N-terminal domains. Docking calculations based on the information obtained from NMR provided a structural basis for the protein-protein contact site. The high sequence and structural similarity within these protein classes suggests a generality of the interaction mode among homologous proteins. The implications of these results on the molecular details of the urease activation process are considered and analyzed.