Dynamics of the protein structure of T169S pyranose 2‐oxidase in solution: Molecular dynamics simulation

Pyranose 2‐oxidase (P2O) from Trametes multicolor contains FAD as cofactor, and forms a tetramer. The protein structure of a mutated P2O, T169S (Thr169 is replaced by Ser), in solution was studied by means of molecular dynamics simulation and analyses of photoinduced electron transfer (ET) from Trp168 to excited isoalloxazine (Iso*), and was compared with wild type (WT) P2O. Hydrogen bonding between Iso and nearby amino acids was very similar as between T169S and WT protein. Distances between Iso and Tyr456 were extremely heterogeneous among the subunits, 1.7 (1.5 in WT) in subunit A (Sub A), 0.97 (2.2 in WT) in Sub B, 1.3 (2.1 in WT) in Sub C, 1.3 nm (2.0 in WT) in Sub D. Mean values of root of mean square fluctuation over all residues were greater by four times than those in WT. This suggests that the protein structure of T169S is much more flexible than that of WT. Electrostatic (ES) energies between Iso anion in one subunit and ionic groups in the entire protein were evaluated. It was found that more than 50% of the total ES energy in each subunit is contributed from other subunits. Reported fluorescence decays were analyzed by a method as WT, previously reported. Electron affinities of Iso* in T169S were appreciably higher than those in WT. Static dielectric constants near Iso and Trp168 were also quite higher in T169S than those in WT.     

Temperature-mediated switching of protectant-denaturant behavior of trimethylamine-N-oxide and consequences on protein stability from a replica exchange molecular dynamics simulation study

The detailed mechanism of protein folding–unfolding processes with the aid of osmolytes has been a leading topic of discussion over many decades. We have used replica-exchange molecular dynamics simulation to propose the molecular mechanism of interaction of a 20-residue mini-protein with urea and trimethylamine N-oxide (TMAO) that act as denaturing and protecting osmolyte, respectively, in binary osmolyte solutions. Urea is found to exert its action by interacting directly with the protein residues. Temperature tolerance of TMAO’s action is particularly emphasised in this study. At lower range of temperature, TMAO acts as a successful protein protectant. Interestingly, the study discloses the tendency of TMAO molecules to prefer self-association at the protein surface at elevated temperature. A greater number of TMAO molecules in the protein hydration shell at higher temperature is also observed. Dihedral angle principal component analysis and free energy landscape plots sampled all possible conformations adopted by the protein that reveal highly folded behaviour of the protein in pure water and binary TMAO solutions and highly unfolded behaviour in presence of urea.     

Comparison of the monomer structure of the FMN-binding protein from Desulfovibrio vulgaris obtained by NMR and molecular dynamics simulation approaches

Flavin mononucleotide (FMN)-binding proteins (FBPs) play an important role in the electron transport process in bacteria. In this study, the structures of the FBP from Desulfovibrio vulgaris (DvFBP) (Miyazaki F) were compared between those obtained experimentally by nuclear magnetic resonance (NMR) spectroscopy and those derived from molecular dynamics simulations (MDSs). A high-residue root of mean square deviation (RMSD) was observed in residues located at both sides of the wings (Gly22, Glu23, Asp24, Ala59, Arg60, Asp61, Glu62, Gly75, Arg76, Asn77, Gly78 and Pro79), while a low-residue RMSD was found in residues located in a hollow of the structure (Asn12, Glu13, Gly14, Val15, Val16, Asn30, Thr31, Trp32, Asn33, Ser34, Gly69, Ser70, Arg71 and Lys72). Inter-planar angles between the Phe7 and Iso and between the Phe7 and Trp106 residues were remarkably different between the MDS- and NMR-derived DvFBP structures. Distribution of the torsion angles around the covalent bonds in the aliphatic chain of FMN was similar in the MDS- and NMR-derived structures, except for those around the C1′–C2′ and C5′–O5′ bonds. Hydrogen bond formation between IsoO2 and the Gly49 or Gly50 peptide NH was formed in both the NMR- and MDS-derived structures. Overall, the MDS-derived structures were found to be considerably different from the NMR-derived structures, which must be considered when the photoinduced electron transfer in flavoproteins is analysed with MDS-derived structures.     

