Molecular dynamics simulations of cyclohexyl modified peptide nucleic acids (PNA)

Parameterization of the miniPEG-Modified γPNA Backbone: Toward Induced γPNA Duplex Dissociation

γ-Modified peptide nucleic acids (γPNAs) serve as potential therapeutic agents against genetic diseases. Miniature poly(ethylene glycol) (miniPEG) has been reported to increase solubility and binding affinity toward genetic targets, yet details of γPNA structure and dynamics are not understood. Within our work, we parameterized missing torsional and electrostatic terms for the miniPEG substituent on the γ-carbon atom of the γPNA backbone in the CHARMM force field. Microsecond timescale molecular dynamics simulations were carried out on six miniPEG-modified γPNA duplexes from NMR structures (PDB ID: 2KVJ). Three NMR models for the γPNA duplex (PDB ID: 2KVJ) were simulated as a reference for structural and dynamic changes captured for the miniPEG-modified γPNA duplex. Principal component analysis performed on the γPNA backbone atoms identified a single isotropic conformational substate (CS) for the NMR simulations, whereas four anisotropic CSs were identified for the ensemble of miniPEG-modified γPNA simulations. The NMR structures were found to have a 23° helical bend toward the major groove, consistent with our simulated CS structure of 19.0°. However, a significant difference between simulated methyl- and miniPEG-modified γPNAs involved the opportunistic invasion of miniPEG through the minor and major groves. Specifically, hydrogen bond fractional analysis showed that the invasion was particularly prone to affect the second G-C base pair, reducing the Watson-Crick base pair hydrogen bond by 60% over the six simulations, whereas the A-T base pairs decreased by only 20%. Ultimately, the invasion led to base stack reshuffling, where the well-ordered base stacking was reduced to segmented nucleobase stacking interactions. Our 6 μs timescale simulations indicate that duplex dissociation suggests the onset toward γPNA single strands, consistent with the experimental observation of decreased aggregation. To complement the insight of miniPEG-modified γPNA structure and dynamics, the new miniPEG force field parameters allow for further exploration of such modified γPNA single strands as potential therapeutic agents against genetic diseases.

Perspectives on conformationally constrained peptide nucleic acid (PNA): insights into the structural design, properties and applications

Peptide nucleic acid or PNA is a synthetic DNA mimic that contains a sequence of nucleobases attached to a peptide-like backbone derived from N-2-aminoethylglycine. The semi-rigid PNA backbone acts as a scaffold that arranges the nucleobases in a proper orientation and spacing so that they can pair with their complementary bases on another DNA, RNA, or even PNA strand perfectly well through the standard Watson-Crick base-pairing. The electrostatically neutral backbone of PNA contributes to its many unique properties that make PNA an outstanding member of the xeno-nucleic acid family. Not only PNA can recognize its complementary nucleic acid strand with high affinity, but it does so with excellent specificity that surpasses the specificity of natural nucleic acids and their analogs. Nevertheless, there is still room for further improvements of the original PNA in terms of stability and specificity of base-pairing, direction of binding, and selectivity for different types of nucleic acids, among others. This review focuses on attempts towards the rational design of new generation PNAs with superior performance by introducing conformational constraints such as a ring or a chiral substituent in the PNA backbone. A large collection of conformationally rigid PNAs developed during the past three decades are analyzed and compared in terms of molecular design and properties in relation to structural data if available. Applications of selected modified PNA in various areas such as targeting of structured nucleic acid targets, supramolecular scaffold, biosensing and bioimaging, and gene regulation will be highlighted to demonstrate how the conformation constraint can improve the performance of the PNA. Challenges and future of the research in the area of constrained PNA will also be discussed.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.

Improved Force Fields for Peptide Nucleic Acids with Optimized Backbone Torsion Parameters

Peptide nucleic acids are promising nucleic acid analogs for antisense therapies as they can form stable duplex and triplex structures with DNA and RNA. Computational studies of PNA-containing duplexes and triplexes are an important component for guiding their design, yet existing force fields have not been well validated and parametrized with modern computational capabilities. We present updated CHARMM and Amber force fields for PNA that greatly improve the stability of simulated PNA-containing duplexes and triplexes in comparison with experimental structures and allow such systems to be studied on microsecond time scales. The force field modifications focus on reparametrized PNA backbone torsion angles to match high-level quantum mechanics reference energies for a model compound. The microsecond simulations of PNA-PNA, PNA-DNA, PNA-RNA, and PNA-DNA-PNA complexes also allowed a comprehensive analysis of hydration and ion interactions with such systems.

