Comparative simulation studies of native and single-site mutant human beta-defensin-1 peptides

Human defensins play important roles in a broad range of biological functions, such as microbial defense and immunity. Yet, little is known about their molecular properties, i.e. secondary structure stability, structural variability, important side chain interactions, surface charge distribution, and resistance to thermal fluctuations, and how these properties are related to their functions. To assess these factors, we studied the native human β-defensin-1 monomer and dimer as well as several single-site mutants using molecular dynamics simulations. The results showed that disulfide bonds are important determinants in maintaining the defensins’ structural integrity, as no structural transitions were observed at 300?K and only minor structural unfolding was detected upon heating to 500?K. The α-helix was less thermally stable than the core β-sheet structure held together by hydrogen bonds and hydrophobic interactions. The monomer α-helix stability was directly correlated, whereas the end-to-end distance was inversely correlated to the experimentally measured β-defensin-1 chemotactic activity, in the order: mutant 2 (Gln24Glu)?>?mutant 3 (Lys31Ala)?=?wild type?>?mutant 1 (Asn4Ala). The structural stability of the β-defensin-1 dimer species exhibited an inverse correlation to their chemotactic activity. In dimers formed by mutants 2 and 3, we observed sliding of one monomer upon the surface of the other in the absence of unbinding. This dynamic sliding feature may enhance the molecular oligomerization of β-defensin-1 peptides contributing to their antibacterial activity. It could also help these peptides orient correctly in the CC chemokine receptor 6 binding site, thereby initiating their chemotactic activity. In agreement with this notion, the remarkable sliding behavior was observed only for the mutants with the highest chemotactic activity.     

Molecular dynamics simulation of chitinase I from Thermomyces lanuginosus SSBP to ensure optimal activity

Abstract

The fungal chitinase I obtained from Thermomyces lanuginosus SSBP, a thermophilic deuteromycete, has an optimum growth temperature and pH of 323.15 K and 6.5, respectively. This enzyme plays an important task in the defence mechanism of organisms against chitin-containing parasites by hydrolysing β-1, 4-linkages in chitin. It acts as both anti-fungal and biofouling agents, with some being thermostable and suitable for the industrial applications. Three-dimensional model of chitinase I enzyme was predicted and analysed using various bioinformatics tools. The structure of chitinase I exhibited a well-defined TIM barrel topology with an eight-stranded α/β domain. Structural analysis and folding studies at temperatures ranging from 300 to 375 K using 10 ns molecular dynamics simulations clearly showed the stability of the protein was evenly distributed even at higher temperatures, in accordance with the experimental results. We also carried out a number of 20 ns constant pH molecular dynamics simulations of chitinase I at a pH range 2–6 in a solvent. This work was aimed at establishing the optimum activity and stability profiles of chitinase I. We observed a strong conformational pH dependence of chitinase I and the enzyme retained their characteristic TIM barrel topology at low pH.     

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg²⁺ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg²⁺ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na⁺ ions. The determining factor of Mg²⁺ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg²⁺ to DNA non electrostatic contributions play a major role.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5     

Transport diffusivities of fluids in nanopores by non-equilibrium molecular dynamics simulation

We present a method to study fluid transport through nanoporous materials using highly efficient non-equilibrium molecular dynamics simulations. A steady flow is induced by applying an external field to the fluid particles within a small slab of the simulation cell. The external field generates a density gradient between both sides of the porous material, which in turn triggers a convective flux through the porous medium. The heat dissipated by the fluid flow is released by a Gaussian thermostat applied to the wall particles. This method is effective for studying diffusivities in a slit pore as well as more natural, complex wall geometries. The dependence of the diffusive flux on the external field sheds light on the transport diffusivities and allows a direct calculation of effective diffusivities. Both pore and fluid particle interactions are represented by coarse-grained molecular models in order to present a proof-of-concept and to retain computational efficiency in the simulations. The application of the method is demonstrated in two different scenarios, namely the effective mass transport through a slit pore and the calculation of the effective self-diffusion through this system. The method allows for a distinction between diffusive and convective contributions of the mass transport.     

