Understanding diversity of human innate immunity receptors: analysis of surface features of leucine-rich repeat domains in NLRs and TLRs.

BACKGROUND: The human innate immune system uses a system of extracellular Toll-like receptors (TLRs) and intracellular Nod-like receptors (NLRs) to match the appropriate level of immune response to the level of threat from the current environment. Almost all NLRs and TLRs have a domain consisting of multiple leucine-rich repeats (LRRs), which is believed to be involved in ligand binding. LRRs, found also in thousands of other proteins, form a well-defined "horseshoe"-shaped structural scaffold that can be used for a variety of functions, from binding specific ligands to performing a general structural role. The specific functional roles of LRR domains in NLRs and TLRs are thus defined by their detailed surface features. While experimental crystal structures of four human TLRs have been solved, no structure data are available for NLRs. RESULTS: We report a quantitative, comparative analysis of the surface features of LRR domains in human NLRs and TLRs, using predicted three-dimensional structures for NLRs. Specifically, we calculated amino acid hydrophobicity, charge, and glycosylation distributions within LRR domain surfaces and assessed their similarity by clustering. Despite differences in structural and genomic organization, comparison of LRR surface features in NLRs and TLRs allowed us to hypothesize about their possible functional similarities. We find agreement between predicted surface similarities and similar functional roles in NLRs and TLRs with known agonists, and suggest possible binding partners for uncharacterized NLRs. CONCLUSION: Despite its low resolution, our approach permits comparison of molecular surface features in the absence of crystal structure data. Our results illustrate diversity of surface features of innate immunity receptors and provide hints for function of NLRs whose specific role in innate immunity is yet unknown.

New insight into long-range nonadditivity within protein double-mutant cycles.

Additivity principles in chemistry, biochemistry, and biophysics have been used extensively for decades. Nevertheless, it is well known that additivity frequently breaks down in complex biomacromolecules. Nonadditivity within protein double mutant free energy cycles of spatially close residue pairs is a generally well-understood phenomenon, whereas a robust description of nonadditivity extending over large distances remains to be developed. Here, we test the hypothesis that the mutational effects tend to be nonadditive if two structurally well-separated mutated residues belong to the same rigid cluster within the wild type protein, and additive if they are located within different clusters. We find the hypothesis to be statistically significant with P-values that range from 10(-5) to 10(-6). To the best of our knowledge, this result represents the first demonstration of a statistically significant preponderance for nonadditivity over long distances. These findings provide new insight into the origins of long-range nonadditivity in double mutant cycles, which complements the conventional wisdom that nonadditivity arises in double mutations involving contacting residues. Consequently, these results should have far-reaching implications for a proper understanding of protein stability, structure/function analyses, and protein design.

On the role of structural class of a protein with two-state folding kinetics in determining correlations between its size, topology, and folding rate.

The time it takes for proteins to fold into their native states varies over several orders of magnitude depending on their native-state topology, size, and amino acid composition. In a number of previous studies, it was found that there is strong correlation between logarithmic folding rates and contact order for proteins that fold with two-state kinetics, while such correlation is absent for three-state proteins. Conversely, strong correlations between folding rates and chain length occur within three-state proteins, but not in two-state proteins. Here, we demonstrate that chain lengths and folding rates of two-state proteins are not correlated with each other only when all-alpha, all-beta, and mixed-class proteins are considered together, which is typically the case. However, when considering all-alpha and all-beta two-state proteins separately, there is significant linear correlation between folding rate and size. Moreover, the sets of data points for the all-alpha and all-beta classes define asymptotes of lower and upper limits on folding rates of mixed-class proteins. By analyzing correlation of other topological parameters with folding rates of two-state proteins, we find that only the long-range order exhibits correlation with folding rates that is uniform over all three classes. It is also the only descriptor to provide statistically significant correlations for each of the three structural classes. We give an interpretation of this observation in terms of Makarov and Plaxco's diffusion-based topomer-search model.

Ensemble properties of network rigidity reveal allosteric mechanisms.

The distance constraint model (DCM) is a unique computational modeling paradigm that integrates mechanical and thermodynamic descriptions of macromolecular structure. That is, network rigidity calculations are used to account for nonadditivity within entropy components, thus restoring the utility of free-energy decomposition. The DCM outputs a large number of structural characterizations that collectively allow for quantified stability-flexibility relationships (QSFR) to be identified. In this review, we describe the theoretical underpinnings of the DCM and introduce several common QSFR metrics. Application of the DCM across protein families highlights the sensitivity within the set of protein structure residue-to-residue couplings. Further, we have developed a perturbation method to identify putative allosteric sites, where large changes in QSFR upon rigidification (mimicking ligand-binding) detect sites likely to invoke allosteric changes.

