Protein-protein and ligand-protein interactions studied by analytical ultracentrifugation

Label-free methods for optical in vitro characterization of protein-protein interactions

Protein-protein interactions are involved in the regulation and function of the majority of cellular processes. As a result, much effort has been aimed at the development of methodologies capable of quantifying protein-protein interactions, with label-free methods being of particular interest due to the associated simplified workflows and minimisation of label-induced perturbations. Here, we review recent advances in optical technologies providing label-free in vitro measurements of affinities and kinetics. We provide an overview and comparison of existing techniques and their principles, discussing advantages, limitations, and recent applications.

Computationally predicting binding affinity in protein-ligand complexes: free energy-based simulations and machine learning-based scoring functions

Accurately predicting protein-ligand binding affinities can substantially facilitate the drug discovery process, but it remains as a difficult problem. To tackle the challenge, many computational methods have been proposed. Among these methods, free energy-based simulations and machine learning-based scoring functions can potentially provide accurate predictions. In this paper, we review these two classes of methods, following a number of thermodynamic cycles for the free energy-based simulations and a feature-representation taxonomy for the machine learning-based scoring functions. More recent deep learning-based predictions, where a hierarchy of feature representations are generally extracted, are also reviewed. Strengths and weaknesses of the two classes of methods, coupled with future directions for improvements, are comparatively discussed.

Resuspending samples in analytical ultracentrifugation

In analytical ultracentrifugation it is often very useful to resuspend samples in situ after sedimentation experiments for further investigation. This can be achieved by manually subjecting the entire sample cell assembly to gentle motion that causes the air bubble in the sample compartment to repeatedly move through the solution and thereby cause convection. Here we describe a cell mixing device that can accomplish the same through axial rotation and slow rocking motion. This cell mixer is low-cost, open-source, and can be easily assembled from readily available components. It can efficiently mix multiple sample cells side-by-side and may be used with various centerpiece designs.

Elucidating Complicated Assembling Systems in Biology Using Size-and-Shape Analysis of Sedimentation Velocity Data

Sedimentation velocity analytical ultracentrifugation (SV-AUC) has seen a resurgence in popularity as a technique for characterizing macromolecules and complexes in solution. SV-AUC is a particularly powerful tool for studying protein conformation, complex stoichiometry, and interacting systems in general. Deconvoluting velocity data to determine a sedimentation coefficient distribution c(s) allows for the study of either individual proteins or multicomponent mixtures. The standard c(s) approach estimates molar masses of the sedimenting species based on determination of the frictional ratio (f/f0) from boundary shapes. The frictional ratio in this case is a weight-averaged parameter, which can lead to distortion of mass estimates and loss of information when attempting to analyze mixtures of macromolecules with different shapes. A two-dimensional extension of the c(s) analysis approach provides size-and-shape distributions that describe the data in terms of a sedimentation coefficient and frictional ratio grid. This allows for better resolution of species with very distinct shapes that may co-sediment and provides better molar mass determinations for multicomponent mixtures. An example case is illustrated using globular and nonglobular proteins of different masses with nearly identical sedimentation coefficients that could only be resolved using the size-and-shape distribution. Other applications of this analytical approach to complex biological systems are presented, focusing on proteins involved in the innate immune response to cytosolic microbial DNA.

A histogram approach to the quality of fit in sedimentation velocity analyses

The quality of fit of sedimentation velocity data is critical to judge the veracity of the sedimentation model and accuracy of the derived macromolecular parameters. Absolute statistical measures are usually complicated by the presence of characteristic systematic errors and run-to-run variation in the stochastic noise of data acquisition. We present a new graphical approach to visualize systematic deviations between data and model in the form of a histogram of residuals. In comparison with the ideally expected Gaussian distribution, it can provide a robust measure of fit quality and be used to flag poor models.

Structural and hydrodynamic analysis of a novel drug delivery vector: ELP[V5G3A2-150]

The therapeutic potential of elastin-like polypeptide (ELP) conjugated to therapeutic compounds is currently being investigated as an approach to target drugs to solid tumors. ELPs are hydrophobic polymers that are soluble at low temperatures and cooperatively aggregate above a transition temperature (TT), allowing for thermal targeting of covalently attached drugs. They have been shown to cooperatively transition from a disordered structure to a repeating type II β-turn structure, forming a β-spiral above the TT. Here we present biophysical measurements of the structural, thermodynamic, and hydrodynamic properties of a specific ELP being investigated for drug delivery, ELP[V5G3A2-150]. We examine the biophysical properties below and above the TT to understand and predict the therapeutic potential of ELP-drug conjugates. We observed that below the TT, ELP[V5G3A2-150] is soluble, with an extended conformation consisting of both random coil and heterogeneous β structures. Sedimentation velocity experiments indicate that ELP[V5G3A2-150] undergoes weak self-association with increasing temperature, and above the TT the hydrophobic effect drives aggregation entropically. These experiments also reveal a previously unreported temperature-dependent critical concentration (Cc) that resembles a solubility constant. Labeling ELP[V5G3A2-150] with fluorescein lowers the TT by 3.5°C at 20 μM, whereas ELP[V5G3A2-150] dissolution in physiological media (fetal bovine serum) increases the TT by ∼2.2°C.

