From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

©2009 Gao, Skolnick. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.doi:10.1371/journal.pcbi.1000341DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Ca deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein

The Effect of Macromolecular Crowding, Ionic Strength and Calcium Binding on Calmodulin Dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.Comment: Accepted to PLoS Comp Biol, 201

Using SMOG 2 to simulate complex biomolecular assemblies

Over the last 20 years, the application of structure-based (Gō-like) models has ranged from protein folding with coarse-grained models to all-atom representations of large-scale molecular assemblies. While there are many variants that may be employed, the common feature of these models is that some (or all) of the stabilizing energetic interactions are defined based on the knowledge of a particular experimentally obtained conformation. With the generality of this approach, there was a need for a versatile computational platform for designing and implementing this class of models. To this end, the SMOG 2 software package provides an easy-to-use interface, where the user has full control of the model parameters. This software allows the user to edit XML-formatted files in order to provide definitions of new structure-based models. SMOG 2 reads these "template" files and maps the interactions onto specific structures, which are provided in PDB format. The force field files produced by SMOG 2 may then be used to perform simulations with a variety of popular molecular dynamics suites. In this chapter, we describe some of the key features of the SMOG 2 package, while providing examples and strategies for applying these techniques to complex (often large-scale) molecular assemblies, such as the ribosome