Molecular dynamics simulation studies on influenza A virus H5N1 complexed with sialic acid and fluorinated sialic acid

Abbreviations HA Hemagglutinin

MD Molecular Dynamics

MM-PBSA Molecular Mechanics Poisson–Boltzmann Surface Area

NA Neuraminidase

NAMD Nanoscale Molecular Dynamic Simulation

PMEMD Particle Mesh Ewald Molecular Dynamics

RMSD Root-Mean-Square Deviation

RMSF Root-Mean-Square Fluctuation

SIA sialic acid

VMD Visual Molecular Dynamics

Communicated by Ramaswamy H. Sarma     

Exploration of interaction mechanism of tyrosol as a potent anti-inflammatory agent

Abstract

Drug discovery for a vigorous and feasible lead candidate is a challenging scientific mission as it requires expertise, experience, and huge investment. Natural products and their derivatives having structural diversity are renowned source of therapeutic agents since many years. Tyrosol (a natural phenylethanoid) has been extracted from olive oil, and its structure was confirmed by elemental analysis, FT-IR, FT-NMR, and single crystal X-ray crystallography. The conformational analysis for tyrosol geometry was performed by Gaussian 09 in terms of density functional theory. Validation of bond lengths and bond angles obtained experimentally as well as theoretically were performed with the help of curve fitting analysis, and values of correlation coefficient (R) obtained as 0.988 and 0.984, respectively. The charge transfer within the tyrosol molecule was confirmed by analysis of HOMO→LUMO molecular orbitals. In molecular docking with COX-2 (PDB ID: 5F1A), tyrosol was found to possess satisfactory binding affinity as compared to other NSAIDs (Aspirin, Ibuprofen, and Naproxen) and a COX-2 selective drug (Celecoxib). ADMET prediction, drug-likeness and bioactivity score altogether confirm the lead/drug like potential of tyrosol. Further investigation of simulation quality plot, RMSD and RMSF plots, ligands behavior plot as well as post simulation analysis manifest the consistency of 5F1A-tyrosol complex throughout the 20?ns molecular simulation process that signifies its compactness and stability within the receptor pocket. Abbreviations ADMET Absorption, Distribution, Metabolism, Excretion and Toxicity

Å Angstrom

COX-2 Cyclooxygenase-2

DFT Density Functional Theory

DMF Dimethylformamide

FMO Frontier Molecular Orbital

FT-IR Fourier-transform Infrared Spectroscopy

FT-NMR Nuclear Magnetic Resonance Spectroscopy

HOMO Highest Occupied Molecular Orbital

LUMO Lowest Unoccupied Molecular Orbital

MD Molecular Dynamics

NS Nanosecond

NSAIDs Non-steroidal anti-inflammatory drugs

OPE Osiris Property Explorer

RMSD Root-Mean-Square Deviation

RMSF Root Sean Square Fluctuation

Communicated by Ramaswamy H. Sarma     

Virtual identification of novel PPARα/γ dual agonists by 3D-QSAR,molecule docking and molecular dynamics studies

Abstract

Peroxisome proliferator-activated receptors (PPARs) are considered important targets for the treatment of Type 2 diabetes (T2DM). To accelerate the discovery of PPAR α/γ dual agonists, the comparative molecular field analysis (CoMFA) were performed for PPARα and PPARγ, respectively. Based on the molecular alignment, highly predictive CoMFA model for PPARα was obtained with a cross-validated q² value of 0.741 and a conventional r² of 0.975 in the non-cross-validated partial least-squares (PLS) analysis, while the CoMFA model for PPARγ with a better predictive ability was shown with q² and r² values of 0.557 and 0.996, respectively. Contour maps derived from the 3D-QSAR models provided information on main factors towards the activity. Then, we carried out structural optimization and designed several new compounds to improve the predicted biological activity. To investigate the binding modes of the predicted compounds in the active site of PPARα/γ, a molecular docking simulation was carried out. Molecular dynamic (MD) simulations indicated that the predicted ligands were stable in the active site of PPARα/γ. Therefore, combination of the CoMFA and structure-based drug design results could be used for further structural alteration and synthesis and development of novel and potent dual agonists. Abbreviations DM diabetes mellitus

T2DM type 2 diabetes

PPARs peroxisome proliferator-activated receptors

LBDD ligand based drug design

3D-QSAR three-dimensional quantitative structure activity relationship

CoMFA comparative molecular field analysis

PLS partial least square

LOO leave-one-out

q² cross-validated correlation coefficient

