Molecular Dynamics Study of Supercoiled DNA Minicircles Tightly Bent and Supercoiled DNA in Atomistic Resolution

Towards the complete understanding of the DNA response to superhelical stress, sequence dependence structural disruptions on the ~100 base pairs supercoiled DNA minicircles were examined through a series of atomistic MD simulations. The results showed the effects from some subtle structural characteristics of DNA on defect formation, including flexibility at base pair step level and anisotropy, whose dynamic information are available only from atomistic MD simulations. For longer supercoiled DNA minicircles (240-340 bp), the molecules adapt into their writhed conformations. Writhe can be calculated by a Gauss’ integral performed along the DNA central axis path. A new mathematical definition for the DNA central axis path was developed for the more exact writhe calculation. Finally, atomistic representation of supercoiled 336 base pairs minicircles was provided by fitting the DNA structure obtained by explicitly solvated MD simulations into the density maps from electron cryo-tomography. Structural data were analysed and provided a decent explanation for the mechanism of the sequence specific binding of the enzyme topoisomerase 1B onto the negatively supercoiled DNA

Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks deter- mine the levels of transcription of a given gene at a given time. Using three complementary computa- tional techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are con- strained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any require- ment for external signalling factors. Our models com- bine atomistic molecular dynamics (MD) with coarse- grained simulations and statistical mechanical calcu- lations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denatura- tion in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informa- tional transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology

Roles of Tryptophan and Charged Residues on the Polymorphisms of Amyloids Formed by K-Peptides of Hen Egg White Lysozyme Investigated through Molecular Dynamics Simulations

Atomistic molecular dynamics simulations of amyloid models, consisting of the previously reported STDY-K-peptides and K-peptides from the hen egg white lysozyme (HEWL), were performed to address the effects of charged residues and pH observed in an in vitro study. Simulation results showed that amyloid models with antiparallel configurations possessed greater stability and compactness than those with parallel configurations. Then, peptide chain stretching and ordering were measured through the end-to-end distance and the order parameter, for which the amyloid models consisting of K-peptides and the STDY-K-peptides at pH 2 displayed a higher level of chain stretching and ordering. After that, the molecular mechanics energy decomposition and the radial distribution function (RDF) clearly displayed the importance of Trp62 to the K-peptide and the STDY-K-peptide models at pH 2. Moreover, the results also displayed how the negatively charged Asp52 disrupted the interaction networks and prevented the amyloid formation from STDY-K-peptide at pH 7. Finally, this study provided an insight into the interplay between pH conditions and molecular interactions underlying the formation of amyloid fibrils from short peptides contained within the HEWL. This served as a basis of understanding towards the design of other amyloids for biomaterial applications

<i>In Silico</i> Identification of Potential Sites for a Plastic-Degrading Enzyme by a Reverse Screening through the Protein Sequence Space and Molecular Dynamics Simulations

The accumulation of polyethylene terephthalate (PET) seriously harms the environment because of its high resistance to degradation. The recent discovery of the bacteria-secreted biodegradation enzyme, PETase, sheds light on PET recycling; however, the degradation efficiency is far from practical use. Here, in silico alanine scanning mutagenesis (ASM) and site-saturation mutagenesis (SSM) were employed to construct the protein sequence space from binding energy of the PETase–PET interaction to identify the number and position of mutation sites and their appropriate side-chain properties that could improve the PETase–PET interaction. The binding mechanisms of the potential PETase variant were investigated through atomistic molecular dynamics simulations. The results show that up to two mutation sites of PETase are preferable for use in protein engineering to enhance the PETase activity, and the proper side chain property depends on the mutation sites. The predicted variants agree well with prior experimental studies. Particularly, the PETase variants with S238C or Q119F could be a potential candidate for improving PETase. Our combination of in silico ASM and SSM could serve as an alternative protocol for protein engineering because of its simplicity and reliability. In addition, our findings could lead to PETase improvement, offering an important contribution towards a sustainable future

Molecular Mechanisms on the Selectivity Enhancement of Ascorbic Acid, Dopamine, and Uric Acid by Serine Oligomers Decoration on Graphene Oxide: A Molecular Dynamics Study

