Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

Rosetta:MSF: a modular framework for multi-state computational protein design

Computational protein design (CPD) is a powerful technique to engineer existing proteins or to design novel ones that display desired properties. Rosetta is a software suite including algorithms for computational modeling and analysis of protein structures and offers many elaborate protocols created to solve highly specific tasks of protein engineering. Most of Rosetta’s protocols optimize sequences based on a single conformation (i. e. design state). However, challenging CPD objectives like multi-specificity design or the concurrent consideration of positive and negative design goals demand the simultaneous assessment of multiple states. This is why we have developed the multi-state framework MSF that facilitates the implementation of Rosetta’s single-state protocols in a multi-state environment and made available two frequently used protocols. Utilizing MSF, we demonstrated for one of these protocols that multi-state design yields a 15% higher performance than single-state design on a ligand-binding benchmark consisting of structural conformations. With this protocol, we designed de novo nine retro-aldolases on a conformational ensemble deduced from a (βα)₈-barrel protein. All variants displayed measurable catalytic activity, testifying to a high success rate for this concept of multi-state enzyme design.

Protein engineering, i. e. the targeted modification or design of proteins has tremendous potential for medical and industrial applications. One generally applicable strategy for protein engineering is rational protein design: based on detailed knowledge of structure and function, computer programs like Rosetta propose the sequence of a protein possessing the desired properties. So far, most computer protocols have used rigid structures for design, which is a simplification because a protein’s structure is more accurately specified by a conformational ensemble. We have now implemented a framework for computational protein design that allows certain design protocols of Rosetta to make use of multiple design states like structural ensembles. An in silico assessment simulating ligand-binding design showed that this new approach generates more reliably native-like sequences than a single-state approach. As a proof-of-concept, we introduced de novo retro-aldolase activity into a scaffold protein and characterized nine variants experimentally, all of which were catalytically active.

Release of Cytochrome C from Bax Pores at the Mitochondrial Membrane

How cytochrome C is released from the mitochondria to the cytosol via Bax oligomeric pores, a process which is required for apoptosis, is still a mystery. Based on experimentally measured residue-residue distances, we recently solved the first atomic model for Bax oligomeric pores at the membranes using computational approaches. Here, we investigate the mechanism at the microsecond time- and nanometer space- scale using MD simulations. Our free energy landscape depicts a low barrier for the permeation of cytochrome C into the Bax C-terminal mouth, with the pathway proceeding to the inner cavity and exiting via the N-terminal mouth. Release is guided by organized charged/hydrophilic surfaces. The hydrophilicity and negative charge of the pore surface gradually increase along the release pathway from the pore entry to the exit opening. Rather than inert passing of the cytochrome C through a rigid pore, the flexible pore may selectively aid the cytochrome C passage. Once the Bax pore is formed in the membrane, with a low energy barrier, the release of cytochrome C may be readily achieved through energy fluctuations. Collectively, our work provides mechanistic insight in atomic detail into the release of cytochrome C through Bax oligomeric pores.

Preferential binding effects on protein structure and dynamics revealed by coarse-grained Monte Carlo simulation

The effect of preferential binding of solute molecules within an aqueous solution on the structure and dynamics of the histone H3.1 protein is examined by a coarse-grained Monte Carlo simulation. The knowledge-based residue-residue and hydropathy-index-based residue-solvent interactions are used as input to analyze a number of local and global physical quantities as a function of the residue-solvent interaction strength (f). Results from simulations that treat the aqueous solution as a homogeneous effective solvent medium are compared to when positional fluctuations of the solute molecules are explicitly considered. While the radius of gyration (Rg) of the protein exhibits a non-monotonic dependence on solvent interaction over a wide range of f within an effective medium, an abrupt collapse in Rg occurs in a narrow range of f when solute molecules rapidly bind to a preferential set of sites on the protein. The structure factor S(q) of the protein with wave vector (q) becomes oscillatory in the collapsed state, which reflects segmental correlations caused by spatial fluctuations in solute-protein binding. Spatial fluctuations in solute binding also modify the effective dimension (D) of the protein in fibrous (D ∼ 1.3), random-coil (D ∼ 1.75), and globular (D ∼ 3) conformational ensembles as the interaction strength increases, which differ from an effective medium with respect to the magnitude of D and the length scale.

