P53 dynamics upon response element recognition explored by molecular simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Terakawa T
PROVIDER: S-EPMC4656996 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Because of high stability and slow unfolding rates of G-quadruplexes (G4), cells have evolved specialized helicases that disrupt these non-canonical DNA and RNA structures in an ATP-dependent manner. One example is DHX36, a DEAH-box helicase, which participates in gene expression and replication by recognizing and unwinding parallel G4s. Here, we studied the molecular basis for the high affinity and specificity of DHX36 for parallel-type G4s using all-atom molecular dynamics simulations. By computing binding free energies, we found that the two main G4-interacting subdomains of DHX36, DSM and OB, separately exhibit high G4 affinity but they act cooperatively to recognize two distinctive features of parallel G4s: the exposed planar face of a guanine tetrad and the unique backbone conformation of a continuous guanine tract, respectively. Our results also show that DSM-mediated interactions are the main contributor to the binding free energy and rely on making extensive van der Waals contacts between the GXXXG motifs and hydrophobic residues of DSM and a flat guanine plane. Accordingly, the sterically more accessible 5'-G-tetrad allows for more favorable van der Waals and hydrophobic interactions which leads to the preferential binding of DSM to the 5'-side. In contrast to DSM, OB binds to G4 mostly through polar interactions by flexibly adapting to the 5'-terminal guanine tract to form a number of strong hydrogen bonds with the backbone phosphate groups. We also identified a third DHX36/G4 interaction site formed by the flexible loop missing in the crystal structure.
| S-EPMC8114077 | biostudies-literature
Project description:We have studied room-temperature structural and dynamic properties of the p53 DNA-binding domain in both DNA-bound and DNA-free states. A cumulative 55 ns of explicit solvent molecular dynamics simulations with the particle mesh Ewald treatment of electrostatics was performed. It was found that the mean structures in the production portions of the trajectories agree well with the crystal structure: backbone root-mean-square deviations are in the range of 1.6 and 2.0 A. In both simulations, noticeable backbone deviations from the crystal structure are observed only in loop L6, due to the lack of crystal packing in the simulations. More deviations are observed in the DNA-free simulation, apparently due to the absence of DNA. Computed backbone B-factor is also in qualitative agreement with the crystal structure. Interestingly, little backbone structural change is observed between the mean simulated DNA-bound and DNA-free structures. A notable difference is observed only at the DNA-binding interface. The correlation between native contacts and inactivation mechanisms of tumor mutations is also discussed. In the H2 region, tumor mutations at sites D281, R282, E285, and E286 may weaken five key interactions that stabilize H2, indicating that their inactivation mechanisms may be related to the loss of local structure around H2, which in turn may reduce the overall stability to a measurable amount. In the L2 region, tumor mutations at sites Y163, K164, E171, V173, L194, R249, I251, and E271 are likely to be responsible for the loss of stability in the protein. In addition to apparent DNA contacts that are related to DNA binding, interactions R175/S183, S183/R196, and E198/N235 are highly occupied only in the DNA-bound form, indicating that they are more likely to be responsible for DNA binding.
| S-EPMC2522240 | biostudies-literature
Project description:Neuraminidase (NA) of influenza is a key target for antiviral inhibitors, and the 150-cavity in group-1 NA provides new insight in treating this disease. However, NA of 2009 pandemic influenza (09N1) was found lacking this cavity in a crystal structure. To address the issue of flexibility of the 150-loop, Hamiltonian replica exchange molecular dynamics simulations were performed on different groups of NAs. Free energy landscape calculated based on the volume of 150-cavity indicates that 09N1 prefers open forms of 150-loop. The turn A (residues 147-150) of the 150-loop is discovered as the most dynamical motif which induces the inter-conversion of this loop among different conformations. In the turn A, the backbone dynamic of residue 149 is highly related with the shape of 150-loop, thus can function as a marker for the conformation of 150-loop. As a contrast, the closed conformation of 150-loop is more energetically favorable in N2, one of group-2 NAs. The D147-H150 salt bridge is found having no correlation with the conformation of 150-loop. Instead the intimate salt bridge interaction between the 150 and 430 loops in N2 variant contributes the stabilizing factor for the closed form of 150-loop. The clustering analysis elaborates the structural plasticity of the loop. This enhanced sampling simulation provides more information in further structural-based drug discovery on influenza virus.
| S-EPMC3622661 | biostudies-literature
Project description:Bcl-X(L), an antiapoptotic Bcl-2 family protein, plays a central role in the regulation of the apoptotic pathway. Heterodimerization of the antiapoptotic Bcl-2 family proteins with the proapoptotic family members such as Bad, Bak, Bim and Bid is a crucial step in the apoptotic regulation. In addition to these conventional binding partners, recent evidences reveal that the Bcl-2 family proteins also interact with noncanonical binding partners such as p53. Our previous NMR studies showed that Bcl-X(L): BH3 peptide and Bcl-X(L): SN15 peptide (a peptide derived from residues S15-N29 of p53) complex structures share similar modes of bindings. To further elucidate the molecular basis of the interactions, here we have employed molecular dynamics simulations coupled with MM/PBSA approach. Bcl-X(L) and other Bcl-2 family proteins have 4 hydrophobic pockets (p1-p4), which are occupied by four systematically spaced hydrophobic residues (h1-h4) of the proapoptotic Bad and Bak BH3 peptides. We observed that three conserved hydrophobic residues (F19, W23 and L26) of p53 (SN15) peptide anchor into three hydrophobic pockets (p2-p4) of Bcl-X(L) in a similar manner as BH3 peptide. Our results provide insights into the novel molecular recognition by Bcl-X(L) with p53.
| S-EPMC3198449 | biostudies-literature
Project description: Not available
| S-BSST248 | biostudies-other
Project description:Plants synthesize small molecule diterpenes composed of 20 carbons from precursor isopentenyl diphosphate and dimethylallyl disphosphate, manufacturing diverse compounds used for defense, signaling, and other functions. Industrially, diterpenes are used as natural aromas and flavoring, as pharmaceuticals, and as natural insecticides or repellents. Despite diterpene ubiquity in plant systems, it remains unknown how plants control diterpene localization and transport. For many other small molecules, plant cells maintain transport proteins that control compound compartmentalization. However, for most diterpene compounds, specific transport proteins have not been identified, and so it has been hypothesized that diterpenes may cross biological membranes passively. Through molecular simulation, we study membrane transport for three complex diterpenes from among the many made by members of the Lamiaceae family to determine their permeability coefficient across plasma membrane models. To facilitate accurate simulation, the intermolecular interactions for leubethanol, abietic acid, and sclareol were parametrized through the standard CHARMM methodology for incorporation into molecular simulations. To evaluate the effect of membrane composition on permeability, we simulate the three diterpenes in two membrane models derived from sorghum and yeast lipidomics data. We track permeation events within our unbiased simulations, and compare implied permeation coefficients with those calculated from Replica Exchange Umbrella Sampling calculations using the inhomogeneous solubility diffusion model. The diterpenes are observed to permeate freely through these membranes, indicating that a transport protein may not be needed to export these small molecules from plant cells. Moreover, the permeability is observed to be greater for plant-like membrane compositions when compared against animal-like membrane models. Increased permeability for diterpene molecules in plant membranes suggest that plants have tailored their membranes to facilitate low-energy transport processes for signaling molecules.
| S-EPMC9923751 | biostudies-literature
Project description:The tumor suppressor p53 is a master transcriptional regulator that affects a diverse range of cellular events. Surprisingly, even with >100 validated p53 response element (RE) sequences available, the effect of p53 binding on transcriptional behavior is seldom predictable and no functional rules have been described. Here, we report a systematic study on the role of specific nucleotides within the p53RE by using p21, a well-known target for p53 activation and contrasting it with Lasp1, a gene recently identified to be repressed by p53. Functional assays revealed a specific dinucleotide core combination within the CWWG motif of the p53RE to be the key factor that determines whether p53 transcriptionally activates or represses a target gene. The triplet RRR and YYY sequences flanking the core CWWG motif were also shown to play an important role in modulating the transcriptional behavior of p53. With the establishment of a set of predictive rules, we were able to reassess 162 published p53REs and showed that the attributed function for 20/162 p53REs studied were in fact erroneous. A significant proportion of p53REs (39/162) were found to be repressive, which is substantially higher than what is currently thought. Hence this clearer definition of the transcriptional behavior of p53 interaction with its RE will provide better insight toward the understanding of its fundamental role in cellular networks.
| S-EPMC2709670 | biostudies-literature
Project description:Chlorotoxin (CTX) is a 36⁻amino acid peptide with eight Cys residues that forms four disulfide bonds. It has high affinity for the glioma-specific chloride channel and matrix metalloprotease-2. Structural and binding properties of CTX analogs with various Cys residue substitutions with l-α-aminobutyric acid (Abu) have been previously reported. Using 4.2 µs molecular dynamics, we compared the conformational and essential space sampling of CTX and analogs with selective substitution of the Cys residues and associated disulfide bonds with either Abu or Ser. The native and substituted peptides maintained a high degree of α-helix propensity from residues 8 through 21, with the exception of substitution of the Cys⁵⁻Cys28 residues with Ser and the Cys16⁻Cys33 residues with Abu. In agreement with previous circular dichroism spectropolarimetry results, the C-terminal β-sheet content varied less from residues 25 through 29 and 32 through 36 and was well conserved in most analogs. The Cys16⁻Cys33 and Cys20⁻Cys35 disulfide-bonded residues appear to be required to maintain the αβ motif of CTX. Selective substitution with the hydrophilic Ser, may mitigate the destabilizing effect of Cys16⁻Cys33 substitution through the formation of an inter residue H-bond from Ser16:OγH to Ser33:OγH bridged by a water molecule. All peptides shared considerable sampled conformational space, which explains the retained receptor binding of the non-native analogs.
| S-EPMC6470725 | biostudies-literature
Project description:The combination of several closely spaced DNA lesions, which can be induced by a single radical hit, constitutes a hallmark in the DNA damage landscape and radiation chemistry. The occurrence of such a tandem base lesion gives rise to a strong coupling with the double helix degrees of freedom and induces important structural deformations, in contrast to DNA strands containing a single oxidized nucleobase. Although such complex lesions are known to be refractory to repair by DNA glycosylases, there is still a lack of structural evidence to rationalize these phenomena. In this contribution, we explore, by numerical modeling and molecular simulations, the behavior of the bacterial glycosylase responsible for base excision repair (MutM), specialized in excising oxidatively-damaged defects such as 7,8-dihydro-8-oxoguanine (8-oxoG). The difference in lesion recognition between a simple damage and a tandem lesion featuring an additional abasic site is assessed at atomistic resolution owing to microsecond molecular dynamics simulations and machine learning postprocessing, allowing to extensively pinpoint crucial differences in the interaction patterns of the damaged bases. Our results reveal substantial changes in the interaction network surrounding the 8-oxoG upon addition of an adjacent abasic site, leading to the perturbation of the intercalation triad which is crucial for lesion recognition and processing. The recognition process might also be impacted by a more constrained MutM-DNA binding upon tandem damage, as shown by the machine learning post-processing. This work advocates for the use of such high throughput numerical simulations for exploring the complex combinatorial chemistry of tandem DNA lesions repair and more generally local multiple damaged sites of the utmost significance in radiation chemistry.
| S-EPMC8141532 | biostudies-literature
Project description:The conformational landscapes of p53 peptide variants and phage derived peptide (12/1) variants, all known to bind to MDM2, are studied using hamiltonian replica exchange molecular dynamics simulations. Complementing earlier observations, the current study suggests that the p53 peptides largely follow the 'conformational selection' paradigm in their recognition of and complexation by MDM2 while the 12/1 peptides likely undergo some element of conformational selection but are mostly driven by 'binding induced folding'. This hypothesis is further supported by pulling simulations that pull the peptides away from their bound states with MDM2. This data extends the earlier mechanisms proposed to rationalize the entropically driven binding of the p53 set and the enthalpically driven binding of the 12/1 set. Using our hypothesis, we suggest mutations to the 12/1 peptide that increase its helicity in simulations and may, in turn, shift the binding towards conformational selection. In summary, understanding the conformational landscapes of the MDM2-binding peptides may suggest new peptide designs with bespoke binding mechanisms.
| S-EPMC5688104 | biostudies-literature