Evaluation of mutations on O6-methylguanine methyl transferase structure and its interactions: molecular dynamics simulation study.

O6-methylguanine DNA methyl transferase (MGMT) is a significant vehicle for the cellular clearance of alkyl lesions, particularly the methyl group of the O-6 and O-4 positions of guanine and thymine, respectively. Many publications have studied the correlation between polymorphisms in MGMT and susceptibility to various cancers. In the present study, we investigated the consequence of L84F, common single-nucleotide polymorphism, K125E, site-specific mutagenesis, and L84F/K125E on conformation, stability, and behavior of MGMT in the free form and interaction with proliferating cell nuclear antigen (PCNA) and DNA as partners in the biochemical network by using molecular dynamics simulation method. Our results showed that all free variants of MGMT differed from the native form. However, among all free variants of MGMT, the L84F/K125E variant exhibited similar properties compared with the wild-type. In contrast, in complex modes, only amino acid residues of the L84F variant are involved in the interactions with PCNA and DNA somewhat differently relative to the wild-type. Furthermore, L84F SNP showed the highest binding free energy compared to other variants and native forms. These alterations in the amino acids and binding free energy of L84F relative to the native are the reasons for changing its region connection compared to the native form. Therefore, we propose conducting further investigations into the impact of inhibitors or chemotherapeutic agents to assess their effectiveness on MGMT variants compared to the wild-type, aiming to reduce the cost of cancer treatment that will depend on inhibiting native MGMT protein.Communicated by Ramaswamy H. Sarma.

Structural view of insulin adsorption on the multiple sizes of Cu nanoparticles; molecular dynamics simulation study.

Copper nanoparticles with different sizes cause significant changes in humans, and their interaction with biological molecules such as proteins is one of molecular biology's most critical and difficult challenges. One of the valuable tools for understating the mechanisms of different protein-nanoparticle interactions are molecular dynamics simulations. The present work attempts to understand how insulin interaction mechanisms Cu nanoparticles of different sizes (2, 6, and 10 nm) may occur with the molecular dynamics simulations methods. According to the results, in comparison to 2 and 6 nm nanoparticles, insulin interaction with the Cu 10 nm nanoparticle surface exhibits higher relative stability. The Van der Waals Forces (VDW) interactions show that by increasing the nanoparticle size from 2 nm to 6 nm (-673.32 and -847.83Kj/mol, respectively), the VDW energy leads to a significant increase. Although in the CU 10 nm, a noticeable decrease in VDW energy interaction was demonstrated (-357.21Kj/mol) due to present of three disulfide bond which act as a node that limits the excessive opening of insulin and another reason is the decline of surface electron density with increasing Cu-NP size. The secondary structure changes of insulin in the interaction with different sizes of Cu nanoparticles show the B-chain of insulin in the 2 and 6 nm nanoparticles systems had an essential role in the formation of interaction energy, in return for Cu 10 nm nanoparticle with the lowest energy level the chain-A has the bolder role. In the insulin-Cu 6 nm nanoparticle ß-sheet structure formation in the GlnB4, HisB5, LeuB6, GlyB20, GluB21, and ArgB22 residues was observed around 30 ns. Also, in these amino acids, a partial ß-sheet formation in B-chain has been apparent in the insulin-Cu 2 nm nanoparticle system. In the presence of the largest Cu nanoparticle in this study, a rise in helix, random coil, and B-bridge in secondary insulin structure was seen. It is considerable that in this system the most significant structural change occurring in the A-chain. The results of the hydrophilic-hydrophobicity properties and solvent accessible surface analyses also indicated that the insulin is unfolded during its interaction with different sizes of Cu nanoparticles.

Tracking the interaction between single-wall carbon nanotube and SARS-Cov-2 spike glycoprotein: A molecular dynamics simulations study.

COVID-19, a newly discovered type of coronavirus, is the cause of the pandemic infection that was first reported in Wuhan, China, in December 2019. One of the most critical problems in this regard is to identify innovative drugs that may reduce or manage this global health concern. Nanoparticles have shown a pivotal role in drug delivery systems in recent decades. The surface of nanoparticles could be covered by a layer composed of different biomolecules (e.g., proteins and macromolecules) following the incubation with a biological fluid. This protein-rich layer is called "Protein Corona." In this study, an all-atom molecular dynamics simulation was used for investigating the monomeric B domain of the spike glycoprotein due to its role in the accessibility of the spike glycoprotein to single-wall carbon nanotubes (SWCNTs). The interaction energy values between the carbon nanotube and B domain of the viral spike glycoprotein were evaluated. The obtained results, based on Lennard-Jones potentials, demonstrated that SWCNTs had an affinity to the B domain of the S1 subunit in the spike glycoprotein. The adsorption of SWCNTs on the B domain surface led to a significant change in solvent-accessible surface, internal hydrogen bonds, and finally in the tertiary structure, which could provide a reasonable method to impede the interaction between the angiotensin-converting enzyme II and SARS-CoV-2 spike glycoprotein. A decrease in the mean square displacement of the B domain was shown after the adsorption of SWCNTs as a result of increasing the hydrophobic-hydrophilic properties of the B domain. The arrangement of SWCNTs on the B domain surface and their interaction using the 2-acetamido-2-deoxy-ß-d-glucopyranose group (988, 991, and 992) demonstrated that a change in the affinity of the S1 subunit could be used as a barrier to viral replication. The analysis of the SWCNT-B domain complex indicated that the presence of SWCNTs is able to cause alterations in the S1 subunit of the spike protein, and these nanotubes could be employed for further in-vitro and in-vivo antiviral studies. Also, SWCNTs are able to be utilized in drug delivery systems.

