The Mechanism of p53 Rescue by SUSP4.

p53 is an important tumor-suppressor protein deactivation of which by mdm2 results in cancers. A SUMO-specific protease 4 (SUSP4) was shown to rescue p53 from mdm2-mediated deactivation, but the mechanism is unknown. The discovery by NMR spectroscopy of a "p53 rescue motif" in SUSP4 that disrupts p53-mdm2 binding is presented. This 29-residue motif is pre-populated with two transient helices connected by a hydrophobic linker. The helix at the C-terminus binds to the well-known p53-binding pocket in mdm2 whereas the N-terminal helix serves as an affinity enhancer. The hydrophobic linker binds to a previously unidentified hydrophobic crevice in mdm2. Overall, SUSP4 appears to use two synergizing modules, the p53 rescue motif described here and a globular-structured SUMO-binding catalytic domain, to stabilize p53. A p53 rescue motif peptide exhibits an anti-tumor activity in cancer cell lines expressing wild-type p53. A pre-structures motif in the intrinsically disordered proteins is thus important for target recognition.

Contribution of proline to the pre-structuring tendency of transient helical secondary structure elements in intrinsically disordered proteins.

BACKGROUND: IDPs function without relying on three-dimensional structures. No clear rationale for such a behavior is available yet. PreSMos are transient secondary structures observed in the target-free IDPs and serve as the target-binding "active" motifs in IDPs. Prolines are frequently found in the flanking regions of PreSMos. Contribution of prolines to the conformational stability of the helical PreSMos in IDPs is investigated. METHODS: MD simulations are performed for several IDP segments containing a helical PreSMo and the flanking prolines. To measure the influence of flanking-prolines on the structural content of a helical PreSMo calculations were done for wild type as well as for mutant segments with Pro→Asp, His, Lys, or Ala. The change in the helicity due to removal of a proline was measured both for the PreSMo region and for the flanking regions. RESULTS: The α-helical content in ~70% of the helical PreSMos at the early stage of simulation decreases due to replacement of an N-terminal flanking proline by other residues whereas the helix content in nearly all PreSMos increases when the same replacements occur at the C-terminal flanking region. The helix destabilizing/terminating role of the C-terminal flanking prolines is more pronounced than the helix promoting effect of the N-terminal flanking prolines. GENERAL SIGNIFICANCE: This work represents a novel example demonstrating that a proline is encoded in an IDP with a defined purpose. The helical PreSMos presage their target-bound conformations. As they most likely mediate IDP-target binding via conformational selection their helical content can be an important feature for IDP function.

The functional and structural characterization of a novel oncogene GIG47 involved in the breast tumorigenesis.

BACKGROUND: A candidate oncogene GIG47, previously known as a neudesin with a neurotrophic activity, was identified by applying the differential expression analysis method. METHODS: As a first step to understand the molecular role of GIG47, we analyzed the expression profile of GIG47 in multiple human cancers including the breast cancer and characterized its function related to human carcinogenesis. Based on this oncogenic role of GIG47, we then embarked on determining the high-resolution structure of GIG47. We have applied multidimensional heteronuclear NMR methods to GIG47. RESULTS: GIG47 was over-expressed in primary breast tumors as well as other human tumors including carcinomas of the uterine cervix, malignant lymphoma, colon, lung, skin, and leukemia. To establish its role in the pathogenesis of breast cancer in humans, we generated stable transfectants of MCF7 cells. The ectopic expression of GIG47 in MCF7 cells promoted the invasiveness in the presence of 50% serum. In addition, it also resulted in the increased tumorigenicity in in vivo tumor formation assay. The tumorigenesis mechanism involving GIG47 might be mediated by the activation of MAPK and PI3K pathways. These results indicate that GIG47 plays a role in the breast tumorigenesis, thus representing a novel target for the treatment of breast cancer. To facilitate the development of GIG47-targeted therapeutics, we determined the structural configuration of GIG47. The high-resolution structure of GIG47 was obtained by combination of NMR and homology modeling. The overall structure of GIG47 has four α-helices and 6 ß-strands, arranged in a ß1-α1-ß2-ß3-α2-ß4-α3-α4-ß5-ß6 topology. There is a potential heme/steroid binding pocket formed between two helices α2 and α3. CONCLUSION: The determined three-dimensional structure of GIG47 may facilitate the development of potential anti-cancer agents.

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis.

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

Characterizing amyloid-beta protein misfolding from molecular dynamics simulations with explicit water.

