Coarse-grained models for proteins

Coarse-grained models for proteins and biomolecular aggregates have recently enjoyed renewed interest. Coarse-grained representations combined with enhanced computer power currently allow the simulation of systems of biologically relevant size (submicrometric) and timescale (microsecond or millisecond). Although these techniques still cannot be considered as predictive as all-atom simulations, noticeable advances have recently been achieved, mainly concerning the use of more rigorous parameterization techniques and novel algorithms for sampling configurational space. Moreover, the simulation size scales and timescales coincide with those that can be reached with the most advanced spectroscopic techniques, making it possible to directly compare simulation and experiment.

Analysis of protein complexes using mass spectrometry

The versatile combination of affinity purification and mass spectrometry (AP-MS) has recently been applied to the detailed characterization of many protein complexes and large protein-interaction networks. The combination of AP-MS with other techniques, such as biochemical fractionation, intact mass measurement and chemical crosslinking, can help to decipher the supramolecular organization of protein complexes. AP-MS can also be combined with quantitative proteomics approaches to better understand the dynamics of protein-complex assembly.

Design and implementation of bimolecular fluorescence complementation (BiFC) assays for the visualization of protein interactions in living cells

Bimolecular fluorescence complementation (BiFC) analysis enables direct visualization of protein interactions in living cells. The BiFC assay is based on the discoveries that two non-fluorescent fragments of a fluorescent protein can form a fluorescent complex and that the association of the fragments can be facilitated when they are fused to two proteins that interact with each other. BiFC must be confirmed by parallel analysis of proteins in which the interaction interface has been mutated. It is not necessary for the interaction partners to juxtapose the fragments within a specific distance of each other because they can associate when they are tethered to a complex with flexible linkers. It is also not necessary for the interaction partners to form a complex with a long half-life or a high occupancy since the fragments can associate in a transient complex and un-associated fusion proteins do not interfere with detection of the complex. Many interactions can be visualized when the fusion proteins are expressed at levels comparable to their endogenous counterparts. The BiFC assay has been used for the visualization of interactions between many types of proteins in different subcellular locations and in different cell types and organisms. It is technically straightforward and can be performed using a regular fluorescence microscope and standard molecular biology and cell culture reagents.

A perspective on the use of iTRAQ reagent technology for protein complex and profiling studies

Proteomic research includes the characterization of protein mixtures in order to understand complex biological systems and determine relationships between proteins, their function, and protein-protein interactions. Often the goal of such research is to monitor changes of proteins in perturbed systems, a type of study referred to as differential expression analysis. To perform these studies requires the ability to execute some type of differential comparison of a given protein state in reference to some type of a control. The iTRAQ reagents are a set of isobaric reagents which are amine specific and allow for the identification and quantitation of up to four different samples simultaneously. The amine specificity of these reagents makes most peptides in a sample amenable to this labeling strategy with no loss of information from samples involving post-translational modifications, such as the scrutiny of signal transduction pathways that often involve phosphorylation phenomena. In addition, the multiplexing capacity of these reagents allows for information replication within certain LC-MS/MS experimental regimes, providing additional statistical validation within any given experiment. The results presented herein demonstrate a few examples of the wide variety of quantitative information that can be realized when undertaking such experimental approaches. These include temporal analysis of drug-induced-protein expression, discovery and elucidation of disease markers, and protein-protein interactions in multi-protein complexes.

Probing cellular protein complexes using single-molecule pull-down

Proteins perform most cellular functions in macromolecular complexes. The same protein often participates in different complexes to exhibit diverse functionality. Current ensemble approaches of identifying cellular protein interactions cannot reveal physiological permutations of these interactions. Here we describe a single-molecule pull-down (SiMPull) assay that combines the principles of a conventional pull-down assay with single-molecule fluorescence microscopy and enables direct visualization of individual cellular protein complexes. SiMPull can reveal how many proteins and of which kinds are present in the in vivo complex, as we show using protein kinase A. We then demonstrate a wide applicability to various signalling proteins found in the cytosol, membrane and cellular organelles, and to endogenous protein complexes from animal tissue extracts. The pulled-down proteins are functional and are used, without further processing, for single-molecule biochemical studies. SiMPull should provide a rapid, sensitive and robust platform for analysing protein assemblies in biological pathways.

Molecular characterization and comparison of the components and multiprotein complexes in the postsynaptic proteome