Peptide nucleic acids: Advanced tools for biomedical applications

•

Peptide Nucleic Acids − DNA/RNA analogues.

•

Different Modifications on PNA backbone and their effects.

•

Neutral backbone − remarkable hybridization properties.

•

PNA based biosensors and their diverse biomedical applications.

•

Potential antigene and antisense agents.

Peptide Nucleic Acids (PNAs) are the DNA/RNA analogues in which sugar-phosphate backbone is replaced by N-2-aminoethylglycine repeating units. PNA contains neutral backbone hence due to the absence of electrostatic repulsion, its hybridization shows remarkable stability towards complementary oligonucleotides. PNAs are highly resistant to cleavage by chemicals and enzymes due to the substrate specific nature of enzymes and therefore not degraded inside the cells. PNAs are emerging as new tools in the market due to their applications in antisense and antigene therapies by inhibiting translation and transcription respectively. Hence, several methods based on PNAs have been developed for designing various anticancer and antigene drugs, detection of mutations or modulation of PCR reactions. The duplex homopurine sequence of DNA may also be recognized by PNA, forming firm PNA/DNA/PNA triplex through strand invasion with a looped-out DNA strand. PNAs have also been found to replace DNA probes in varied investigative purposes. There are several disadvantages regarding cellular uptake of PNA, so modifications in PNA backbone or covalent coupling with cell penetrating peptides is necessary to improve its delivery inside the cells. In this review, hybridization properties along with potential applications of PNA in the field of diagnostics and pharmaceuticals are elaborated.

Conformational Dynamics of Cyanocobalamin and Its Conjugates with Peptide Nucleic Acids

Vitamin B₁₂ also called cobalamin (Cbl) is an important enzymatic cofactor taken up by mammalian and also by many bacterial cells. Peptide nucleic acid (PNA) is a synthetic DNA analogue that has the ability to bind in a complementary manner to natural nucleic acids. Provided that PNA is efficiently delivered to cells, it could act as a steric blocker of functional DNA or RNA and regulate gene expression at the level of transcription or translation. Recently, Cbl has been examined as a transporter of various molecules to cells. Also, PNA, if covalently linked with Cbl, can be delivered to bacterial cells, but it is crucial to verify that Cbl does not change the desired PNA biological properties. We have analyzed the structure and conformational dynamics of conjugates of Cbl with a PNA monomer and oligomer. We synthesized a cyanocobalamin derivative with a PNA monomer C connected via the triazole linker and determined its NMR spectra. Using microsecond-long molecular dynamics simulations, we examined the internal dynamics of cyanocobalamin-C, its conjugate with a 14-mer PNA, and free PNA. The results suggest that all compounds acquire rather compact structures but the PNA oligomer conformations vary. For the Cbl-C conjugate the cross-peaks from the ROESY spectrum corroborated with the clusters from molecular dynamics trajectories. Within PNA the dominant interaction is stacking but the stacking bases are not necessarily neighboring in the PNA sequence. More bases stack in free PNA than in PNA of the conjugate, but stacking is less stable in free PNA. PNA in the conjugate is slightly more exposed to solvent. Overall, cyanocobalamin attached to a PNA oligomer increases the flexibility of PNA in a way that could be beneficial for its hybridization with natural nucleic acid oligomers.

Effect of backbone chemistry on hybridization thermodynamics of oligonucleic acids: a coarse-grained molecular dynamics simulation study

