12 research outputs found
DataSheet1.PDF
<p>The computational modeling of peptide inhibitors to target protein-protein binding interfaces is growing in interest as these are often too large, too shallow, and too feature-less for conventional small molecule compounds. Here, we present a rare successful application of an alchemical binding free energy method for the calculation of converged absolute binding free energies of a series of protein-peptide complexes. Specifically, we report the binding free energies of a series of cyclic peptides derived from the LEDGF/p75 protein to the integrase receptor of the HIV1 virus. The simulations recapitulate the effect of mutations relative to the wild-type binding motif of LEDGF/p75, providing structural, energetic and dynamical interpretations of the observed trends. The equilibration and convergence of the calculations are carefully analyzed. Convergence is aided by the adoption of a single-decoupling alchemical approach with implicit solvation, which circumvents the convergence difficulties of conventional double-decoupling protocols. We hereby present the single-decoupling methodology and critically evaluate its advantages and limitations. We also discuss some of the challenges and potential pitfalls of binding free energy calculations for complex molecular systems which have generally limited their applicability to the quantitative study of protein-peptide binding equilibria.</p
Large Scale Affinity Calculations of Cyclodextrin HostâGuest Complexes: Understanding the Role of Reorganization in the Molecular Recognition Process
Hostâguest
inclusion complexes are useful models for understanding
the structural and energetic aspects of molecular recognition. Due
to their small size relative to much larger proteinâligand
complexes, converged results can be obtained rapidly for these systems
thus offering the opportunity to more reliably study fundamental aspects
of the thermodynamics of binding. In this work, we have performed
a large scale binding affinity survey of 57 ÎČ-cyclodextrin (CD)
hostâguest systems using the binding energy distribution analysis
method (BEDAM) with implicit solvation (OPLS-AA/AGBNP2). Converged
estimates of the standard binding free energies are obtained for these
systems by employing techniques such as parallel Hamiltonian replica
exchange molecular dynamics, conformational reservoirs, and multistate
free energy estimators. Good agreement with experimental measurements
is obtained in terms of both numerical accuracy and affinity rankings.
Overall, average effective binding energies reproduce affinity rank
ordering better than the calculated binding affinities, even though
calculated binding free energies, which account for effects such as
conformational strain and entropy loss upon binding, provide lower
root-mean-square errors when compared to measurements. Interestingly,
we find that binding free energies are superior rank order predictors
for a large subset containing the most flexible guests. The results
indicate that, while challenging, accurate modeling of reorganization
effects can lead to ligand design models of superior predictive power
for rank ordering relative to models based only on ligandâreceptor
interaction energies
Free Energy-Based Virtual Screening and Optimization of RNase H Inhibitors of HIVâ1 Reverse Transcriptase
We
report the results of a binding free energy-based virtual screening
campaign of a library of 77 α-hydroxytropolone derivatives against
the challenging RNase H active site of the reverse transcriptase (RT)
enzyme of human immunodeficiency virus-1. Multiple protonation states,
rotamer states, and binding modalities of each compound were individually
evaluated. The work involved more than 300 individual absolute alchemical
binding free energy parallel molecular dynamics calculations and over
1 million CPU hours on national computing clusters and a local campus
computational grid. The thermodynamic and structural measures obtained
in this work rationalize a series of characteristics of this system
useful for guiding future synthetic and biochemical efforts. The free
energy model identified key ligand-dependent entropic and conformational
reorganization processes difficult to capture using standard docking
and scoring approaches. Binding free energy-based optimization of
the lead compounds emerging from the virtual screen has yielded four
compounds with very favorable binding properties, which will be the
subject of further experimental investigations. This work is one of
the few reported applications of advanced-binding free energy models
to large-scale virtual screening and optimization projects. It further
demonstrates that, with suitable algorithms and automation, advanced-binding
free energy models can have a useful role in early-stage drug-discovery
programs
Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series
The accurate prediction of proteinâligand binding
affinities
is crucial for drug discovery. Alchemical free energy calculations
have become a popular tool for this purpose. However, the accuracy
and reliability of these methods can vary depending on the methodology.
In this study, we evaluate the performance of a relative binding free
energy protocol based on the alchemical transfer method (ATM), a novel
approach based on a coordinate transformation that swaps the positions
of two ligands. The results show that ATM matches the performance
of more complex free energy perturbation (FEP) methods in terms of
Pearson correlation but with marginally higher mean absolute errors.
This study shows that the ATM method is competitive compared to more
traditional methods in speed and accuracy and offers the advantage
of being applicable with any potential energy function
Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation
Molecular
dynamics modeling of complex biological systems is limited by finite
simulation time. The simulations are often trapped close to local
energy minima separated by high energy barriers. Here, we introduce
Hamiltonian replica exchange (H-REMD) with torsional flattening in
the Binding Energy Distribution Analysis Method (BEDAM), to reduce
energy barriers along torsional degrees of freedom and accelerate
sampling of intramolecular degrees of freedom relevant to proteinâligand
binding. The method is tested on a standard benchmark (T4 Lysozyme/L99A/<i>p</i>-xylene complex) and on a library of HIV-1 integrase complexes
derived from the SAMPL4 blind challenge. We applied the torsional
flattening strategy to 26 of the 53 known binders to the HIV Integrase
LEDGF site found to have a binding energy landscape funneled toward
the crystal structure. We show that our approach samples the conformational
space more efficiently than the original method without flattening
when starting from a poorly docked pose with incorrect ligand dihedral
angle conformations. In these unfavorable cases convergence to a binding
pose within 2â3 Ă
from the crystallographic pose is obtained
within a few nanoseconds of the Hamiltonian replica exchange simulation.
We found that torsional flattening is insufficient in cases where
trapping is due to factors other than torsional energy, such as the
formation of incorrect intramolecular hydrogen bonds and stacking.
Work is in progress to generalize the approach to handle these cases
and thereby make it more widely applicable
Simulating Replica Exchange: Markov State Models, Proposal Schemes, and the Infinite Swapping Limit
Replica
exchange molecular dynamics is a multicanonical simulation
technique commonly used to enhance the sampling of solvated biomolecules
on rugged free energy landscapes. While replica exchange is relatively
easy to implement, there are many unanswered questions about how to
use this technique most efficiently, especially because it is frequently
the case in practice that replica exchange simulations are not fully
converged. A replica exchange cycle consists of a series of molecular
dynamics steps of a set of replicas moving under different Hamiltonians
or at different thermodynamic states followed by one or more replica
exchange attempts to swap replicas among the different states. How
the replica exchange cycle is constructed affects how rapidly the
system equilibrates. We have constructed a Markov state model of replica
exchange (MSMRE) using long molecular dynamics simulations of a hostâguest
binding system as an example, in order to study how different implementations
of the replica exchange cycle can affect the sampling efficiency.
We analyze how the number of replica exchange attempts per cycle,
the number of MD steps per cycle, and the interaction between the
two parameters affects the largest implied time scale of the MSMRE
simulation. The infinite swapping limit is an important concept in
replica exchange. We show how to estimate the infinite swapping limit
from the diagonal elements of the exchange transition matrix constructed
from MSMRE âsimulations of simulationsâ as well as from
relatively short runs of the actual replica exchange simulations
Immunoselection scheme used with 2F5 Library
<p><b>I</b>. The 2F5 Library I was subjected to selection using varying amounts of competitive peptide Ac-EQELLELDKWASSLW-NH<sub>2</sub> in hopes of finding chimeras that were well recognized by 2F5. The chimeras used for immunization studies were subjected to a single round of immunoselection.</p
Representative neutralization curves reflecting the ability of antiserum GFA-278, derived from immunization with chimeric virus 12B1, to neutralize 12 HIV pseudoviruses from six subtypes.
<p>Neutralizing titers corresponding to 50% inhibition of the reporter luciferase activity (compared with inhibition by pooled normal guinea pig serum values) is denoted by the dashed lines.</p
Antigenicity assays reflecting differences among individual chimeric viruses.
<div><p><b>A</b>. Direct ELISA using 4E10 to capture the chimeras.</p>
<p><b>B</b>. Competitive ELISA using 4E10 to capture chimeras that could bind in the presence of competitive peptide, NK-15.</p>
<p><b>C</b>. MTT cell-killing assay using 4E10 to bind to and neutralize the chimeric viruses, preventing HeLa cell death.</p></div
Docking of HRV14 (shown in yellow) with a crystallographic structure (1TZG from the Protein Data Bank) of a complex of the 4E10 Fab (orange) complexed with the cognate Ac-KGWNWFDITNWGK-NH<sub>2</sub> peptide (blue ribbon and gold stick-and-balls model).
<p>Docking provides ideas for how to connect the 4E10 epitope to the surface of HRV.</p