Comparison of the structural characteristics of Cu2+-bound and unbound α-syn12 peptide obtained in simulations using different force fields

The effects of Cu²⁺ binding and the utilization of different force fields when modeling the structural characteristics of α-syn12 peptide were investigated. To this end, we performed extensive temperature replica exchange molecular dynamics (T-REMD) simulations on Cu²⁺-bound and unbound α-syn12 peptide using the GROMOS 43A1, OPLS-AA, and AMBER03 force fields. Each replica was run for 300 ns. The structural characteristics of α-syn12 peptide were studied based on backbone dihedral angle distributions, free-energy surfaces obtained with different reaction coordinates, favored conformations, the formation of different Turn structures, and the solvent exposure of the hydrophobic residues. The findings show that AMBER03 prefers to sample helical structures for the unbound α-syn12 peptide and does not sample any β-hairpin structure for the Cu²⁺-bound α-syn12 peptide. In contrast, the central structure of the major conformational clusters for the Cu²⁺-bound and unbound α-syn12 peptide according to simulations performed using the GROMOS 43A1 and OPLS-AA force fields is a β-hairpin with Turn_9-6. Cu²⁺ can also promote the formation of the β-hairpin and increase the solvent exposure of hydrophobic residues, which promotes the aggregation of α-syn12 peptide. This study can help us to understand the mechanisms through which Cu²⁺ participates in the fibrillation of α-syn12 peptide at the atomic level, which in turn represents a step towards elucidating the nosogenesis of Parkinson’s disease.

Investigation of the binding network of IGF-I on the cavity surface of IGFBP4

Insulin-like growth factor-binding proteins (IGFBPs) control bioactivity and distribution of insulin-like growth factors (IGFs) through high-affinity complex of IGFBP and IGF. To get more insight into the binding interaction of IGF system, the site-directed mutagenesis and force-driving desorption methods were employed to study the interaction mechanism of IGFBP4 and IGF-I by molecular dynamics (MD) simulation. In IGF-I, residues Gly7 to Asp12 were found to be the hot spots and they mainly anchored on the N-domain of IGFBP4. The contact area, the shape and size of protein, the surroundings of the binding site, the hydrophobic and electrostatic interaction between the two proteins worked as a complex network to regulate the protein-protein interaction. It was also found that the unfolding of the helix was not inevitable in the mutant, and it could be regulated by careful selection of the substituted amino acid.

Effects of protein binding on a lipid bilayer containing local anesthetic articaine, and the potential of mean force calculation: a molecular dynamics simulation approach

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer’s order parameter and the diffusion coefficient of drug are discussed.

Protein dynamics of a β-sheet protein

Rhodnius prolixus Nitrophorin 4 (abbreviated NP4) is an almost pure β-sheet heme protein. Its dynamics is investigated by X-ray structure determination at eight different temperatures from 122 to 304 K and by means of Mössbauer spectroscopy. A comparison of this β-sheet protein with the pure α-helical protein myoglobin (abbreviated Mbmet) is performed. The mean square displacement derived from the Mössbauer spectra increases linearly with temperature below a characteristic temperature T _c. It is about 10 K larger than that of myoglobin. Above T _c the mean square displacements increase dramatically. The Mössbauer spectra are analyzed by a two state model. The increased mean square displacements are caused by very slow motions occurring on a time scale faster than 140 ns. With respect to these motions NP4 shows the same protein specific modes as Mbmet. There is, however, a difference in the fast vibration regime. The B values found in the X-ray structures vary linearly over the entire temperature range. The mean square displacements in NP4 increase with slopes which are 60% larger than those observed for Mbmet. This indicates that nitrophorin has a larger structural distribution which makes it more flexible than myoglobin.     

Comparative analysis of thermal unfolding simulations of RNA recognition motifs (RRMs) of TAR DNA-binding protein 43 (TDP-43)

