Calculations of the Free Energy of Interaction of the c-Fos−c-Jun Coiled Coil: Effects of the Solvation Model and the Inclusion of Polarization Effects

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Evidence for Homodimerization of the c-Fos Transcription Factor in Live Cells Revealed by Fluorescence Microscopy and Computer Modeling

The c-Fos and c-Jun transcription factors, members of the activator protein 1 (AP-1) complex, form heterodimers and bind to DNA via a basic leucine zipper and regulate the cell cycle, apoptosis, differentiation, etc. Purified c-Jun leucine zipper fragments could also form stable homodimers, whereas c-Fos leucine zipper homodimers were found to be much less stable in earlier in vitro studies. The importance of c-Fos overexpression in tumors and the controversy in the literature concerning c-Fos homodimerization prompted us to investigate Fos homodimerization. Förster resonance energy transfer (FRET) and molecular brightness analysis of fluorescence correlation spectroscopy data from live HeLa cells transfected with fluorescent-protein-tagged c-Fos indicated that c-Fos formed homodimers. We developed a method to determine the absolute concentrations of transfected and endogenous c-Fos and c-Jun, which allowed us to determine dissociation constants of c-Fos homodimers (Kd = 6.7 ± 1.7 μM) and c-Fos-c-Jun heterodimers (on the order of 10 to 100 nM) from FRET titrations. Imaging fluorescence cross-correlation spectroscopy (SPIM-FCCS) and molecular dynamics modeling confirmed that c-Fos homodimers were stably associated and could bind to the chromatin. Our results establish c-Fos homodimers as a novel form of the AP-1 complex that may be an autonomous transcription factor in c-Fos-overexpressing tissues and could contribute to tumor development.

Cloning and serotonergic regulation of RING finger protein38 (rnf38) in the brain of medaka (Oryzias latipes)

Serotonin (5-HT) is a key regulator of mood and sexual behaviors. 5-HT reuptake inhibitors have been used as antidepressants. Really interesting new gene (RING) finger proteins have been associated with 5-HT regulation but their role remains largely unknown. Some RING finger proteins are involved in the serotonergic system, therefore, we speculate that the gene expression of RING finger protein38 (rnf38) is regulated by the serotonergic system. In the present study, we aimed to identify the full length sequence of medaka (Oryzias latipes) rnf38 mRNA and investigate its association with the serotonergic system using an antidepressant, citalopram (CIT). We identified the full length rnf38 cDNA, which consisted of 2726 nucleotides spanning 12 exons and the deduced protein sequence consisting of 518 amino acid residues including a RING finger domain, a KIT motif and a coiled-coil domain. Medaka exposed to 10(-7)M of CIT showed anxiety-like behavior. The expressions of 5-HT-related genes, pet1, solute carrier family 6, member 4A (slc6a4) and tryptophan hydroxylase (tph2) were significantly low (P<0.05) in the hindbrain. On the other hand, rnf38 gene was significantly high (P<0.05) in the telencephalon and the hypothalamus. This shows that 5-HT synthesis and transport in the hindbrain is suppressed by CIT, which induces rnf38 gene expression in the forebrain where 5-HT neurons project. Thus, the expression of rnf38 is negatively regulated by the serotonergic system.

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC₅₀ values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC₅₀) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆G_solv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.

Identifying MurI uncompetitive inhibitors by correlating decomposed binding energies with bioactivity

MurI uncompetitive inhibitors can be virtually identified by a new method that correlates decomposed binding free energies with the bioactivity.

Free energy of binding of coiled-coil complexes with different electrostatic environments: the influence of force field polarisation and capping

Coiled-coils are well known protein-protein interaction motifs, with the leucine zipper region of activator protein-1 (AP-1) consisting of the c-Jun and c-Fos proteins being a typical example. Molecular dynamics (MD) simulations using the MM/GBSA method have been used to predict the free energy of interaction of these proteins. The influence of force field polarisation and capping on the predicted free energy of binding of complexes with different electrostatic environments (net charge) were investigated. Although both force field polarisation and peptide capping are important for the prediction of the absolute free energy of binding, peptide capping has the largest influence on the predicted free energy of binding. Polarisable simulations appear better suited to determine structural properties of the complexes of these proteins while non-polarisable simulations seem to give better predictions of the associated free energies of binding.

Molecular mechanism of the effects of salt and pH on the affinity between protein A and human immunoglobulin G1 revealed by molecular simulations

Protein A from the bacterium Staphylococcus aureus (SpA) has been widely used as an affinity ligand for purification of immunoglobulin G (IgG). The affinity between SpA and IgG is affected differently by salt and pH, but their molecular mechanisms still remain unclear. In this work, molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area analysis were performed to investigate the salt (NaCl) and pH effects on the affinity between SpA and human IgG1 (hIgG1). It is found that salt and pH affect the interactions of the hot spots of SpA by different mechanisms. In the salt solution, the compensations between helices I and II of SpA as well as between the nonpolar and electrostatic energies make the binding free energy independent of salt concentration. At pH 3.0, the unfavorable electrostatic interactions increase greatly and become the driving force for dissociation of the SpA-hIgG1 complex. They mainly come from the strong electrostatic repulsions between positively charged residues (H137, R146, and K154) of SpA and the positively charged residues of hIgG1. It is considered to be the molecular basis for hIgG1 elution from SpA-based affinity adsorbents at pH 3.0. The dissociation mechanism is then used to refine the binding model of SpA to hIgG1. The model is expected to help design high-affinity peptide ligands of IgG.

Bound water at protein-protein interfaces: partners, roles and hydrophobic bubbles as a conserved motif

Background

There is a great interest in understanding and exploiting protein-protein associations as new routes for treating human disease. However, these associations are difficult to structurally characterize or model although the number of X-ray structures for protein-protein complexes is expanding. One feature of these complexes that has received little attention is the role of water molecules in the interfacial region.

Methodology

A data set of 4741 water molecules abstracted from 179 high-resolution (≤ 2.30 Å) X-ray crystal structures of protein-protein complexes was analyzed with a suite of modeling tools based on the HINT forcefield and hydrogen-bonding geometry. A metric termed Relevance was used to classify the general roles of the water molecules.

Results

The water molecules were found to be involved in: a) (bridging) interactions with both proteins (21%), b) favorable interactions with only one protein (53%), and c) no interactions with either protein (26%). This trend is shown to be independent of the crystallographic resolution. Interactions with residue backbones are consistent for all classes and account for 21.5% of all interactions. Interactions with polar residues are significantly more common for the first group and interactions with non-polar residues dominate the last group. Waters interacting with both proteins stabilize on average the proteins' interaction (−0.46 kcal mol⁻¹), but the overall average contribution of a single water to the protein-protein interaction energy is unfavorable (+0.03 kcal mol⁻¹). Analysis of the waters without favorable interactions with either protein suggests that this is a conserved phenomenon: 42% of these waters have SASA ≤ 10 Ų and are thus largely buried, and 69% of these are within predominantly hydrophobic environments or “hydrophobic bubbles”. Such water molecules may have an important biological purpose in mediating protein-protein interactions.