All-atom Molecular Dynamics Simulations Research Articles

Molecular simulations of doxorubicin complexed with native and modified cyclodextrins in water

In this study, all-atom molecular dynamics simulations have been performed to investigate from a nanoscopic point of view the interactions in supramolecular complexes formed between doxorubicin and cyclodextrins in an aqueous phase at a temperature of 310 K. All the simulations are on the time scales of microseconds and are carried out on High Performance Computing resources provided by the Greek Supercomputer ARIS. Doxorubicin, which is a well-known anticancer drug belonging to the class of anthracyclines, is simulated with two cyclodextrins of broad utilization in drug delivery, i.e.: γ-cyclodextrin and 2-hydroxypropyl-β-cyclodextrin. In both cases, the complexation process takes place spontaneously during the course of the simulations and a detailed analysis of the formed complexes is conducted in terms of structural, dynamical and thermodynamic properties. In addition, we present a quantitative analysis of how doxorubicin affects, upon complexation, the geometry, the hydrogen bonding network and the hydration of the selected cyclodextrins. Special attention is paid to the thermodynamic description of the simulated systems. More specifically, the standard binding free energy is calculated with an extended version of the linear interaction energy method. The presented in silico approach reveals that doxorubicin forms stable complexes with the chosen cyclodextrin molecules, there is a water solubility enhancement of the drug with respect to its unbound state and also cyclodextrins act protectively to the drug preventing undesired effects, such as the glycosidic degradation. Excellent agreement is found between the presented results and published in silico and experimental findings.

Growth kinetics of amyloid-like fibrils: An integrated atomistic simulation and continuum theory approach.

Amyloid fibrils have long been associated with many neurodegenerative diseases. The conventional picture of the formation and proliferation of fibrils from unfolded proteins comprises primary and secondary nucleation of oligomers followed by elongation and fragmentation thereof. In this work, we first employ extensive all-atom molecular dynamics (MD) simulations of short peptides to investigate the governing processes of fibril growth at the molecular scale. We observe that the peptides in the bulk solution can bind onto and subsequently diffuse along the fibril surface, which leads to fibril elongation via either bulk- or surface-mediated docking mechanisms. Then, to guide the quantitative interpretation of these observations and to provide a more comprehensive picture of the growth kinetics of single fibrils, a continuum model which incorporates the key processes observed in the MD simulations is formulated. The model is employed to investigate how relevant physical parameters affect the kinetics of fibril growth and identify distinct growth regimes. In particular, it is shown that fibrils which strongly bind peptides may undergo a transient exponential growth phase in which the entire fibril surface effectively acts as a sink for peptides. We also demonstrate how the relevant model parameters can be estimated from the MD trajectories. Our results provide compelling evidence that the overall fibril growth rates are determined by both bulk and surface peptide fluxes, thereby contributing to a more fundamental understanding of the growth kinetics of amyloid-like fibrils.

Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-fluororouracil targeting human pyruvate kinase m2

This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (k q) and the effective binding constants (K b) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be −118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2. Communicated by Ramaswamy H. Sarma

Molecular Mechanisms of the Generation and Accumulation of Gas at the Interface.

Gas-evolving reactions are widespread in chemical and energy fields. However, the generated gas will accumulate at the interface, which reduces the rate of gas generation. Understanding the microscopic processes of the generation and accumulation of gas at the interface is crucial for improving the efficiency of gas generation. Here, we develop an algorithm to reproduce the process of catalytic gas generation at the molecular scale based on the all-atom molecular dynamics simulations and obtain the quantitative evolution of the gas generation, which agrees well with the experimental results. In addition, we demonstrate that under an external electric field, the generated gas molecules do not accumulate at the electrode surface, which implies that the electric field can significantly increase the rate of the gas generation. The results suggest that the external electric field changes the structure of the water molecules near the electrode surface, making it difficult for gas molecules to accumulate on the electrode surface. Furthermore, it is found that gas desorption from the electrode surface is an entropy-driven process, and its accumulation at the electrode surface depends mainly on the competition between the entropy and the enthalpy of the water molecules under the influence of the electric field. These results provide deep insight into gas generation and inhibition of gas accumulation.

All-Atom Molecular Dynamics Simulations of Uncharged Linear Polymer Bottlebrushes: Effect of the Brush Sizes and the Number of Side-Chain Monomers.