Extracellular deposition of amyloid-beta (Aß) protein, a fragment of membrane glycoprotein called ß-amyloid precursor transmembrane protein (ßAPP), is the major characteristic for the Alzheimer's disease (AD). However, the structural and mechanistic information of forming Aß protein aggregates in a lag phase in cell exterior has been still limited. Here, we have performed multiple all-atom molecular dynamics simulations for physiological 42-residue amyloid-beta protein (Aß42) in explicit water to characterize most plausible aggregation-prone structure (APS) for the monomer and the very early conformational transitions for Aß42 protein misfolding process in a lag phase. Monitoring the early sequential conformational transitions of Aß42 misfolding in water, the APS for Aß42 monomer is characterized by the observed correlation between the nonlocal backbone H-bond formation and the hydrophobic side-chain exposure. Characteristics on the nature of the APS of Aß42 allow us to provide new insight into the higher aggregation propensity of Aß42 over Aß40, which is in agreement with the experiments. On the basis of the structural features of APS, we propose a plausible aggregation mechanism from APS of Aß42 to form fibril. The structural and mechanistic observations based on these simulations agree with the recent NMR experiments and provide the driving force and structural origin for the Aß42 aggregation process to cause AD.

Rescuing p53 from mdm2 by a pre-structured motif in intrinsically unfolded SUMO specific protease 4.

Many intrinsically unstructured/unfolded proteins (IUPs) contain transient local secondary structures even though they are "unstructured" in a tertiary sense. These local secondary structures are named "pre-structured motifs (PreSMos)" and in fact are the specificity determinants for IUP-target binding, i.e., the active sites in IUPs. Using high-resolution NMR we have delineated a PreSMo active site in the intrinsically unfolded mid-domain (residues 201-300) of SUMO-specific protease 4 (SUSP4). This 29-residue motif which we termed a p53 rescue motif can protect p53 from mdm2 quenching by binding to the p53-helix binding pocket in mdm2(3-109). Our work demonstrates that the PreSMo approach is quite effective in providing a structural rationale for interactions of p53-mdm2- SUSP4 and opens a novel avenue for designing mdm2- inhibiting anticancer compounds. [BMB Reports 2017; 50(10): 485-486].

Structural investigation on the intrinsically disordered N-terminal region of HPV16 E7 protein.

Human papillomavirus (HPV) is the major cause of cervical cancer, a deadly threat to millions of females. The early oncogene product (E7) of the high-risk HPV16 is the primary agent associated with HPV-related cervical cancers. In order to understand how E7 contributes to the transforming activity, we investigated the structural features of the flexible N-terminal region (46 residues) of E7 by carrying out N-15 heteronuclear NMR experiments and replica exchange molecular dynamics simulations. Several NMR parameters as well as simulation ensemble structures indicate that this intrinsically disordered region of E7 contains two transient (10-20% populated) helical pre-structured motifs that overlap with important target binding moieties such as an E2F-mimic motif and a pRb-binding LXCXE segment. Presence of such target-binding motifs in HPV16 E7 provides a reasonable explanation for its promiscuous target-binding behavior associated with its transforming activity. [BMB Reports 2016; 49(8): 431-436].

Wide-line NMR and DSC studies on intrinsically disordered p53 transactivation domain and its helically pre-structured segment.

Wide-line 1H NMR intensity and differential scanning calorimetry measurements were carried out on the intrinsically disordered 73-residue full transactivation domain (TAD) of the p53 tumor suppressor protein and two peptides: one a wild type p53 TAD peptide with a helix pre-structuring property, and a mutant peptide with a disabled helix-forming propensity. Measurements were carried out in order to characterize their water and ion binding characteristics. By quantifying the number of hydrate water molecules, we provide a microscopic description for the interactions of water with a wild-type p53 TAD and two p53 TAD peptides. The results provide direct evidence that intrinsically disordered proteins (IDPs) and a less structured peptide not only have a higher hydration capacity than globular proteins, but are also able to bind a larger amount of charged solute ions. [BMB Reports 2016; 49(9): 497-501].

A pre-structured helix in the intrinsically disordered 4EBP1.

The eIF4E-binding protein 1 (4EBP1) has long been known to be completely unstructured without any secondary structures, which contributed significantly to the proposal of the induced fit mechanism for target binding of intrinsically disordered proteins. We show here that 4EBP1 is not completely unstructured, but contains a pre-structured helix.

Elucidating the Molecular Origin of Hydrolysis Energy of Pyrophosphate in Water.