In silico profiling and structural insights of zinc metal ion on O6-methylguanine methyl transferase and its interactions using molecular dynamics approach.

O6-methylguanine DNA methyl transferase (MGMT) is a metalloenzyme participating in the repair of alkylated DNA. In this research, we performed a comparative study for evaluating the impact of zinc metal ion on the behavior and interactions of MGMT in the both enzymatic forms of apo MGMT and holo MGMT. DNA and proliferating cell nuclear antigen (PCNA), as partners of MGMT, were utilized to evaluate molecular interactions by virtual microscopy of molecular dynamics simulation. The stability and conformational alterations of each forms (apo and holo) MGMT-PCNA, and (apo and holo) MGMT-DNA complexes were calculated by MM/PBSA method. A total of seven systems including apo MGMT, holo MGMT, free PCNA, apo MGMT-PCNA, holo MGMT-PCNA, apo MGMT-DNA, and holo MGMT-DNA complexes were simulated. In this study, we found that holo MGMT was more stable and had better folding and functional properties than that of apo MGMT. Simulation analysis of (apo and holo) MGMT-PCNA complexes displayed that the sequences of the amino acids involved in the interactions were different in the two forms of MGMT. The important amino acids of holo MGMT involved in its interaction with PCNA included E92, K101, A119, G122, N123, P124, and K125, whereas the important amino acids of apo MGMT included R128, R135, S152, N157, Y158, and L162. Virtual microscopy of molecular dynamics simulation showed that the R128 and its surrounding residues were important amino acids involved in the interaction of holo MGMT with DNA that was exactly consistent with X-ray crystallography structure. In the apo form of the protein, the N157 and its surrounding residues were important amino acids involved in the interaction with DNA. The binding free energies of - 387.976, - 396.226, - 622.227, and - 617.333 kcal/mol were obtained for holo MGMT-PCNA, apo MGMT-PCNA, holo MGMT-DNA, and apo MGMT-DNA complexes, respectively. The principle result of this research was that the area of molecular interactions differed between the two states of MGMT. Therefore, in investigations of metalloproteins, the metal ion must be preserved in their structures. Finally, it is recommended to use the holo form of metalloproteins in in vitro and in silico researches.

Molecular dynamics simulations of nano-confined methanol and methanol-water mixtures between infinite graphite plates: Structure and dynamics.

Molecular dynamics simulations are used to investigate microscopic structures and dynamics of methanol and methanol-water binary mixture films confined between hydrophobic infinite parallel graphite plate slits with widths, H, in the range of 7-20 Å at 300 K. The initial geometric densities of the liquids were chosen to be the same as bulk methanol at the same temperature. For the two narrowest slit widths, two smaller initial densities were also considered. For the nano-confined system with H = 7 Å and high pressure, a solid-like hexagonal arrangement of methanol molecules arranged perpendicular to the plates is observed which reflects the closest packing of the molecules and partially mirrors the structure of the underlying graphite structure. At lower pressures and for larger slit widths, in the contact layer, the methanol molecules prefer having the C-O bond oriented parallel to the walls. Layered structures of methanol parallel to the wall were observed, with contact layers and additional numbers of central layers depending on the particular slit width. For methanol-water mixtures, simulations of solutions with different composition were performed between infinite graphite slits with H = 10 and 20 Å at 300 K. For the nanoslit with H = 10 Å, in the solution mixtures, three layers of molecules form, but for all mole fractions of methanol, methanol molecules are excluded from the central fluid layer. In the nanopore with H = 20 Å, more than three fluid layers are formed and methanol concentrations are enhanced near the confining plates walls compared to the average solution stoichiometry. The self-diffusion coefficients of methanol and water molecules in the solution show strong dependence on the solution concentration. The solution mole fractions with minimal diffusivity are the same in confined and non-confined bulk methanol-water mixtures.

Simulations of structural and dynamic anisotropy in nano-confined water between parallel graphite plates.

We use molecular dynamics simulations to study the structure, dynamics, and transport properties of nano-confined water between parallel graphite plates with separation distances (H) from 7 to 20 Å at different water densities with an emphasis on anisotropies generated by confinement. The behavior of the confined water phase is compared to non-confined bulk water under similar pressure and temperature conditions. Our simulations show anisotropic structure and dynamics of the confined water phase in directions parallel and perpendicular to the graphite plate. The magnitude of these anisotropies depends on the slit width H. Confined water shows "solid-like" structure and slow dynamics for the water layers near the plates. The mean square displacements (MSDs) and velocity autocorrelation functions (VACFs) for directions parallel and perpendicular to the graphite plates are calculated. By increasing the confinement distance from H = 7 Å to H = 20 Å, the MSD increases and the behavior of the VACF indicates that the confined water changes from solid-like to liquid-like dynamics. If the initial density of the water phase is set up using geometric criteria (i.e., distance between the graphite plates), large pressures (in the order of ~10 katm), and large pressure anisotropies are established within the water. By decreasing the density of the water between the confined plates to about 0.9 g cm(-3), bubble formation and restructuring of the water layers are observed.