Bottlebrush polymers (BBPs), characterized by grafted polymer side chains on linear backbone polymer chain, have emerged as a unique and versatile class of macromolecules with extensive applications in the fields of material science, electronics, battery materials, self-healing technology, etc. In this paper, we employ all-atom molecular dynamics (MD) simulations to present a comprehensive study of poly(methyl methacrylate)-g-poly(2-ethyl-2-oxazoline) (PMMA-g-PEtOx) BBP and its structural and hydration properties for varying number of backbone monomers (NBB) and side chain monomers (NSC), as well as properties of water molecules supported by the BBP. We find that the radius of gyration follows a scaling of Rg ∼NSC0.36 for smaller grafts and Rg ∼ NSC0.52-0.58 for longer grafts. We also find that the overall shape of the bottlebrush goes from a rod to sphere-like shape with the increase in NSC. Both the hydration per side chain monomer and hydrogen bonds (HBs) per oxygen and nitrogen of the side chain monomer reduce with an increase in NSC, caused by a corresponding enhancement in localization of the side chain monomers in the interior of the BBP. Furthermore, steric influences ensure the number of water-oxygen HBs is much more than the number of water-nitrogen HBs (with oxygen and nitrogen atoms belonging to the monomer side chains). Also, the BBP-supported water molecules demonstrate two distinctly ordered domains with one more structured and one less structured. The more structured domain disappears with an increase in NSC that causes more side chain monomers to localize in the interior of the BBPs. Finally, we observe that despite the highly negative partial charges of the oxygen and nitrogen atoms (of the side chain monomers), the dipole orientation distributions of water molecules around these atoms exhibit the presence of a neutral environment rather than an anionic environment. Overall, we anticipate that our study will generate significant interest in probing the various BBP systems in greater atomistic detail.

Iron uptake pathway of Escherichia coli as an entry route for peptide nucleic acids conjugated with a siderophore mimic.

Bacteria secrete various iron-chelators (siderophores), which scavenge Fe3+ from the environment, bind it with high affinity, and retrieve it inside the cell. After the Fe3+ uptake, bacteria extract the soluble iron(II) from the siderophore. Ferric siderophores are transported inside the cell via the TonB-dependent receptor system. Importantly, siderophore uptake paths have been also used by sideromycins, natural antibiotics. Our goal is to hijack the transport system for hydroxamate-type siderophores to deliver peptide nucleic acid oligomers into Escherichia coli cells. As siderophore mimics we designed and synthesized linear and cyclic Nδ-acetyl-Nδ-hydroxy-l-ornithine based peptides. Using circular dichroism spectroscopy, we found that iron(III) is coordinated by the linear trimer with hydroxamate groups but not by the cyclic peptide. The internal flexibility of the linear siderophore oxygen atoms and their interactions with Fe3+ were confirmed by all-atom molecular dynamics simulations. Using flow cytometry we found that the designed hydroxamate trimer transports PNA oligomers inside the E. coli cells. Growth recovery assays on various E. coli mutants suggest the pathway of this transport through the FhuE outer-membrane receptor, which is responsible for the uptake of the natural iron chelator, ferric-coprogen. This pathway also involves the FhuD periplasmic binding protein. Docking of the siderophores to the FhuE and FhuD receptor structures showed that binding of the hydroxamate trimer is energetically favorable corroborating the experimentally suggested uptake path. Therefore, this siderophore mimic, as well as its conjugate with PNA, is most probably internalized through the hydroxamate pathway.

Comparative analysis of structural dynamics and allosteric mechanisms of RecA/Rad51 family proteins: Integrated atomistic MD simulation and network-based analysis

