SiCoDEA: A Simple, Fast and Complete App for Analyzing the Effect of Individual Drugs and Their Combinations.

The administration of combinations of drugs is a method widely used in the treatment of different pathologies as it can lead to an increase in the therapeutic effect and a reduction in the dose compared to the administration of single drugs. For these reasons, it is of interest to study combinations of drugs and to determine whether a specific combination has a synergistic, antagonistic or additive effect. Various mathematical models have been developed, which use different methods to evaluate the synergy of a combination of drugs. We have developed an open access and easy to use app that allows different models to be explored and the most fitting to be chosen for the specific experimental data: SiCoDEA (Single and Combined Drug Effect Analysis). Despite the existence of other tools for drug combination analysis, SiCoDEA remains the most complete and flexible since it offers options such as outlier removal or the ability to choose between different models for analysis. SiCoDEA is an easy to use tool for analyzing drug combination data and to have a view of the various steps and offer different results based on the model chosen.

Understanding and Controlling Food Protein Structure and Function in Foods: Perspectives from Experiments and Computer Simulations.

The structure and interactions of proteins play a critical role in determining the quality attributes of many foods, beverages, and pharmaceutical products. Incorporating a multiscale understanding of the structure-function relationships of proteins can provide greater insight into, and control of, the relevant processes at play. Combining data from experimental measurements, human sensory panels, and computer simulations through machine learning allows the construction of statistical models relating nanoscale properties of proteins to the physicochemical properties, physiological outcomes, and tastes of foods. This review highlights several examples of advanced computer simulations at molecular, mesoscale, and multiscale levels that shed light on the mechanisms at play in foods, thereby facilitating their control. It includes a practical simulation toolbox for those new to in silico modeling.

Protonation of the Biliverdin IXα Chromophore in the Red and Far-Red Photoactive States of a Bacteriophytochrome.

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Chromophore-Protein Interplay during the Phytochrome Photocycle Revealed by Step-Scan FTIR Spectroscopy.

Phytochrome proteins regulate many photoresponses of plants and microorganisms. Light absorption causes isomerization of the biliverdin chromophore, which triggers a series of structural changes to activate the signaling domains of the protein. However, the structural changes are elusive, and therefore the molecular mechanism of signal transduction remains poorly understood. Here, we apply two-color step-scan infrared spectroscopy to the bacteriophytochrome from Deinococcus radiodurans. We show by recordings in H2O and D2O that the hydrogen bonds to the biliverdin D-ring carbonyl become disordered in the first intermediate (Lumi-R) forming a dynamic microenvironment, then completely detach in the second intermediate (Meta-R), and finally reform in the signaling state (Pfr). The spectra reveal via isotope labeling that the refolding of the conserved "PHY-tongue" region occurs with the last transition between Meta-R and Pfr. Additional changes in the protein backbone are detected already within microseconds in Lumi-R. Aided by molecular dynamics simulations, we find that a strictly conserved salt bridge between an arginine of the PHY tongue and an aspartate of the chromophore binding domains is broken in Lumi-R and the arginine is recruited to the D-ring C═O. This rationalizes how isomerization of the chromophore is linked to the global structural rearrangement in the sensory receptor. Our findings advance the structural understanding of phytochrome photoactivation.

Accurate Three States Model for Amino Acids with Two Chemically Coupled Titrating Sites in Explicit Solvent Atomistic Constant pH Simulations and pK(a) Calculations.

Correct protonation of titratable groups in biomolecules is crucial for their accurate description by molecular dynamics simulations. In the context of constant pH simulations, an additional protonation degree of freedom is introduced for each titratable site, allowing the protonation state to change dynamically with changing structure or electrostatics. Here, we extend previous approaches for an accurate description of chemically coupled titrating sites. A second reaction coordinate is used to switch between two tautomeric states of an amino acid with chemically coupled titratable sites, such as aspartate (Asp), glutamate (Glu), and histidine (His). To this aim, we test a scheme involving three protonation states. To facilitate charge neutrality as required for periodic boundary conditions and Particle Mesh Ewald (PME) electrostatics, titration of each respective amino acid is coupled to a "water" molecule that is charged in the opposite direction. Additionally, a force field modification for Amber99sb is introduced and tested for the description of carboxyl group protonation. Our three states model is tested by titration simulations of Asp, Glu, and His, yielding a good agreement, reproducing the correct geometry of the groups in their different protonation forms. We further show that the ion concentration change due to the neutralizing "water" molecules does not significantly affect the protonation free energies of the titratable groups, suggesting that the three states model provides a good description of biomolecular dynamics at constant pH.

Charge-Neutral Constant pH Molecular Dynamics Simulations Using a Parsimonious Proton Buffer.

In constant pH molecular dynamics simulations, the protonation states of titratable sites can respond to changes of the pH and of their electrostatic environment. Consequently, the number of protons bound to the biomolecule, and therefore the overall charge of the system, fluctuates during the simulation. To avoid artifacts associated with a non-neutral simulation system, we introduce an approach to maintain neutrality of the simulation box in constant pH molecular dynamics simulations, while maintaining an accurate description of all protonation fluctuations. Specifically, we introduce a proton buffer that, like a buffer in experiment, can exchange protons with the biomolecule enabling its charge to fluctuate. To keep the total charge of the system constant, the uptake and release of protons by the buffer are coupled to the titration of the biomolecule with a constraint. We find that, because the fluctuation of the total charge (number of protons) of a typical biomolecule is much smaller than the number of titratable sites of the biomolecule, the number of buffer sites required to maintain overall charge neutrality without compromising the charge fluctuations of the biomolecule, is typically much smaller than the number of titratable sites, implying markedly enhanced simulation and sampling efficiency.

Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics.