Visualizing single-stranded nucleic acids in solution

Single-stranded nucleic acids (ssNAs) are ubiquitous in many key cellular functions. Their flexibility limits both the number of high-resolution structures available, leaving only a small number of protein-ssNA crystal structures, while forcing solution investigations to report ensemble averages. A description of the conformational distributions of ssNAs is essential to more fully characterize biologically relevant interactions. We combine small angle X-ray scattering (SAXS) with ensemble-optimization methods (EOM) to dynamically build and refine sets of ssNA structures. By constructing candidate chains in representative dinucleotide steps and refining the models against SAXS data, a broad array of structures can be obtained to match varying solution conditions and strand sequences. In addition to the distribution of large scale structural parameters, this approach reveals, for the first time, intricate details of the phosphate backbone and underlying strand conformations. Such information on unperturbed strands will critically inform a detailed understanding of an array of problems including protein-ssNA binding, RNA folding and the polymer nature of NAs. In addition, this scheme, which couples EOM selection with an iteratively refining pool to give confidence in the underlying structures, is likely extendable to the study of other flexible systems.

Elevated μs-ms timescale backbone dynamics in the transition state analog form of arginine kinase

Arginine kinase catalyzes reversible phosphoryl transfer between arginine and ATP. Crystal structures of arginine kinase in an open, substrate-free form and closed, transition state analog (TSA) complex indicate that the enzyme undergoes substantial domain and loop rearrangements required for substrate binding, catalysis, and product release. Nuclear magnetic resonance (NMR) has shown that substrate-free arginine kinase is rigid on the ps-ns timescale (average S²=0.84±0.08) yet quite dynamic on the µs-ms timescale (35 residues with R_ex, 12%), and that movements of the N-terminal domain and the loop comprising residues I182-G209 are rate-limiting on catalysis. Here, NMR of the TSA-bound enzyme shows similar rigidity on the ps-ns timescale (average S²=0.91±0.05) and substantially increased μs-ms timescale dynamics (77 residues; 22%). Many of the residues displaying μs-ms dynamics in NMR Carr-Purcell-Meiboom-Gill (CPMG) ¹⁵N backbone relaxation dispersion experiments of the TSA complex are also dynamic in substrate-free enzyme. However, the presence of additional dynamic residues in the TSA-bound form suggests that dynamics extend through much of the C-terminal domain, which indicates that in the closed form, a larger fraction of the protein takes part in conformational transitions to the excited state(s). Conformational exchange rate constants (k_ex) of the TSA complex are all approximately 2500s^-1, higher than any observed in the substrate-free enzyme (800-1900s^-1). Elevated μs-ms timescale protein dynamics in the TSA-bound enzyme is more consistent with recently postulated catalytic networks involving multiple interconnected states at each step of the reaction, rather than a classical single stabilized transition state.

Applications of NMR and computational methodologies to study protein dynamics

Overwhelming evidence now illustrates the defining role of atomic-scale protein flexibility in biological events such as allostery, cell signaling, and enzyme catalysis. Over the years, spin relaxation nuclear magnetic resonance (NMR) has provided significant insights on the structural motions occurring on multiple time frames over the course of a protein life span. The present review article aims to illustrate to the broader community how this technique continues to shape many areas of protein science and engineering, in addition to being an indispensable tool for studying atomic-scale motions and functional characterization. Continuing developments in underlying NMR technology alongside software and hardware developments for complementary computational approaches now enable methodologies to routinely provide spatial directionality and structural representations traditionally harder to achieve solely using NMR spectroscopy. In addition to its well-established role in structural elucidation, we present recent examples that illustrate the combined power of selective isotope labeling, relaxation dispersion experiments, chemical shift analyses, and computational approaches for the characterization of conformational sub-states in proteins and enzymes.