Homologous recombination plays a key role in double-strand break repair, stalled replication fork repair, and meiosis. The RecA/Rad51 family recombinases catalyze the DNA strand invasion reaction that occurs during homologous recombination. However, the high sequence differences between homologous groups have hindered the thoroughly studies of this ancient protein family. The dynamic mechanisms of the family, particularly at the residual level, remain poorly understood. In this work, five representative RecA/Rad51 recombinase family members from all major kingdoms of living organisms: prokaryotes, eukaryotes, archaea, and viruses, were selected to explore the molecular mechanisms behind their conserved biological significance. A variety of techniques, including all-atom molecular dynamics simulation, perturbation response scanning, and protein structure network analysis, were used to examine the flexibility and correlation of protein domains, distribution of sensors and effectors and conserved hub residues. Furthermore, the potential communication routes between the ATP-binding region and the DNA-binding region of each recombinase were identified. Our results demonstrate the conserved molecular dynamics of these recombinases in the early stage of homologous recombination, including cooperative motions between regions, conserved sensing and effecting functional residue distribution, and conserved hub residues. Meanwhile, the unique ATP-DNA communication routes of each recombinase was also revealed. These results provide new insights into the mechanism of RecA/Rad51 family proteins, and provide new theoretical guidance for the development of allosteric inhibitors and the application of RecA/Rad51 family proteins.

Computational investigation on the conformational dynamics of C-terminal truncated α-synuclein bound to membrane

Accelerated progression rates in Parkinson’s disease (PD) have been linked to C-terminal domain (CTD) truncations of monomeric α-Synuclein (α-Syn), which have been suggested to increase amyloid aggregation in vivo and in vitro. In the brain of PD patients, CTD truncated α-Syn was found to have lower cell viability and tends to increase in the formation of fibrils. The CTD of α-Syn acts as a guard for regulating the normal functioning of α-Syn. The absence of the CTD may allow the N-terminal of α-Syn to interact with the membrane thereby affecting the normal functioning of α-Syn, and all of which will affect the etiology of PD. In this study, the conformational dynamics of CTD truncated α-Syn (1–99 and 1–108) monomers and their effect on the protein–membrane interactions were demonstrated using the all-atom molecular dynamics (MD) simulation method. From the MD analyses, it was noticed that among the two truncated monomers, α-Syn (1–108) was found to be more stable, shows rigidness at the N-terminal region and contains a significant number of intermolecular hydrogen bonds between the non-amyloid β-component (NAC) region and membrane, and lesser number of extended strands. Further, the bending angle in the N-terminal domain was found to be lesser in the α-Syn (1–108) in comparison with the α-Syn (1–99). Our findings suggest that the truncation on the CTD of α-Syn affects its interaction with the membrane and subsequently has an impact on the aggregation. Communicated by Ramaswamy H. Sarma

Ionomer structure and component transport in the cathode catalyst layer of PEM fuel cells: A molecular dynamics study.

The transport of water and protons in the cathode catalyst layer (CCL) of proton exchange membrane (PEM) fuel cells is critical for cell performance, but the underlying mechanism is still unclear. Herein, the ionomer structure and the distribution/transport characteristics of water and protons in CCLs are investigated via all-atom molecular dynamics simulations. The results show that at low water contents, isolated water clusters form in ionomer pores, while proton transport is mainly via the charged sites of the ionomer side chains and the Grotthuss mechanism. Moreover, with increasing water content, water clusters are interconnected to form continuous water channels, which provide effective paths for proton transfer via the vehicular and Grotthuss mechanisms. Increasing the ionomer mass content can enhance the dense arrangement of the ionomer, which, in turn, increases the density of charge sites and improves the proton transport efficiency. When the ionomer mass content is high, the clustering effect reduces the space for water diffusion, increases the proton transport path, and finally decreases the proton transport efficiency. By providing physics insights into the proton transport mechanism, this study is helpful for the structural design and performance improvement of CCLs of PEM fuel cells.

A complex network of interdomain interactions underlies the conformational ensemble of monomeric TDP-43 and modulates its phase behavior.

TAR DNA-binding protein 43 (TDP-43) is a multidomain protein involved in the regulation of RNA metabolism, and its aggregates have been observed in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Numerous studies indicate TDP-43 can undergo liquid-liquid phase separation (LLPS) in vitro and is a component of biological condensates. Homo-oligomerization via the folded N-terminal domain (aa:1-77) and the conserved helical region (aa:319-341) of the disordered, C-terminal domain is found to be an important driver of TDP-43 phase separation. However, a comprehensive molecular view of TDP-43 phase separation, particularly regarding the nature of heterodomain interactions, is lacking due to the challenges associated with its stability and purification. Here, we utilize all-atom and coarse-grained (CG) molecular dynamics (MD) simulations to uncover the network of interdomain interactions implicated in TDP-43 phase separation. All-atom simulations uncovered the presence of transient, interdomain interactions involving flexible linkers, RNA-recognition motif (RRM) domains and a charged segment of disordered C-terminal domain (CTD). CG simulations indicate these inter-domain interactions which affect the conformational landscape of TDP-43 in the dilute phase are also prevalent in the condensed phase. Finally, sequence and surface charge distribution analysis coupled with all-atom simulations (at high salt) confirmed that the transient interdomain contacts are predominantly electrostatic in nature. Overall, our findings from multiscale simulations lead to a greater appreciation of the complex interaction network underlying the structural landscape and phase separation of TDP-43.