ONC optimal number of principal components

r² non-cross-validated correlation coefficient

SEE standard error of estimate

F the Fischer ratio

r²_pred predictive correlation coefficient

DBD DNA binding domain

MD molecular dynamics

RMSD root-mean-square deviation

RMSF root mean square fluctuations

Communicated by Ramaswamy H. Sarma     

Impact of Gln94Glu mutation on the structure and function of protection of telomere 1, a cause of cutaneous familial melanoma

Abstract

Protection of telomere 1 (POT1) is a key component of shelterin complex, essential for maintaining telomere length and its regulation. It consists of N-terminal domain (residues 1–299), which interacts with telomeric ssDNA, and the C-terminal domain (residues 320–634) that binds to the tripeptidyl-peptidase I (TPP1). A large number of naturally occurring mutations in the POT1 gene are associated with glioma, cardiac angiosarcoma and cutaneous familial melanoma (FM). In particular, Q94E mutation disrupts the interaction of POT1 with telomeric DNA which subsequently enhances telomere uncapping and elongation and promotes the development of cutaneous familial melanoma. To understand the underlying mechanism of familial melanoma developed by Q94E-mutation, we have performed extensive structure analysis of WT and mutant protein followed by molecular dynamics simulations. Q94E mutation causes a dramatic change in the structure and stability of POT1 protein. A considerable decrease in the flexibility, fluctuation and solvent accessibility of Q94E was observed in comparison to the WT, indicating overall destabilization of protein. Essential dynamics and Anisotropic Network Mode analysis have quantified a significant change in direction and magnitude of conformational motion in Q94E mutant compared to WT. A significant loss of frustration due to Q94E mutation was also observed. Our findings indicate the loss of protein stability and dynamics of POT1 protein by Q94E mutation may be associated with the familial melanoma. Abbreviations ANM anisotropic network mode

ED essential dynamics

FM familial melanoma

MD molecular dynamics

POT1 protection of telomere 1

Rg radius of gyration

RMSD root-mean-square deviation

RMSF root-mean-square fluctuations

SASA solvent accessible surface area

SIFT sorting Intolerant from Tolerant

TPP1 tripeptidyl-peptidase I

WT wild type

Communicated by Ramaswamy H. Sarma     

Scottish Sarcodon imbricatus under scrutiny

Summary

Molecular techniques applied to Scottish collections of Sarcodon imbricatus showed that these collections did not belong to this species but to the closely related S. squamosus. This latter species is accepted as a member of the British mycota.     

A combination of bioactive and nonbioactive alkyl-peptides form a more stable nanofiber structure for differentiating neural stem cells: a molecular dynamics simulation survey

Self-assembling alkyl-peptides are important molecules due to their ability to construct nano-level structures such as nanofibers to be utilized as tissue engineering scaffolds. The bioactive epitope of FAQRVPP which acts as neural stem cells (NSCs) outgrowth inducing factor is used in nanofiber structures. Based on previous experimental studies the density and distribution pattern of the epitopes on the surface of the nanofibers plays an important role in the differentiation function efficiency. We decided to survey and compare the stability of two pre-constructed fiber structures in the forms of all-functionalized nanofiber (containing only bioactive alkyl-peptides) and distributed functionalized nanofiber (a combination of nonbioactive and bioactive alkyl-peptides with ratio 2:1). Our findings reveal that the all-functionalized fiber shows an unstable structure and is split into intermediate micelle-like structures to reduce compactness and steric hindrance of functional epitopes whereas the distributed functionalized fiber shows an integrated stable nanofiber with a more amount of beta sheets that are well-organized and oriented around the hydrophobic core. The hydrogen bonds and energy profiles of the structures indicate the role of hydrophobic interactions during the alkyl-chain core formation and the important role of electrostatic interactions and hydrogen bond network in the stability of the final structures. Finally, it seems that the possibility of the presence of intermediate structure is increased in the all-functionalized nanofiber environment, and it can reduce functional efficiency of the scaffolds. These findings can help to design more efficient nanofiber structures with different goals in scaffolds for tissue engineering. Abbreviations MD Molecular Dynamics

NSCs Neural Stem Cells

PME Particle mesh Ewald

RDF Radial Distribution Function

RG Radius of gyration

RASA Relative Accessible Surface Area

RMSD Root Mean Square Deviations

SASA Solvent Accessible Surface Area.

Communicated by Ramaswamy H. Sarma     

iTRAQ-based proteomics reveals novel biomarkers of osteoarthritis