Dynamics of coregulator-induced conformational perturbations in androgen receptor ligand binding domain.

Androgen receptor (AR) coregulators modulate ligand-induced gene expression in a tissue specific manner. The molecular events that follow coactivator binding to AR and the mechanisms that govern the sequence-specific effects of AR coregulators are poorly understood. Using consensus coactivator sequence D11-FxxLF and biophysical techniques, we show that coactivator association is followed by conformational rearrangement in AR ligand binding domain (AR-LBD) that is enthalpically and entropically favorable with activation energy of 29.8±4.2 kJ/mol. Further characterization of ARA70 and SRC3-1 based consensus sequences reveal that each coactivator induces a distinct conformational state in the dihydrotestosterone:AR-LBD:coactivator complex. Complementary computational modeling revealed that coactivator induced specific alterations in the backbone flexibility of AR-LBD distant from the site of coactivator binding and that the intramolecular rearrangements in AR-LBD backbone induced by the two coactivator peptides were different. These data suggest that coactivators may impart specificity in the transcriptional machinery by changing the steady-state conformation of AR-LBD. These data provide direct evidence that even in the presence of same ligand, AR-LBD can occupy distinct conformational states depending on its interactions with specific coactivators in the tissues. We posit that this coactivator-specific conformational gating may then dictate subsequent binding partners and interaction/affinity for the DNA-response elements.

Kinetic and thermodynamic characterization of dihydrotestosterone-induced conformational perturbations in androgen receptor ligand-binding domain.

Ligand-induced conformational perturbations in androgen receptor (AR) are important in coactivator recruitment and transactivation. However, molecular rearrangements in AR ligand-binding domain (AR-LBD) associated with agonist binding and their kinetic and thermodynamic parameters are poorly understood. We used steady-state second-derivative absorption and emission spectroscopy, pressure and temperature perturbations, and 4,4'-bis-anilinonaphthalene 8-sulfonate (bis-ANS) partitioning to determine the kinetics and thermodynamics of the conformational changes in AR-LBD after dihydrotestosterone (DHT) binding. In presence of DHT, the second-derivative absorption spectrum showed a red shift and a change in peak-to-peak distance. Emission intensity increased upon DHT binding, and center of spectral mass was blue shifted, denoting conformational changes resulting in more hydrophobic environment for tyrosines and tryptophans within a more compact DHT-bound receptor. In pressure perturbation calorimetry, DHT-induced energetic stabilization increased the Gibbs free energy of unfolding to 8.4 +/- 1.3 kcal/mol from 3.5 +/- 1.6 kcal/mol. Bis-ANS partitioning studies revealed that upon DHT binding, AR-LBD underwent biphasic rearrangement with a high activation energy (13.4 kcal/mol). An initial, molten globule-like burst phase (k approximately 30 sec(-1)) with greater solvent accessibility was followed by rearrangement (k approximately 0.01 sec(-1)), leading to a more compact conformation than apo-AR-LBD. Molecular simulations demonstrated unique sensitivity of tyrosine and tryptophan residues during pressure unfolding with rearrangement of residues in the coactivator recruitment surfaces distant from the ligand-binding pocket. In conclusion, DHT binding leads to energetic stabilization of AR-LBD domain and substantial rearrangement of residues distant from the ligand-binding pocket. DHT binding to AR-LBD involves biphasic receptor rearrangement including population of a molten globule-like intermediate state.

Elliptic and circular dichroism effects in two-photon double ionization of atoms.

An ab initio parametrization of the two-photon double ionization amplitude from an s2 subshell of an atom in a 1S state is presented and used to predict two light polarization effects on photoelectron angular distributions that do not exist in single-photon double ionization: (i) elliptic dichroism and (ii) circular dichroism at equal energy sharing. Estimates for He show large magnitudes for these effects, which provide a means for polarization control of double ionization by vacuum ultraviolet light.

Nondipole effects in photo-double-ionization of He by a vuv photon.

We present a dynamical model-independent, ab initio parametrization of the quadrupole transition amplitude for photo-double-ionization of He. An asymmetry of the triple differential cross section induced by the nondipole corrections is discussed and shown to be significant even for an excess energy as low as 80 eV. We provide predictions for two different kinds of experiments in which nondipole effects should be observable with current experimental capabilities.