A molecular dynamics study of solvent behavior around a protein

The solvent structure and behavior around a protein were examined by analyzing a trajectory of molecular dynamics simulation of thetrp-holorepressor in a periodic box of water. The calculated selfdiffusion coefficient indicated that the solvent within 10 Å of the protein had lower mobility. Examination of the solvent diffusion around different atoms of different kinds of residues showed no general tendency. Thisfact suggested that the solvent mobility is not influenced significantly bythe kind of the atom or residue they solvated. Distribution analysis aroundthe protein revealed two peaks of water oxygen: a sharp one at 2.8 Å around polar and charged atoms and a broad one at ～3.4 Å aroundapolar atoms. The former was stabilized by water–protein hydrogen bonds, and the latter was stabilized by water-lwater hydrogen bonds, suggesting the existence of a hydrophobic shell. An analysis of protein atom–water radial distribution functions confirmed these shell structures around polar or charged atoms and apolar ones. © 1993 Wiley-Liss, Inc.     

Prediction of protein mutation effects based on dehydration and hydrogen bonding – A large‐scale study

Reliable computational prediction of protein side chain conformations and the energetic impact of amino acid mutations are the key aspects for the optimization of biotechnologically relevant enzymatic reactions using structure‐based design. By improving the protein stability, higher yields can be achieved. In addition, tuning the substrate selectivity of an enzymatic reaction by directed mutagenesis can lead to higher turnover rates. This work presents a novel approach to predict the conformation of a side chain mutation along with the energetic effect on the protein structure. The HYDE scoring concept applied here describes the molecular interactions primarily by evaluating the effect of dehydration and hydrogen bonding on molecular structures in aqueous solution. Here, we evaluate its capability of side‐chain conformation prediction in classic remutation experiments. Furthermore, we present a new data set for evaluating “cross‐mutations,” a new experiment that resembles real‐world application scenarios more closely. This data set consists of protein pairs with up to five point mutations. Thus, structural changes are attributed to point mutations only. In the cross‐mutation experiment, the original protein structure is mutated with the aim to predict the structure of the side chain as in the paired mutated structure. The comparison of side chain conformation prediction (“remutation”) showed that the performance of HYDE_protein is qualitatively comparable to state‐of‐the art methods. The ability of HYDE_protein to predict the energetic effect of a mutation is evaluated in the third experiment. Herein, the effect on protein stability is predicted correctly in 70% of the evaluated cases. Proteins 2017; 85:1550–1566. © 2017 Wiley Periodicals, Inc.     

A molecular dynamics study of acylphosphatase in aggregation-promoting conditions: the influence of trifluoroethanol/water solvent

The 98-residue protein acylphosphatase exhibits a high propensity for aggregation under certain conditions. Aggregates formed from wild-type acylphosphatase in the presence of 2,2,2-trifluoroethanol and from highly destabilized mutants are essentially identical in structure. Furthermore, it has been shown by mutational studies that different regions of the protein are important for aggregation and folding. In the present molecular dynamics study, we compare the behavior of the protein in aqueous solution and in a 25% (v/v) 2,2,2-trifluoroethanol/water environment mimicking the experimental conditions. The 2,2,2-trifluoroethanol surrounding affects the structure of the protein mostly in the regions important for aggregation, in good agreement with experimental data. This suggests that the early step of (partly) unfolding, which precedes the aggregation process, has been observed.     

N-glycosylation induced changes in tau protein dynamics reveal its role in tau misfolding and aggregation: A microsecond long molecular dynamics study

Various posttranslational modifications like hyperphosphorylation, O-GlcNAcylation, and acetylation have been attributed to induce the abnormal folding in tau protein. Recent in vitro studies revealed the possible involvement of N-glycosylation of tau protein in the abnormal folding and tau aggregation. Hence, in this study, we performed a microsecond long all atom molecular dynamics simulation to gain insights into the effects of N-glycosylation on Asn-359 residue which forms part of the microtubule binding region. Trajectory analysis of the stimulations coupled with essential dynamics and free energy landscape analysis suggested that tau, in its N-glycosylated form tends to exist in a largely folded conformation having high beta sheet propensity as compared to unmodified tau which exists in a large extended form with very less beta sheet propensity. Residue interaction network analysis of the lowest energy conformations further revealed that Phe378 and Lys353 are the functionally important residues in the peptide which helped in initiating the folding process and Phe378, Lys347, and Lys370 helped to maintain the stability of the protein in the folded state.     

