A historical perspective on the application of molecular dynamics (MD) to biological macromolecules is presented. Recent developments combining state-of-the-art force fields with continuum solvation calculations have allowed us to reach the fourth era of MD applications in which one can often derive both accurate structure and accurate relative free energies from molecular dynamics trajectories. We illustrate such applications on nucleic acid duplexes, RNA hairpins, protein folding trajectories, and proteinligand, protein-protein, and protein-nucleic acid interactions.

No abstract

Fluorinated compounds are synthesized in pharmaceutical research on a routine basis and many marketed compounds contain fluorine. The present review summarizes some of the most frequently employed strategies for using fluorine substituents in medicinal chemistry. Quite often, fluorine is introduced to improve the metabolic stability by blocking metabolically labile sites. However, fluorine can also be used to modulate the physicochemical properties, such as lipophilicity or basicity. It may exert a substantial effect on the conformation of a molecule. Increasingly, fluorine is used to enhance the binding affinity to the target protein. Recent 3D-structure determinations of protein complexes with bound fluorinated ligands have led to an improved understanding of the nonbonding protein-ligand interactions that involve fluorine.

We report calculations of free energies of binding, DeltaG(bind), between a diverse set of nine ligands and avidin as well as between a peptide and streptavidin using the recently developed MM/PBSA approach. This method makes use of a molecular dynamics simulation of the ligand-protein complex to generate a thermally averaged ensemble of conformations of the molecules that are involved in the complex formation. Based on this set of structures, a free energy of binding is calculated using molecular mechanical and continuum solvent energies as well as including estimates of the nonpolar solvation free energy and solute entropy. We compare in our simulations different classes of ligands, involving biotin derivatives, the dye 2-(4'-hydroxyazobenzene)benzoic acid (HABA), and a cyclic hexapeptide, which cover a large range of binding free energies from -5 to -20 kcal/mol. Our calculations are able to reproduce experimental DeltaG(bind) values with a very good correlation coefficient of r(2) = 0.92. This agreement is considerably better than the results obtained with an alternate approach, the linear interaction energy approximation, for this system (r(2) = 0.55).

The formation of intramolecular hydrogen bonds has a very pronounced effect on molecular structure and properties. We study both aspects in detail with the aim of enabling a more rational use of this class of interactions in medicinal chemistry. On the basis of exhaustive searches in crystal structure databases, we derive propensities for intramolecular hydrogen bond formation of five- to eight-membered ring systems of relevance in drug discovery. A number of motifs, several of which are clearly underutilized in drug discovery, are analyzed in more detail by comparing small molecule and protein-ligand X-ray structures. To investigate effects on physicochemical properties, sets of closely related structures with and without the ability to form intramolecular hydrogen bonds were designed, synthesized, and characterized with respect to membrane permeability, water solubility, and lipophilicity. We find that changes in these properties depend on a subtle balance between the strength of the hydrogen bond interaction, geometry of the newly formed ring system, and the relative energies of the open and closed conformations in polar and unpolar environments. A number of general guidelines for medicinal chemists emerge from this study.

Halogen bonding triggers activity: Increasing binding affinity was observed for a series of covalent human Cathepsin L inhibitors by exchanging an aryl ring H atom with Cl, Br, and I, which undergo halogen bonding with the CO group of Gly61 in the S3 pocket of the enzyme. Fluorine, in contrast, strongly avoids halogen bonding (see scheme). The strong distance and angle dependence of halogen bonding was confirmed for biological systems.

The MM-PBSA approach has become a popular method for calculating binding affinities of biomolecular complexes. Published application examples focus on small test sets and few proteins and, hence, are of limited relevance in assessing the general validity of this method. To further characterize MM-PBSA, we report on a more extensive study involving a large number of ligands and eight different proteins. Our results show that applying the MM-PBSA energy function to a single, relaxed complex structure is an adequate and sometimes more accurate approach than the standard free energy averaging over molecular dynamics snapshots. The use of MM-PBSA on a single structure is shown to be valuable (a) as a postdocking filter in further enriching virtual screening results, (b) as a helpful tool to prioritize de novo design solutions, and (c) for distinguishing between good and weak binders (DeltapIC(50) > or = 2-3), but rarely to reproduce smaller free energy differences.

Based on torsion angle distributions of frequently occurring substructures, conformation preferences of druglike molecules are presented, accompanied by a review of the relevant literature. First, the relevance of the Cambridge Structural Database (CSD) for drug design is demonstrated by comparing substructures present in compounds entering clinical trials with those found in the CSD and protein-bound ligands in the Protein Data Bank (PDB). Next, we briefly highlight preferred conformations of elementary acyclic systems, followed by a discussion of sulfonamide conformations. Due to their central role in medicinal chemistry, we discuss properties of aryl ring substituents in depth, including biaryl systems and systems of two aryl rings connected by two acyclic bonds. For a subset of torsion motifs, we also compare torsion angle histograms derived from CSD structures with those derived from ligands in the PDB. Furthermore, selected properties of some six-and seven-membered ring systems are discussed. The article closes with a section on attractive sulfuroxygen contacts.