Characterization of the composition of the postsynaptic proteome (PSP) provides a framework for understanding the overall organization and function of the synapse in normal and pathological conditions. We have identified 698 proteins from the postsynaptic terminal of mouse CNS synapses using a series of purification strategies and analysis by liquid chromatography tandem mass spectrometry and large-scale immunoblotting. Some 620 proteins were found in purified postsynaptic densities (PSDs), nine in AMPA-receptor immuno-purifications, 100 in isolates using an antibody against the NMDA receptor subunit NR1, and 170 by peptide-affinity purification of complexes with the C-terminus of NR2B. Together, the NR1 and NR2B complexes contain 186 proteins, collectively referred to as membrane-associated guanylate kinase-associated signalling complexes. We extracted data from six other synapse proteome experiments and combined these with our data to provide a consensus on the composition of the PSP. In total, 1124 proteins are present in the PSP, of which 466 were validated by their detection in two or more studies, forming what we have designated the Consensus PSD. These synapse proteome data sets offer a basis for future research in synaptic biology and will provide useful information in brain disease and mental disorder studies.

A role for the MLL fusion partner ENL in transcriptional elongation and chromatin modification

Chimeric proteins joining the histone methyltransferase MLL with various fusion partners trigger distinctive lymphoid and myeloid leukemias. Here, we immunopurified proteins associated with ENL, a protein commonly fused to MLL. Identification of these ENL-associated proteins (EAPs) by mass spectrometry revealed enzymes with a known role in transcriptional elongation (RNA polymerase II C-terminal domain kinase [RNAPolII CTD] positive transcription elongation factor b [pTEFb]), and in chromatin modification (histone-H3 methyltransferase DOT1L) as well as other frequent MLL partners (AF4, AF5q31, and LAF4), and polycomb group members (RING1, CBX8, and BCoR). The composition of EAP was further verified by coimmunoprecipitation, 2-hybrid analysis, pull-down, and colocalization experiments. Purified EAP showed a histone H3 lysine 79-specific methylase activity, displayed a robust RNAPolII CTD kinase function, and counteracted the effect of the pTEFb inhibitor 5,6-dichloro-benzimidazole-riboside. In vivo, an ENL knock-down diminished genome-wide as well as gene-specific H3K79 dimethylation, reduced global run-on elongation, and inhibited transient transcriptional reporter activity. According to structure-function data, DOT1L recruitment was important for transformation by the MLL-ENL fusion derivative. These results suggest a function of ENL in histone modification and transcriptional elongation.

Real-time measurements of actin filament polymerization by total internal reflection fluorescence microscopy

Understanding the mechanism of actin polymerization and its regulation by associated proteins requires an assay to monitor polymerization dynamics and filament topology simultaneously. The only assay meeting these criteria is total internal reflection fluorescence microscopy (Amann and Pollard, 2001; Fujiwara et al., 2002). The fluorescence signal is fourfold stronger with actin labeled on Cys-374 with Oregon green rather than rhodamine. To distinguish growth at barbed and pointed ends we used image drift correction and maximum intensity projections to reveal points where single N-ethylmaleimide inactivated myosins attach filaments to the glass coverslip. We estimated association rates at high actin concentrations and dissociation rates near and below the critical actin concentration. At the barbed end, the association rate constant for Mg-ATP-actin is 7.4 microM(-1) s(-1) and the dissociation rate constant is 0.89 s(-1). At the pointed end the association and dissociation rate constants are 0.56 microM(-1) s(-1) and 0.19 s(-1). When vitamin D binding protein sequesters all free monomers, ADP-actin dissociates from barbed ends at 1.4 s(-1) and from pointed ends at 0.16 s(-1) regardless of buffer nucleotide.

A network representation of protein structures: implications for protein stability

This study views each protein structure as a network of noncovalent connections between amino acid side chains. Each amino acid in a protein structure is a node, and the strength of the noncovalent interactions between two amino acids is evaluated for edge determination. The protein structure graphs (PSGs) for 232 proteins have been constructed as a function of the cutoff of the amino acid interaction strength at a few carefully chosen values. Analysis of such PSGs constructed on the basis of edge weights has shown the following: 1), The PSGs exhibit a complex topological network behavior, which is dependent on the interaction cutoff chosen for PSG construction. 2), A transition is observed at a critical interaction cutoff, in all the proteins, as monitored by the size of the largest cluster (giant component) in the graph. Amazingly, this transition occurs within a narrow range of interaction cutoff for all the proteins, irrespective of the size or the fold topology. And 3), the amino acid preferences to be highly connected (hub frequency) have been evaluated as a function of the interaction cutoff. We observe that the aromatic residues along with arginine, histidine, and methionine act as strong hubs at high interaction cutoffs, whereas the hydrophobic leucine and isoleucine residues get added to these hubs at low interaction cutoffs, forming weak hubs. The hubs identified are found to play a role in bringing together different secondary structural elements in the tertiary structure of the proteins. They are also found to contribute to the additional stability of the thermophilic proteins when compared to their mesophilic counterparts and hence could be crucial for the folding and stability of the unique three-dimensional structure of proteins. Based on these results, we also predict a few residues in the thermophilic and mesophilic proteins that can be mutated to alter their thermal stability.

