scholarly journals Thermodynamic Integration in 3n Dimensions without Biases or Alchemy for Protein Interactions

2017 ◽  
Author(s):  
Liao Y Chen

ABSTRACTThermodynamic integration (TI), a powerful formalism for computing the Gibbs free energy, has been implemented for many biophysical processes characterized by one-dimensional order parameters with alchemical schemes that require delicate human efforts to choose/design biasing potentials for sampling the desired biophysical events and to remove their artifactitious consequences afterwards. Theoretically, an alchemical scheme is exact but practically, it causes error amplification. Small relative errors in the interaction parameters can be amplified many times in their propagation into the computed free energy [due to subtraction of similar numbers such as (105 ± 5) − (100 ± 5) = 5 ± 7], which would render the results significantly less accurate than the input interaction parameters. In this paper, we present an unsophisticated implementation of TI in 3n dimensions (3nD) (n=1,2,3…) without alchemy or biasing potentials. In TI3nD, the errors in the interaction parameters will not be amplified and human efforts are not required to design biasing potentials that generate unphysical consequences. Using TI3nD, we computed the standard free energies of three protein complexes: trometamol in Salmonella effector SpvD (n=1), biotin in avidin (n=2), and Colicin E9 endonuclease with cognate immunity protein Im9 (n=3) and the hydration energies of ten biologically relevant compounds (n=1 for water, acetamide, urea, glycerol, trometamol, ammonium and n=2 for erythritol, 1,3-propanediol, xylitol, biotin). The computed results all agree with available experimental data. Each of the 13 computations is accomplishable within two (for a hydration problem) to ten (for the protein-recognition problem) days on an inexpensive workstation (two Xeon E5-2665 2.4GHz CPUs and one nVidia P5000 GPU).

2021 ◽  
Author(s):  
Alexander Wade ◽  
Agastya Bhati ◽  
Shunzhou Wan ◽  
Peter Coveney

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat myriad diseases. In this work we examine the computation of alchemical relative binding free energies with an eye to assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2 and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2 and NAMD3. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between packages of 0.5 $kcal/mol$ The correlation between packages is very good with the lowest Spearman's, Pearson's and Kendall's tau correlation coefficient between two packages being 0.91, 0.89 and 0.74 respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.


Biochemistry ◽  
1995 ◽  
Vol 34 (42) ◽  
pp. 13751-13759 ◽  
Author(s):  
Russell Wallis ◽  
Kit-Yi Leung ◽  
Ansgar J. Pommer ◽  
Hortense Videler ◽  
Geoffrey R. Moore ◽  
...  

2021 ◽  
Vol 9 ◽  
Author(s):  
Elisa Martino ◽  
Sara Chiarugi ◽  
Francesco Margheriti ◽  
Gianpiero Garau

Because of the key relevance of protein–protein interactions (PPI) in diseases, the modulation of protein-protein complexes is of relevant clinical significance. The successful design of binding compounds modulating PPI requires a detailed knowledge of the involved protein-protein system at molecular level, and investigation of the structural motifs that drive the association of the proteins at the recognition interface. These elements represent hot spots of the protein binding free energy, define the complex lifetime and possible modulation strategies. Here, we review the advanced technologies used to map the PPI involved in human diseases, to investigate the structure-function features of protein complexes, and to discover effective ligands that modulate the PPI for therapeutic intervention.


Author(s):  
Sherlyn Jemimah ◽  
Masakazu Sekijima ◽  
M Michael Gromiha

Abstract Motivation Protein–protein interactions are essential for the cell and mediate various functions. However, mutations can disrupt these interactions and may cause diseases. Currently available computational methods require a complex structure as input for predicting the change in binding affinity. Further, they have not included the functional class information for the protein–protein complex. To address this, we have developed a method, ProAffiMuSeq, which predicts the change in binding free energy using sequence-based features and functional class. Results Our method shows an average correlation between predicted and experimentally determined ΔΔG of 0.73 and mean absolute error (MAE) of 0.86 kcal/mol in 10-fold cross-validation and correlation of 0.75 with MAE of 0.94 kcal/mol in the test dataset. ProAffiMuSeq was also tested on an external validation set and showed results comparable to structure-based methods. Our method can be used for large-scale analysis of disease-causing mutations in protein–protein complexes without structural information. Availability and implementation Users can access the method at https://web.iitm.ac.in/bioinfo2/proaffimuseq/. Supplementary information Supplementary data are available at Bioinformatics online.


2019 ◽  
Vol 21 (25) ◽  
pp. 13826-13834 ◽  
Author(s):  
Piero Procacci

Free energy perturbation (FEP) approaches with stratification have seen widespread and increasing use in computational studies of biologically relevant molecules.


