Thomas Steinbrecher

Towards accurate free energy calculations in ligand protein-binding studies.

Cells contain a multitude of different chemical reaction paths running simultaneously and quite independently next to each other. This amazing feat is enabled by molecular recognition, the ability of biomolecules to form stable and specific complexes with each other and with their substrates. A better understanding of this process, i.e. of the kinetics, structures and thermodynamic properties of biomolecule binding, would be invaluable in the study of biological systems. In addition, as the mode of action of many pharmaceuticals is based upon their inhibition or activation of biomolecule targets, predictive models of small molecule receptor binding are very helpful tools in rational drug design. Since the goal here is normally to design a new compound with a high inhibition strength, one of the most important thermodynamic properties is the binding free energy DeltaG(0). The prediction of binding constants has always been one of the major goals in the field of computational chemistry, because the ability to reliably assess a hypothetical compounds binding properties without having to synthesize it first would save a tremendous amount of work. The different approaches to this question range from fast and simple empirical descriptor methods to elaborate simulation protocols aimed at putting the computation of free energies onto a solid foundation of statistical thermodynamics. While the later methods are still not suited for the screenings of thousands of compounds that are routinely performed in computational drug design studies, they are increasingly put to use for the detailed study of protein ligand interactions. This review will focus on molecular mechanics force field based free energy calculations and their application to the study of protein ligand interactions. After a brief overview of other popular methods for the calculation of free energies, we will describe recent advances in methodology and a variety of exemplary studies of molecular dynamics simulation based free energy calculations.

Soft‐core potentials in thermodynamic integration: Comparing one‐ and two‐step transformations

Molecular dynamics‐based free energy calculations allow the determination of a variety of thermodynamic quantities from computer simulations of small molecules. Thermodynamic integration (TI) calculations can suffer from instabilities during the creation or annihilation of particles. This “singularity” problem can be addressed with “soft‐core” potential functions which keep pairwise interaction energies finite for all configurations and provide smooth free energy curves. “One‐step” transformations, in which electrostatic and van der Waals forces are simultaneously modified, can be simpler and less expensive than “two‐step” transformations in which these properties are changed in separate calculations. Here, we study solvation free energies for molecules of different hydrophobicity using both models. We provide recommended values for the two parameters αLJ and βC controlling the behavior of the soft‐core Lennard–Jones and Coulomb potentials and compare one‐ and two‐step transformations with regard to their suitability for numerical integration. For many types of transformations, the one‐step procedure offers a convenient and accurate approach to free energy estimates.

Controlling biological activity with light: diarylethene-containing cyclic peptidomimetics.

Photobiological processes in nature are usually triggered by nonpeptidic chromophores or by modified side chains. A system is presented in which the polypeptide backbone itself can be conformationally switched by light. An amino acid analogue was designed and synthesized based on a reversibly photoisomerizable diarylethene scaffold. This analogue was incorporated into the cyclic backbone of the antimicrobial peptide gramicidin S at several sites. The biological activity of the resulting peptidomimetics could then be effectively controlled by ultraviolet/visible light within strictly defined spatial and temporal limits.

Nonadiabatic QM/MM Simulations of Fast Charge Transfer in Escherichia coli DNA Photolyase

In this report, we study the photoactivation process in Escherichia coli DNA photolyase, involving long-range electron transport along a conserved chain of Trp residues between the protein surface and the flavin adenine dinucleotide (FAD) cofactor. Fully coupled nonadiabatic (Ehrenfest) quantum mechanics/molecular mechanics (QM/MM) simulations allow us to follow the time evolution of charge distributions over the natural time scale of multiple charge transfer events and conduct rigorous statistical analysis. Charge transfer rates in excellent agreement with experimental data are obtained without the need for any system-specific parametrization. The simulations are shown to provide a more detailed picture of electron transfer than a classical analysis of Marcus parameters. The protein and solvent both strongly influence the localization and transport properties of a positive charge, but the directionality of the process is mainly caused by solvent polarization. The time scales of charge movement, delocalization, protein relaxation and solvent reorganization overlap and lead to nonequilibrium reaction conditions. All these contributions are explicitly considered and fully resolved in the model used and provide an intricate picture of multistep biochemical electron transfer in a flexible, heterogeneous environment.

Binding of hydroxyquinoline probes to human serum albumin: combining molecular modeling and Förster's resonance energy transfer spectroscopy to understand flexible ligand binding.

Human serum albumin (HSA) is the most abundant protein in blood plasma. It has high relevance for the lipid metabolism, and its ability to bind a large variety of natural and pharmaceutical compounds makes it a crucial determinant of drug pharmaco-kinetics and -dynamics. The drug binding properties of HSA can be characterized by spectroscopic analysis of bound probe molecules. We have recently characterized the subdomain IIA binding site of HSA using three hydroxyquinoline derivatives. In this work, we extend our study by combining data from energy transfer experiments, ligand docking, and long molecular dynamics (MD) simulations. Multiple possible binding locations are found within the subdomain IIA site, and their solvent accessibility and interactions with ligands are analyzed in detail. Binding pockets appear well hydrated during simulations, with ligands in direct contact to water molecules at all times. Binding free energies in good agreement to experiment are calculated. The HSA apoprotein is found to exhibit significant conformational flexibility over 250 ns of simulation time, but individual domains remain structurally stable. Two rotamers of Trp214 were observed on a time scale longer than 50 ns in the MD simulations, supporting the experimental observation of two fluorescence lifetime components. The flexible protein structure and heterogeneous nature of its binding sites explain the ability of HSA to act as a versatile molecular transporter. The combination of experimental and computational molecular distance information allows the conclusion that hydroxyquinoline probes bind in a binding mode similar to the anticoagulant drug warfarin.

