Structural basis of specific binding between Aurora A and TPX2 by molecular dynamics simulations

Peptide probes derived from pertuzumab by molecular dynamics modeling for HER2 positive tumor imaging

A high level of HER2 expression in breast cancer correlates with a higher tumor growth rate, high metastatic potential, and a poor long-term patient survival rate. Pertuzumab, a human monoclonal antibody, can reduce the effect of HER2 overexpression by preventing HER2 dimerization. In this study, a combination protocol of molecular dynamics modeling and MM/GBSA binding free energy calculations was applied to design peptides that interact with HER2 based on the HER2/pertuzumab crystal structure. Based on a β hairpin in pertuzumab from Glu46 to Lys65—which plays a key role in interacting with HER2—mutations were carried out in silico to improve the binding free energy of the hairpin that interacts with the Phe256-Lys314 of the HER2 protein. Combined the use of one-bead-one-compound library screening, among all the mutations, a peptide (58F63Y) with the lowest binding free energy was confirmed experimentally to have the highest affinity, and it may be used as a new probe in diagnosing and treating HER2-positive breast cancer.

Many therapeutic approaches, including the human monoclonal antibodies trastuzumab and pertuzumab, target the human epidermal growth factor receptor 2 (HER2) of any breast cancer that features HER2 overexpression. Compared to these antibodies, peptides have many advantages, including lower cost, easier synthesis, high affinity, and lower toxicity. Here, we first designed peptides that interact with HER2 protein based on the HER2/pertuzumab crystal structure (PDB entry: 1S78), using a combination protocol of molecular dynamics modeling, molecular mechanics/generalized Born solvent-accessible surface area (MM/GBSA) binding free energy calculations. Then, combined with the peptide library screening, six peptides were selected for further analysis and experimental validations. The results of ex vivo and in vivo experiments confirmed that one peptide (58F63Y) in particular has a strong affinity and high specificity to HER2-overexpressing tumors. This may due to more paired residues and lower binding free energy in peptide 58F63Y and HER2 complex based on free energy decomposition analysis and distances calculation. While both in silico and in vitro screenings point to the same high-affinity peptide, the findings suggest that in silico screening based on calculated binding free energies is rather reliable. Additionally, based on the calculation of binding free energies among mutants, we can reduce the library capacity of one-bead-one-compound screening. In summary, we present a rather simple and rapid means of deriving a peptide with a clear binding site to its target protein.

The Application of Ligand-Mapping Molecular Dynamics Simulations to the Rational Design of Peptidic Modulators of Protein-Protein Interactions

A computational ligand-mapping approach to detect protein surface pockets that interact with hydrophobic moieties is presented. In this method, we incorporated benzene molecules into explicit solvent molecular dynamics simulations of various protein targets. The benzene molecules successfully identified the binding locations of hydrophobic hot-spot residues and all-hydrocarbon cross-links from known peptidic ligands. They also unveiled cryptic binding sites that are occluded by side chains and the protein backbone. Our results demonstrate that ligand-mapping molecular dynamics simulations hold immense promise to guide the rational design of peptidic modulators of protein-protein interactions, including that of stapled peptides, which show promise as an exciting new class of cell-penetrating therapeutic molecules.

MD and QM/MM studies on long-chain L-α-hydroxy acid oxidase: substrate binding features and oxidation mechanism

Long-chain L-α-hydroxy acid oxidase (LCHAO) is a flavin mononucleotide (FMN)-dependent oxidase that dehydrogenates l-α-hydroxy acids to keto acids. There were two different mechanisms, named as hydride transfer (HT) mechanism and carbanion (CA) mechanism, respectively, proposed about the catalytic process for the FMN-dependent L-α-hydroxy acid oxidases on the basis of biochemical data. However, crystallographic and kinetic studies could not provide enough evidence to prove one of the mechanisms or eliminate the alternative. In the present studies, theoretical computations were carried out to study the molecular mechanism for LCHAO-catalyzed dehydrogenation of L-lactate. Our molecular dynamics (MD) simulations indicated that L-lactate prefers to bind with LCHAO in a hydride transfer mode rather than a carbanion mode. Quantum mechanics/molecular mechanics (QM/MM) calculations were further carried out to obtain the optimized structures of reactants, transition states, and products at the level of ONIOM-EE (B3LYP/6-311++G(d,p)//B3LYP/6-31G(d,p):AMBER). Quantum chemical studies indicated that LCHAO-catalyzed dehydrogenation of L-lactate would be a stepwise catalytic reaction in a hydride transfer mechanism but not a carbanion mechanism. MD simulations, binding free energy calculations, and QM/MM computations were also implemented on the complex between L-lactate and Y129F mutant LCHAO. By comparing the Y129F mutant system with the wild-type system, it was further confirmed that the key residue Tyr129 in the active site of LCHAO would not affect L-lactate's binding to LCHAO but play an important role on the catalytic reaction process through an H-bond interaction.