Analysis of the binding energies of testosterone, 5alpha-dihydrotestosterone, androstenedione and dehydroepiandrosterone sulfate with an antitestosterone antibody

Probing the structural and energetic basis of kinesin-microtubule binding using computational alanine-scanning mutagenesis

Kinesin-microtubule (MT) binding plays a critical role in facilitating and regulating the motor function of kinesins. To obtain a detailed structural and energetic picture of kinesin-MT binding, we performed large-scale computational alanine-scanning mutagenesis based on long-time molecular dynamics (MD) simulations of the kinesin-MT complex in both ADP and ATP states. First, we built three all-atom kinesin-MT models: human conventional kinesin bound to ADP and mouse KIF1A bound to ADP and ATP. Then, we performed 30 ns MD simulations followed by kinesin-MT binding free energy calculations for both the wild type and mutants obtained after substitution of each charged residue of kinesin with alanine. We found that the kinesin-MT binding free energy is dominated by van der Waals interactions and further enhanced by electrostatic interactions. The calculated mutational changes in kinesin-MT binding free energy are in excellent agreement with results of an experimental alanine-scanning study with a root-mean-square error of ~0.32 kcal/mol [Woehlke, G., et al. (1997) Cell 90, 207-216]. We identified a set of important charged residues involved in the tuning of kinesin-MT binding, which are clustered on several secondary structural elements of kinesin (including well-studied loops L7, L8, L11, and L12, helices α4, α5, and α6, and less-explored loop L2). In particular, we found several key residues that make different contributions to kinesin-MT binding in ADP and ATP states. The mutations of these residues are predicted to fine-tune the motility of kinesin by modulating the conformational transition between the ADP state and the ATP state of kinesin.

Reducing CDK4/6-p16(INK4a) interface. Computational alanine scanning of a peptide bound to CDK6 protein

The tumor suppressor gene p16INK4a is commonly found altered in numerous and different types of cancer. The encoded protein arrests cell cycle in G1 phase by binding to CDK4 and CDK6, inhibiting their kinase function. In 1995, a 20-residue peptide, extracted from p16INK4a protein sequence, was discovered that retains the cell cycle inhibition properties of the endogenous tumor suppressor. However, its structure has not been determined yet. In this article, the features of a theoretical structure of the peptide bound to CDK6 are reported. The complex was modeled from CDK6-p16INK4a X-ray crystal structure and through molecular dynamics. Final structure was assessed by comparing computed binding free energy changes, when single-alanine substitutions were brought about on the peptide, to experimental data. Better concordance was obtained when including a high level of solvation effects. Solute-solvent vdW energy and electrostatic energy between solute and first shells of water, computed through a force field and considering explicit waters, were also to be included to achieve reasonably good concordance between theoretical and experimental data.

Stereochemistry and position-dependent effects of carcinogens on TATA/TBP binding

The TATA-box binding protein (TBP) is required by eukaryotic RNA polymerases to bind to the TATA box, an eight-basepair DNA promoter element, to initiate transcription. Carcinogen adducts that bind to the TATA box can hamper this important process. Benzo[a]pyrene (BP) is a representative chemical carcinogen that can be metabolically converted to highly reactive benzo[a]pyrene diol epoxides (BPDE), which in turn can form chemically stereoisomeric BP-DNA adducts. Depending on the TATA-bound adduct's location and stereochemistry, TATA/TBP binding can be decreased or increased. Our previous study interpreted the location-dependent effect in terms of conformational freedom and major-groove space available to BP. Here we further explore specific structural changes of the TATA/TBP complex to help interpret the stereochemical effect in terms of the flexibility of the TATA bases that frame the intercalated adduct. Thermodynamic analyses using molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) yield large standard deviations, which make the computed binding free energies the same within the error bars and point to current limitations of free energy calculations of large and highly charged systems like DNA/protein complexes.

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.

MM-PBSA free energy analysis of endo-1,4-xylanase II (XynII)–substrate complexes: binding of the reactive sugar in a skew boat and chair conformation

The binding of xylotetraose in different conformations to the active site of endo-1,4-beta-xylanase II (XynII) from Trichoderma reesei was studied using molecular dynamics (MD) simulations and free energy analyses employing the MM-PBSA (Molecular Mechanics-Poisson-Boltzmann Surface Area) method. MD simulations of 1 ns were done for the substrate xylotetraose having the reactive sugar, which is bound in the -1 subsite of XynII in the 4C1 (chair) and 2So (skew boat) ground state conformations, and for the transition state of the XynII catalysed hydrolysis of the beta-glycosidic linkage. According to the simulations and free energy analysis, XynII binds the substrate with the -1 sugar in the 2So conformation 59.8 kJ mol(-1) tighter than the substrate with the sugar in the 4C1 conformation. The reactive 2So conformation resembles closely the reaction transition state and has the breaking glycosidic bond in a pseudo-axial orientation ready for facile bond cleavage. The transition state was calculated to be bound 77.1 kJ mol(-1) tighter than the 4C1 ground state conformation. The molecular mechanical interaction energy between the enzyme and the reactive pyranoside unit at the -1 subsite was 75.7 kJ mol(-1) more favorable for the binding of the 2So conformation than the 2C1 conformation, explaining the clearly tighter binding of the reactive structure The results of this study indicate that in the Michaelis complex XynII, a member of the family 11 xylanase, the substrate is bound in a skew boat conformation and in the catalytic reaction, the -1 sugar proceeds from the 4C1 conformation through 2So to the transition state with the -1 sugar in the 2,5B conformation.