Entropy redistribution controls allostery in a metalloregulatory protein

Allosteric communication between two ligand-binding sites in a protein is a central aspect of biological regulation that remains mechanistically unclear. Here we show that perturbations in equilibrium picosecond-nanosecond motions impact zinc (Zn)-induced allosteric inhibition of DNA binding by the Zn efflux repressor CzrA (chromosomal zinc-regulated repressor). DNA binding leads to an unanticipated increase in methyl side-chain flexibility and thus stabilizes the complex entropically; Zn binding redistributes these motions, inhibiting formation of the DNA complex by restricting coupled fast motions and concerted slower motions. Allosterically impaired CzrA mutants are characterized by distinct nonnative fast internal dynamics "fingerprints" upon Zn binding, and DNA binding is weakly regulated. We demonstrate the predictive power of the wild-type dynamics fingerprint to identify key residues in dynamics-driven allostery. We propose that driving forces arising from dynamics can be harnessed by nature to evolve new allosteric ligand specificities in a compact molecular scaffold.

Predicting Protein Dynamics and Allostery Using Multi-Protein Atomic Distance Constraints

The related concepts of protein dynamics, conformational ensembles and allostery are often difficult to study with molecular dynamics (MD) due to the timescales involved. We present ExProSE (Exploration of Protein Structural Ensembles), a distance geometry-based method that generates an ensemble of protein structures from two input structures. ExProSE provides a unified framework for the exploration of protein structure and dynamics in a fast and accessible way. Using a dataset of apo/holo pairs it is shown that existing coarse-grained methods often cannot span large conformational changes. For T4-lysozyme, ExProSE is able to generate ensembles that are more native-like than tCONCOORD and NMSim, and comparable with targeted MD. By adding additional constraints representing potential modulators, ExProSE can predict allosteric sites. ExProSE ranks an allosteric pocket first or second for 27 out of 58 allosteric proteins, which is similar and complementary to existing methods. The ExProSE source code is freely available.

Understanding the connection between conformational changes of peptides and equilibrium thermal fluctuations

Despite the increasing evidence that conformational transitions in peptides and proteins are driven by specific vibrational energy pathways along the molecule, the current experimental techniques of analysis do as yet not allow to study these biophysical processes in terms of anisotropic energy flows. Computational methods offer a complementary approach to obtain a more detailed understanding of the vibrational and conformational dynamics of these systems. Accordingly, in this work we investigate jointly the vibrational energy distribution and the conformational dynamics of trialanine peptide in water solution at room temperature by applying the Instantaneous Normal Mode analysis to the results derived from equilibrium molecular dynamics simulations. It is shown that conformational changes in trialanine are triggered by the vibrational energy accumulated in the low-frequency modes of the molecule, and that excitation is caused exclusively by thermal fluctuations of the solute-solvent system, thus excluding the possibility of an intramolecular vibrational energy redistribution process.

Native State Volume Fluctuations in Proteins as a Mechanism for Dynamic Allostery

Allostery enables tight regulation of protein function in the cellular environment. Although existing models of allostery are firmly rooted in the current structure-function paradigm, the mechanistic basis for allostery in the absence of structural change remains unclear. In this study, we show that a typical globular protein is able to undergo significant changes in volume under native conditions while exhibiting no additional changes in protein structure. These native state volume fluctuations were found to correlate with changes in internal motions that were previously recognized as a source of allosteric entropy. This finding offers a novel mechanistic basis for allostery in the absence of canonical structural change. The unexpected observation that function can be derived from expanded, low density protein states has broad implications for our understanding of allostery and suggests that the general concept of the native state be expanded to allow for more variable physical dimensions with looser packing.

Guest-Induced Folding and Self-Assembly of Conformationally Adaptive Macrocycles into Nanosheets and Nanotubes