TAR DNA-binding protein 43 (TDP-43) inclusions have been found in Amyotrophic lateral sclerosis (ALS) and several other neurodegenerative diseases. Many studies suggest the involvement of RNA recognition motifs (RRMs) in TDP-43 proteinopathy. To elucidate the structural stability and the unfolding dynamics of RRMs, we have carried out atomistic molecular dynamics simulations at two different temperatures (300 and 500 K). The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two domains and RRM1 unfolds faster than RRM2 in accordance with the lower thermal stability found experimentally. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet and structural rearrangements of β-sheets results in formation of additional α-helices. At higher temperature, RRM1 exhibit increased overall flexibility and unfolding than RRM2. The temperature-dependent free energy landscapes consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds in RRM2, thus rendering the RRM2 more prone to misfolding. The structural rearrangements of RRM2 could lead to aberrant protein–protein interactions that may account for enhanced aggregation and toxicity of TDP-43. Our analysis, thus identify the structural and thermodynamic characteristics of the RRMs of TDP-43, which will serve to uncover molecular mechanisms and driving forces in TDP-43 misfolding and aggregation.     

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

PrP assemblies: Spotting the responsible regions in prion propagation

The “protein only” hypothesis states that the key phenomenon in prion pathogenesis is the conversion of the host protein (PrP^C) into a β-sheet enriched polymeric and pathogenic conformer (PrP^Sc). However the region of PrP bearing the information for structural transfer is still controversial. In a recent report, we highlighted the role of the C terminal part i.e., the helixes H2 and H3, using mutation approaches on recombinant PrP. The H2H3 was shown to be the minimal region necessary to reproduce the oligomerization pattern of the full-length protein. The oligomers produced from isolated H2H3 domain presented the same structural characteristics as the oligomers formed from the full-length PrP. Combining other groups'' results, this paper further discusses the relative, direct or indirect role of different PrP regions in assembly. The H2H3 region represents the core of PrP oligomers and fibrils, whereas the N terminus could explain divergences among different aggregates. Finally this review evocates the possibility to separate the domain involved in prion information transference (i.e., prion replication) from the domain bearing the cytotoxicity properties.Key words: prion, H2H3, amyloid, domain of replication, unfolding, strain, polymer, fibersTransmissible spongiform encephalopathies (TSE), fatal neurodegenerative diseases affecting humans and other mammalians, induce in most cases loss of motor control and dementia. PrP is a protein physiologically present in parts of the animal kingdom (in mammals, birds, reptiles and fishes). According to the “protein-only” hypothesis,^1,2 the key phenomenon in the pathogenesis is the conversion of the α-helix rich host-encoded PrP form (PrP^C) into a pathogenic conformer (PrP^Sc) characterized by a higher content in β-sheet and a polymeric state. The conversion to an enriched β-sheet structure is supposed to be due to the modification—induced only by a PrP^Sc-like state acting as a template- of PrP^C into the PrP^Sc conformer. This hypothesis was first proposed by Griffith in 1967³ and revisited by Lansbury et al. in 1993.⁴ The prion hypothesis has now found increasing support from experimental evidence based on the synthetic production of β-sheeted recombinant PrP which shows pathogenic properties in a wide variety of physico-chemical conditions.^5–7 However, the molecular basis of prion conversion remains unclear, especially the various structural landscape of the PrP^Sc, which is the basis of the strain phenomenon.⁸To understand the mechanisms of transfer of the structural information, two mains issues have to be addressed: (1) we need to understand which region(s) of the protein act as template for conversion and (2) what is the “pathogenic” state of this domain. In this review, we shall assume that the region bearing the infectious information for replication and the region responsible for polymerization are identical. However, the link between the propensity of a domain to form aggregates and the ability to contain the necessary information for prion replication is far from being trivial. Generally the formation of amyloid assemblies results from the aggregation of disordered peptides or in some cases from disordered regions of a folded protein.⁸ If we consider that prion replication is only supported by the globular part of PrP⁹ the currently available model involves the folded domain. Since all structural transitions need at least a partial unfolding and refolding process, pre-required structural events should be considered prior to the conversion process.     