Solvent sensitivity of protein aggregation in Cu,Zn superoxide dismutase: a molecular dynamics simulation study

Misfolding and aggregation of Cu, Zn Superoxide Dismutase (SOD1) is often found in amyotrophic lateral sclerosis (ALS) patients. The central apo SOD1 barrel was involved in protein maturation and pathological aggregation in ALS. In this work, we employed atomistic molecular dynamics (MD) simulations to study the conformational dynamics of SOD1^barrel monomer in different concentrations of trifluoroethanol (TFE). We find concentration dependence unusual structural and dynamical features, characterized by the local unfolding of SOD1^barrel. This partially unfolded structure is characterized by the exposure of hydrophobic core, is highly dynamic in nature, and is the precursor of aggregation seen in SOD1^barrel. Our computational studies supports the hypothesis of the formation of aggregation ‘building blocks’ by means of local unfolding of apo monomer as the mechanism of SOD1 fibrillar aggregation. The non-monotonic TFE concentration dependence of protein conformational changes was explored through simulation studies. Our results suggest that altered protein conformation and dynamics within its structure may underlie the aggregation of SOD1 in ALS.     

A molecular dynamics simulation study of segment B1 of protein G

The immunoglobulin binding protein, segment B1 of protein G, has been studied experimentally as a paradigm for protein folding. This protein consists of 56 residues, includes both β sheet and α helix and contains neither disulfide bonds nor proline residues. We report an all-atom molecular dynamics study of the native manifold of the protein in explicit solvent. A 2-ns simulation starting from the nuclear magnetic resonance (NMR) structure and a 1-ns control simulation starting from the x-ray structure were performed. The difference between average structures calculated over the equilibrium portion of trajectories is smaller than the difference between their starting conformations. These simulation averages are structurally similar to the x-ray structure and differ in systematic ways from the NMR-determined structure. Partitioning of the fluctuations into fast (<20 ps) and slow (<20 ps) components indicates that the β sheet displays greater long-time mobility than does the α helix. Clore and Gronenborn [J. Mol. Biol. 223:853–856, 1992] detected two long-residence water molecules by NMR in a solution structure of segment B1 of protein G. Both molecules were found in the fully exposed regions and were proposed to be stabilized by bifurcated hydrogen bonds to the protein backbone. One of these long-residence water molecules, found near an exposed loop region, is identified in both of our simulations, and is seen to be involved in the formation of a stable water-mediated hydrogen bond bridge. The second water molecule, located near the middle of the α helix, is not seen with an exceptional residence time in either as a result of the conformation being closer to the x-ray structure in this region of the protein. Proteins 29:193–202, 1997. © 1997 Wiley-Liss, Inc.     

REMO: A new protocol to refine full atomic protein models from C‐alpha traces by optimizing hydrogen‐bonding networks

Protein structure prediction approaches usually perform modeling simulations based on reduced representation of protein structures. For biological utilizations, it is an important step to construct full atomic models from the reduced structure decoys. Most of the current full atomic model reconstruction procedures have defects which either could not completely remove the steric clashes among backbone atoms or generate final atomic models with worse topology similarity relative to the native structures than the reduced models. In this work, we develop a new protocol, called REMO, to generate full atomic protein models by optimizing the hydrogen‐bonding network with basic fragments matched from a newly constructed backbone isomer library of solved protein structures. The algorithm is benchmarked on 230 nonhomologous proteins with reduced structure decoys generated by I‐TASSER simulations. The results show that REMO has a significant ability to remove steric clashes, and meanwhile retains good topology of the reduced model. The hydrogen‐bonding network of the final models is dramatically improved during the procedure. The REMO algorithm has been exploited in the recent CASP8 experiment which demonstrated significant improvements of the I‐TASSER models in both atomic‐level structural refinement and hydrogen‐bonding network construction. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Molecular basis of the structural stability of hemochromatosis factor E: A combined molecular dynamic simulation and GdmCl‐induced denaturation study