The selectivity in the simultaneous detection of ascorbic acid (AA), dopamine (DA), and uric acid (UA) has been an open problem in the biosensing field. Many surface modification methods were carried out for glassy carbon electrodes (GCE), including the use of graphene oxide and amino acids as a selective layer. In this work, molecular dynamics (MD) simulations were performed to investigate the role of serine oligomers on the selectivity of the AA, DA, and UA analytes. Our models consisted of a graphene oxide (GO) sheet under a solvent environment. Serine tetramers were added into the simulation box and were adsorbed on the GO surface. Then, the adsorption of each analyte on the mixed surface was monitored from MD trajectories. It was found that the adsorption of AA was preferred by serine oligomers due to the largest number of hydrogen-bond forming functional groups of AA, causing a 10-fold increase of hydrogen bonds by the tetraserine adsorption layer. UA was the least preferred due to its highest aromaticity. Finally, the role of hydrogen bonds on the electron transfer selectivity of biosensors was discussed with some previous studies. AA radicals received electrons from serine through hydrogen bonds that promoted oxidation reaction and caused the negative shifts and separation of the oxidation potential in experiments, as DA and UA were less affected by serine. Agreement of the in vitro and in silico results could lead to other in silico designs of selective layers to detect other types of analyte molecules

Atomistic Molecular Dynamics Simulations of DNA Minicircle Topoisomers: A Practical Guide to Setup, Performance and Analysis

While DNA supercoiling is ubiquitous in vivo, the structure of supercoiled DNA is more challenging to study experimentally than simple linear sequences because the DNA must have a closed topology in order to sustain superhelical stress. DNA minicircles, which are closed circular double-stranded DNA sequences typically containing between 60 and 500 base pairs, have proven to be useful biochemical tools for the study of supercoiled DNA mechanics. We present detailed protocols for constructing models of DNA minicircles in silico, for performing atomistic molecular dynamics (MD) simulations of supercoiled minicircle DNA, and for analyzing the results of the calculations. These simulations are computationally challenging due to the large system sizes. However, improvements in parallel computing software and hardware promise access to improve conformational sampling and simulation timescales. Given the concurrent improvements in the resolution of experimental techniques such as atomic force microscopy (AFM) and cryo-electron microscopy, the study of DNA minicircles will provide a more complete understanding of both the structure and the mechanics of supercoiled DNA

Effects of Site-Directed Mutations on the Communicability between Local Segments and Binding Pocket Distortion of Engineered GH11 Xylanases Visualized through Network Topology Analysis

Mutations occurred within the binding pocket of enzymes directly modified the interaction network between an enzyme and its substrate. However, some mutations affecting the catalytic efficiency occurred far from the binding pocket and the explanation regarding mechanisms underlying the transmission of the mechanical signal from the mutated site to the binding pocket was lacking. In this study, network topology analysis was used to characterize and visualize the changes of interaction networks caused by site-directed mutations on a GH11 xylanase from our previous study. For each structure, coordinates from molecular dynamics (MD) trajectory were obtained to create networks of representative atoms from all protein and xylooligosaccharide substrate residues, in which edges were defined between pairs of residues within a cutoff distance. Then, communicability matrices were extracted from the network to provide information on the mechanical signal transmission from the number of possible paths between any residue pairs or local protein segments. The analysis of subgraph centrality and communicability clearly showed that site-direct mutagenesis at non-reducing or reducing ends caused binding pocket distortion close to the opposite ends and created denser interaction networks. However, site-direct mutagenesis at both ends cancelled the binding pocket distortion, while enhancing the thermostability. Therefore, the network topology analysis tool on the atomistic simulations of engineered proteins could play some roles in protein design for the minimization to the correction of binding pocket tilting, which could affect the functionality and efficacy of enzymes

Effects of helix and fingertip mutations on the thermostability of xyn11A investigated by molecular dynamics simulations and enzyme activity assays

<p>Local conformational changes and global unfolding pathways of wildtype xyn11A recombinant and its mutated structures were studied through a series of atomistic molecular dynamics (MD) simulations, along with enzyme activity assays at three incubation temperatures to investigate the effects of mutations at three different sites to the thermostability. The first mutation was to replace an unstable negatively charged residue at a surface beta turn near the active site (D32G) by a hydrophobic residue. The second mutation was to create a disulphide bond (S100C/N147C) establishing a strong connection between an alpha helix and a distal beta hairpin associated with the thermally sensitive Thumb loop, and the third mutation add an extra hydrogen bond (A155S) to the same alpha helix. From the MD simulations performed, MM/PBSA energy calculations of the unfolding energy were in a good agreement with the enzyme activities measured from the experiment, as all mutated structures demonstrated the improved thermostability, especially the S100C/N147C proved to be the most stable mutant both by the simulations and the experiment. Local conformational analysis at the catalytic sites and the xylan access region also suggested that mutated xyn11A structures could accommodate xylan binding. However, the analysis of global unfolding pathways showed that structural disruptions at the beta sheet regions near the N-terminal were still imminent. These findings could provide the insight on the molecular mechanisms underlying the enhanced thermostability due to mutagenesis and changes in the protein unfolding pathways for further protein engineering of the GH11 family xylanase enzymes.</p