A conformationally adaptive macrocycle is presented, namely zorb[4]arene, which exists in multiple conformations in the uncomplexed state. The binding cavity of zorb[4]arene is concealed, either due to a collapsed conformation or by self-inclusion. The zorb[4]arene with long alkyl chains manifests itself with surprisingly low melting point and thus exist as an oil at room temperature. Binding of a guest molecule induces the folding and conformational rigidity of zorb[4]arene and leads to well-defined three-dimensional structures, which can further self-assemble into nanosheets or nanotubes upon solvent evaporation, depending on guest molecules and the conformations they can induce.

R, C. S. Expediting dynamics approach to understand the influence of 14-3-3ζ causing metastatic cancer through the interaction of YAP1 and β-TRCP

The 14-3-3ζ protein acts as a molecular switch in regulating the TGF-β pathway, which alters from a tumor suppressor in the early stage of breast cancer to a promoter of metastasis in the late stage.

Allosteric control of antibody-prion recognition through oxidation of a disulfide bond between the CH and CL chains

Molecular details of the recognition of disordered antigens by their cognate antibodies have not been studied as extensively as folded protein antigens and much is still unknown. To follow the conformational changes in the antibody and cross-talk between its subunits and with antigens, we performed molecular dynamics (MD) simulations of the complex of Fab and prion-associated peptide in the apo and bound forms. We observed that the inter-chain disulfide bond in constant domains restrains the conformational changes of Fab, especially the loops in the CH1 domain, resulting in inhibition of the cross-talk between Fab subdomains that thereby may prevent prion peptide binding. We further identified several negative and positive correlations of motions between the peptide and Fab constant domains, which suggested structural cross-talks between the constant domains and the antigen. The cross-talk was influenced by the inter-chain disulfide bond, which reduced the number of paths between them. Importantly, network analysis of the complex and its bound water molecules observed that those water molecules form an integral part of the Fab/peptide complex network and potential allosteric pathways. On-going work focuses on developing strategies aimed to incorporate these new network communications-including the associated water molecules-toward the grand challenge of antibody design.

Controlled self-sorting in self-assembled cage complexes

The construction of complex molecular machines and devices requires control of sub-component self-assembly. Varying the size or angle, and introducing secondary interactions into the backbone allows selective self-sorting between similar sub-components in metal–ligand cage complexes.

Effects of side chains of oxatub[4]arene on its conformational interconversion, molecular recognition and macroscopic self-assembly

The side-chain length of oxatub[4]arenes was shown to affect its conformational interconversion, molecular recognition and macroscopic self-assembly behavior.

Unexpectedly Enhanced Solubility of Aromatic Amino Acids and Peptides in an Aqueous Solution of Divalent Transition-Metal Cations

We experimentally observed considerable solubility of tryptophan (Trp) in a CuCl_{2} aqueous solution, which could reach 2-5 times the solubility of Trp in pure water. Theoretical studies show that the strong cation-π interaction between Cu^{2+} and the aromatic ring in Trp modifies the electronic distribution of the aromatic ring to enhance significantly the water affinity of Trp. Similar solubility enhancement has also been observed for other divalent transition-metal cations (e.g., Zn^{2+} and Ni^{2+}), another aromatic amino acid (phenylalanine), and three aromatic peptides (Trp-Phe, Phe-Phe, and Trp-Ala-Phe).

Induced Fit in Protein Multimerization: The HFBI Case

Hydrophobins, produced by filamentous fungi, are small amphipathic proteins whose biological functions rely on their unique surface-activity properties. Understanding the mechanistic details of the multimerization process is of primary importance to clarify the interfacial activity of hydrophobins. We used free energy calculations to study the role of a flexible β-hairpin in the multimerization process in hydrophobin II from Trichoderma reesei (HFBI). We characterized how the displacement of this β-hairpin controls the stability of the monomers/dimers/tetramers in solution. The regulation of the oligomerization equilibrium of HFBI will necessarily affect its interfacial properties, fundamental for its biological function and for technological applications. Moreover, we propose possible routes for the multimerization process of HFBI in solution. This is the first case where a mechanism by which a flexible loop flanking a rigid patch controls the protein-protein binding equilibrium, already known for proteins with charged binding hot-spots, is described within a hydrophobic patch.

