In order to study protein-protein nonbonded interactions, we present the development of a new reduced protein model that represents each amino acid residue with one to three coarse grains, whose physical properties are derived in a consistent bottom-up procedure from the higher-resolution all-atom AMBER force field. The resulting potential energy function is pairwise additive and includes distinct van-der-Waals and Coulombic terms. The van-der-Waals effective interactions are deduced from preliminary molecular dynamics simulations of all possible amino acid homodimers. They are best represented by a soft 1/r6 repulsion and a Gaussian attraction, with parameters obeying Lorentz-Berthelot mixing rules. For the Coulombic interaction, coarse grain charges are optimized for each separate protein in order to best represent the all-atom electrostatic potential outside the protein core. This approach leaves the possibility of using any implicit solvent model to describe solvation effects and electrostatic screening. The coarse-grained force field is tested carefully for a small homodimeric complex, the magainin. It is shown to reproduce satisfactorily the specificity of the all-atom underlying potential, in particular within a PB/SA solvation model. The coarse-grained potential is applied to the redocking prediction of three different protein-protein complexes: the magainin dimer, the barnase-barstar, and the trypsin-BPTI complexes. It is shown to provide per se an efficient and discriminating scoring energy function for the protein-protein docking problem that remains pertinent at both the global and refinement stage.

The structurally correlated dihedral angles epsilon and zeta are known for their large variability within the B-DNA backbone. We have used molecular modelling to study both energetic and mechanical features of these variations which can produce BI/BII transitions. Calculations were carried out on DNA oligomers containing either YpR or RpY dinucleotides steps within various sequence environments. The results indicate that CpA and CpG steps favour the BI/BII transition more than TpA or any RpY step. The stacking energy and its intra- and inter-strand components explain these effects. Analysis of neighbouring base pairs reveals that BI/BII transitions of CpG and CpA are easiest within (Y)n(R)n sequences. These can also induce a large vibrational amplitude for TpA steps within the BI conformation.