Molecular complexes at a glance: automated generation of two-dimensional complex diagrams

MOTIVATION

In this paper a new algorithmic approach is presented, which automatically generates structure diagrams of molecular complexes. A complex diagram contains the ligand, the amino acids of the protein interacting with the ligand and the hydrophilic interactions schematized as dashed lines between the corresponding atoms. The algorithm is based on a combinatorial optimization strategy which solves parts of the layout problem non-heuristically. The depicted molecules are represented as structure diagrams according to the chemical nomenclature. Due to the frequent usage of complex diagrams in the scientific literature as well as in text books dealing with structural biology, biochemistry and medicinal chemistry, the new algorithm is a key element for computer applications in these areas.

RESULTS

The method was implemented in the new software tool PoseView. It was tested on a representative dataset containing 305 protein-ligand complexes in total from the Brookhaven Protein Data Bank. PoseView was able to find collision-free layouts for more than three quarters of all complexes. In the following the layout generation algorithm is presented and, additional to the statistical results, representative test cases demonstrating the challenges of the layout generation will be discussed.

AVAILABILITY

The method is available as a webservice at http://www.zbh.uni-hamburg.de/poseview.

Models for excluded volume interaction between an unfolded protein and rigid macromolecular cosolutes: macromolecular crowding and protein stability revisited

Statistical-thermodynamic models for the excluded volume interaction between an unfolded polypeptide chain and a hard sphere or hard rod cosolute are presented, permitting estimation of the free energy of transfer of a polypeptide chain with fixed radius of gyration from a dilute (ideal) solution to a solution containing volume fraction of either cosolute. Also presented is a general thermodynamic description of the equilibrium between a unique native state and a manifold of unfolded or partially unfolded states of a protein distinguished by their respective radii of gyration. Together with results of a Monte Carlo calculation of the distribution of radii of gyration of four different unfolded proteins published by Goldenberg in 2003, these models are used to estimate the effect of intermolecular excluded volume upon an experimentally measurable apparent two-state constant for equilibrium between native and nonnative conformations of each of the four proteins, and upon the experimentally measurable root mean-square radius of gyration of the unfolded protein. Model calculations predict that addition of inert cosolutes at volume fractions exceeding 0.1 stabilizes the native state relative to unfolded states by an amount that increases strongly with and with the size of the native protein relative to the size of inert cosolute, and results in significant compaction of the manifold of unfolded states. Predicted effects are in qualitative and/or semiquantitative accord with the results of several published experimental studies.

The cellular EJC interactome reveals higher-order mRNP structure and an EJC-SR protein nexus

In addition to sculpting eukaryotic transcripts by removing introns, pre-mRNA splicing greatly impacts protein composition of the emerging mRNP. The exon junction complex (EJC), deposited upstream of exon-exon junctions after splicing, is a major constituent of spliced mRNPs. Here, we report comprehensive analysis of the endogenous human EJC protein and RNA interactomes. We confirm that the major "canonical" EJC occupancy site in vivo lies 24 nucleotides upstream of exon junctions and that the majority of exon junctions carry an EJC. Unexpectedly, we find that endogenous EJCs multimerize with one another and with numerous SR proteins to form megadalton sized complexes in which SR proteins are super-stoichiometric to EJC core factors. This tight physical association may explain known functional parallels between EJCs and SR proteins. Further, their protection of long mRNA stretches from nuclease digestion suggests that endogenous EJCs and SR proteins cooperate to promote mRNA packaging and compaction.

Self-assembled peptide nanotubes are uniquely rigid bioinspired supramolecular structures

