Performance of four different force fields for simulations of dipeptide conformations: GlyGly, GlyGly-, GlyGly · Cl-, GlyGly · Na+ and GlyGly · (H2O)2

Comparing the Performances of Force Fields in Conformational Searching of Hydrogen-Bond-Donating Catalysts

Here, we compare the relative performances of different force fields for conformational searching of hydrogen-bond-donating catalyst-like molecules. We assess the force fields by their predictions of conformer energies, geometries, low-energy, nonredundant conformers, and the maximum numbers of possible conformers. Overall, MM3, MMFFs, and OPLS3e had consistently strong performances and are recommended for conformationally searching molecules structurally similar to those in this study.

Specific DNA sequences allosterically enhance protein-protein interaction in a transcription factor through modulation of protein dynamics: implications for specificity of gene regulation

Most genes are regulated by multiple transcription factors, often assembling into multi-protein complexes in the gene regulatory region. Understanding of the molecular origin of specificity of gene regulatory complex formation in the context of the whole genome is currently inadequate. A phage transcription factor λ-CI forms repressive multi-protein complexes by binding to multiple binding sites in the genome to regulate the lifecycle of the phage. The protein-protein interaction between two DNA-bound λ-CI molecules is stronger when they are bound to the correct pair of binding sites, suggesting allosteric transmission of recognition of correct DNA sequences to the protein-protein interaction interface. Exploration of conformation and dynamics by time-resolved fluorescence anisotropy decay and molecular dynamics suggests a change in protein dynamics to be a crucial factor in mediating allostery. A lattice-based model suggests that DNA-sequence induced allosteric effects could be crucial underlying factors in differentially stabilizing the correct site-specific gene regulatory complexes. We conclude that transcription factors have evolved multiple mechanisms to augment the specificity of DNA-protein interactions in order to achieve an extraordinarily high degree of spatial and temporal specificities of gene regulatory complexes, and DNA-sequence induced allostery plays an important role in the formation of sequence-specific gene regulatory complexes.

Searching for bioactive conformations of drug-like ligands with current force fields: how good are we?

Drug-like ligands obtained from protein-ligand complexes deposited in the Protein Databank were subjected to conformational searching using various force fields and solvation settings. For each ligand, the resulting conformer pool was examined for the presence of the bioactive (crystal pose-like) conformation. Similarity of conformers toward the crystal-pose was quantified as the best achievable root mean squared deviation (RMSD, heavy atoms only). Analyzing the conformer pools generated by various force fields revealed only small differences in the likelihood of finding a crystal pose-like conformation. However, employing different solvents in the conformational search was found to be very important for achieving RMSDs below 1.0 Å. The best statistical values of likelihood were observed with a recently released force field covering a large portion of dihedral angles occurring in drug-like compounds in combination with the water as solvent. In order to enable computational chemists and modelers to efficiently use available software tools, we have additionally performed several focused analyses on ligands, grouped according to descriptors most relevant for the rational drug design.