A class V chitinase from <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>: gene responses,enzymatic properties,and crystallographic analysis

Expression of a class V chitinase gene (At4g19810, AtChiC) in Arabidopsis thaliana was examined by quantitative real-time PCR and by analyzing microarray data available at Genevestigator. The gene expression was induced by the plant stress-related hormones abscisic acid (ABA) and jasmonic acid (JA) and by the stress resulting from the elicitor flagellin, NaCl, and osmosis. The recombinant AtChiC protein was produced in E. coli, purified, and characterized with respect to the structure and function. The recombinant AtChiC hydrolyzed N-acetylglucosamine oligomers producing dimers from the non-reducing end of the substrates. The crystal structure of AtChiC was determined by the molecular replacement method at 2.0 Å resolution. AtChiC was found to adopt an (β/α)₈ fold with a small insertion domain composed of an α-helix and a five-stranded β-sheet. From docking simulation of AtChiC with pentameric substrate, the amino acid residues responsible for substrate binding were found to be well conserved when compared with those of the class V chitinase from Nicotiana tabacum (NtChiV). All of the structural and functional properties of AtChiC are quite similar to those obtained for NtChiV, and seem to be common to class V chitinases from higher plants.     

Secondary structure perturbations in salt-induced protein precipitates

The secondary structure implications of precipitation induced by a chaotropic salt, KSCN, and a structure stabilizing salt, Na₂SO₄, were studied for twelve different proteins. α-helix and β-sheet content of precipitate and native structures were estimated from the analysis of amide I band Raman spectra. A statistical analysis of the estimated perturbations in the secondary structure contents indicated that the most significant event is the formation of β-sheet structures with a concomitant loss of α-helix on precipitation with KSCN. The conformational changes for each protein were also analyzed with respect to elements of primary, secondary and tertiary structure existing in the native protein; primary structure was quantified by the fractions of hydrophobic and charged amino acids, secondary structure by x-ray estimates of α-helix and β-sheet contents of native proteins and tertiary structure by the dipole moment and solvent-accessible surface area. For the KSCN precipitates, factors affecting β-sheet content included the fraction of charged amino acids in the primary sequence and the surface area. Changes in α-helix content were influenced by the initial helical content and the dipole moment. The enhanced β-sheet contents of precipitates observed in this work parallel protein structural changes occurring in other aggregative phenomena.     

MD simulation of self-diffusion and structure in some n-alkanes over a wide temperature range at high pressures

Self-diffusion and structural properties of n-alkanes have been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 600 K at pressures up to 300 MPa. The simulated results of lower n-alkanes are in good agreement with the existing experimental data, and support the reliability of results of the simulations of self-diffusion coefficients obtained at the extreme conditions. We predict the self-diffusion coefficients for methane, ethane, propane and n-butane at the similar reduced temperatures and pressures to draw a comparison between them. Then the correlation between self-diffusion and structural properties are further investigated by calculating the coordination numbers. Moreover, we define four distances and their corresponding relative deviations to characterize the flexibility of long-chain n-alkanes. The simulated results show that the self-diffusion of n-alkane molecules is mainly affected by the close packing, and the flexibility has a strong impact on the self-diffusion of longer n-alkane molecules.

In-silico screening of cancer associated mutation on PLK1 protein and its structural consequences