We recently presented a novel class of self-assembled diphenylalanine-based peptide nanotubes. Here, for the first time, we present their mechanical properties, which we directly measured through indentation type experiments using atomic force microscopy. We find that the averaged point stiffness of the nanotubes is 160 N/m, and that they have a correspondingly high Young's modulus of approximately 19 GPa, as calculated by finite element analysis. This high value places these peptide nanotubes among the stiffest biological materials presently known, making them attractive building blocks for the design and assembly of biocompatible nanodevices.

Molecular architecture of a kinetochore-microtubule attachment site

Kinetochore attachment to spindle microtubule plus-ends is necessary for accurate chromosome segregation during cell division in all eukaryotes. The centromeric DNA of each chromosome is linked to microtubule plus-ends by eight structural-protein complexes. Knowing the copy number of each of these complexes at one kinetochore-microtubule attachment site is necessary to understand the molecular architecture of the complex, and to elucidate the mechanisms underlying kinetochore function. We have counted, with molecular accuracy, the number of structural protein complexes in a single kinetochore-microtubule attachment using quantitative fluorescence microscopy of GFP-tagged kinetochore proteins in the budding yeast Saccharomyces cerevisiae. We find that relative to the two Cse4p molecules in the centromeric histone, the copy number ranges from one or two for inner kinetochore proteins such as Mif2p, to 16 for the DAM-DASH complex at the kinetochore-microtubule interface. These counts allow us to visualize the overall arrangement of a kinetochore-microtubule attachment. As most of the budding yeast kinetochore proteins have homologues in higher eukaryotes, including humans, this molecular arrangement is likely to be replicated in more complex kinetochores that have multiple microtubule attachments.

RADIOLYTIC PROTEIN FOOTPRINTING WITH MASS SPECTROMETRY TO PROBE THE STRUCTURE OF MACROMOLECULAR COMPLEXES

Structural proteomics approaches using mass spectrometry are increasingly used in biology to examine the composition and structure of macromolecules. Hydroxyl radical-mediated protein footprinting using mass spectrometry has recently been developed to define structure, assembly, and conformational changes of macromolecules in solution based on measurements of reactivity of amino acid side chain groups with covalent modification reagents. Accurate measurements of side chain reactivity are achieved using quantitative liquid-chromatography-coupled mass spectrometry, whereas the side chain modification sites are identified using tandem mass spectrometry. In addition, the use of footprinting data in conjunction with computational modeling approaches is a powerful new method for testing and refining structural models of macromolecules and their complexes. In this review, we discuss the basic chemistry of hydroxyl radical reactions with peptides and proteins, highlight various approaches to map protein structure using radical oxidation methods, and describe state-of-the-art approaches to combine computational and footprinting data.

Protein structural perturbation and aggregation on homogeneous surfaces

We have demonstrated that globular proteins, such as hen egg lysozyme in phosphate buffered saline at room temperature, lose native structural stability and activity when adsorbed onto well-defined homogeneous solid surfaces. This structural loss is evident by alpha-helix to turns/random during the first 30 min and followed by a slow alpha-helix to beta-sheet transition. Increase in intramolecular and intermolecular beta-sheet content suggests conformational rearrangement and aggregation between different protein molecules, respectively. Amide I band attenuated total reflection/Fourier transformed infrared (ATR/FTIR) spectroscopy was used to quantify the secondary structure content of lysozyme adsorbed on six different self-assembled alkanethiol monolayer surfaces with -CH3, -OPh, -CF3, -CN, -OCH3, and -OH exposed functional end groups. Activity measurements of adsorbed lysozyme were in good agreement with the structural perturbations. Both surface chemistry (type of functional groups, wettability) and adsorbate concentration (i.e., lateral interactions) are responsible for the observed structural changes during adsorption. A kinetic model is proposed to describe secondary structural changes that occur in two dynamic phases. The results presented in this article demonstrate the utility of the ATR/FTIR spectroscopic technique for in situ characterization of protein secondary structures during adsorption on flat surfaces.

Ferulic acid destabilizes preformed beta-amyloid fibrils in vitro

Inhibition of the formation of beta-amyloid fibrils (fAbeta), as well as the destabilization of preformed fAbeta in the CNS, would be attractive therapeutic targets for the treatment of Alzheimer's disease (AD). We reported previously that curcumin (Cur) inhibits fAbeta formation from Abeta and destabilizes preformed fAbeta in vitro. Using fluorescence spectroscopic analysis with thioflavin T and electron microscopic studies, we examined the effects of ferulic acid (FA) on the formation, extension, and destabilization of fAbeta at pH 7.5 at 37 degrees C in vitro. We next compared the anti-amyloidogenic activities of FA with Cur, rifampicin, and tetracycline. Ferulic acid dose-dependently inhibited fAbeta formation from amyloid beta-peptide, as well as their extension. Moreover, it destabilized preformed fAbetas. The overall activity of the molecules examined was in the order of: Cur > FA > rifampicin = tetracycline. FA could be a key molecule for the development of therapeutics for AD.