Abstract

Objective: We performed comprehensive proteomic analyses of articular cartilage by using the isobaric tags for relative and absolute quantitation (iTRAQ) method, and searched for candidate biomarkers for osteoarthritis (OA).

Methods: Articular cartilage was collected from patients with OA or femoral neck fracture for the control group. Molecular variations were detected by the iTRAQ method, and quantitative analyses were performed by western blot.

Results: Using the iTRAQ method, we identified 76 proteins with different expression levels in OA patients and the control group. Among these proteins, we selected LECT2 (leukocyte cell-derived chemotaxin-2), BAALC (brain and acute leukemia, cytoplasmic), and PRDX6 (peroxiredoxin-6), which had not been reported as biomarkers for OA.

Conclusions: Use of these proteins in combination with conventional OA biomarkers may better reflect the grade and prognosis of OA.     

How alike are the shapes of two random chains?

Two completely flexible random chains containing N atoms joined by links of length L can be superimposed by rigid-body motions to match their structures as closely as possible. The statistics of the root-mean-square best-fit distance D are discussed. For long chains, D² tends of 0.125L²N. Thus, for a protein of 100 residues, where the C^α-C^α distance is L = 3.8 Å, the expected root-mean-square atomic distance is 13.4 Å. This result gives a surprisingly good fit to the observed results for unrelated structures.     

Water velocity dependence of phosphate uptake on a coral-dominated fringing reef flat,La Réunion Island,Indian Ocean

Phosphate uptake (P-uptake) into coral reef communities has been hypothesized to be mass-transfer limited. One method of demonstrating mass-transfer limitation of P-uptake is to show dependence of P-uptake on water velocity. Water velocity across reef flats varies with tides and swell; thus, we measured P-uptake over the entire reef flat on eight different days, representing a range in water velocities. P-uptake was calculated from changes in P concentration of the water column. Changes in P concentration were measured by water sampling at six sites along a 300-m cross-reef transect while simultaneously measuring water velocity. To smooth the variability in phosphate concentrations, peristaltic pumps were used to get time-integrated water samples for 4–6 h at each site. Water velocities were measured in the middle of the transect using an acoustic Doppler current profiler and were averaged to match the time-integrated water sampling. Depth-averaged cross-reef water velocities were 0.031 ± 0.013 m s⁻¹ (mean ± SD), while the root-mean-square water velocities, accounting for oscillatory flow, averaged 3.3 times higher, 0.101 ± 0.021 m s⁻¹ (mean ± SD). Phosphate decreased along all transects. The first-order rate constant for P-uptake (S) was 8.5 ± 2.4 m d⁻¹ (mean ± SD) and increased linearly with root-mean-square water velocity. The Stanton number derived from oscillatory flow, the ratio of the first-order rate constant for P-uptake to the root-mean-square water velocity (S/U _rms), was (9.4 ± 1.2) × 10⁻⁴ (mean ± SD). P-uptake ranged from 0.2 to 1.1 mmol P m⁻² d⁻¹, demonstrating that P-uptake is variable on short time scales and is directly related to P concentration and water velocity.

     

Light scattering from wormlike chains with excluded volume effects

This paper reports a calculation of the angular dependence of light scattering from wormlike chains with excluded volume effects. The Daniels distribution function, modified for excluded volume effects, is used to compute averages for scattering elements separated by contour lengths which are long with respect to the persistence length of the chain. An expansion in terms of exactly known moments of the distribution for the wormlike coil without excluded volume is used for short contour lengths. The results are applied to scattering from calf thymus (M = 18.1 × 10⁶) and T7 (M = 25.4 × 10⁶) DNA. It is demonstrated that the same values of excluded volume parameter (ε = 0.11) and statistical segment length (1/λ′ = 900 Å) which explain the sedimentation and viscosity behavior of DNA also account satisfactorily for the scattering behavior. Molecular weights and root-mean-square radii estimated by extrapolation from scattering data obtained in the angular region from 10° to 25° will be 5–10% too large for DNA of molecular weight 20 × 10⁶–30 × 10⁶.     

Utility of AFLP markers for the assessment of molecular relationships in Ceratozamia Brongn. (Zamiaceae)

Abstract