Spermine as a possible endogenous allosteric activator of carboxypeptidase A: multispectroscopic and molecular simulation studies

Abstract

Carboxypeptidase A (EC.3.4.17.1) is a zinc-containing proteolytic enzyme that removes the C-terminal amino acid from a peptide chain with the free carboxylate-terminal. In this study, the effect of spermine interaction on the structure and thermal stability of Carboxypeptidase A was investigated by ultraviolet???visible spectroscopy, fluorescence spectroscopy, circular dichroism, Kinetic measurement, molecular docking and simulation studies have also been followed at the pH of 7.5. The transition temperature of Carboxypeptidase A, as a criterion of protein thermal stability, in the presence of spermine was enhanced by increasing the concentration of spermine. The results of fluorescence intensity changes, at two temperatures of 308 and 318?K, also suggested that spermine had a great ability to quench the fluorescence of Carboxypeptidase A through the static quenching procedure. The thermodynamic parameters changes, including standard Gibbs free-energy, entropy and enthalpy, showed that the binding of spermine to Carboxypeptidase A was spontaneous and the hydrogen bonding and van der Waals interactions played a major role in stabilizing the Carboxypeptidase A–spermine complex. The changes in the content of the α-helix and the β-sheet of the Carboxypeptidase A with binding to spermine were shown by the CD spectra method. Further, kinetic studies revealed that by increasing concentration of spermine, the activity of Carboxypeptidase A was enhanced. Also, the docking study revealed that the hydrogen bonding and van der Waals interactions played a major role in stabilizing the Carboxypeptidase A–spermine complex. As a result, spermine could be considered as an activator and a stabilizer for Carboxypeptidase A.

Communicated by Ramaswamy H. Sarma     

Understanding the effects of two bound glucose in Sudlow site I on structure and function of human serum albumin: theoretical studies

Human serum albumin (HSA) is the most abundant protein found in blood serum. It carries essential metabolites and many drugs. The glycation of HSA causes abnormal biological effects. Importantly, glycated HSA (GHSA) is of interest as a biomarker for diabetes. Recently, the first HSA structure with bound pyranose (GLC) and open-chain (GLO) glucose at Sudlow site I has been crystallised. We therefore employed molecular dynamics (MD) simulations and ONIOM calculations to study the dynamic nature of two bound glucose in a pre-glycated HSA (pGHSA) and observe how those sugars alter a protein structure comparing to wild type (Apo) and fatty acid-bound HSA (FA). Our analyses show that the overall structural stability of pGHSA is similar to Apo and FA, except Sudlow site II. Having glucose induces large protein flexibility at Sudlow site II. Besides, the presence of glucose causes W214 to reorient resulting in a change in W214 microenvironment. Considering sugars, both sugars are exposed to water, but GLO is more solvent-accessible. ONIOM results show that glucose binding is favoured for HSA (?115.04 kcal/mol) and GLO (?85.10 kcal/mol) is more preferable for Sudlow site I over GLC (?29.94 kcal/mol). GLO can strongly react with K195 and K199, whereas K195 and K199 provide slightly repulsive forces for GLC. This can confirm that an open-chain GLO is more favourable inside a pocket.     

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.     

Inhibitor binding studies of Mycobacterium tuberculosis MraY (Rv2156c): Insights from molecular modeling,docking, and simulation studies

Abstract