Direct simulation of electron transfer reactions in DNA radical cations

The electron transfer properties of DNA radical cations are important in DNA damage and repair processes. Fast long-range charge transfer has been demonstrated experimentally, but the subtle influences that experimental conditions as well as DNA sequences and geometries have on the details of electron transfer parameters are still poorly understood. In this work, we employ an atomistic QM/MM approach, based on a one-electron tight binding Hamiltonian and a classical molecular mechanics forcefield, to conduct nanosecond length MD simulations of electron holes in DNA oligomers. Multiple spontaneous electron transfer events were observed in 100 ns simulations with neighboring adenine or guanine bases. Marcus parameters of charge transfer could be extracted directly from the simulations. The reorganization energy lambda for hopping between neighboring bases was found to be ca. 25 kcal/mol and charge transfer rates of 4.1 x 10(9) s(-1) for AA hopping and 1.3 x 10(9) s(-1) for GG hopping were obtained.

Peptide-Lipid Interactions of the Stress-Response Peptide TisB That Induces Bacterial Persistence

The bacterial stress-response peptide TisB in Escherichia coli has been suggested to dissipate the transmembrane potential, such that the depletion of ATP levels induces the formation of dormant persister cells which can eventually form biofilms. We studied the structure and membrane interactions of TisB to find out whether it forms pores or other proton-selective channels. Circular dichroism revealed an amphiphilic α-helical structure when reconstituted in lipid vesicles, and oriented circular dichroism showed that the helix assumes a transmembrane alignment. The addition of TisB to dye-loaded vesicles caused leakage only at very high peptide concentration, notably with a Hill coefficient of 2, which suggests that dimers must be involved. Coarse-grained molecular dynamics simulations showed that membrane binding of monomeric TisB is rapid and spontaneous, and transmembrane insertion is energetically feasible. When TisB oligomers are assembled as transmembrane pores, these channels collapse during the simulations, but transmembrane dimers are found to be stable. Given the pattern of charges on the amphiphilic TisB helix, we postulate that antiparallel dimers could be assembled via a ladder of salt bridges. This electrostatic charge-zipper could enable protons to pass along a wire of trapped water molecules across the hydrophobic membrane.

Charge transfer in model peptides: obtaining Marcus parameters from molecular simulation

Charge transfer within and between biomolecules remains a highly active field of biophysics. Due to the complexities of real systems, model compounds are a useful alternative to study the mechanistic fundamentals of charge transfer. In recent years, such model experiments have been underpinned by molecular simulation methods as well. In this work, we study electron hole transfer in helical model peptides by means of molecular dynamics simulations. A theoretical framework to extract Marcus parameters of charge transfer from simulations is presented. We find that the peptides form stable helical structures with sequence dependent small deviations from ideal PPII helices. We identify direct exposure of charged side chains to solvent as a cause of high reorganization energies, significantly larger than typical for electron transfer in proteins. This, together with small direct couplings, makes long-range superexchange electron transport in this system very slow. In good agreement with experiment, direct transfer between the terminal amino acid side chains can be dicounted in favor of a two-step hopping process if appropriate bridging groups exist.

FASCAPLYSIN as a Specific Inhibitor for CDK4: Insights from Molecular Modelling

Cyclin-dependent kinases (CDKs) play a key role in the cell cycle and are important anti-cancer drug targets. The natural product fascaplysin inhibits CDK4 with surprising selectivity (IC50 = 0.4 µM) compared to the close homolog CDK2 (IC50 = 500 µM). Free energy calculations of the positively charged fascaplysin and an uncharged iso-electronic derivative in the CDK2 and CDK4 inhibitor complexes indicate that the positive charge of fascaplysin is crucial for selectivity. This finding will guide further improvements in the design of fascaplysin-based selective inhibitors for CDK4.

Charge Transfer in E. coli DNA Photolyase: Understanding Polarization and Stabilization Effects via QM/MM Simulations

We study fast hole transfer events in E. coli DNA photolyase, a key step in the photoactivation process, using a multiscale computational method that combines nonadiabatic propagation schemes and linear-scaling quantum chemical methods with molecular mechanics force fields. This scheme allows us to follow the time-dependent evolution of the electron hole in an unbiased fashion; that is, no assumptions about hole wave function localization, time scale separation, or adiabaticity of the process have to be made beforehand. DNA photolyase facilitates an efficient long-range charge transport between its flavin adenine dinucleotide (FAD) cofactor and the protein surface via a chain of evolutionary conserved Trp residues on the sub-nanosecond time scale despite the existence of multiple potential trap states. By including a large number of aromatic residues along the charge transfer pathway into the quantum description, we are able to identify the main pathway among alternative possible routes. The simulations show that charge transfer, which is extremely fast in this protein, occurs on the same time scale as the protein response to the electrostatic changes; that is, time-scale separation as often presupposed in charge transfer studies seems to be inappropriate for this system. Therefore, coupled equations of motion, which propagate electrons and nuclei simultaneously, appear to be necessary. The applied computational model is shown to capture the essentials of the reaction kinetics and thermodynamics while allowing direct simulations of charge transfer events on their natural time scale.