The genus Ceratozamia has a distribution from Mexico to Central America. Molecular relationships were explored among Mexican Ceratozamia species using variable molecular markers, namely Amplified Fragment Length Polymorphisms (AFLP). Three AFLP primer-combinations generated 77 loci for a total of 1,684 fragments. Molecular variation was 75.72% among species and 24.28% within species, and a low correlation was present between genetic and geographical distances. Cluster analysis produced a phenogram in which two distinct clusters are clearly recognizable, one including C. alvarezii, C. sabatoi, C. zaragozae, C. kuesteriana, C. mexicana and C. robusta, and another one with C. zoquorum, C. miqueliana, C. latifolia, C. microstrobila, C. hildae and C. morettii. The AFLP proved to be suitable to study the relationships among species of Ceratozamia, and the results corroborated those from previous morphology-based studies.     

Ewald Summation in the Molecular Dynamics Simulation of Large Ionic Systems

Abstract

The accuracy and efficiency of the direct Ewald summation are discussed in terms of the size of a Molecular Dynamics (MD) ionic system and the ranges of the r-space and q-space summations. The dependence of the convergence parameter α on the size of the system and on the choice of cut-off radius for the short-range potential is given. The possibility of neglecting the q-space term for large ionic systems is discussed in terms of the accuracy and efficiency of the simulation.     

Molecular Dynamic Simulations of the N-Terminal Receiver Domain of NtrC Reveal Intrinsic Conformational Flexibility in the Inactive State

Abstract

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop β3-α3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop β3-α3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.     

Gibbs-Ensemble Molecular Dynamics: Liquid-Gas Equilibria for Lennard-Jones Spheres and n-Hexane

Abstract

We present a novel method to simulate phase equilibria in atomic and molecular systems. The method is a Molecular Dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles (or molecules) in a thermodynamically consistent fashion. Here we pres the derivation of the generalized equations of motion and discuss the relation of the resulting trajectory averages to the relevant ensemble. We test this Gibbs-Ensemble Molecular Dynamics algorithm by applying it to an atomic and a molecular system, i.e. to the liquid-gas coexistence in a Lennard-Jones fluid and in n-hexane. In both cases our results are in good accord with previous mean field and Gibbs-Ensemble Monte Carlo results as well as with the experimental data in the case of hexane. We also show that our Gibbs-Ensemble Molecular Dynamics algorithm like other Molecular Dynamics techniques can be used to study the dynamics of the system. Self-diffusion coefficients calculated with this method are in agreement with the result of conventional constant temperature Molecular Dynamics.     

On the Efficiency of Vectorized Molecular Dynamics Algorithms of Order N

Abstract

We discuss in this paper several factors which decide the efficiency of a vectorized Molecular Dynamics algorithm of order N, handling short-range interactions. General rules have been formulated and a fulfilment of these ought to ensure a very high speed of vector processing. The principles described make it possible to perform on presently available supercomputers the simulation of a 3-D system as large as 10⁵ (or even more) particles in rather modest cpu time.     

New heat shock protein (Hsp90) inhibitors,designed by pharmacophore modeling and virtual screening: synthesis,biological evaluation and molecular dynamics studies

Abstract

Inhibition of heat shock protein 90 (Hsp90) is known to be a significantly effective strategy in cancer therapy. Here, pyrazolopyranopyrimidine derivatives were characterized as new Hsp90 inhibitors. The molecules’ key structure (ZINC02819805) was determined by utilizing a pharmacophore model virtual screening workflow. Structural optimization was then carried out on compound ZINC02819805, pyrazolopyranopyrimidine derivatives were designed and six chosen derivatives were synthesized. The inhibition of Hsp90 ATPase activity of synthesized compounds revealed that para methylphenyl derivative of pyrazolopyranopyrimidine (HM3) was the most potent inhibitor (IC₅₀ = 5.5?µM). The anti-proliferative activity of this compound was evaluated against a panel of cell lines including MCF-7, HeLa and HUVEC (IC₅₀ = 1.28?µM, IC₅₀ = 1.74?µM and IC₅₀ = 61.48?µM respectively) by MTT method. The western blot analysis of treated MCF-7 cells with compound HM3 showed that the expression level of Hsp70 and Her2 proteins changed. The high level of Hsp70 expression and low level of Her2 expression suggest that compound HM3 exhibits inhibitory effect on Hsp90. Finally, the key interactions between HM3 and Hsp90 protein were studied by molecular dynamics simulation and showed that compound HM3 was stable in Hsp90 active cite during 200?ns simulation. Abbreviations Hsp90 Heat shock protein 90

MTT 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide

ATP adenosine triphosphate

MD molecular dynamics simulation

RMSD root-mean-square deviation

RMSF root-mean-square fluctuation

Rg gyration radius