Enhancing protein dynamics analysis with hydrophilic polyethylene glycol cross-linkers.

Cross-linkers play a critical role in capturing protein dynamics in chemical cross-linking mass spectrometry techniques. Various types of cross-linkers with different backbone features are widely used in the study of proteins. However, it is still not clear how the cross-linkers' backbone affect their own structure and their interactions with proteins. In this study, we systematically characterized and compared methylene backbone and polyethylene glycol (PEG) backbone cross-linkers in terms of capturing protein structure and dynamics. The results indicate the cross-linker with PEG backbone have a better ability to capture the inter-domain dynamics of calmodulin, adenylate kinase, maltodextrin binding protein and dual-specificity protein phosphatase. We further conducted quantum chemical calculations and all-atom molecular dynamics simulations to analyze thermodynamic and kinetic properties of PEG backbone and methylene backbone cross-linkers. Solution nuclear magnetic resonance was employed to validate the interaction interface between proteins and cross-linkers. Our findings suggest that the polarity distribution of PEG backbone enhances the accessibility of the cross-linker to the protein surface, facilitating the capture of sites located in dynamic regions. By comprehensively benchmarking with disuccinimidyl suberate (DSS)/bis-sulfosuccinimidyl-suberate(BS3), bis-succinimidyl-(PEG)2 revealed superior advantages in protein dynamic conformation analysis in vitro and in vivo, enabling the capture of a greater number of cross-linking sites and better modeling of protein dynamics. Furthermore, our study provides valuable guidance for the development and application of PEG backbone cross-linkers.

Water Exchange from the Buried Binding Sites of Cytochrome P450 Enzymes 1A2, 2D6, and 3A4 Correlates with Conformational Fluctuations.

Human cytochrome P450 enzymes (CYPs) are critical for the metabolism of small-molecule pharmaceuticals (drugs). As such, the prediction of drug metabolism by and drug inhibition of CYP activity is an important component of the drug discovery and design process. Relative to the availability of a wide range of experimental atomic-resolution CYP structures, the development of structure-based CYP activity models has been limited. To better characterize the role of CYP conformational fluctuations in CYP activity, we perform multiple microsecond-scale all-atom explicit-solvent molecular dynamics (MD) simulations on three CYP isoforms, 1A2, 2D6, and 3A4, which together account for the majority of CYP-mediated drug metabolism. The MD simulations employ a variety of positional restraints, ranging from keeping all CYP atoms close to their experimentally determined coordinates to allowing full flexibility. We find that, with full flexibility, large fluctuations in the CYP binding sites correlate with efficient water exchange from these buried binding sites. This is especially true for 1A2, which, when restrained to its crystallographic conformation, is unable to exchange water between the binding site and bulk solvent. These findings imply that, in addition to crystal structures, a representative ensemble of conformational states ought to be included when developing structure-based CYP activity models.

Interaction between Permeation Enhancers and Lipid Bilayers.

Permeation enhancers (PEs) are a class of molecules that interact with the epithelial membrane and transiently increase its transcellular permeability. Although there have been few clinical trials of PE coformulated drugs, the mechanism of action of PEs remains elusive. In this paper, the interaction between two archetypes of PEs [salcaprozate sodium (SNAC) and sodium caprate (C10)] and membranes is investigated with extensive all-atom molecular dynamics simulations. The simulations show that (1) the association between the neutral PEs and membranes is favored in free energy, (2) the propensity of neutral PE aggregation is larger in aqueous solution than in lipid bilayers, (3) the equilibrium distribution of neutral PEs in membranes is fast, e.g., accessible with unbiased MD simulations, and (4) the micelle of neutral PEs formed in aqueous solution does not rupture the membranes (e.g., not forming pores or breaking up the membrane) under simulation conditions. All results combined, this study indicates that PEs insert into the membranes in an equilibrium or near equilibrium process. This study lays the foundation for future investigations of how PEs impact the free energy of permeation for small molecules.

