Search

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 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].

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].

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.

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.

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.

A case report. This study discusses a case of spinal segmental myoclonus caused by thoracic myelopathy, mimicking hiccup spasms. Spinal myoclonus caused by thoracic myelopathy is extremely rare. It can be misdiagnosed as chronic intractable hiccups due to similar clinical manifestations. Korea University Anam Hospital, Seoul, Republic of Korea. A 42-year-old man presented with a history of involuntary jerky movement of the upper abdominal wall muscles that had been continuing for over 3 years. A neurological examination, brain computed tomography and electroencephalogram did not reveal a cause of the symptoms. Electromyography was performed on the abdominal muscles and the findings revealed were compatible with spinal myoclonus. The spinal myoclonus had started in the abdominal muscles, with a spinal magnetic resonance imaging revealing a disc protrusion compressing the anterior spinal cord. The cause of the spinal myoclonus was determined to be spinal myelopathy due to mild T7 disc protrusion. The patient refused surgical or invasive interventions and was conservatively treated with clonazepam. The symptoms were reported to be less frequent following the treatment. Compressive myelopathy developed from disc protrusion may cause spinal myoclonus mimicking as hiccup spasms.

) have been theoretically investigated for the first time using molecular dynamic (MD) simulations and den-sity functional theory (MPWB1K/6-31G(d)//B3LYP/6-31G(d)) methods. The formation of the metal ion complex by the host is mainly driven by the electrostatic attraction between crown-5 oxygens and a cation together with the minor contribution of the cation-π interaction between two facing phenyl rings around the cation. The computed binding energies and the atomic charge distribution analysis for the metal binding complexes indicate the selectivity toward a potassium ion. The theoretical results herein explain the experimentally observed extractability order by this host towards various alkali metal ions. The physical nature and the driving forces for cation recognition by this host are discussed in detail.

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.

Amide I local mode frequencies and vibrational coupling constants in various multiple-stranded antiparallel β-sheet polyalanines are calculated by using the semiempirical calculation method and Hessian matrix reconstruction methods. The amide I local mode frequency is strongly dependent on the nature and number of hydrogen bonds. Vibrational couplings among amide I local modes in the multiple-stranded β-sheets are shown to be fully characterized by eight different coupling constants. The intrastrand coupling constants are found to be much smaller than the interstrand ones. Introducing newly defined inverse participation ratios and phase-correlation factors, the extent of two-dimensional delocalization and the vibrational phase relationship of amide I normal modes are elucidated. The A−E1 frequency splitting magnitude is shown to be strongly dependent on the number of strands but not on the length of each strand. A reduced one-dimensional Frenkel exciton model is used to describe the observed A−E1 frequenc...

Chain length and site dependencies of amide I local mode frequencies of α-helical polyalanines are theoretically studied by carrying out semiempirical quantum chemistry calculations. A theoretical model that can be used to quantitatively predict both the local amide I mode frequencies and coupling constants between two different local amide I modes is developed. Using this theoretical model and performing molecular dynamics simulation of an α-helical polyalanine in liquid water, we investigate conformational fluctuation and hydrogen-bonding dynamics by monitoring amide I frequency fluctuations. The instantaneous normal-mode analysis method is used to obtain densities of states of the one- and two-exciton bands and to quantitatively investigate the extent of delocalization of the instantaneous amide I normal modes. Also, by introducing a novel concept of the so-called weighted phase-correlation factor, the symmetric natures of the delocalized amide I normal modes are elucidated, and it is also shown that t...

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.

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.

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. alpha-Conotoxin neuropeptides are highly selective probes capable of discriminating different subtypes of nAChRs. In this study, we performed homology modeling to generate the neuronal alpha 3, beta 2 and beta 4 subunits using the x-ray structure of the alpha 1 subunit as a template. The structures of the extracellular domains containing ligand binding sites in the alpha 3 beta 2 and alpha 3 beta 4 nAChR subtypes were constructed using MD simulations and ligand docking processes in their free and ligand-bound states using alpha-conotoxin GIC, which exhibited the highest alpha 3 beta 2 vs. alpha 3 beta 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. [BMB Reports 2012; 45(5): 275-280]

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.

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.

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.

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.

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.

The amyloid-β (Aβ) proteins are responsible for amyloid plaques in Alzheimer's disease and have been the most widely studied subject in the process of fibril growth. Although much progress has been made to elucidate amyloid fibril properties at a molecular level, the full identification and characterization of all the conformational states and oligomeric structures in the aggregation process and all the conformational changes that link between those different states are still needed to be revealed. Here, we present the results of targeted molecular dynamics (TMD) simulations with explicit water to investigate the structural and mechanistic aspects of the association and the dissociation of the Aβ42 dimer. We will discuss the reversibility and the driving forces of the Aβ42 dimerization process with several order parameters along the protein aggregation pathway.

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.

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.

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.