Tuberculosis (TB) is a contagious disease caused by Mycobacterium tuberculosis (M.tb) or tubercule bacillus, and H37Rv is the most studied clinical strain. The recent development of resistance to existing drugs is a global health-care challenge to control and cure TB. Hence, there is a critical need to discover new drug targets in M.tb. The members of peptidoglycan biosynthesis pathway are attractive target proteins for antibacterial drug development. We have performed in silico analysis of M.tb MraY (Rv2156c) integral membrane protein and constructed the three-dimensional (3D) structure model of M.tb MraY based on homology modeling method. The validated model was complexed with antibiotic muraymycin D2 (MD2) and was used to generate structure-based pharmacophore model (e-pharmacophore). High-throughput virtual screening (HTVS) of Asinex database and molecular docking of hits was performed to identify the potential inhibitors based on their mode of interactions with the key residues involved in M.tb MraY–MD2 binding. The validation of these molecules was performed using molecular dynamics (MD) simulations for two best identified hit molecules complexed with M.tb MraY in the lipid bilayer, dipalmitoylphosphatidyl-choline (DPPC) membrane. The results indicated the stability of the complexes formed and retained non-bonding interactions similar to MD2. These findings may help in the design of new inhibitors to M.tb MraY involved in peptidoglycan biosynthesis.     

Hydrogen-bond dynamics of water in presence of an amphiphile,tetramethylurea: signature of confinement-induced effects

Various experimental and simulation studies have suggested that the presence of amphiphilic molecules in aqueous solutions substantially perturbs the tetrahedral hydrogen-bond (H-bond) network of neat liquid water. Such structural perturbation is expected to impact H-bond lifetime of liquid water. Tetramethylurea (TMU) is an example of an amphiphile because it possesses both hydrophobic and hydrophilic moieties. Molecular dynamics simulations of (water+TMU) binary mixtures at various compositions have been performed in order to investigate the microscopic mechanism through which the amphiphiles influence the H-bond dynamics of liquid water at room temperature. Present simulations indicate lengthening of both water–water H-bond lifetime and H-bond structural relaxation time upon addition of TMU in aqueous solution. At the highest TMU mole fraction studied, H-bond lifetime and structural relaxation time are, respectively, ～4 and ～8 times longer than those in neat water. This is comparable with the slowing down of H-bond dynamics for water molecules confined in cyclodextrin cavities. Simulated relaxation profiles are multi-exponential in character at all mixture compositions, and simulated radial distribution functions suggest enhanced water–water and water–TMU interactions upon addition of TMU. No evidence for complete encapsulation of TMU by water H-bond network has been found.     

The adsorption of glycated human serum albumin-selective aptamer onto a graphene sheet: simulation studies

The adsorption of glycated human serum albumin (gHSA)-selective DNA aptamer, where gHSA is a diabetes biomarker, on a mobile graphene in electrolyte and salt-free solutions has been investigated using molecular dynamics (MD) simulations. This work was done to benefit the development of diabetes graphene-based fluorescent aptasensor. For aptasensor, an aptamer is quenched upon GRA binding and becomes fluorescent when leaving and binding to its analyte (gHSA). Different DNA conformations are obtained due to the GRA flexibility. Apart from a single-side deposition, the clipping DNA conformation is obtained. Our results show the faster DNA adsorption under a salt condition. Landing onto GRA with different aptamer parts affects adhered DNA conformations. The clipping conformation is found when the mid of aptamer is firstly deposited on a GRA surface. In contrast, if both 3′ and 5′ termini arrive together, they act as legs to stand on a GRA sheet generating the U-loop like structure before collapsing. Like other studies, laying flat onto a GRA sheet is the most favourable. The high nucleobase-GRA contacts and highly water-exposed phosphate backbone are observed. This implies key roles of nucleobases for GRA adsorption and a phosphate backbone for binding analytes and triggering further desorption.     

Conformational flexibility of N-glycans in solution studied by REMD simulations

Protein–glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of “conformer selection” in protein–glycan recognition.     

Comparative study of stability and transport of molecules through cyclic peptide nanotube and aquaporin: a molecular dynamics simulation approach

Abstract