Hemochromatosis factor E (HFE) is a member of class I MHC family and plays a significant role in the iron homeostasis. Denaturation of HFE induced by guanidinium chloride (GdmCl) was measured by monitoring changes in [θ]₂₂₂ (mean residue ellipticity at 222 nm), intrinsic fluorescence emission intensity at 346 nm (F₃₄₆) and the difference absorption coefficient at 287 nm (Δε₂₈₇) at pH 8.0 and 25°C. Coincidence of denaturation curves of these optical properties suggests that GdmCl‐induced denaturation (native (N) state ? denatured (D) state) is a two‐state process. The GdmCl‐induced denaturation was found reversible in the entire concentration range of the denaturant. All denaturation curves were analyzed for , Gibbs free energy change associated with the denaturation equilibrium (N state ? D state) in the absence of GdmCl, which is a measure of HFE stability. We further performed molecular dynamics simulation for 40 ns to see the effect of GdmCl on the structural stability of HFE. A well defined correlation was established between in vitro and in silico studies. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 133–142, 2016.     

3Drefine: Consistent protein structure refinement by optimizing hydrogen bonding network and atomic‐level energy minimization

One of the major limitations of computational protein structure prediction is the deviation of predicted models from their experimentally derived true, native structures. The limitations often hinder the possibility of applying computational protein structure prediction methods in biochemical assignment and drug design that are very sensitive to structural details. Refinement of these low‐resolution predicted models to high‐resolution structures close to the native state, however, has proven to be extremely challenging. Thus, protein structure refinement remains a largely unsolved problem. Critical assessment of techniques for protein structure prediction (CASP) specifically indicated that most predictors participating in the refinement category still did not consistently improve model quality. Here, we propose a two‐step refinement protocol, called 3Drefine, to consistently bring the initial model closer to the native structure. The first step is based on optimization of hydrogen bonding (HB) network and the second step applies atomic‐level energy minimization on the optimized model using a composite physics and knowledge‐based force fields. The approach has been evaluated on the CASP benchmark data and it exhibits consistent improvement over the initial structure in both global and local structural quality measures. 3Drefine method is also computationally inexpensive, consuming only few minutes of CPU time to refine a protein of typical length (300 residues). 3Drefine web server is freely available at http://sysbio.rnet.missouri.edu/3Drefine/ . Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

Cumulative deamidations of the major lens protein γS‐crystallin increase its aggregation during unfolding and oxidation

Age‐related lens cataract is the major cause of blindness worldwide. The mechanisms whereby crystallins, the predominant lens proteins, assemble into large aggregates that scatter light within the lens, and cause cataract, are poorly understood. Due to the lack of protein turnover in the lens, crystallins are long‐lived. A major crystallin, γS, is heavily modified by deamidation, in particular at surface‐exposed N14, N76, and N143 to introduce negative charges. In this present study, deamidated γS was mimicked by mutation with aspartate at these sites and the effect on biophysical properties of γS was assessed via dynamic light scattering, chemical and thermal denaturation, hydrogen‐deuterium exchange, and susceptibility to disulfide cross‐linking. Compared with wild type γS, a small population of each deamidated mutant aggregated rapidly into large, light‐scattering species that contributed significantly to the total scattering. Under partially denaturing conditions in guanidine hydrochloride or elevated temperature, deamidation led to more rapid unfolding and aggregation and increased susceptibility to oxidation. The triple mutant was further destabilized, suggesting that the effects of deamidation were cumulative. Molecular dynamics simulations predicted that deamidation augments the conformational dynamics of γS. We suggest that these perturbations disrupt the native disulfide arrangement of γS and promote the formation of disulfide‐linked aggregates. The lens‐specific chaperone αA‐crystallin was poor at preventing the aggregation of the triple mutant. It is concluded that surface deamidations cause minimal structural disruption individually, but cumulatively they progressively destabilize γS‐crystallin leading to unfolding and aggregation, as occurs in aged and cataractous lenses.     

Influence of pressure on the low‐frequency vibrational modes of lysozyme and water: A complementary inelastic neutron scattering and molecular dynamics simulation study