Author(s):  
Shunzhou Wan ◽  
Peter V Coveney ◽  
Darren R Flower

The binding to the T cell receptor of wild-type and variant HTLV-1 Tax peptide complexed to the major histocompatibility complex has been investigated by means of molecular dynamics simulations. The binding free energy difference is calculated using the molecular mechanics Poisson–Boltzmann surface area and linear interaction energy methods. These methods extract useful information on the binding energetics from simulations of the physical states of the ligands, which are more computationally expedient than the commonly used thermodynamic integration method. The successful reproduction of the relative binding free energies shows that these methods can be useful for free energy calculations and the rational design of drugs and vaccines.


2019 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Esteban Vöhringer-Martinez

In molecular modeling the description of the interactions between molecules forms the basis for a correct prediction of macroscopic observables. Here, we derive atomic charges from the implicitly polarized electron density of eleven molecules in the SAMPL6 challenge using the Hirshfeld-I and Minimal Basis Set Iterative Stockholder(MBIS) partitioning method. These atomic charges combined with other parameters in the GAFF force field and different water/octanol models were then used in alchemical free energy calculations to obtain hydration and solvation free energies, which after correction for the polarization cost, result in the blind prediction of the partition coefficient. From the tested partitioning methods and water models the S-MBIS atomic charges with the TIP3P water model presented the smallest deviation from the experiment. Conformational dependence of the free energies and the energetic cost associated with the polarization of the electron density are discussed.


2019 ◽  
Vol 26 (21) ◽  
pp. 3890-3910 ◽  
Author(s):  
Branislava Gemovic ◽  
Neven Sumonja ◽  
Radoslav Davidovic ◽  
Vladimir Perovic ◽  
Nevena Veljkovic

Background: The significant number of protein-protein interactions (PPIs) discovered by harnessing concomitant advances in the fields of sequencing, crystallography, spectrometry and two-hybrid screening suggests astonishing prospects for remodelling drug discovery. The PPI space which includes up to 650 000 entities is a remarkable reservoir of potential therapeutic targets for every human disease. In order to allow modern drug discovery programs to leverage this, we should be able to discern complete PPI maps associated with a specific disorder and corresponding normal physiology. Objective: Here, we will review community available computational programs for predicting PPIs and web-based resources for storing experimentally annotated interactions. Methods: We compared the capacities of prediction tools: iLoops, Struck2Net, HOMCOS, COTH, PrePPI, InterPreTS and PRISM to predict recently discovered protein interactions. Results: We described sequence-based and structure-based PPI prediction tools and addressed their peculiarities. Additionally, since the usefulness of prediction algorithms critically depends on the quality and quantity of the experimental data they are built on; we extensively discussed community resources for protein interactions. We focused on the active and recently updated primary and secondary PPI databases, repositories specialized to the subject or species, as well as databases that include both experimental and predicted PPIs. Conclusion: PPI complexes are the basis of important physiological processes and therefore, possible targets for cell-penetrating ligands. Reliable computational PPI predictions can speed up new target discoveries through prioritization of therapeutically relevant protein–protein complexes for experimental studies.


2020 ◽  
Vol 27 (37) ◽  
pp. 6306-6355 ◽  
Author(s):  
Marian Vincenzi ◽  
Flavia Anna Mercurio ◽  
Marilisa Leone

Background:: Many pathways regarding healthy cells and/or linked to diseases onset and progression depend on large assemblies including multi-protein complexes. Protein-protein interactions may occur through a vast array of modules known as protein interaction domains (PIDs). Objective:: This review concerns with PIDs recognizing post-translationally modified peptide sequences and intends to provide the scientific community with state of art knowledge on their 3D structures, binding topologies and potential applications in the drug discovery field. Method:: Several databases, such as the Pfam (Protein family), the SMART (Simple Modular Architecture Research Tool) and the PDB (Protein Data Bank), were searched to look for different domain families and gain structural information on protein complexes in which particular PIDs are involved. Recent literature on PIDs and related drug discovery campaigns was retrieved through Pubmed and analyzed. Results and Conclusion:: PIDs are rather versatile as concerning their binding preferences. Many of them recognize specifically only determined amino acid stretches with post-translational modifications, a few others are able to interact with several post-translationally modified sequences or with unmodified ones. Many PIDs can be linked to different diseases including cancer. The tremendous amount of available structural data led to the structure-based design of several molecules targeting protein-protein interactions mediated by PIDs, including peptides, peptidomimetics and small compounds. More studies are needed to fully role out, among different families, PIDs that can be considered reliable therapeutic targets, however, attacking PIDs rather than catalytic domains of a particular protein may represent a route to obtain selective inhibitors.


Sign in / Sign up

Export Citation Format

Share Document