The structural stability and transport properties of the cyclic peptide nanotube (CPN) 8?×?[Cys–Gly–Met–Gly]₂ in different phospholipid bilayers such as POPA (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) and POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine) with water have been investigated using molecular dynamics (MD) simulation. The hydrogen bonds and non-bonded interaction energies were calculated to study the stability in different bilayers. One µs MD simulation in POPA lipid membrane reveals the stability of the cyclic peptide nanotube, and the simulations at various temperatures manifest the higher stability of 8?×?[Cys–Gly–Met–Gly]₂. We demonstrated that the presence of sulphur-containing amino acids in CPN enhances the stability through disulphide bonds between the adjacent rings. Further, the water permeation coefficient of the CPN is calculated and compared with human aquaporin-2 (AQP2) channel protein. It is found that the coefficients are highly comparable to the AQP2 channel though the mechanism of water transport is not similar to AQP 2; the flow of water in the CPN is taking place as a two-line 1–2–1–2 file fashion. In addition to that, the transport behavior of Na⁺ and K⁺ ions, single water molecule, urea and anti-cancer drug fluorouracil were investigated using pulling simulation and potential of mean force calculation. The above transport behavior shows that Na⁺ is trapped in CPN for a longer time than other molecules. Also, the interactions of the ions and molecules in Cα and mid-Cα plane were studied to understand the transport behavior of the CPN. Abbreviations AQP2 Aquaporin-2

CPN Cyclic peptide nanotube

MD Molecular dynamics

POPA 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid

POPE 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine

POPG 1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol

POPS 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine

Communicated by Ramaswamy H. Sarma     

The partition and transport behavior of cytotoxic ionic liquids (ILs) through the DPPC bilayer: Insights from molecular dynamics simulation

A molecular dynamics (MD) simulation with atomistic details was performed to examine the partitioning and transport behavior of moderately cytotoxic ionic liquids (ILs), namely choline bis(2-ethylhexyl) phosphate (CBEH), choline bis(2,4,4-trimethylpentyl) phosphinate (CTMP) and choline O,O-diethyl dithiophosphate (CDEP) in a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer in the fluid phase at 323?K. The structure of ILs was so selected to understand if the role of dipole and dispersion forces in the ILs distribution in the membrane can be possible. Several analyses including mass density, electrostatic potential, order parameter, diffusion coefficients and hydrogen bond formation, was carried out to determine the precise location of the anionic species inside the membrane. Moreover, the potential of the mean force (PMF) method was used to calculate free energy profile for transferring anionic species from the DPPC membrane into the bulk water. While less cytotoxic DEP is located within the bulk water, more cytotoxic TMP and BEH ILs were found to remain in the membrane and the energy barrier for crossing through the bilayer center of BEH was higher. Various ILs have no significant effect on P–N vector. The thickness of lipid bilayer decreased in all systems comprising ILs, while area per lipid increased.     

Molecular dynamics simulation of the size-dependent morphological stability of cubic shape silver nanoparticles

The morphological stability of sharp-edged silver nanoparticles is examined by the classical molecular dynamics (MD) simulations. The crystalline structure and the perfect fcc atom packing of a series of silver nanocubes (AgNC) of different sizes varying from 63 up to 1099 atoms are compared against quasi-spherical nanoparticles of the same sizes at temperature 303 K. Our MD simulations demonstrate that starting from the preformed perfect crystalline structures the cubic shape is preserved for AgNCs composed of 365–1099 atoms. Surprisingly, the rapid loss of the cubic shape morphology and transformation into the non-fcc-structure are found for smaller AgNCs composed of less than ~256 atoms. No such loss of the preformed crystalline structure is seen for quasi-spherical nanoparticles composed of 38–1007 atoms. The analysis of the temperature dependence and the binding energy of outermost Ag surface atoms suggests that the loss of the perfect cubic shape, rounding and smoothing of sharp edges and corners are driven by the tendency towards the increase in their coordination number. In addition, we revealed that AgNC₁₀₉₉ partially loses its sharp edges and corners in the aqueous environment; however, the polymer coating with poly(vinyl alcohol) (PVA) was able to preserve the well-defined cubic morphology. Finally, these results help improve the understanding of the role of surface capping agents in solution phase synthesis of Ag nanocubes.     

Computational studies of proton transport through the M2 channel