m-SABNPs boehmite nanoparticles-supported sulfamic acid

Communicated by Ramaswamy H. Sarma     

An Extended Phosphate Linkage: Synthesis,Hybridization and Modeling Studies of Modified Oligonucleotides

Abstract

Novel stretched oligonucleotides (A-D) containing a 3′-α-C-methylene phosphodiester bridge (5-atoms long) have been synthesized on an automated synthesizer utilizing phosphoramidite chemistry. The key building-block 1-[3′-O-β-cyanoethyldiisopropylaminophosphiryl-2,3-dideoxy-5-O-dimethoxytriphenylmethyl-3-C-(hydroxymethyl)-β-D-erythro-pentofuranosyl]thymine (21) was prepared in a stereoselective manner from thymidine. Hybridization studies indicated a drop (1.8–3.0°C/mod.) in affinity for the complementary RNA and DNA targets. Molecular modeling results indicated that the 5-atom modified backbone had a different geometry around the phosphodiester linkage compared to the natural phosphodiester linkage. The stretched backbone may not be useful for antisense or triplex constructs, however it may find applications in biochemical/enzyme studies.

     

MD Simulations of Liquids with Td,Oh Molecular Symmetry

Abstract

A new method is proposed for the calculation of intermolecular interactions in Molecular Dynamics simulations of liquids with T_d, O_h molecular symmetry. The new algorithm is based on the separation of the pair potential into a short-range and a long-range contribution described by a site-site and a spherical centre-centre potential model respectively using an additional cutoff distance. Test calculations for the Lennard-Jones fluids CCl₄ and SF₆ show significant savings in CPU time. We compare thermodynamic properties, pair correlation functions and a few dynamic autocorrelation functions obtained with the novel strategy with results of the commonly used algorithm for systems containing 864 molecules. Since no significant differences appear the new algorithm may be suggested as a useful contribution to the area of Molecular Dynamics simulation of liquids with these rather high molecular symmetries.     

Location of F1 ATPase-like Genes on the Physical Map of the Bacillus subtilis 168 Chromosome

A water-soluble polysaccharide, FI, extracted from the mycelium of Granoderma tsugae, was fractionated and purified by ion-exchange chromatography, gel filtration, and affinity chromatography. Sixteen polysaccharides obtained were examined for antitumor effects on Sarcoma 180 in mice.

The three active polysaccharides obtained were as follows:

FI₀-a: A glycan-protein complex containing 9.3% protein, and having a hetero-glyco-chain of mannose and xylose.

FI₀-b-α: Molecular weight 10,000, glucan-protein complex containing 25.8% protein. The inhibition ratio was 61.8% against the solid cancer Sarcoma 180/mice; the survival ratio was more than 194% of the control group (100).

FA-1-b-α: Molcular weight 16,000, a complex of glycan: protein = 42:58 w/w, consisting of glucose as a main component, and associated with arabinose, mannose, xylose, and galactose. This had a tumor inhibition ratio of 56% and a survival ratio of more than 182%.

Comparison of active glycan with the fruiting body and mycelium: Among water-soluble polysaccharides of fruiting body, FI₀-a and FA-1, with antitumor activity, were both glucogalactan-protein complexes of molecular weight 10,000, but that of mycelium was a homoglucan protein complex in FI₀-b-α and heteroglucan protein in both FA-l-a and FA-l-b-α. The heteroglucan had a low tumor inhibition ratio, but caused a high survival ratio in mice.     

A New Species of Ctenomys (Rodentia: Ctenomyidae) from Patagonia Related to C. sociabilis

The genus Ctenomys includes a high number of taxa, with at least ten species from Patagonia and three recently described species for northeastern Chubut Province (Argentina). Ctenomys sociabilis is a social species of the genus Ctenomys and is currently distributed in the surrounding area of Sierra Cuyin Manzano (Neuquén Province), with a recently extinct population that occurred in Laguna Nahuelquir (Cushamen, Chubut Province). Molecular analyses have placed C. sociabilis at the base of Ctenomys clade, as the sister species to all other Ctenomys. Based on a morphological assessment (qualitative and quantitative) and DNA sequencing, we describe a new species of Ctenomys from Esquel, Chubut Province. Phylogenetic analysis shows the new species to be closely related to C. sociabilis, with evidence of solitary behavior. This new species is the first reported to be closely related phylogenetically to Ctenomys sociabilis at the base of the Ctenomys phylogeny. We provide anatomical comparisons between the new species and other species of Ctenomys from Patagonia, especially C. sociabilis.