pH is an important parameter in condensed-phase systems, because it determines the protonation state of titratable groups and thus influences the structure, dynamics, and function of molecules in solution. In most force field simulation protocols, however, the protonation state of a system (rather than its pH) is kept fixed and cannot adapt to changes of the local environment. Here, we present a method, implemented within the MD package GROMACS, for constant pH molecular dynamics simulations in explicit solvent that is based on the λ-dynamics approach. In the latter, the dynamics of the titration coordinate λ, which interpolates between the protonated and deprotonated states, is driven by generalized forces between the protonated and deprotonated states. The hydration free energy, as a function of pH, is included to facilitate constant pH simulations. The protonation states of titratable groups are allowed to change dynamically during a simulation, thus reproducing average protonation probabilities at a certain pH. The accuracy of the method is tested against titration curves of single amino acids and a dipeptide in explicit solvent.

Inclusion of ionization states of ligands in affinity calculations.

When estimating binding affinities of a ligand, which can exists in multiple forms, for a target molecule, one must consider all possible competing equilibria. Here, a method is presented that estimates the contribution of the protonation equilibria of a ligand in solution to the measured or calculated binding affinity. The method yields a correction to binding constants that are based on the total concentration of inhibitor (the sum of all ionized forms of the inhibitor in solution) to account for the complexed form of the inhibitor only. The method is applied to the calculation of the difference in binding affinity of two inhibitors, 2-phosphoglycolate (PGA) and its phoshonate analog 3-phosphonopropionate (3PP), for the glycolytic enzyme triosephosphate isomerase. Both inhibitors have three titrating sites and exist in solution as a mixture of different forms. In this case the form that actually binds to the enzyme is present at relative low concentrations. The contributions of the alternative forms to the difference in binding energies is estimated by means of molecular dynamics simulations and corrections. The inhibitors undergo a pK(a) shift upon binding that is estimated by ab initio calculations. An interesting finding is that the affinity difference of the two inhibitors is not due to different interactions in the active site of the enzyme, but rather due to the difference in the solvation properties of the inhibitors.

Conformational changes upon calcium binding and phosphorylation in a synthetic fragment of calmodulin.

We have recently investigated by far-UV circular dichroism (CD) the effects of Ca(2+) binding and the phosphorylation of Ser 81 for the synthetic peptide CaM [54-106] encompassing the Ca(2+)-binding loops II and III and the central alpha helix of calmodulin (CaM) (Arrigoni et al., Biochemistry 2004, 43, 12788-12798). Using computational methods, we studied the changes in the secondary structure implied by these spectra with the aim to investigate the effect of Ca(2+) binding and the functional role of the phosphorylation of Ser 81 in the action of the full-length CaM. Ca(2+) binding induces the nucleation of helical structure by inducing side chain stacking of hydrophobic residues. We further investigated the effect of Ca(2+) binding by using near-UV CD spectroscopy. Molecular dynamics simulations of different fragments containing the central alpha-helix of CaM using various experimentally determined structures of CaM with bound Ca(2+) disclose the structural effects provided by the phosphorylation of Ser 81. This post-translational modification is predicted to alter the secondary structure in its surrounding and also to hinder the physiological bending of the central helix of CaM through an alteration of the hydrogen bond network established by the side chain of residue 81. Using quantum mechanical methods to predict the CD spectra for the frames obtained during the MD simulations, we are able to reproduce the relative experimental intensities in the far-UV CD spectra for our peptides. Similar conformational changes that take place in CaM [54-106] upon Ca(2+) binding and phosphorylation may occur in the full-length CaM.

The planar conformation of a strained proline ring: a QM/MM study.

QM and QM/MM energy calculations have been carried out on an atomic resolution structure of liganded triosephosphate isomerase (TIM) that has an active site proline (Pro168) in a planar conformation. The origin of the planarity of this proline has been identified. Steric interactions between the atoms of the proline ring and a tyrosine ring (Tyr166) on one side of the proline prevent the ring from adopting the up pucker (chi1 is approximately -30 degrees), while the side chain of a nearby alanine (Ala171) forbids the down pucker (chi1 is approximately +30 degrees). To obtain a proline conformation that is in agreement with the experimentally observed planar state, a quantum system of sufficient size is required and should at least include the nearby side chains of Tyr166, Ala171, and Glu129 to provide enough stabilization. It is argued that the current force fields for structure optimization do not describe strained protein fragments correctly. The proline is part of a catalytic loop that closes upon ligand binding. Comparison of the proline conformation in different TIM X-ray structures, indicates that in the closed conformation of TIM the proline is planar or nearly planar, while in the open conformation it is down puckered. This suggests that the planarity possibly plays a role in the overall catalytic cycle of TIM, presumable acting as a reservoir of energy that becomes available upon loop opening.

Incorporating the effect of ionic strength in free energy calculations using explicit ions.

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution.

Calculation of affinities of peptides for proteins.

Several methodologies were employed to calculate the Gibbs standard free energy of binding for a collection of protein-ligand complexes, where the ligand is a peptide and the protein is representative for various protein families. Almost 40 protein-ligand complexes were employed for a continuum approach, which considers the protein and the peptide at the atomic level, but includes solvent as a polarizable continuum. Five protein-ligand complexes were employed for an all-atom approach that relies on a combination of the double decoupling method with thermodynamic integration and molecular dynamics. These affinities were also computed by means of the linear interaction energy method. Although it generally proved rather difficult to predict the absolute free energies correctly, for some protein families the experimental ranking order was correctly reproduced by the continuum and all-atom approach. Considerable attention has also been given to correctly analyze the affinities of charged peptides, where it is required to judge the effect of one or more ions that are being decoupled in an all-atom approach to preserve electroneutrality. The various methods are further judged upon their merits.