The M2 ion channel is an essential component of the influenza A virus. This low-pH gated channel has a high selectivity for protons. Evidence from various experimental data has indicated that the essential structure responsible for the channel is a parallel homo-tetrameric alpha-helix bundle having a left-handed twist with each helix tilted with respect to the membrane normal. A backbone structure has been determined by solid state nuclear magnetic resonance (NMR). Though detailed structures for the side chains are not available yet, evidence has indicated that His37 and Trp41 in the alpha-helix are implicated in the local molecular structure responsible for the selectivity and channel gate. More definitive conformations for the two residues were recently suggested based on the known backbone structure and recently obtained NMR data. While two competitive proton-conductance mechanisms have been proposed, the actual proton-conductance mechanism remains an unsolved problem. Computer simulations of an excess proton in the channel and computational studies of the His37/Trp41 conformations have provided insights into these structural and mechanism issues.     

Exploring the role of elongation Factor-Like 1 (EFL1) in Shwachman-Diamond syndrome through molecular dynamics

Abstract

Shwachman-Diamond Syndrome (SDS) is an autosomal recessive disorder whose patients present mutations in two ribosome assembly proteins, the Shwachman-Bodian-Diamond Syndrome protein (SBDS) and the Elongation Factor-Like 1 (EFL1). Due to the lack of knowledge of the molecular mechanisms responsible for SDS pathogenesis, current therapy is nonspecific and focuses only at alleviating the symptoms. Building on the recent observation that EFL1 single-point mutations clinically manifest as SDS-like phenotype, we carried out comparative Molecular Dynamics (MD) simulations on three mutants, T127A, M882K and R1095Q and wild type EFL1. As supported by small angle X-ray scattering experiments, the obtained data improve the static EFL1 model resulting from the Cryo-electron microscopy and clearly show that all the mutants experience a peculiar rotation, around the hinge region, of domain IV with respect to domains I and II leading to a different conformation respect to that of wild type protein. This study supports the notion that EFL1 function is governed by an allosteric mechanism involving the concerted action of GTPase domain (domain I) and the domain IV and can help point towards new approaches to SDS treatment.

Communicated by Ramaswamy H. Sarma     

Deciphering the mechanism behind the varied binding activities of COXIBs through Molecular Dynamic Simulations,MM-PBSA binding energy calculations and per-residue energy decomposition studies

COX-2 is a well-known drug target in inflammatory disorders. COX-1/COX-2 selectivity of NSAIDs is crucial in assessing the gastrointestinal side effects associated with COX-1 inhibition. Celecoxib, rofecoxib, and valdecoxib are well-known specific COX-2 inhibiting drugs. Recently, polmacoxib, a COX-2/CA-II dual inhibitor has been approved by the Korean FDA. These COXIBs have similar structure with diverse activity range. Present study focuses on unraveling the mechanism behind the 10-fold difference in the activities of these sulfonamide-containing COXIBs. In order to obtain insights into their binding with COX-2 at molecular level, molecular dynamics simulations studies, and MM-PBSA approaches were employed. Further, per-residue decomposition of these energies led to the identification of crucial amino acids and interactions contributing to the differential binding of COXIBs. The results clearly indicated that Leu338, Ser339, Arg499, Ile503, Phe504, Val509, and Ser516 (Leu352, Ser353, Arg513, Ile517, Phe518, Val523, and Ser530 in PGHS-1 numbering) were imperative in determining the activity of these COXIBs. The binding energies and energy contribution of various residues were similar in all the three simulations. The results suggest that hydrogen bond interaction between the hydroxyl group of Ser516 and five-membered ring of diarylheterocycles augments the affinity in COXIBs. The SAR of the inhibitors studied and the per-residue energy decomposition values suggested the importance of Ser516. Additionally, the positive binding energy obtained with Arg106 explains the binding of COXIBs in hydrophobic channel deep in the COX-2 active site. The findings of the present work would aid in the development of potent COX-2 inhibitors.