Tracking conformational transitions of the gonadotropin hormone receptors in a bilayer of (SDPC) poly-unsaturated lipids from all-atom molecular dynamics simulations.

Glycoprotein hormone receptors [thyrotropin (TSHR), luteinizing hormone/chorionic gonadotropin (LHCGR), and follicle stimulating hormone (FSHR) receptors] are rhodopsin-like G protein-coupled receptors. These receptors display common structural features including a prominent extracellular domain with leucine-rich repeats (LRR) stabilized by β-sheets and a long and flexible loop known as the hinge region (HR), and a transmembrane (TM) domain with seven α-helices interconnected by intra- and extracellular loops. Binding of the ligand to the LRR resembles a hand coupling transversally to the α- and β-subunits of the hormone, with the thumb being the HR. The structure of the FSH-FSHR complex suggests an activation mechanism in which Y335 at the HR binds into a pocket between the α- and β-chains of the hormone, leading to an adjustment of the extracellular loops. In this study, we performed molecular dynamics (MD) simulations to identify the conformational changes of the FSHR and LHCGR. We set up a FSHR structure as predicted by AlphaFold (AF-P23945); for the LHCGR structure we took the cryo-electron microscopy structure for the active state (PDB:7FII) as initial coordinates. Specifically, the flexibility of the HR domain and the correlated motions of the LRR and TM domain were analyzed. From the conformational changes of the LRR, TM domain, and HR we explored the conformational landscape by means of MD trajectories in all-atom approximation, including a membrane of polyunsaturated phospholipids. The distances and procedures here defined may be useful to propose reaction coordinates to describe diverse processes, such as the active-to-inactive transition, and to identify intermediaries suited for allosteric regulation and biased binding to cellular transducers in a selective activation strategy.

Elucidation of temperature-induced water structuring on cellulose surfaces for environmental and energy sustainability

Optimizing drying energy in the forest products industry is critical for integrating lignocellulosic feedstocks across all manufacturing sectors. Despite substantial efforts to reduce thermal energy consumption during drying, further enhancements are possible. Cellulose, the main component of forest products, is Earth's most abundant biopolymer and a promising renewable feedstock. This study employs all-atom molecular dynamics (MD) simulations to explore the structural dynamics of a small Iβ-cellulose microcrystallite and surrounding water layers during drying. Molecular and atomistic profiles revealed localized water near the cellulose surface, with water structuring extending beyond 8 Å into the water bulk, influencing solvent-accessible surface area and solvation energy. With increasing temperature, there was a ∼20 % reduction in the cellulose surface available for interaction with water molecules, and a ∼22 % reduction in solvation energy. The number of hydrogen bonds increased with thicker water layers, facilitated by a “bridging” effect. Electrostatic interactions dominated the intermolecular interactions at all temperatures, creating an energetic barrier that hinders water removal, slowing the drying processes. Understanding temperature-dependent cellulose-water interactions at the molecular level will help in designing novel strategies to address drying energy consumption, advancing the adoption of lignocellulosics as viable manufacturing feedstocks.

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations.

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

Mechanistic Insights into the Adsorption of Monoclonal Antibodies at the Water/Vapor Interface.

Monoclonal antibodies (mAbs) are active components of therapeutic formulations that interact with the water-vapor interface during manufacturing, storage, and administration. Surface adsorption has been demonstrated to mediate antibody aggregation, which leads to a loss of therapeutic efficacy. Controlling mAb adsorption at interfaces requires a deep understanding of the microscopic processes that lead to adsorption and identification of the protein regions that drive mAb surface activity. Here, we report all-atom molecular dynamics (MD) simulations of the adsorption behavior of a full IgG1-type antibody at the water/vapor interface. We demonstrate that small local changes in the protein structure play a crucial role in promoting adsorption. Also, interfacial adsorption triggers structural changes in the antibody, potentially contributing to the further enhancement of surface activity. Moreover, we identify key amino acid sequences that determine the adsorption of antibodies at the water-air interface and outline strategies to control the surface activity of these important therapeutic proteins.