In this paper we study how varying oligonucleic acid backbone chemistry affects the hybridization/melting thermodynamics of oligonucleic acids. We first describe the coarse-grained (CG) model with tunable parameters that we developed to enable the study of both naturally occurring oligonucleic acids, such as DNA, and their chemically-modified analogues, such as peptide nucleic acids (PNAs) and locked nucleic acids (LNAs). The DNA melting curves obtained using such a CG model and molecular dynamics simulations in an implicit solvent and with explicit ions match with the melting curves obtained using the empirical nearest-neighbor models. We use these CG simulations to then elucidate the effect of backbone flexibility, charge, and nucleobase spacing along the backbone on the melting curves, potential energy and conformational entropy change upon hybridization and base-pair hydrogen bond residence time. We find that increasing backbone flexibility decreases duplex thermal stability and melting temperature mainly due to increased conformational entropy loss upon hybridization. Removing charges from the backbone enhances duplex thermal stability due to the elimination of electrostatic repulsion and as a result a larger energetic gain upon hybridization. Lastly, increasing nucleobase spacing decreases duplex thermal stability due to decreasing stacking interactions that are important for duplex stability.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

Molecular dynamics simulations of PNA-PNA and PNA-DNA duplexes by the use of new parameters implemented in the GROMACS package: a conformational and dynamics study

Peptide Nucleic Acids (PNAs) still represent a growing research area thanks to their potential applications in many fields of science from chemistry and biology to medicine. In these years, structural investigations by means of either experimental or computational techniques have proved to be very useful for the understanding of the structural organization and the binding properties of PNA. In this context, we here report an all-atoms Molecular Dynamics (MD) study of a PNA-PNA duplex and a PNA-DNA hetero-duplex with the well known GROMACS simulation package, by using new force field parameters properly derived for PNA molecules. The good agreement of our results with the crystallographic and NMR data, available for both the systems under investigation, confirms the validity of our approach. Moreover, our simulations reveal new interesting features related to the conformational-dynamic behavior of the studied systems, thus demonstrating the ability of MD simulations to gain insights into the dynamic properties of biologically relevant systems. This force field parametrization represents a good starting point for the implementation of a computational platform, based on the GROMACS package, useful for the rational design of modified PNA molecules with improved conformational features for selective binding toward DNA or RNA.

Conformational studies of chiral D-Lys-PNA and achiral PNA system in binding with DNA or RNA through a molecular dynamics approach

The growing interest in peptide nucleic acid (PNA) oligomers has led to the development of a very wide variety of PNA derivatives. Among others, the introduction of charged chiral groups on a PNA oligomer has proven effective in improving DNA binding ability, complexation direction and cellular uptake. In particular, the introduction of three adjacent chiral monomers based on D-Lys in the middle of the PNA sequence (D-Lys-PNA) has produced noteworthy results in modulating the directionality of the binding with the DNA complementary strand and in mismatch detection. Here, through a molecular dynamics approach, a comparative study has been carried out to investigate the structural properties that drive the interaction of the chiral D-Lys-PNA and the corresponding achiral PNA system with DNA as well as RNA complementary strands, starting from the crystal structure of D-Lys-PNA in complex with DNA. The results obtained complement experimental data and indicate that the binding with the RNA molecule, compared to DNA, is differently affected by the addition of three D-Lys groups on the PNA backbone, suggesting that this modification could be taken into account for the development of new PNA-based molecules able to discriminate between DNA and RNA.

Influence of PNA containing 8-aza-7-deazaadenine on structure stability and binding affinity of PNA·DNA duplex: insights from thermodynamics, counter ion, hydration and molecular dynamics analysis

This paper describes the synthesis of a novel 8-aza-7-deazapurin-2,6-diamine (DPP)-containing peptide nucleic acid (PNA) monomer and Boc protecting group-based oligomerization of PNA, replacing adenine (A) with DPP monomers in the PNA strand. The PNA oligomers were synthesized against the biologically relevant SV40 promoter region (2494-AATTTTTTTTATTTA-2508) of pEGFP-N3 plasmid. The DPP-PNA·DNA duplex showed enhanced stability as compared to normal duplex (A-PNA·DNA). The electronic distribution of DPP monomer suggested that DPP had better electron donor properties over 2,6-diamino purine. UV melting and thermodynamic analysis revealed that the PNA oligomer containing a diaminopyrazolo(3,4-d)pyrimidine moiety (DPP) stabilized the PNA·DNA hybrids compared to A-PNA·DNA. DPP-PNA·DNA duplex showed higher water activity (Δnw = 38.5) in comparison to A-PNA·DNA duplex (Δnw = 14.5). The 50 ns molecular dynamics simulations of PNA·DNA duplex containing DPP or unmodified nucleobase-A showed average H-bond distances in the DPP-dT base pair of 2.90 Å (OH-N bond) and 2.91 Å (NH-N bond), which were comparably shorter than in the A-dT base pair, in which the average distances were 3.18 Å (OH-N bond) and 2.97 Å (NH-N bond), and there was one additional H-bond in the DPP-dT base pair of around 2.98 Å (O2H-N2 bond), supporting the higher stability of DPP-PNA·DNA. The analysis of molecular dynamics simulation data showed that the system binding free energy increased at a rate of approximately -4.5 kcal mol(-1) per DPP base of the PNA·DNA duplex. In summary, increased thermal stability, stronger hydrogen bonding and more stable conformation in the DPP-PNA·DNA duplex make it a better candidate as antisense/antigene therapeutic agents.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.

