Calculation of Relative Binding Free Energy in the Water-Filled Active Site of Oligopeptide-Binding Protein A

The periplasmic oligopeptide binding protein A (OppA) represents a well-known example of water-mediated protein-ligand interactions. Here, we perform free-energy calculations for three different ligands binding to OppA, using a thermodynamic integration approach. The tripeptide ligands share a high structural similarity (all have the sequence KXK), but their experimentally-determined binding free energies differ remarkably. Thermodynamic cycles were constructed for the ligands, and simulations conducted in the bound and (freely solvated) unbound states. In the unbound state, it was observed that the difference in conformational freedom between alanine and glycine leads to a surprisingly slow convergence, despite their chemical similarity. This could be overcome by increasing the softness parameter during alchemical transformations. Discrepancies remained in the bound state however, when comparing independent simulations of the three ligands. These difficulties could be traced to a slow relaxation of the water network within the active site. Fluctuations in the number of water molecules residing in the binding cavity occur mostly on a timescale larger than the simulation time along the alchemical path. After extensive simulations, relative binding free energies that were converged to within thermal noise could be obtained, which agree well with available experimental data

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.Fil: Hernández González, Jorge Enrique. Universidade Estadual Paulista Julio de Mesquita Filho; Brasil. Universidade Federal do Rio de Janeiro; Brasil. University of Natural Resources and Life Sciences; AustriaFil: Salas Sarduy, Emir. Universidad Nacional de San Martín. Instituto de Investigaciones Biotecnológicas. - Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Instituto de Investigaciones Biotecnológicas; ArgentinaFil: Hernández Alvarez, Lilian. Universidade Estadual Paulista Julio de Mesquita Filho; BrasilFil: Barreto Gomes, Diego Enry. Universidade Federal do Rio de Janeiro; Brasil. Universidade Federal de Juiz de Fora; BrasilFil: Pascutti, Pedro Geraldo. Universidade Federal do Rio de Janeiro; BrasilFil: Oostenbrink, Chris. University Of Natural Resources And Life Sciences; AustriaFil: Leite, Vitor B. P.. Universidade Estadual Paulista Julio de Mesquita Filho; Brasi

An automated force field topology builder (ATB) and repository: Version 1.0

The Automated force field Topology Builder (ATB, http://compbio.biosci.uq.edu.au/atb) is a Web-accessible server that can provide topologies and parameters for a wide range of molecules appropriate for use in molecular simulations, computational drug design, and X-ray refinement. The ATB has three primary functions: (1) to act as a repository for molecules that have been parametrized as part of the GROMOS family of force fields, (2) to act as a repository for pre-equilibrated systems for use as starting configurations in molecular dynamics simulations (solvent mixtures, lipid systems pre-equilibrated to adopt a specific phase, etc.), and (3) to generate force field descriptions of novel molecules compatible with the GROMOS family of force fields in a variety of formats (GROMOS, GROMACS, and CNS). Force field descriptions of novel molecules are derived using a multistep process in which results from quantum mechanical (QM) calculations are combined with a knowledge-based approach to ensure compatibility (as far as possible) with a specific parameter set of the GROMOS force field. The ATB has several unique features: (1) It requires that the user stipulate the protonation and tautomeric states of the molecule. (2) The symmetry of the molecule is analyzed to ensure that equivalent atoms are assigned identical parameters. (3) Charge groups are assigned automatically. (4) Where the assignment of a given parameter is ambiguous, a range of possible alternatives is provided. The ATB also provides several validation tools to assist the user to assess the degree to which the topology generated may be appropriate for a given task. In addition to detailing the steps involved in generating a force field topology compatible with a specific GROMOS parameter set (GROMOS 53A6), the challenges involved in the automatic generation of force field parameters for atomic simulations in general are discussed