Unfolding of protein using MoS2/SnS2 heterostructure for nanopore-based sequencing

Protein sequencing is crucial for understanding the complex mechanisms driving biological functions. However, proteins are usually folded in their native state and the mechanism of fast protein conformation transitions still remains unclear, which make protein sequencing challenging. Molecular dynamics simulations with accurate force field are now able to observe the entire folding/unfolding process, providing valuable insights into protein folding mechanisms. Given that proteins can be unfolded, nanopore technology shows great potential for protein sequencing. In this study, we proposed to use MoS2/SnS2 heterostructures to firstly unfold proteins and then detect them by a nanopore in the heterostructural membrane. All-atom molecular dynamics simulations performed in this work provided rich atomic-level information for a comprehensive understanding of protein unfolding process and mechanism on the MoS2/SnS2 heterostructure, it was found that the strong binding of protein to SnS2 nanostripe and hydrogen bond breaking were the main reasons for unfolding the protein on the heterostructure. After the protein was fully unfolded, it was restrained on the nanostripe because of the affinity of protein to the SnS2 nanostripe. Thus by integrating the proposed unfolding technique with nanopore technology, detection of linear unfolded peptide was realized in this work, allowing for the identification of protein components, which is essential for sequencing proteins in the near future.

Molecular Dynamics Simulations on Spike Protein Mutants Binding with Human β Defensin Type 2.

Human β defensin type 2 (hBD-2), a cationic cysteine-rich peptide secreted from the human innate immune system, can bind Spike-RBD at the same site as receptor ACE2, thus blocking viral entry into ACE2-expressing cells. In order to find out the impact of CoV-2 mutations on hBD-2's antiviral activity, it is important to investigate the binding and interaction of hBD-2 with RBD mutants. All-atom molecular dynamics simulations were conducted on typical RBD mutants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N, and N501Y-E484K-K417N, binding with hBD-2. Starting from the stable binding structure of hBD-2 and wt-RBD and ClusPro and HADDOCK docking-predicted initial structures, the RBD variants bound with hBD-2 simulations were set up, and NAMD simulations were conducted. Based on the structure and dynamics analysis, it was found that most RBD variants can still form a similar number of hydrogen bonds with hBD-2, in addition to having a similar-sized buried surface area (BSA) and a similar binding interface to the RBD wildtype. However, the RBD triple mutant formed a less stable binding structure with hBD-2 than other variants. Additionally, the free energy perturbation (FEP) method was applied to calculate the contribution of key mutant residues to the binding and the free energy change caused by the mutations. The result shows that N439K, K417N, and the trimutation increase the binding free energy of RBD with hBD-2; thus, RBD should bind less stably with hBD-2. E484K decreases the binding free energy, thus it should bind more stably with hBD-2, while N501Y, S477N, T478I, and P479S almost do not change the binding free energy with hBD-2. The MM-GBSA method predicted the binding interaction energy which shows that the trimutant should be able to escape the binding with hBD-2 but N501Y should not. The result can provide insight into understanding the functional mechanism of hBD-2 combating SARS-CoV-2 mutants.

Mechanism of water transport through the lipid membrane with trichogin GA IV. Molecular dynamics study

We perform all-atom molecular dynamics simulations of the antimicrobial peptide trichogin GA IV (TG) inside a model DOPC lipid membrane to clarify the mechanism of water transport increase. We considered the behavior of individual TG molecules inside the membrane as well as their aggregates: dimers, tetramers, and octamers. Our calculations show that a TG molecule moves into the membrane from the surrounding water and is localized in the lipid heads region. In this case, TG does not promote the permeation of water molecules through the membrane. The penetration of water into the hydrophobic region of the membrane occurred along the TG molecules that have formed transmembrane clusters (starting from dimers). In this case, the water molecules moved along the peptides by thermally induced hydrogen bond hopping allowing them to pass through the membrane. We show that the transmembrane clusters formed from non-helical TG molecules are more stable and provide increased permeability, since they can form additional hydrogen bonds, both with each other and with water, compared to the helical ones. We note that all considered TG clusters are compact and embedded in the membrane. They do not lead to a noticeable disruption of the membrane structure in their environment, and there are no pores or channels for the movement of water molecules, even in the case of an octamer. The movement of water molecules along all clusters occurs exclusively due to the hydrogen bond hopping mechanism.