Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

Molecular Mechanisms of ZnO Nanoparticle Dispersion in Solution: Modeling of Surfactant Association,Electrostatic Shielding and Counter Ion Dynamics

Molecular models of 5 nm sized ZnO/Zn(OH)2 core-shell nanoparticles in ethanolic solution were derived as scale-up models (based on an earlier model created from ion-by-ion aggregation and self-organization) and subjected to mechanistic analyses of surface stabilization by block-copolymers. The latter comprise a poly-methacrylate chain accounting for strong surfactant association to the nanoparticle by hydrogen bonding and salt-bridges. While dangling poly-ethylene oxide chains provide only a limited degree of sterical hindering to nanoparticle agglomeration, the key mechanism of surface stabilization is electrostatic shielding arising from the acrylates and a halo of Na⁺ counter ions associated to the nanoparticle. Molecular dynamics simulations reveal different solvent shells and distance-dependent mobility of ions and solvent molecules. From this, we provide a molecular rationale of effective particle size, net charge and polarizability of the nanoparticles in solution.     

Study of interactions between polymer nanoparticles and cell membranes at atomistic levels

Knowledge of how the structure of nanoparticles and the interactions with biological cell membranes is important not only for understanding nanotoxicological effects on human, animal health and the environment, but also for better understanding of nanoparticle fabrication for biomedical applications. In this work, we use molecular modelling techniques, namely molecular dynamics (MD) simulations, to explore how polymer nanoparticles interact with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid cell membranes. Two different polymers have been considered: 100 monomer units of polyethylene (approx. 2.83 kDa) and polystyrene (approx. 10.4 kDa), both of which have wide industrial applications. We found that, despite the polar lipid head groups acting as an effective barrier to prevent the nanoparticles from interacting with the membrane surface, irreversible adhesion can be initiated by insertion of dangling chain ends from the polymer into the hydrophobic interior of the membrane. In addition, alignment of chain segments from the polymers with that of hydrocarbon chains in the interior of the membrane facilitates the complete immersion of the nanoparticles into the cell membrane. These findings highlight the importance of the surface and the topological structures of the polymer particles that dictate the absorption behaviour into the membrane and, subsequently, induce the possible translocation into the cell.     

The effect of surface coverage and conformation of poly(ethylene oxide) (PEO) chains of poloxamer 407 on the biological fate of model colloidal drug carriers

Poloxamer 407 was adsorbed onto the surface of model colloidal drug carriers, polystyrene nanoparticles of 40, 70 and 137 nm in diameter, and the effect of the degree of surface coverage and the conformation of the poly(ethylene oxide) (PEO) chains on biological fate was studied. The relationship between the physicochemical and the biological properties of the nanoparticle systems was also investigated. The adsorbed layer of poloxamer 407 was characterised in terms of percentage surface coverage, thickness of the adsorbed layer and average surface area per PEO chain. Computer modelling of the adsorbed layer was performed (applying the self-consistent field technique), to obtain the structural information of the PEO chains in the layer. The in vitro interaction of the nanoparticles with different degrees of poloxamer 407 surface coverage with serum components and the in vivo biodistribution in the rat model were assessed. The results demonstrated that an increase in the surface coverage with poloxamer 407 resulted in an increased volume fraction of the PEO in the adsorbed layer, further extension of the PEO chains from the surface and closer packing of the chains at the surface. With regard to the interaction with the serum components, an increased surface coverage resulted in a reduction of the amount of serum proteins adsorbed, and, importantly, affected the type of proteins adsorbed. High molecular weight proteins were not adsorbed onto the nanoparticles with a surface coverage above approx. 25%. Following the intravenous administration to rats, even the nanoparticles with the lowest degree of surface coverage (approx. 5%) showed improved circulation profiles relative to the uncoated nanoparticles. The effect was more pronounced for the 40 nm nanoparticles. A further increase in the surface coverage to approx. 25% resulted in a significant increase in circulation time, as compared to uncoated and 5% coated systems, for all sizes of nanoparticles. Importantly, it was found that a long in vivo blood circulation time could be achieved for nanoparticles with a relatively low degree of surface coverage with PEO chains.     

Molecular dynamics studies of the antimicrobial peptides piscidin 1 and its mutants with a DOPC lipid bilayer