Fungi proliferate by creating a complex hyphal network growing within a wet environment. However, for most fungi to colonize new territories, they must produce spores carried by aerial hyphae and spread them into the air. Aerial structures need to overcome the surface tension of the surrounding water in order to grow into the air. This process requires hydrophobins, a remarkable class of self-associating fungal proteins which lower the surface tension at the air/water interface by creating a thin amphipathic layer. In solution they form multimers in equilibrium with the interfacial layer. Due to their unique surface-activity properties, hydrophobins have been used for a variety of biotechnical applications. We used enhanced sampling molecular dynamics simulations methods to study the multimerization process in solution of a hydrophobin from Trichoderma reesei (HFBI). We clarified the fundamental role of a small flexible region within the HFBI monomer involved in the formation of multimers. A flexible loop flanking a rigid interaction patch is able to fine-tune the interaction energy. This mechanism, already known for charged binding patches, is described here for hydrophobic hot-spots. This result is remarkably important in order to clarify the mechanism of arranging at the interface and enhancing hydrophobin-based technological applications.

Conformational selection in amyloid-based immunotherapy: Survey of crystal structures of antibody-amyloid complexes

BACKGROUND

The dominant feature in neurodegenerative diseases is protein aggregations that lead to neuronal loss. Immunotherapies using antibodies or antibody fragments to target the aggregations are a highly perused approach. The molecular mechanisms underlying the amyloid-based immunotherapy are complex. Deciphering the properties of amyloidogenic proteins responsible for these diseases is essential to obtain insights into antibody recognition of the amyloid antigens.

SCOPE OF REVIEW

We systematically explore all available crystal structures of antibody-amyloid complexes related to neurodegenerative diseases, including antibodies that recognize the Aβ peptide, tau protein, prion protein, alpha-synuclein, huntingtin protein (mHTT), and polyglutamine.

MAJOR CONCLUSIONS

We found that antibodies mostly use the conformational selection mechanism to recognize the highly flexible amyloid antigens. In particular, solanezumab bound to Aβ12-28 tripeptide motif conformation (F19F20A21), which is shared with the Aβ42 fibril. This motif, which is trapped by the antibody, may provide the missing link in amyloid formation. Water molecules often bridge between the antibody and amyloid, contributing to the recognition.

GENERAL SIGNIFICANCE

This paper provides the structural basis for antibody recognition of amyloidogenic proteins. The analysis and discussion of known structures are expected to help in the design and optimization of antibodies in neurodegenerative diseases. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

TRAF3 signaling: Competitive binding and evolvability of adaptive viral molecular mimicry

BACKGROUND

The tumor necrosis factor receptor (TNFR) associated factor 3 (TRAF3) is a key node in innate and adaptive immune signaling pathways. TRAF3 negatively regulates the activation of the canonical and non-canonical NF-κB pathways and is one of the key proteins in antiviral immunity.

SCOPE OF REVIEW

Here we provide a structural overview of TRAF3 signaling in terms of its competitive binding and consequences to the cellular network. For completion, we also include molecular mimicry of TRAF3 physiological partners by some viral proteins.

MAJOR CONCLUSIONS

By out-competing host partners, viral proteins aim to subvert TRAF3 antiviral action. Mechanistically, dynamic, competitive binding by the organism's own proteins and same-site adaptive pathogen mimicry follow the same conformational selection principles.

GENERAL SIGNIFICANCE

Our premise is that irrespective of the eliciting event - physiological or acquired pathogenic trait - pathway activation (or suppression) may embrace similar conformational principles. However, even though here we largely focus on competitive binding at a shared site, similar to physiological signaling other pathogen subversion mechanisms can also be at play. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

Protein dynamics: Conformational footprints

The mechanisms by which proteins evolve new functions can be slow and mysterious. Comprehensive structural analysis of enzyme variants reveal how gradual enrichments of pre-existing populations with the right productive dynamics for new functions can accomplish this aim.

