We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates.

MM-PBSA is a post-processing end-state method to calculate free energies of molecules in solution. MMPBSA.py is a program written in Python for streamlining end-state free energy calculations using ensembles derived from molecular dynamics (MD) or Monte Carlo (MC) simulations. Several implicit solvation models are available with MMPBSA.py, including the Poisson-Boltzmann Model, the Generalized Born Model, and the Reference Interaction Site Model. Vibrational frequencies may be calculated using normal mode or quasi-harmonic analysis to approximate the solute entropy. Specific interactions can also be dissected using free energy decomposition or alanine scanning. A parallel implementation significantly speeds up the calculation by dividing frames evenly across available processors. MMPBSA.py is an efficient, user-friendly program with the flexibility to accommodate the needs of users performing end-state free energy calculations. The source code can be downloaded at http://ambermd.org/ with AmberTools, released under the GNU General Public License.

No abstract

No abstract

Estimating protein-protein interaction energies is a very challenging task for current simulation protocols. Here, absolute binding free energies are reported for the complex H-Ras/C-Raf1 using the MM-PB(GB)SA approach, testing the internal consistency and model dependence of the results. Averaging gas-phase energies (MM), solvation free energies as determined by Generalized Born models (GB/SA), and entropic contributions calculated by normal mode analysis for snapshots obtained from 10 ns explicit-solvent molecular dynamics in general results in an overestimation of the binding affinity when a solvent-accessible surface area-dependent model is used to estimate the nonpolar solvation contribution. Applying the sum of a cavity solvation free energy and explicitly modeled solute-solvent van der Waals interaction energies instead provides less negative estimates for the nonpolar solvation contribution. When the polar contribution to the solvation free energy is determined by solving the Poisson-Boltzmann equation (PB) instead, the calculated binding affinity strongly depends on the atomic radii set chosen. For three GB models investigated, different absolute deviations from PB energies were found for the unbound proteins and the complex. As an alternative to normal-mode calculations, quasiharmonic analyses have been performed to estimate entropic contributions due to changes of solute flexibility upon binding. However, such entropy estimates do not converge after 10 ns of simulation time, indicating that sampling issues may limit the applicability of this approach. Finally, binding free energies estimated from snapshots of the unbound proteins extracted from the complex trajectory result in an underestimate of binding affinity. This points to the need to exercise caution in applying the computationally cheaper "one-trajectory-alternative" to systems where there may be significant changes in flexibility and structure due to binding. The best estimate for the binding free energy of Ras-Raf obtained in this study of -8.3 kcal mol(-1) is in good agreement with the experimental result of -9.6 kcal mol(-1), however, further probing the transferability of the applied protocol that led to this result is necessary.

The influence of a xenobiotic compound on an organism is usually summarized by the expression biological activity. If a controlled, therapeutically relevant, and regulatory action is observed the compound has potential as a drug, otherwise its toxicity on the biological system is of interest. However, what do we understand by the biological activity? In principle, the overall effect on an organism has to be considered. However, because of the complexity of the interrelated processes involved, as a simplification primarily the "main action" on the organism is taken into consideration. On the molecular level, biological activity corresponds to the binding of a (low-molecular weight) compound to a macromolecular receptor, usually a protein. Enzymatic reactions or signal-transduction cascades are thereby influenced with respect to their function for the organism. We regard this binding as a process under equilibrium conditions; thus, binding can be described as an association or dissociation process. Accordingly, biological activity is expressed as the affinity of both partners for each other, as a thermodynamic equilibrium quantity. How well do we understand these terms and how well are they theoretically predictable today? The holy grail of rational drug design is the prediction of the biological activity of a compound. The processes involving ligand binding are extremely complicated, both ligand and protein are flexible molecules, and the energy inventory between the bound and unbound states must be considered in aqueous solution. How sophisticated and reliable are our experimental approaches to obtaining the necessary insight? The present review summarizes our current understanding of the binding affinity of a small-molecule ligand to a protein. Both theoretical and empirical approaches for predicting binding affinity, starting from the three-dimensional structure of a protein-ligand complex, will be described and compared. Experimental methods, primarily microcalorimetry, will be discussed. As a perspective, our own knowledge-based approach towards affinity prediction and experimental data on factorizing binding contributions to protein-ligand binding will be presented.

Detailed knowledge of how molecules recognize interaction partners and of the conformational preferences of biomacromolecules is pivotal for understanding biochemical processes. Such knowledge also provides the foundation for the design of novel molecules, as undertaken in pharmaceutical research. Computer-based free energy calculations enable a detailed investigation of the energetic factors that are responsible for molecular stability or binding affinity. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach is an efficient method for the calculation of free energies of diverse molecular systems. Here we describe the concepts of this approach and outline the practical proceeding. Furthermore we give an overview of the wide spectrum of problems that have been addressed with this method and of successful analyses carried out, thereby focussing on ambitious and recent studies. Limits of the approach in terms of accuracy and applicability are discussed. Despite these limitations MM-PBSA is a method with great potential that allows comparative free energy analyses for various molecular systems at low computational cost.

We present a comprehensive toolkit for Förster resonance energy transfer (FRET)-restrained modeling of biomolecules and their complexes for quantitative applications in structural biology. A dramatic improvement in the precision of FRET-derived structures is achieved by explicitly considering spatial distributions of dye positions, which greatly reduces uncertainties due to flexible dye linkers. The precision and confidence levels of the models are calculated by rigorous error estimation. The accuracy of this approach is demonstrated by docking a DNA primer-template to HIV-1 reverse transcriptase. The derived model agrees with the known X-ray structure with an r.m.s. deviation of 0.5 Å. Furthermore, we introduce FRET-guided 'screening' of a large structural ensemble created by molecular dynamics simulations. We used this hybrid approach to determine the formerly unknown configuration of the flexible single-strand template overhang.