Piscidin 1 (Pis‐1) has a high broad‐spectrum activity against bacteria, fungi, and viruses but it also has a moderate hemolytic activities. To improve the antibacterial activity and to reduce toxicity, mutants Pis‐1AA (G8A/G13A double mutant) and Pis‐1PG (G8P mutant) have been designed based on the crystal structure of Pis‐1. Eighteen independent molecular dynamics (MD) simulations of Pis‐1 and its mutants with membranes are conducted in this article. Furthermore, 60 independent MD simulations of three peptides in water box have also been discussed for comparison. The results indicate that the unfolding process starts at the middle of the peptide. Pis‐1 disrupts easily in the region of Val10‐Lys14. Pis‐1PG has a flexible N‐terminal region, and the interaction between N‐terminal and C‐terminal is very weak. Pis‐1AA has the most stable helical structure. In addition, percentage of native contacts and hydrogen bonds analysis are also performed. Lipid‐peptide interaction analysis suggests that Pis‐1 and Pis‐1AA has a stronger interaction with the zwitterionic dioleoylphosphatidylcholine (DOPC) lipid bilayer than Pis‐1PG. When compared with the results of peptide with membrane, peptides are unstable and unfolding quickly in water solution. Our results are applicable in examining diversities on hemolytic, antibacterial, and selectivity of antimicrobial peptides. © 2012 Wiley Periodicals, Inc. Biopolymers 97:998–1009, 2012.     

Structure and energy landscape of a photoswitchable peptide: a replica exchange molecular dynamics study

A replica exchange molecular dynamics (REMD) simulation of a bicyclic azobenzene peptide in explicit dimethyl sulfoxide solution is presented in order to characterize the conformational structures and energy landscape of a photoswitchable peptide. It is shown that an enhanced-sampling technique such as the REMD method is essential to obtain a converged conformational sampling of the peptide at room temperature. This is because conventional MD simulations of less than approximately 100-ns length are either trapped in local minima (at 295 K) or-if run at high temperature-do not resemble the room-temperature REMD results. Calculating various nuclear Overhauser effects (NOEs) and (3)J-couplings, a good overall agreement between the REMD simulations and the NMR experiments of Renner et al. (Biopolymers 2000;54:501-514) is found. In particular, the REMD study confirms the general picture drawn by Renner et al. that the trans-isomer of the azobenzene peptide exhibits a well-defined structure, while the cis-isomer is a conformational heterogeneous system; that is, the trans-isomer occurs in 2 well-defined conformers, while the cis-isomer represents an energetically frustrated system that leads to an ensemble of conformational structures. Employing a principal component analysis of the REMD data, the free energy landscape of the systems is studied at various temperatures. The implications for the folding and unfolding pathways of the system are discussed.     

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Mimicking the action of folding chaperones by Hamiltonian replica-exchange molecular dynamics simulations: application in the refinement of de novo models

The efficiency of using a variant of Hamiltonian replica‐exchange molecular dynamics (Chaperone H‐replica‐exchange molecular dynamics [CH‐REMD]) for the refinement of protein structural models generated de novo is investigated. In CH‐REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH‐REMD approach sampled structures in which the root‐mean‐square deviation (RMSD) of secondary structure elements (SSE‐RMSD) with respect to the experimental structure was more than 1.0 Å lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 Å. The ability of three different statistical potentials to identify near‐native conformations was also examined. Little correlation between the SSE‐RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE‐RMSD of the best scoring structures was on average 0.3 Å lower than the initial model. Overall the work demonstrates that targeted enhanced‐sampling techniques such as CH‐REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near‐native structures are still needed. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.     

The Selective Interaction between Silica Nanoparticles and Enzymes from Molecular Dynamics Simulations