Exploiting computationally derived out-of-the-box protein conformations for drug design

Structural plasticity is an intrinsic property of proteins that allows each gene product to accomplish its tasks in a strictly regulated manner at a precise time and cellular location. Moreover, protein motions allow protein-ligand and protein-protein recognition. The knowledge of the conformational ensemble that a drug target populates may be crucial for the design of small molecules that can differently modulate its function. X-ray crystallography and NMR have endlessly provided snapshots of protein states. However, experimental structure determination is not always straightforward. Therefore, attempts have been made to depict protein conformational landscapes through molecular dynamics and enhanced sampling methods. Here, we review how accounting for protein dynamics through in silico generated out-of-the-box protein conformations has started to impact on drug discovery.

Conformational dynamics of cancer-associated MyD88-TIR domain mutant L252P (L265P) allosterically tilts the landscape toward homo-dimerization

MyD88 is an essential adaptor protein, which mediates the signaling of the toll-like and interleukin-1 receptors' superfamily. The MyD88 L252P (L265P) mutation has been identified in diffuse large B-cell lymphoma. The identification of this mutation has been a major advance in the diagnosis of patients with aldenstrom macroglobulinemia and related lymphoid neoplasms. Here we used computational methods to characterize the conformational effects of the mutation. Our molecular dynamics simulations revealed that the mutation allosterically quenched the global conformational dynamics of the toll/IL-1R (TIR) domain, and readjusted its salt bridges and dynamic community network. Specifically, the mutation changed the orientation and reduced the fluctuation of α-helix 3, possibly through eliminating/weakening ~8 salt bridges and enhancing the salt bridge D225-K258. Using the energy landscape of the TIR domains of MyD88, we identified two dynamic conformational basins, which correspond to the binding sites used in homo- and hetero-oligomerization, respectively. Our results indicate that the mutation stabilizes the core of the homo-dimer interface of the MyD88-TIR domain, and increases the population of homo-dimer-compatible conformational states in MyD88 family proteins. However, the dampened motion restricts its ability to heterodimerize with other TIR domains, thereby curtailing physiological signaling. In conclusion, the L252P both shifts the landscape toward homo-dimerization and restrains the dynamics of the MyD88-TIR domain, which disfavors its hetero-dimerization with other TIR domains. We further put these observations within the framework of MyD88-mediated cell signaling.

Conformational selection underpins recognition of multiple DNA sequences by proteins and consequent functional actions

Recognition of multiple functional DNA sequences by a DNA-binding protein occurs widely in nature. The physico-chemical basis of this phenomenon is not well-understood. The E. coli gal repressor, a gene regulatory protein, binds two homologous but non-identical sixteen basepair sequences in the gal operon and interacts by protein-protein interaction to regulate gene expression. The two sites have nearly equal affinities for the Gal repressor. Spectroscopic studies of the Gal repressor bound to these two different DNA sequences detected significant conformational differences between them. Comprehensive single base-substitution and binding measurements were carried out on the two sequences to understand the nature of the two protein-DNA interfaces. Magnitudes of basepair-protein interaction energy show significant variation between homologous positions of the two DNA sequences. Magnitudes of variation are such that when summed over the whole sequence they largely cancel each other out, thus producing nearly equal net affinity. Modeling suggests significant alterations in the protein-DNA interface in the two complexes, which are consistent with conformational adaptation of the protein to different DNA sequences. The functional role of the two sequences was studied by substitution of one site by the other and vice versa. In both cases, substitution reduces repression in vivo. This suggests that naturally occurring DNA sequence variations play functional roles beyond merely acting as high-affinity anchoring points. We propose that two different pre-existing conformations in the conformational ensemble of the free protein are selected by two different DNA sequences for efficient sequence read-out and the conformational difference of the bound proteins leads to different functional roles.

The principle of conformational signaling

Specific conformations of signaling proteins can serve as “signals” in signal transduction by being recognized by receptors.