The molecular origin of the energy produced by the ATP hydrolysis has been one of the long-standing fundamental issues. A classical view is that the negative hydrolysis free energy of ATP originates from intramolecular effects connected with the backbone P-O bond, so called "high-energy bond". On the other hand, it has also been recognized that solvation effects are essential in determining the hydrolysis free energy. Here, using the 3D-RISM-SCF (three-dimensional reference interaction site model self-consistent field) theory that integrates the ab initio quantum chemistry method and the statistical mechanical theory of liquids, we investigate the molecular origin of hydrolysis free energy of pyrophosphate, an ATP analogue, in water. We demonstrate that our theory quantitatively reproduces the experimental results without the use of empirical parameters. We clarify the crucial role of water in converting the hydrolysis free energy in the gas phase determined solely by intramolecular effects, which ranges from endothermic, thermoneutral, to highly exothermic depending on the charged state of pyrophosphate, into moderately exothermic in the aqueous phase irrespective of the charged state as observed in experimental data. We elucidate that this is brought about by different natures of solute-water interactions depending on the charged state of solute species: the hydration free energy of low-charged state is mainly subjected to short-range hydrogen-bonds, while that of high-charged state is dominated by long-range electrostatic interactions. We thus provide unambiguous evidence on the critical role of water in determining the ATP hydrolysis free energy.

Molecular docking study on the α3ß2 neuronal nicotinic acetylcholine receptor complexed with α-conotoxin GIC.

Nicotinic acetylcholine receptors (nAChRs) are a diverse family of homo- or heteropentameric ligand-gated ion channels. Understanding the physiological role of each nAChR subtype and the key residues responsible for normal and pathological states is important. α-Conotoxin neuropeptides are highly selective probes capable of discriminating different subtypes of nAChRs. In this study, we performed homology modeling to generate the neuronal α3, ß2 and ß4 subunits using the x-ray structure of the α1 subunit as a template. The structures of the extracellular domains containing ligand binding sites in the α3ß2 and α3ß4 nAChR subtypes were constructed using MD simulations and ligand docking processes in their free and ligand-bound states using α-conotoxin GIC, which exhibited the highest α3ß2 vs. α3ß4 discrimination ratio. The results provide a reasonable structural basis for such a discriminatory ability, supporting the idea that the present strategy can be used for future investigations on nAChR-ligand complexes.

Understanding pre-structured motifs (PreSMos) in intrinsically unfolded proteins.

Intrinsically unfolded proteins (IUPs) do not obey the golden rule of structural biology, 3D structure = function, as they manifest their inherent functions without resorting to three-dimensional structures. Absence of a compact globular topology in these proteins strongly implies that their ligand recognition processes should involve factors other than spatially well-defined binding pockets. Heteronuclear multidimensional (HetMulD) NMR spectroscopy assisted with a stable isotope labeling technology is a powerful tool for quantitatively investigating detailed structural features in IUPs. In particular, it allows us to delineate the presence and locations of pre-structured motifs (PreSMos) on a per-residue basis. PreSMos are the transient local structural elements that presage target-bound conformations and act as specificity determinants for IUP recognition by target proteins. Here, we present a brief chronicle of HetMulD NMR studies on IUPs carried out over the past two decades along with a discussion on the functional significance of PreSMos in IUPs.

Cu2+ ion-induced self-assembly of pyrenylquinoline with a pyrenyl excimer formation.

Synthesis of a novel pyrene derivative sensor (1) and its intermolecular binding pattern to Cu(2+) in CH(3)CN were investigated. Upon Cu(2+) binding, the sensor exhibited a strong static excimer emission at 460 nm, along with a weak monomer emission at 388 nm. The excimer emission intensity induced by the Cu(2+) ion declined as the spacer length between the pyrene and quinolinylamide unit increased. The Cu(2+) ion-induced self-assembled pyrenyl excimer formation is rationalized by fluorescence experiments and theoretical DFT calculations.

Nitrile and thiocyanate IR probes: quantum chemistry calculation studies and multivariate least-square fitting analysis.

Hydration effects on the C[Triple Bond]N stretching mode frequencies of MeCN and MeSCN are investigated by carrying out ab initio calculations for a number of MeCN-water and MeSCN-water complexes with varying number of water molecules. It is found that the CN frequency shift induced by the hydrogen-bonding interactions with water molecules originate from two different ways to form hydrogen bonds with the nitrogen atom of the CN group. Considering the MeCN- and MeSCN-water cluster calculation results as databases, we first examined the validity of vibrational Stark effect relationship between the CN frequency and the electric field component parallel to the CN bond and found no strong correlation between the two. However, taking into account of additional electric field vector components is a simple way to generalize the vibrational Stark theory for the nitrile chromophore. Also, the electrostatic potential calculation method has been proposed and examined in detail. It turned out that the interactions of water molecules with nitrogen atom's lone pair orbital and with nitrile pi orbitals can be well described by the electrostatic potential calculation method. The present computational results will be of use to quantitatively simulate various linear and nonlinear vibrational spectra of nitrile compounds in solutions.

Vibrational dynamics of DNA: IV. Vibrational spectroscopic characteristics of A-, B-, and Z-form DNA's.