Nanoscale particles have become promising materials in many fields, such as cancer therapeutics, diagnosis, imaging, drug delivery, catalysis, as well as biosensors. In order to stimulate and facilitate these applications, there is an urgent need for the understanding of the interaction mode between the nano-particles and proteins. In this study, we investigate the orientation and adsorption between several enzymes (cytochrome c, RNase A, lysozyme) and 4 nm/11 nm silica nanoparticles (SNPs) by using molecular dynamics (MD) simulation. Our results show that three enzymes are adsorbed onto the surfaces of both 4 nm and 11 nm SNPs during our MD simulations and the small SNPs induce greater structural stabilization. The active site of cytochrome c is far away from the surface of 4 nm SNPs, while it is adsorbed onto the surface of 11 nm SNPs. We also explore the influences of different groups (-OH, -COOH, -NH₂ and CH₃) coated onto silica nanoparticles, which show significantly different impacts. Our molecular dynamics results indicate the selective interaction between silicon nanoparticles and enzymes, which is consistent with experimental results. Our study provides useful guides for designing/modifying nanomaterials to interact with proteins for their bio-applications.     

Evaluation of conformational changes in diabetes‐associated mutation in insulin a chain: A molecular dynamics study

Insulin plays a central role in the regulation of metabolism in humans. Mutations in the insulin gene can impair the folding of its precursor protein, proinsulin, and cause permanent neonatal‐onset diabetes mellitus known as Mutant INS‐gene induced Diabetes of Youth (MIDY) with insulin deficiency. To gain insights into the molecular basis of this diabetes‐associated mutation, we perform molecular dynamics simulations in wild‐type and mutant (Cys^A7 to Tyr or C(A7)Y) insulin A chain in aqueous solutions. The C(A7)Y mutation is one of the identified mutations that impairs the protein folding by substituting the cysteine residue which is required for the disulfide bond formation. A comparative analysis reveals structural differences between the wild‐type and the mutant conformations. The analyzed mutant insulin A chain forms a metastable state with major effects on its N‐terminal region. This suggests that MIDY mutant involves formation of a partially folded intermediate with conformational change in N‐terminal region in A chain that generates flexible N‐terminal domain. This may lead to the abnormal interactions with other proinsulins in the aggregation process. Proteins 2015; 83:662–669. © 2015 Wiley Periodicals, Inc.     

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.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.     

Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone

The molecular interaction between common polymer chains and the cell membrane is unknown. Molecular dynamics simulations offer an emerging tool to characterise the nature of the interaction between common degradable polymer chains used in biomedical applications, such as polycaprolactone, and model cell membranes. Herein we characterise with all-atomistic and coarse-grained molecular dynamics simulations the interaction between single polycaprolactone chains of varying chain lengths with a phospholipid membrane. We find that the length of the polymer chain greatly affects the nature of interaction with the membrane, as well as the membrane properties. Furthermore, we next utilise advanced sampling techniques in molecular dynamics to characterise the two-dimensional free energy surface for the interaction of varying polymer chain lengths (short, intermediate, and long) with model cell membranes. We find that the free energy minimum shifts from the membrane-water interface to the hydrophobic core of the phospholipid membrane as a function of chain length. Finally, we perform coarse-grained molecular dynamics simulations of slightly larger membranes with polymers of the same length and characterise the results as compared with all-atomistic molecular dynamics simulations. These results can be used to design polymer chain lengths and chemistries to optimise their interaction with cell membranes at the molecular level.     

Conformational fluctuations versus constraints in amino acid side chains: the evolution of information content from free amino acids to proteins

Like all other complex biological systems, proteins exhibit properties not found in free amino acids (i.e., emergent properties). Here, we explore top-down constraints experienced by the residue side chains in proteins compared to amino acids in increasingly complex molecular environments: free amino acids, end-capped amino acids, and the central residue in an alpha-helical nonapeptide. The crystalline structure of the contractile protein profilin Ib and the enzyme trypsin were chosen as objects of study, and submitted to 10 ns molecular dynamics (MD) simulations. The results revealed increased conformational constraints on the side chains when going from the simpler to the more complex compounds. A Shannon entropy (SE) analysis of the conformational behavior of the side chains showed in most cases a progressive and marked decrease in the SE of the chi1 and chi2 dihedral angles. This is equivalent to stating that conformational constraints on the side chain of residues increase their information content and, hence, recognition specificity compared to free amino acids. In other words, the vastly increased information content of a protein relative to its free monomers is embedded not only in the tertiary structure of the backbone, but also in the conformational behavior of the side chains. The postulated implication is that both backbone and side chains, by virtue of being conformationally constrained, contribute to the protein's recognition specificity toward other macromolecules and ligands.     

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.