SUSP1 antagonizes formation of highly SUMO2/3-conjugated species

Small ubiquitin-related modifier (SUMO) processing and deconjugation are mediated by sentrin-specific proteases/ubiquitin-like proteases (SENP/Ulps). We show that SUMO-specific protease 1 (SUSP1), a mammalian SENP/Ulp, localizes within the nucleoplasm. SUSP1 depletion within cell lines expressing enhanced green fluorescent protein (EGFP) fusions to individual SUMO paralogues caused redistribution of EGFP-SUMO2 and -SUMO3, particularly into promyelocytic leukemia (PML) bodies. Further analysis suggested that this change resulted primarily from a deficit of SUMO2/3-deconjugation activity. Under these circumstances, PML bodies became enlarged and increased in number. We did not observe a comparable redistribution of EGFP-SUMO1. We have investigated the specificity of SUSP1 using vinyl sulfone inhibitors and model substrates. We found that SUSP1 has a strong paralogue bias toward SUMO2/3 and that it acts preferentially on substrates containing three or more SUMO2/3 moieties. Together, our findings argue that SUSP1 may play a specialized role in dismantling highly conjugated SUMO2 and -3 species that is critical for PML body maintenance.

How aggregation and conformational scrambling of unfolded states govern fluorescence emission spectra

In a case study on five homologous alpha-amylases we analyzed the properties of unfolded states as obtained from treatments with GndHCl and with elevated temperatures. In particular the wavelength of the tryptophan fluorescence emission peak (lambda(max)) is a valuable parameter to characterize properties of the unfolded state. In all cases with a typical red shift of the emission spectrum occurring during structural unfolding we observed a larger magnitude of this shift for GndHCl-induced unfolding as compared to thermal unfolding. Although a quantitative relation between aggregation and reduction of the unfolding induced red shifts cannot be given, our data indicate that protein aggregation contributes significantly to smaller magnitudes of red shifts as observed during thermal unfolding. In addition, other properties of the unfolded states, most probable structural compactness or simply differences in the conformational scrambling, also affect the magnitude of red shifts. For the irreversible unfolding alpha-amylases studied here, transition temperatures and magnitudes of red shifts are strongly depending on heating rates. Lower protein concentrations and smaller heating rates lead to larger red shifts upon thermal unfolding, indicating that under these conditions the protein aggregation is less pronounced.

Cohesin, condensin, and the intramolecular centromere loop together generate the mitotic chromatin spring

Sister chromatid cohesion provides the mechanistic basis, together with spindle microtubules, for generating tension between bioriented chromosomes in metaphase. Pericentric chromatin forms an intramolecular loop that protrudes bidirectionally from the sister chromatid axis. The centromere lies on the surface of the chromosome at the apex of each loop. The cohesin and condensin structural maintenance of chromosomes (SMC) protein complexes are concentrated within the pericentric chromatin, but whether they contribute to tension-generating mechanisms is not known. To understand how pericentric chromatin is packaged and resists tension, we map the position of cohesin (SMC3), condensin (SMC4), and pericentric LacO arrays within the spindle. Condensin lies proximal to the spindle axis and is responsible for axial compaction of pericentric chromatin. Cohesin is radially displaced from the spindle axis and confines pericentric chromatin. Pericentric cohesin and condensin contribute to spindle length regulation and dynamics in metaphase. Together with the intramolecular centromere loop, these SMC complexes constitute a molecular spring that balances spindle microtubule force in metaphase.

Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces

The 16-22 amino-acid fragment of the beta-amyloid peptide associated with the Alzheimer's disease, Abeta, is capable of forming amyloid fibrils. Here we study the aggregation mechanism of Abeta16-22 peptides by unbiased thermodynamic simulations at the atomic level for systems of one, three, and six Abeta16-22 peptides. We find that the isolated Abeta16-22 peptide is mainly a random coil in the sense that both the alpha-helix and beta-strand contents are low, whereas the three- and six-chain systems form aggregated structures with a high beta-sheet content. Furthermore, in agreement with experiments on Abeta16-22 fibrils, we find that large parallel beta-sheets are unlikely to form. For the six-chain system, the aggregated structures can have many different shapes, but certain particularly stable shapes can be identified.