Structural interconversions modulate activity of Escherichia coli ribonucleotide reductase

Essential for DNA biosynthesis and repair, ribonucleotide reductases (RNRs) convert ribonucleotides to deoxyribonucleotides via radical-based chemistry. Although long known that allosteric regulation of RNR activity is vital for cell health, the molecular basis of this regulation has been enigmatic, largely due to a lack of structural information about how the catalytic subunit (α(2)) and the radical-generation subunit (β(2)) interact. Here we present the first structure of a complex between α(2) and β(2) subunits for the prototypic RNR from Escherichia coli. Using four techniques (small-angle X-ray scattering, X-ray crystallography, electron microscopy, and analytical ultracentrifugation), we describe an unprecedented α(4)β(4) ring-like structure in the presence of the negative activity effector dATP and provide structural support for an active α(2)β(2) configuration. We demonstrate that, under physiological conditions, E. coli RNR exists as a mixture of transient α(2)β(2) and α(4)β(4) species whose distributions are modulated by allosteric effectors. We further show that this interconversion between α(2)β(2) and α(4)β(4) entails dramatic subunit rearrangements, providing a stunning molecular explanation for the allosteric regulation of RNR activity in E. coli.

A novel tetrameric PilZ domain structure from xanthomonads

PilZ domain is one of the key receptors for the newly discovered secondary messenger molecule cyclic di-GMP (c-di-GMP). To date, several monomeric PilZ domain proteins have been identified. Some exhibit strong c-di-GMP binding activity, while others have barely detectable c-di-GMP binding activity and require an accessory protein such as FimX to indirectly respond to the c-di-GMP signal. We now report a novel tetrameric PilZ domain structure of XCC6012 from the plant pathogen Xanthomonas campestris pv. campestris (Xcc). It is one of the four PilZ domain proteins essential for Xcc pathogenicity. Although the monomer adopts a structure similar to those of the PilZ domains with very weak c-di-GMP binding activity, it is nevertheless interrupted in the middle by two extra long helices. Four XCC6012 proteins are thus self-assembled into a tetramer via the extra heptad repeat α3 helices to form a parallel four-stranded coiled-coil, which is further enclosed by two sets of inclined α2 and α4 helices. We further generated a series of XCC6012 variants and measured the unfolding temperatures and oligomeric states in order to investigate the nature of this novel tetramer. Discovery of this new PilZ domain architecture increases the complexity of c-di-GMP-mediated regulation.

Recognition of DNA by the helix-turn-helix global regulatory protein Lrp is modulated by the amino terminus

The AsnC/Lrp family of regulatory proteins links bacterial and archaeal transcription patterns to metabolism. In Escherichia coli, Lrp regulates approximately 400 genes, over 200 of them directly. In earlier studies, lrp genes from Vibrio cholerae, Proteus mirabilis, and E. coli were introduced into the same E. coli background and yielded overlapping but significantly different regulons. These differences were seen despite amino acid sequence identities of 92% (Vibrio) and 98% (Proteus) to E. coli Lrp, including complete conservation of the helix-turn-helix motifs. The N-terminal region contains many of the sequence differences among these Lrp orthologs, which led us to investigate its role in Lrp function. Through the generation of hybrid proteins, we found that the N-terminal diversity is responsible for some of the differences between orthologs in terms of DNA binding (as revealed by mobility shift assays) and multimerization (as revealed by gel filtration, dynamic light scattering, and analytical ultracentrifugation). These observations indicate that the N-terminal tail plays a significant role in modulating Lrp function, similar to what is seen for a number of other regulatory proteins.

The use of analytical sedimentation velocity to extract thermodynamic linkage