We performed complementary inelastic neutron scattering (INS) experiments and molecular dynamics (MD) simulations to study the influence of pressure on the low‐frequency vibrational modes of lysozyme in aqueous solution in the 1 atm–6 kbar range. Increasing pressure induces a high‐frequency shift of the low‐frequency part (<10 meV = 80 cm^?1) of the vibrational density of states (VDOS), g(ω), of both lysozyme and water that reveals a stiffening of the interactions ascribed to the reduction of the protein and water volumes. Accordingly, high pressures increase the curvature of the free energy profiles of the protein quasiharmonic vibrational modes. Furthermore, the nonlinear influence of pressure on the g(ω) of lysozyme indicates a change of protein dynamics that reflects the nonlinear pressure dependence of the protein compressibility. An analogous dynamical change is observed for water and stems from the distortion of its tetrahedral structure under pressure. Moreover, our study reveals that the structural, dynamical, and vibrational properties of the hydration water of lysozyme are less sensitive to pressure than those of bulk water, thereby evidencing the strong influence of the protein surface on hydration water. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Mass spectrometry‐based methods to study protein architecture and dynamics

Mass spectrometry is now an indispensable tool in the armamentarium of molecular biophysics, where it is used for tasks ranging from protein sequencing and mapping of post‐translational modifications to studies of higher order structure, conformational dynamics, and interactions of proteins with small molecule ligands and other biopolymers. This mini‐review highlights several popular mass spectrometry‐based tools that are now commonly used for structural studies of proteins beyond their covalent structure with a particular emphasis on hydrogen exchange and direct electrospray ionization mass spectrometry.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Structural analysis of the N‐terminal fragment of the antiangiogenic protein endostatin: A molecular dynamics study

Endostatin is a potent antiangiogenic protein derived from the noncollagenous domain 1 (NC1) of collagen XVIII. The mechanism by which endostatin exerts its antiangiogenic effect is still incompletely understood. It has been shown that the 27 amino acid N‐terminal fragment of murine endostatin has antitumor, antimigration, and antipermeability activities comparable to the full soluble protein. To understand how this peptide can exert such elaborate function, we performed structural analysis using molecular dynamics to evaluate the behavior of this fragment in aqueous environment. Here, we show that the N‐terminal peptide of murine endostatin is able to assume a well‐defined structure, folding into a zinc‐dependent β‐hairpin conformation. Analyzing the folding mechanism, we were able to understand why the N‐terminal peptide of human endostatin with the same length failed to acquire a stable conformation. Conversely, we were able to predict the successful folding of the R4Q mutant and of a shorter form of the human peptide with 25 residues. Finally, we show that the β‐hairpin conformation assumed by the zinc‐bound peptide of murine endostatin has a high structural similarity with fragments of another family of angiogenesis inhibitors: the integrin‐binding portion of the NC1 domain of collagen IV. Indeed, our docking simulations show that arresten, canstatin, and the endostatin peptide bind to the same spot of αVβ3 integrin, suggesting similar interactions via a common binding site on this receptor. Proteins 2011;. © 2011 Wiley‐Liss, Inc.     

Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation study

Substrate phosphorylation by cAMP-dependent-protein kinase A (protein kinase A, PKA) has been studied extensively. Phosphoryl transfer was found to be fast, whereas ADP release was found to be the slow, rate-limiting step. There is also evidence that ADP release may be preceded by a partially rate-limiting conformational change. However, the atomic details of the conformational change and the mode of ADP release are difficult to obtain experimentally. In this work, we studied ADP release from PKA by carrying out molecular dynamics simulations with different pulling forces applied to the ligand. The detailed ADP release pathway and the associated conformational changes were analyzed. The ADP release process was found to involve a swinging motion with the phosphate of ADP anchored to the Gly-rich loop, so that the more buried adenine base and ribose ring came out before the phosphate. In contrast to the common belief that a hinge-bending motion was responsible for the opening of the ligand-binding cleft, our simulations showed that the small lobe exhibited a large amplitude "rocking" motion when the ligand came out. The largest conformational change of the protein was observed at about the first quarter time point along the release pathway. Two prominent intermediate states were observed in the release process.