Protein folding using contact maps

We present the development of the idea to use dynamics in the space of contact maps as a computational approach to the protein folding problem. We first introduce two important technical ingredients, the reconstruction of a three dimensional conformation from a contact map and the Monte Carlo dynamics in contact map space. We then discuss two approximations to the free energy of the contact maps and a method to derive energy parameters based on perceptron learning. Finally we present results, first for predictions based on threading and then for energy minimization of crambin and of a set of 6 immunoglobulins. The main result is that we proved that the two simple approximations we studied for the free energy are not suitable for protein folding. Perspectives are discussed in the last section.Comment: 29 pages, 10 figure

Electrostatic Steering Accelerates C3d:CR2 Association.

Electrostatic effects are ubiquitous in protein interactions and are found to be pervasive in the complement system as well. The interaction between complement fragment C3d and complement receptor 2 (CR2) has evolved to become a link between innate and adaptive immunity. Electrostatic interactions have been suggested to be the driving factor for the association of the C3d:CR2 complex. In this study, we investigate the effects of ionic strength and mutagenesis on the association of C3d:CR2 through Brownian dynamics simulations. We demonstrate that the formation of the C3d:CR2 complex is ionic strength-dependent, suggesting the presence of long-range electrostatic steering that accelerates the complex formation. Electrostatic steering occurs through the interaction of an acidic surface patch in C3d and the positively charged CR2 and is supported by the effects of mutations within the acidic patch of C3d that slow or diminish association. Our data are in agreement with previous experimental mutagenesis and binding studies and computational studies. Although the C3d acidic patch may be locally destabilizing because of unfavorable Coulombic interactions of like charges, it contributes to the acceleration of association. Therefore, acceleration of function through electrostatic steering takes precedence to stability. The site of interaction between C3d and CR2 has been the target for delivery of CR2-bound nanoparticle, antibody, and small molecule biomarkers, as well as potential therapeutics. A detailed knowledge of the physicochemical basis of C3d:CR2 association may be necessary to accelerate biomarker and drug discovery efforts

Simulations of driven overdamped frictionless hard spheres

We introduce an event-driven simulation scheme for overdamped dynamics of frictionless hard spheres subjected to external forces, neglecting hydrodynamic interactions. Our event-driven approach is based on an exact equation of motion which relates the driving force to the resulting velocities through the geometric information characterizing the underlying network of contacts between the hard spheres. Our method allows for a robust extraction of the instantaneous coordination of the particles as well as contact force statistics and dynamics, under any chosen driving force, in addition to shear flow and compression. It can also be used for generating high-precision jammed packings under shear, compression, or both. We present a number of additional applications of our method.Comment: 12 pages, 9 figure

Ligand-biased ensemble receptor docking (LigBEnD): a hybrid ligand/receptor structure-based approach.

Ligand docking to flexible protein molecules can be efficiently carried out through ensemble docking to multiple protein conformations, either from experimental X-ray structures or from in silico simulations. The success of ensemble docking often requires the careful selection of complementary protein conformations, through docking and scoring of known co-crystallized ligands. False positives, in which a ligand in a wrong pose achieves a better docking score than that of native pose, arise as additional protein conformations are added. In the current study, we developed a new ligand-biased ensemble receptor docking method and composite scoring function which combine the use of ligand-based atomic property field (APF) method with receptor structure-based docking. This method helps us to correctly dock 30 out of 36 ligands presented by the D3R docking challenge. For the six mis-docked ligands, the cognate receptor structures prove to be too different from the 40 available experimental Pocketome conformations used for docking and could be identified only by receptor sampling beyond experimentally explored conformational subspace