Structural significance of hypermodified nucleic acid base hydroxywybutine (OHyW) which occur at 37th position in the anticodon loop of yeast tRNA(Phe)

Conformational preferences of hypermodified nucleic acid base hydroxywybutine (OHyW) have been studied using quantum chemical single point semi-empirical PM3 method. Automated geometry optimization using semi-empirical RM1, molecular mechanics force field (MMFF) along with ab-initio HF-SCF (6-31G** basis set) and DFT (B3LYP/6-31G** basis set) calculations have also been made to compare the salient features. Molecular electrostatic potentials (MEPs) depict the polarities of hydroxywybutine (OHyW) side chain. Another conformational study showed that hydroxywybutosine side chain interacts with adjacent bases within the anticodon loop of tRNA(Phe). The solvent accessible surface area (SASA) calculations revealed the structural role of hydroxywybutine in anticodon loop. Explicit molecular dynamics (MD) simulation has been done over the PM3 most stable structure of OHyW. The hydroxywybutine side chain prefers 'distal' conformation i.e. spreads away from the cyclic five membered imidazole moiety of modified tricyclic guanine base. The predicted preferred conformation of hydroxywybutine may prevent extended Watson-Crick base pairing during protein biosynthesis process. This conformation of OHyW stabilized by intramolecular interactions between O(6)⋯HO(16), O(6)⋯HC(15) and O(20)⋯HC(17). Further stabilization is also expected from interactions between O(22)⋯HC(16) and O(23)⋯HC(15). Explicit molecular dynamics (MD) simulation over the PM3 most stable structure of OHyW support the preferred geometry by preserving the 'distal' orientation of hydroxywybutine side chain and intramolecular hydrogen bonding interactions. MD simulation study revealed the role of hydroxyl group of OHyW to avoid fluctuations and prevent multiple iso-energetic conformations of hydroxywybutine side chain as compared to wybutine (yW).

Network models reveal stability and structural rearrangement of signal recognition particle

The signal recognition particle (SRP) and its receptors (SR) mediate the cotranslational targeting of the membrane and secretory proteins in all cells. In Escherichia coli, SRP is composed of the Ffh protein and the 4.5S SRP RNA. Ffh is a multidomain protein comprising a methionine-rich (M) domain, a helical N domain, and a Ras-like guanine triphosphatase (GTPase) (G) domain. The N and G domains are commonly referred to as one structural unit, the NG domain. In this article, the complex structure of SRP and SR is investigated with the Gaussian network model (GNM) and anisotropic network model (ANM). GNM provides the information of structure stability. It is found that the intermolecular interactions between SRP and SR can obviously decrease the fluctuation of NG domains. Nevertheless, the large structural rearrangement will take place during the cotranslational protein targeting cycle. Hence, the moving directions of fluctuation regions are further ascertained by using cross-correlation analysis and the ANM. The NG domain of Ffh undergoes a clockwise rotation around the GM linker and the M domain of Ffh shows an opposite direction to the NG domain. These functional movements will facilitate the SRP structure to transform into the free form and the sequence-bound form. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of protein assembly and supramolecular systems.

Mechanism of antisense oligonucleotide interaction with natural RNAs

Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.

Using computer-aided drug design and medicinal chemistry strategies in the fight against diabetes