Linear and nonlinear IR spectroscopic studies of nucleic acids can provide crucial information on solution conformations of DNA double helix and its complex with other molecules. Carrying out density functional theory calculations of A-, B-, and Z-form DNA's, the authors obtained vibrational spectroscopic properties as well as coupling constants between different basis modes. The vibrational couplings that determine the extent of exciton delocalization are strongly dependent on DNA conformation mainly because the interlayer distance between two neighboring base pairs changes with respect to the DNA conformation. The Z-DNA has comparatively small interlayer vibrational coupling constants so that its vibrational spectrum depends little on the number of base pairs, whereas the A-DNA shows a notable dependency on the size. Furthermore, it is shown that a few distinctively different line shape changes in both IR and two-dimensional IR spectra as the DNA conformation changes from B to A or from B to Z can be used as marker bands and characteristic features distinguishing different DNA conformations.

Vibrational dynamics of DNA. II. Deuterium exchange effects and simulated IR absorption spectra.

In Paper I, we studied vibrational properties of normal bases, base derivatives, Watson-Crick base pairs, and multiple layer base pair stacks in the frequency range of 1400-1800 cm(-1). However, typical IR absorption spectra of single- and double-stranded DNA have been measured in D(2)O solution. Consequently, the more relevant bases and base pairs are those with deuterium atoms in replacement with labile amino hydrogen atoms. Thus, we have carried out density functional theory vibrational analyses of properly deuterated bases, base pairs, and stacked base pair systems. In the frequency range of interest, both aromatic ring deformation modes and carbonyl stretching modes appear to be strongly IR active. Basis mode frequencies and vibrational coupling constants are newly determined and used to numerically simulate IR absorption spectra. It turns out that the hydration effects on vibrational spectra are important. The numerically simulated vibrational spectra are directly compared with experiments. Also, the (18)O-isotope exchange effect on the poly(dG):poly(dC) spectrum is quantitatively described. The present calculation results will be used to further simulate two-dimensional IR photon echo spectra of DNA oligomers in the companion Paper III.

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings.

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra.

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

Phenol-benzene complexation dynamics: quantum chemistry calculation, molecular dynamics simulations, and two dimensional IR spectroscopy.

Molecular dynamics (MD) simulations and quantum mechanical electronic structure calculations are used to investigate the nature and dynamics of the phenol-benzene complex in the mixed solvent, benzene/CCl4. Under thermal equilibrium conditions, the complexes are continuously dissociating and forming. The MD simulations are used to calculate the experimental observables related to the phenol hydroxyl stretching mode, i.e., the two dimensional infrared vibrational echo spectrum as a function of time, which directly displays the formation and dissociation of the complex through the growth of off-diagonal peaks, and the linear absorption spectrum, which displays two hydroxyl stretch peaks, one for the complex and one for the free phenol. The results of the simulations are compared to previously reported experimental data and are found to be in quite reasonable agreement. The electronic structure calculations show that the complex is T shaped. The classical potential used for the phenol-benzene interaction in the MD simulations is in good accord with the highest level of the electronic structure calculations. A variety of other features is extracted from the simulations including the relationship between the structure and the projection of the electric field on the hydroxyl group. The fluctuating electric field is used to determine the hydroxyl stretch frequency-frequency correlation function (FFCF). The simulations are also used to examine the number distribution of benzene and CCl4 molecules in the first solvent shell around the phenol. It is found that the distribution is not that of the solvent mole fraction of benzene. There are substantial probabilities of finding a phenol in either a pure benzene environment or a pure CCl4 environment. A conjecture is made that relates the FFCF to the local number of benzene molecules in phenol's first solvent shell.

Characteristic two-dimensional IR spectroscopic features of antiparallel and parallel beta-sheet polypeptides: simulation studies.

The antiparallel and parallel beta sheets are two of the most abundant secondary structures found in proteins. Although various spectroscopic methods have been used to distinguish these two different structures, the linear spectroscopic measurements could not provide incisive information for distinguishing an antiparallel beta sheet from a parallel beta sheet. After carrying out quantum-chemistry calculations and model simulations, we show that the polarization-controlled two-dimensional (2D) IR photon echo spectroscopy can be of critical use in distinguishing these two different beta sheets. Particularly, the ratio between the diagonal peak and the cross peak is found to be strongly dependent on the quasi-2D array of the amide I local-mode transition dipole vectors. The relative intensities of the cross peaks in the 2D difference spectrum of an antiparallel beta sheet are significantly larger than those of the diagonal peaks, whereas the cross-peak amplitudes in the 2D difference spectrum of a parallel beta sheet are much weaker than the main diagonal-peak amplitudes. A detailed discussion on the origin of the diagonal- and cross-peak intensity distributions of both the antiparallel and parallel beta sheets is presented by examining vibrational exciton delocalization, relative angles between two different normal-mode transition dipoles, and natures of the cross peaks in the 2D difference spectrum.