For 25 years, the Gibbs Conference on Biothermodynamics has focused on the use of thermodynamics to extract information about the mechanism and regulation of biological processes. This includes the determination of equilibrium constants for macromolecular interactions by high precision physical measurements. These approaches further reveal thermodynamic linkages to ligand binding events. Analytical ultracentrifugation has been a fundamental technique in the determination of macromolecular reaction stoichiometry and energetics for 85 years. This approach is highly amenable to the extraction of thermodynamic couplings to small molecule binding in the overall reaction pathway. In the 1980s this approach was extended to the use of sedimentation velocity techniques, primarily by the analysis of tubulin-drug interactions by Na and Timasheff. This transport method necessarily incorporates the complexity of both hydrodynamic and thermodynamic nonideality. The advent of modern computational methods in the last 20 years has subsequently made the analysis of sedimentation velocity data for interacting systems more robust and rigorous. Here we review three examples where sedimentation velocity has been useful at extracting thermodynamic information about reaction stoichiometry and energetics. Approaches to extract linkage to small molecule binding and the influence of hydrodynamic nonideality are emphasized. These methods are shown to also apply to the collection of fluorescence data with the new Aviv FDS.

Analysis of PKR-RNA interactions by sedimentation velocity

PKR is an interferon-induced kinase that plays a pivotal role in the innate immunity pathway for defense against viral infection. PKR is activated to undergo autophosphorylation upon binding to RNAs that contain duplex regions. Some highly structured viral RNAs do not activate and function as PKR inhibitors. In order to define the mechanisms of activation and inhibition of PKR by RNA, it is necessary to characterize the stoichiometries, affinities, and free energy couplings governing the assembly of the relevant complexes. We have found sedimentation velocity analytical ultracentrifugation to be particularly useful in the study of PKR-RNA interactions. Here, we describe protocols for designing and analyzing sedimentation velocity experiments that are generally applicable to studies of protein-nucleic acid interactions. Initially, velocity data obtained at multiple protein:RNA ratios are analyzed using the dc/dt method's to define the association model and to test whether the system is kinetically limited. The sedimentation velocity data obtained at multiple loading concentrations are then globally fitted to this model to determine the relevant association constants. The frictional ratios of the complexes are calculated using the fitted sedimentation coefficients to determine whether the hydrodynamic properties are physically reasonable. We demonstrate the utility of this approach using examples from our studies of PKR interactions with simple dsRNAs, the HIV TAR RNA, and the VAI RNA from adenovirus.

The analysis of macromolecular interactions by sedimentation equilibrium

The study of macromolecular interactions by sedimentation equilibrium is a highly technical method that requires great care in both the experimental design and data analysis. The complexity of the interacting system that can be analyzed is only limited by the ability to deconvolute the exponential contributions of each of the species to the overall concentration gradient. This is achieved in part through the use of multi-signal data collection and the implementation of soft mass conservation. We illustrate the use of these constraints in SEDPHAT through the study of an A+B+B⇌AB+B⇌ABB system and highlight some of the technical challenges that arise. We show that both the multi-signal analysis and mass conservation result in a precise and robust data analysis and discuss improvements that can be obtained through the inclusion of data from other methods such as sedimentation velocity and isothermal titration calorimetry.

Some statistical properties of differencing schemes for baseline correction of sedimentation velocity data

For the detailed analysis of sedimentation velocity data, the consideration of radial-dependent baseline offsets is indispensable. Two main approaches are data differencing ("delta-c" approach) and explicit inclusion of baseline parameters in the model ("direct boundary model" of the raw data). The current work aims to clarify the relationships between the two approaches. To this end, a simple model problem is examined. We show that the explicit consideration of the baseline in the model is equivalent to a differencing scheme where the average value is subtracted from all data points. Pairwise differencing in the delta-c approach always results in higher parameter uncertainty. For equidistant time points, the increase is smallest when the reference points are taken at intervals of 1/3 or 2/3 of the total number of time points. If the difference data are misinterpreted to be statistically independent samples, errors in the calculation of the parameter uncertainties can occur. Contrary to claims in the literature, we observe that there is no distinction in the approaches regarding their "model dependence"; both approaches arise from the integral or differential form of the same model, and both approaches can and should provide explicit estimates of the baseline values in the original data space for optimal discrimination between macromolecular sedimentation models.

Sedimentation velocity ultracentrifugation analysis for hydrodynamic characterization of G-quadruplex structures

Analytical ultracentrifugation (AUC) is a powerful technique for the characterization of hydrodynamic and thermodynamic properties. The intent of this article is to demonstrate the utility of sedimentation velocity (SV) studies to obtain hydrodynamic information for G-quadruplex (GQ) systems and to provide insights into one part of this process, namely, data analysis of existing SV data. An array of data analysis software is available, mostly written and continually developed by established researchers in the AUC field, with particularly rapid advances in the analysis of SV data. Each program has its own learning curve, and this article is intended as a resource in the data analysis process for beginning researchers in the field. We discuss the application of three of the most commonly used data analysis programs, DCDT+, Sedfit, and SedAnal, to the interpretation of SV data obtained in our laboratory on two GQ systems.