The aim of this work is to present a simple, practical and efficient protocol for drug design, in particular Diabetes, which includes selection of the illness, good choice of a target as well as a bioactive ligand and then usage of various computer aided drug design and medicinal chemistry tools to design novel potential drug candidates in different diseases. We have selected the validated target dipeptidyl peptidase IV (DPP-IV), whose inhibition contributes to reduce glucose levels in type 2 diabetes patients. The most active inhibitor with complex X-ray structure reported was initially extracted from the BindingDB database. By using molecular modification strategies widely used in medicinal chemistry, besides current state-of-the-art tools in drug design (including flexible docking, virtual screening, molecular interaction fields, molecular dynamics, ADME and toxicity predictions), we have proposed 4 novel potential DPP-IV inhibitors with drug properties for Diabetes control, which have been supported and validated by all the computational tools used herewith.

Screening from the world's largest TCM database against H1N1 virus

The swine influenza virus (H1N1) 2009 pandemic highlights the importance of having effective anti-viral strategies. Recently, oseltamivir (Tamiflu) resistant influenza viruses are identified; which further emphasizes the urgency in developing new antiviral agents. In influenza virus replication cycle, viral surface glycoprotein, hemagglutinin, is responsible for viral entry into host cells. Hence, a potentially effective antiviral strategy is to inhibit viral entry mechanism. To develop novel antiviral agent that inhibits viral entry, we analyzed 20,000 traditional Chinese medicine (TCM) ingredients in hemagglutinin subtype H1 sialic acid binding site found on H1N1 virus. We then performed molecular dynamics simulations to investigate receptor-ligand interaction of the candidates obtained from docking. Here, we report three TCM derivatives that have high binding affinities to H1 sialic acid binding site residues based on structure-based calculations. The top three derivatives, xylopine_2, rosmaricine_14 and rosmaricine_15, all have an amine group that interact with Glu83 and a pyridinium group that interact with Asp103. Molecular dynamics simulations show that these derivatives form strong hydrogen bonding with Glu83 but interact transiently with Asp103. We therefore suggest that an enhanced hemagglutinin inhibitor, based on our scaffold, should be designed to bind both Glu83 and Asp103 with high affinity.

Thermal stability and unfolding pathways of Sso7d and its mutant F31A: insight from molecular dynamics simulation

The thermo-stability and unfolding behaviors of a small hyperthermophilic protein Sso7d as well as its single-point mutation F31A are studied by molecular dynamics simulation at temperatures of 300 K, 371 K and 500 K. Simulations at 300 K show that the F31A mutant displays a much larger flexibility than the wild type, which implies that the mutation obviously decreases the protein's stability. In the simulations at 371 K, although larger fluctuations were observed, both of these two maintain their stable conformations. High temperature simulations at 500 K suggest that the unfolding of these two proteins evolves along different pathways. For the wild-type protein, the C-terminal alpha-helix is melted at the early unfolding stage, whereas it is destroyed much later in the unfolding process of the F31A mutant. The results also show that the mutant unfolds much faster than its parent protein. The deeply buried aromatic cluster in the F31A mutant dissociates quickly relative to the wild-type protein at high temperature. Besides, it is found that the triple-stranded antiparallel β-sheet in the wild-type protein plays an important role in maintaining the stability of the entire structure.

Potent inhibitor design against H1N1 swine influenza: structure-based and molecular dynamics analysis for M2 inhibitors from traditional Chinese medicine database

The rapid spread of influenza virus subtype H1N1 poses a great threat to million lives worldwide. To search for new anti-influenza compounds, we performed molecular docking and molecular dynamics simulation to identify potential traditional Chinese medicine (TCM) constituents that could block influenza M2 channel activity. Quinic acid, genipin, syringic acid, cucurbitine, fagarine, and methyl isoferulate all have extremely well docking results as compared to control amantadine. Further de novo drug design suggests that derivatives of genipin and methyl isoferulate could have enhanced binding affinity towards M2 channel. Selected molecular dynamics simulations of M2-derivative complexes show stable hydrogen bond interactions between the derivatives and M2 residues, Ser10 and Ala9. To our best knowledge, this is the first study on the anti-viral activity of the above listed TCM compounds.