The Polo-like kinases (Plks) are a conserved subfamily of serine-threonine protein kinases that have significant roles in cell proliferation. The serine/threonine protein kinases or polo-like kinase 1 (PLK1) exist in centrosome during interphase and is an important regulatory enzyme in cell cycle progression during M phase. Mutations in mammalian PLK1 were found to be over expressed in various human cancers and it is disrupting the binding ability of polo box domain with target peptide. In this analysis we implemented a computational approach to filter the most deleterious and cancer associated mutation on PLK1 protein. We found W414F as the most deleterious and cancer associated by Polyphen 2.0, SIFT, I-mutant 3.0, PANTHER, PhD-SNP, SNP&GO, Mutpred and Dr Cancer tools. Molecular docking and molecular dynamics simulation (MDS) approach was used to investigate the structural and functional behavior of PLK1 protein upon mutation. MDS and docking results showed stability loss in mutant PLK1 protein. Due to mutation, PLK1 protein became more flexible and alters the dynamic property of protein which might affect the interaction with target peptide and leads to cell proliferation. Our study provided a well designed computational methodology to examine the cancer associated nsSNPs and their molecular mechanism. It further helps scientists to develop a drug therapy against PLK1 cancer-associated diseases.

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

Molecular dynamics study on the stability of wild-type and the R220K mutant of human prion protein

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases related to the structure transition of α-helix into β-sheet. In order to gain more direct insight into the molecular basis of the disease, the stability of the wild-type human prion protein (hPrPc) and the R220K mutant (m-hPrPc) was studied by molecular dynamics (MD) and flow MD simulation. Both the thermodynamic stability and the mechanical properties of hPrPc were investigated in this work. It was found that β-sheet was more readily to be unfolded in m-hPrPc. In the case of hPrPc, less content of helix was preserved after water turbulence. The H-bond network formed by the mutation-related residue 220 was found to play a key role in the stability of hPrPc.     

Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics

Periplasmic binding proteins are the initial receptors for the transport of various substrates over the inner membrane of gram-negative bacteria. The binding proteins are composed of two domains, and the substrate is entrapped between these domains. For several of the binding proteins it has been established that a closed-up conformation exists even without substrate present, suggesting a highly flexible apo-structure which would compete with the ligand-bound protein for the transporter interaction. For the leucine binding protein (LBP), structures of both open and closed conformations are known, but no closed-up structure without substrate has been reported. Here we present molecular dynamics simulations exploring the conformational flexibility of LBP. Coarse grained models based on the MARTINI force field are used to access the microsecond timescale. We show that a standard MARTINI model cannot maintain the structural stability of the protein whereas the ELNEDIN extension to MARTINI enables simulations showing a stable protein structure and nanosecond dynamics comparable to atomistic simulations, but does not allow the simulation of conformational flexibility. A modification to the MARTINI-ELNEDIN setup, referred to as domELNEDIN, is therefore presented. The domELNEDIN setup allows the protein domains to move independently and thus allows for the simulation of conformational changes. Microsecond domELNEDIN simulations starting from either the open or the closed conformations consistently show that also for LBP, the apo-structure is flexible and can exist in a closed form.

Computational analysis on conformational dynamics of bone morphogenetic protein-2 (BMP-2)

BMP-2 is widely used for bone regeneration because of its ability to induce osteoblast differentiation and proliferation. The pharmaceutical application of BMP-2 as bone implant makes the studies on stability and conformational dynamics very relevant as proteins are functional only in their native three-dimensional state. Knowing the factors affecting BMP-2 structure becomes essential for designing bone implants activated by BMP-2. In order to explore the influence of temperature and hydration on protein conformation, we have performed the molecular dynamics (MD) simulations at the time scale of 100 ns with two different force fields. We have examined the dynamic behaviour of BMP-2 monomer and dimer in aqueous medium as well as in vacuum at four different temperatures (300, 350, 400 and 450 K). MD simulation of BMP-2 monomer and dimer in water and vacuum environments shows the major contribution of water in structure stabilization. Temperature of the system affects the secondary structure differently in case of monomer and dimer simulation and the dynamics also depends on the environment viz. vacuum and aqueous. Vacuum simulations show very early loss of the major secondary structure content. On the other hand, BMP-2 monomer and dimer in aqueous environment show the unfolding of α-helix with increasing temperature. This unfolded α-helix is converted into β-sheet at 400 K in monomer of BMP-2. Contrary to this, we did not observe β-sheet formation in dimer BMP-2 even at 450 K indicating that monomers are more aggregation prone entity as compared to dimers of BMP-2.     