Assessing oligomerization of membrane proteins by four-pulse DEER: pH-dependent dimerization of NhaA Na+/H+ antiporter of E. coli

The pH dependence of the structure of the main Na(+)/H(+) antiporter NhaA of Escherichia coli is studied by continuous-wave (CW) and pulse electron paramagnetic resonance (EPR) techniques on singly spin-labeled mutants. Residues 225 and 254 were selected for site-directed spin labeling, as previous work suggested that they are situated in domains undergoing pH-dependent structural changes. A well-defined distance of 4.4 nm between residues H225R1 in neighboring molecules is detected by a modulation in double electron-electron resonance data. This indicates that NhaA exists as a dimer, as previously suggested by a low-resolution electron density map and cross-linking experiments. The modulation depth decreases reversibly when pH is decreased from 8 to 5.8. A quantitative analysis suggests a dimerization equilibrium, which depends moderately on pH. Furthermore, the mobility and polarity of the environment of a spin label attached to residue 225 change only slightly with changing pH, while no other changes are detected by CW EPR. As antiporter activity of NhaA changes drastically in the studied pH range, residues 225 and 254 are probably located not in the sensor or ion translocation sites themselves but in domains that convey the signal from the pH sensor to the translocation site.

New paradigm for macromolecular crystallography experiments at SSRL: automated crystal screening and remote data collection

Complete automation of the macromolecular crystallography experiment has been achieved at SSRL through the combination of robust mechanized experimental hardware and a flexible control system with an intuitive user interface. These highly reliable systems have enabled crystallography experiments to be carried out from the researchers' home institutions and other remote locations while retaining complete control over even the most challenging systems. A breakthrough component of the system, the Stanford Auto-Mounter (SAM), has enabled the efficient mounting of cryocooled samples without human intervention. Taking advantage of this automation, researchers have successfully screened more than 200 000 samples to select the crystals with the best diffraction quality for data collection as well as to determine optimal crystallization and cryocooling conditions. These systems, which have been deployed on all SSRL macromolecular crystallography beamlines and several beamlines worldwide, are used by more than 80 research groups in remote locations, establishing a new paradigm for macromolecular crystallography experimentation.

Template-synthesized protein nanotubes

A layer-by-layer deposition strategy for preparing protein nanotubes within the pores of a nanopore alumina template membrane is described. This method entails alternately exposing the template membrane to a solution of the desired protein and then to a solution of glutaraldehyde, which acts as cross-linking agent to hold the protein layers together. The number of layers of protein that make up the nanotube walls can be controlled at will by varying the number of alternate protein/glutaraldehyde cycles. After the desired number of layers have been deposited on the pore walls, the alumina template can be dissolved to liberate the protein nanotubes. We show here that glucose oxidase nanotubes prepared in this way catalyze glucose oxidation and that hemoglobin nanotubes retain their heme electroactivity. Furthermore, for the glucose oxidase nanotubes, the enzymatic activity increases with the nanotube wall thickness.

Optimal clustering for detecting near-native conformations in protein docking

Clustering is one of the most powerful tools in computational biology. The conventional wisdom is that events that occur in clusters are probably not random. In protein docking, the underlying principle is that clustering occurs because long-range electrostatic and/or desolvation forces steer the proteins to a low free-energy attractor at the binding region. Something similar occurs in the docking of small molecules, although in this case shorter-range van der Waals forces play a more critical role. Based on the above, we have developed two different clustering strategies to predict docked conformations based on the clustering properties of a uniform sampling of low free-energy protein-protein and protein-small molecule complexes. We report on significant improvements in the automated prediction and discrimination of docked conformations by using the cluster size and consensus as a ranking criterion. We show that the success of clustering depends on identifying the appropriate clustering radius of the system. The clustering radius for protein-protein complexes is consistent with the range of the electrostatics and desolvation free energies (i.e., between 4 and 9 Angstroms); for protein-small molecule docking, the radius is set by van der Waals interactions (i.e., at approximately 2 Angstroms). Without any a priori information, a simple analysis of the histogram of distance separations between the set of docked conformations can evaluate the clustering properties of the data set. Clustering is observed when the histogram is bimodal. Data clustering is optimal if one chooses the clustering radius to be the minimum after the first peak of the bimodal distribution. We show that using this optimal radius further improves the discrimination of near-native complex structures.