The Effects of Lipid Oxidation Product Acrolein on the Structure and Gel Properties of Rabbit Meat Myofibrillar Proteins

We investigated the effect of acrolein, a byproduct of lipid oxidation, on the structure and gel properties of myofibrillar proteins (MPs) isolated from rabbit meat. As the acrolein concentration increased, the protein carbonyl compounds significantly accumulated (p?<?0.05), and the total sulfhydryl content was significantly lost (p?<?0.05). The results of circular dichroism spectra, surface hydrophobicity, UV absorption spectra and intrinsic fluorescence spectra evidenced that acrolein caused the disruption of α-helix structure, the exposure of hydrophobic sites and the unfolding of MPs. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) analysis suggested that medium (0–1 mM) and high (5–10 mM) concentrations of acrolein could induce protein cross-linkage and protein aggregation, respectively. These structural changes could affect gelling properties of MPs involving gel strength and water holding capacity (WHC). The results of Raman spectroscopy indicated that moderate oxidative modification caused protein unfolding as well as the decline of α-helix structure and the increase of β-sheets structure in gels, thereby influencing the gel properties. Moderate oxidative modification (0–1 mM) improved gel strength and WHC, while excessive oxidative modification (5–10 mM) resulted in decreased gel properties.     

Structural elucidation of supramolecular alpha-cyclodextrin dimer/aliphatic monofunctional molecules complexes

The structural elucidation of 2α-cyclodextrin/1-octanethiol, 2α-cyclodextrin/1-octylamine and 2α-cyclodextrin/1-nonanoic acid inclusion complexes by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling has been achieved. The detailed spatial configurations are proposed for the three inclusion complexes based on 2D NMR method. ROESY experiments confirm the inclusion of guest molecules inside the α-cyclodextrin (α-CD) cavity. On the other hand, the host-guest ratio observed was 2:1 for three complexes. The detailed spatial configuration proposed based on 2D NMR methods were further interpreted using molecular modeling studies. The theoretical calculations are in good agreement with the experimental data.

Trifluoroethanol-induced conformational change of tetrameric and monomeric soybean agglutinin: Role of structural organization and implication for protein folding and stability

2,2,2-Trifuoroethanol (TFE)-induced conformational structure change of a β-sheet legume lectin, soybean agglutinin (SBA) has been investigated employing its exclusive structural forms in quaternary (tetramer) and tertiary (monomer) states, by far- and near-UV CD, FTIR, fluorescence, low temperature phosphorescence and chemical modification. Far-UV CD results show that, for SBA tetramer, native atypical β-conformation transforms to a highly α-helical structure, with the helical content reaching 57% in 95% TFE. For SBA monomer, atypical β-sheet first converts to typical β-sheet at low TFE concentration (10%), which then leads to a nonnative α-helix at higher TFE concentration. From temperature-dependent studies (5–60 °C) of TFE perturbation, typical β-sheet structure appears to be less stable than atypical β-sheet and the induced helix entails reduced thermal stability. The heat induced transitions are reversible except for atypical to typical β-sheet conversion. FTIR results reveal a partial α-helix conversion at high protein concentration but with quantitative yield. However, aggregation is detected with FTIR at lower TFE concentration, which disappears in more TFE. Near-UV CD, fluorescence and phosphorescence studies imply the existence of an intermediate with native-like secondary and tertiary structure, which could be related to the dissociation of tetramer to monomer. This has been further supported by concentration dependent far-UV CD studies. Chemical modification with N-bromosuccinimide (NBS) shows that all six tryptophans per monomer are solvent-exposed in the induced α-helical conformation. These results may provide novel and important insights into the perturbed folding problem of SBA in particular, and β-sheet oligomeric proteins in general.