Explicit Solvent Molecular Dynamics Simulations Research Articles

Conformation and dynamics of the C-terminal region in human phosphoglycerate mutase 1.

Phosphoglycerate mutase 1 (PGAM1), an important enzyme in glycolysis, is overexpressed in a number of human cancers, thus has been proposed as a promising metabolic target for cancer treatments. The C-terminal portion of the available crystal structures of PGAM1 and its homologous proteins is partially disordered, as evidenced by weak electron density. In this study, we identified the conformational behavior of the C-terminal region of PGAM1 as well as its role during the catalytic cycle. Using the PONDR-FIT server, we demonstrated that the C-terminal region was intrinsically disordered. We applied the Monte Carlo (MC) method to explore the conformational space of the C-terminus and conducted a series of explicit-solvent molecular dynamics (MD) simulations, and revealed that the C-terminal region is inherently dynamic; large-scale conformational changes in the C-terminal segment led to the structural transition of PGAM1 from the closed state to the open state. Furthermore, the C-terminal segment influenced 2,3-bisphosphoglycerate (2,3-BPG) binding. The proposed swing model illustrated a critical role of the C-terminus in the catalytic cycle through the conformational changes. In conclusion, the C-terminal region induces large movements of PGAM1 from the closed state to the open state and influences cofactor binding during the catalytic cycle. This report describes the dynamic features of the C-terminal region in detail and should aid in design of novel and efficient inhibitors of PGAM1. A swing mechanism of the C-terminal region is proposed, to facilitate further studies of the catalytic mechanism and the physiological functions of its homologues.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

Effect of physicochemical properties of peptides from soy protein on their antimicrobial activity

Antimicrobial peptides (AMPs) kill microbial cells through insertion and damage/permeabilization of the cytoplasmic cell membranes and has applications in food safety and antibiotic replacement. Soy protein is an attractive, abundant natural source for commercial production of AMPs. In this research, explicit solvent molecular dynamics (MD) simulation was employed to investigate the effects of (i) number of total and net charges, (ii) hydrophobicity (iii) hydrophobic moment and (iv) helicity of peptides from soy protein on their ability to bind to lipid bilayer and their transmembrane aggregates to form pores. Interaction of possible AMP segments from soy protein with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPC/POPG) bilayers, a mimic of bacterial cell membrane, was investigated. Pore formation was insensitive to helicity and occurred for hydrophobicity threshold in the range of −0.3–0kcal/mol, hydrophobic moment threshold of 0.3kcal/mol, net charge threshold of 2. Though low hydrophobicity and high number of charges help in the formation of water channel for transmembrane aggregates, insertion of peptides with these properties requires overcome of energy barrier, as shown by potential of mean force calculations, thereby resulting in low antimicrobial activity. Experimental evaluation of antimicrobial activity of these peptides against Gram positive L. monocytogenes and Gram negative E. coli as obtained by spot-on-lawn assay was consistent with simulation results. These results should help in the development of guidelines for selection of peptides with antimicrobial activity based on their physicochemical properties.

Conformational and energy evaluations of novel peptides binding to dengue virus envelope protein

Effective novel peptide inhibitors which targeted the domain III of the dengue envelope (E) protein by blocking dengue virus (DENV) entry into target cells, were identified. The binding affinities of these peptides towards E-protein were evaluated by using a combination of docking and explicit solvent molecular dynamics (MD) simulation methods. The interactions of these complexes were further investigated by using the Molecular Mechanics-Poisson Boltzmann Surface Area (MMPBSA) and Molecular Mechanics Generalized Born Surface Area (MMGBSA) methods. Free energy calculations of the peptides interacting with the E-protein demonstrated that van der Waals (vdW) and electrostatic interactions were the main driving forces stabilizing the complexes. Interestingly, calculated binding free energies showed good agreement with the experimental dissociation constant (Kd) values. Our results also demonstrated that specific residues might play a crucial role in the effective binding interactions. Thus, this study has demonstrated that a combination of docking and molecular dynamics simulations can accelerate the identification process of peptides as potential inhibitors of dengue virus entry into host cells.

Discrimination of Polynucleotide Transport through a Highly Hydrophobic Uncharged Nanopore

The importance of DNA structures in medicine has motivated researchers to develop sequencing methods with low costs and high speeds. The potential use of nanofluidics for the analysis of the translocation of polynucleotides through a nanopore is a very interesting solution. Herein, we demonstrate both by experiments and by explicit solvent molecular dynamics simulations that the translocation of polyadenylic acid and polycytidylic acid can significantly affect the ionic flux behavior. In particular, the nanochannel enhances the velocity difference between polynucleotides during the translocation process. Our results revealed the important role of the strand structures during their translocation, as already observed in biochannels. The folding/unfolding process leads to stronger interactions with the nanopore wall that decrease the velocity of the strands. This behavior could lead to the detection of different polynucleotides in a mixture.

Small Molecules Engage Hot Spots through Cooperative Binding To Inhibit a Tight Protein-Protein Interaction.

Protein-protein interactions drive every aspect of cell signaling, yet only a few small-molecule inhibitors of these interactions exist. Despite our ability to identify critical residues known as hot spots, little is known about how to effectively engage them to disrupt protein-protein interactions. Here, we take advantage of the ease of preparation and stability of pyrrolinone 1, a small-molecule inhibitor of the tight interaction between the urokinase receptor (uPAR) and its binding partner, the urokinase-type plasminogen activator uPA, to synthesize more than 40 derivatives and explore their effect on the protein-protein interaction. We report the crystal structure of uPAR bound to previously discovered pyrazole 3 and to pyrrolinone 12. While both 3 and 12 bind to uPAR and compete with a fluorescently labeled peptide probe, only 12 and its derivatives inhibit the full uPAR·uPA interaction. Compounds 3 and 12 mimic and engage different hot-spot residues on uPA and uPAR, respectively. Interestingly, 12 is involved in a π-cation interaction with Arg-53, which is not considered a hot spot. Explicit-solvent molecular dynamics simulations reveal that 3 and 12 exhibit dramatically different correlations of motion with residues on uPAR. Free energy calculations for the wild-type and mutant uPAR bound to uPA or 12 show that Arg-53 interacts with uPA or with 12 in a highly cooperative manner, thereby altering the contributions of hot spots to uPAR binding. The direct engagement of peripheral residues not considered hot spots through π-cation or salt-bridge interactions could provide new opportunities for enhanced small-molecule engagement of hot spots to disrupt challenging protein-protein interactions.

Free Energy Perturbation Calculations of the Thermodynamics of Protein Side-Chain Mutations

Protein side-chain mutation is fundamental both to natural evolutionary processes and to the engineering of protein therapeutics, which constitute an increasing fraction of important medications. Molecular simulation enables the prediction of the effects of mutation on properties such as binding affinity, secondary and tertiary structure, conformational dynamics, and thermal stability. A number of widely differing approaches have been applied to these predictions, including sequence-based algorithms, knowledge-based potential functions, and all-atom molecular mechanics calculations. Free energy perturbation theory, employing all-atom and explicit-solvent molecular dynamics simulations, is a rigorous physics-based approach for calculating thermodynamic effects of, for example, protein side-chain mutations. Over the past several years, we have initiated an investigation of the ability of our most recent free energy perturbation methodology to model the thermodynamics of protein mutation for two specific problems: protein–protein binding affinities and protein thermal stability. We highlight recent advances in the field and outline current and future challenges.

A Hybrid Molecular Dynamics/Multiconformer Continuum Electrostatics (MD/MCCE) Approach for the Determination of Surface Charge of Nanomaterials

The surface charge of nanomaterials determines their stability in solution and interaction with other molecules and surfaces, yet experimental determination of surface charge of complex nanomaterials is not straightforward. We propose a hybrid approach that iteratively integrates explicit solvent molecular dynamics simulations and a multiconformer continuum electrostatic model (MCCE) to efficiently sample the configurational and titration spaces of surface ligands of nanomaterials. Test calculations of model systems indicate that the iterative approach converges rapidly even for systems that contain hundreds of titratable sites, making the approach complementary to more elaborate methods such as explicit solvent-based constant-pH molecular dynamics. The hybrid method is applied to analyze the pKa distribution of alkylamines attached to a carbon-based nanoparticle as a function of ligand density, nanoparticle surface curvature, and ligand heterogeneity. The results indicate that functionalization strategie...

Structural insight into the antiprion compound inhibition mechanism of native prion folding over misfolding.

Transition of a physiological folded prion (PrPC ) into a pathogenic misfolded prion (PrPSc ) causes lethal neurodegenerative disorders and prion diseases. Antiprion compounds have been developed to prevent this conversion; however, their mechanism of action remains unclear. Recently, we reported two antiprion compounds, BMD29 and BMD35, identified by in silico and in vitro screening. In this study, we used extensive explicit-solvent molecular dynamics simulations to investigate ligand-binding inhibition by antiprion compounds in prion folding over misfolding behavior at acidic pH. The two antiprion compounds and the previously reported GN8 compound resulted in a remarkably stabilized intermediate by binding to the hotspot region of PrPC , whereas free PrPC and the inactive compound BMD01 destabilized the structure of PrPC leading to the misfolded form. The results uncovered a secondary structural transition of free PrPC and transition suppression by the antiprion compounds. One of the major misfolding processes in PrPC , alternation of hydrophobic core residues, disruption of intramolecular interactions, and the increase in residue solvent exposure were significantly inhibited by both antiprion compounds. These findings provide insights into prion misfolding and inhibition by antiprion compounds.

Nucleosome Opening Kinetics and the Influence of Histone Modifications Studied by Single Molecule FRET

Nucleosomes compact the genome and regulate access by wrapping and unwrapping DNA. Here we characterized NaCl-induced opening of mononucleosomes by single-molecule FRET, using donor/acceptor labels on the DNA and various histones. By species-selective fluorescence lifetime and photon distribution analysis we identified new nucleosome opening intermediates and developed a kinetic model [1]. Opening proceeds through a weakening of the H2A-H2B dimer/(H3-H4)2 tetramer interface on a 0.1 ms time scale, then by a slower two-step release of the dimers coupled to DNA unwrapping, extending from several ms to minutes. Nucleosome opening and detachment of histone dimers proceed asymmetrically depending on the DNA sequence.Mutations at the H2A/H3 interface (H2A R81A, R88A, R81A/R88A, R81E/R88E) facilitate the initial opening, confirming the importance of the dimer:tetramer interface for nucleosome stability. This is also supported by molecular dynamics (MD) simulations that show enhanced DNA fluctuations for the mutants. Partially opened states such as described here might be a convenient nucleation point for DNA-recognizing proteins.For characterizing the role of histone tails in nucleosome interaction and stability, we performed MD simulations in explicit and implicit solvent on the H4 tail interacting with the acidic patch of a neighboring nucleosome. The simulations show that H4K16 acetylation decreases the interactions between H4 tail and the acidic patch of the neighboring nucleosome, a process that is central to the regulation of chromatin compaction [2,3].

Glycosaminoglycan Sequence Influence on Iduronate Ring Puckering

Iduronate (IdoA) is a major monosaccharide constituent of the linear polysaccharides, called glycosaminoglycans (GAGs), that form covalent conjugates with core proteins to make proteoglycans (PGs). PGs are known as important structural components of extracellular matrix, and are increasingly recognized for their roles in cell signaling. Toward advancing understanding of the atomic-resolution structural biology of PGs, we have performed all-atom explicit-solvent molecular dynamics (MD) simulations on IdoA as a monosaccharide and in the context of oligosaccharides corresponding to fragments of the GAG dermatan sulfate (DS). The stereochemistry of the carbon atoms in the pyranose ring of IdoA and the resulting spatial distribution of functional groups attached to the ring result in a delicate thermodynamic balance between preferred ring puckering conformation for IdoA. Therefore, Cremer-Pople ring pucker analysis was applied to understand the sequence influence on IdoA in the context of DS oligosaccharide vs. IdoA monosaccharide. These results are compared with analogous simulations for glucuronate (GlcA), which is a single-site epimer (at C6) of IdoA, and makes up the related GAG chondroitin sulfate (CS).

Promiscuous Contacts and Heightened Dynamics Increase Thermostability in an Engineered Variant of the Engrailed Homeodomain

The Engrailed Homeodomain (EnHD) is a three-helix-bundle transcription factor that is popular within the protein folding community due to its ultra-fast un/folding kinetics. In 2007, the Mayo group used EnHD as a backbone template for their repacking algorithm and designed UVF as a proof of principle. UVF not only folded as designed but was extremely stable; its Tm is >99°C whereas the Tm of EnHD is 52°C. Here, we present a structural and dynamical explanation for UVF's increased thermostability based on 1 μs of all-atom, explicit-solvent molecular dynamics simulations of the two proteins at room temperature (25°C) and high temperature (100°C). UVF had heightened dynamics at room temperature, relative to EnHD, as measured by Cα motion and interchangeability of side-chain contacts. UVF maintained these dynamics at 100°C whereas EnHD began to unfold. UVF's binary distribution of amino acids, with only hydrophobic residues at buried positions and only hydrophilic residues on the surface, allowed buried residues to move independently from the surface, effectively acting as an entropy sink. Hydrogen bonds and salt bridges between buried and surface residues of EnHD coupled core motion to surface motion, resulting in EnHD denaturing at 100°C. Based on our simulation data, we hypothesize that UVF benefits from more favorable entropy of folding in addition to its rationally designed, more-favorable enthalpy of folding. Future work will investigate this hypothesis experimentally and determine how the heightened dynamics effect native DNA-binding function.

Novel PARP-1 Inhibitor Scaffolds Disclosed by a Dynamic Structure-Based Pharmacophore Approach.

PARP-1 inhibition has been studied over the last decades for the treatment of various diseases. Despite the fact that several molecules act as PARP-1 inhibitors, a reduced number of compounds are used in clinical practice. To identify new compounds with a discriminatory PARP-1 inhibitory function, explicit-solvent molecular dynamics simulations using different inhibitors bound to the PARP-1 catalytic domain were performed. The representative structures obtained were used to generate structure-based pharmacophores, taking into account the dynamic features of receptor-ligand interactions. Thereafter, a virtual screening of compound databases using the pharmacophore models obtained was performed and the hits retrieved were subjected to molecular docking-based scoring. The drug-like molecules featuring the best ranking were evaluated for their PARP-1 inhibitory activity and IC50 values were calculated for the top scoring docked compounds. Altogether, three new PARP-1 inhibitor chemotypes were identified.

Molecular Dynamics Simulations for Deciphering the Structural Basis of Recognition of Pre-let-7 miRNAs by LIN28.

LIN28 protein inhibits biogenesis of miRNAs belonging to the let-7 family by binding to precursor forms of miRNAs. Overexpression of LIN28 and low levels of let-7 miRNAs are associated with several forms of cancer cells. We have performed multiple explicit solvent molecular dynamics simulations ranging from 200 to 500 ns in length on different isoforms of preE-let-7 in complex with LIN28 and also in isolation to identify structural features and key specificity-determining residues (SDRs) that are important for the inhibitory role of LIN28. Our simulations suggest that a conserved structural feature of the loop regions of preE-let-7 miRNAs is more important for LIN28 recognition than sequence conservation among members of the let-7 family or the presence of the GGAG motif in the 3' region. The loop region consisting of a minimum of five nucleotides helps pre-miRNAs to acquire a conformation ideal for binding to LIN28, but pre-let-7c-2 prefers a conformation with a three-nucleotide loop. Thus, our simulations provide a theoretical rationale for the recent experimental observation of the escape of LIN28-mediated repression by pre-let-7c-2. The essential structural and sequence features highlighted in this study might aid in designing synthetic small molecule inhibitors for modulating LIN28-let-7 interaction in malignant cells. We have also identified crucial SDRs of the LIN28-preE-let-7 complex involving 13 residues of LIN28 and 10 residues of the pre-miRNA. On the basis of the conservation profile of these 13 SDRs, we have identified 10 novel proteins that are not annotated as LIN28 like but are similar in sequence, domain, or fold level to LIN28.

Rigidity and flexibility in the tetrasaccharide linker of proteoglycans from atomic-resolution molecular simulation.

Proteoglycans (PGs) are covalent conjugates between protein and carbohydrate (glycosaminoglycans). Certain classes of glycosaminoglycans such as chondroitin sulfate/dermatan sulfate and heparan sulfate utilize a specific tetrasaccharide linker for attachment to the protein component: GlcAβ1-3Galβ1-3Galβ1-4Xylβ1-O-Ser. Toward understanding the conformational preferences of this linker, the present work used all-atom explicit-solvent molecular dynamics (MD) simulations combined with Adaptive Biasing Force (ABF) sampling to determine high-resolution, high-precision conformational free energy maps ΔG(φ, ψ) for each glycosidic linkage between constituent disaccharides, including the variant where GlcA is substituted with IdoA. These linkages are characterized by single, predominant (> 97% occupancy), and broad (45° × 60° for ΔG(φ, ψ) < 1 kcal/mol) free-energy minima, while the Xyl-Ser linkage has two such minima similar in free-energy, and additional flexibility from the Ser sidechain dihedral. Conformational analysis of microsecond-scale standard MD on the complete tetrasaccharide-O-Ser conjugate is consistent with ABF data, suggesting (φ, ψ) probabilities are independent of the linker context, and that the tetrasaccharide acts as a relatively rigid unit whereas significant conformational heterogeneity exists with respect to rotation about bonds connecting Xyl to Ser. © 2017 Wiley Periodicals, Inc.

Molecular dynamics simulations of β2-microglobulin interaction with hydrophobic surfaces.

Hydrophobic surfaces are known to adsorb and unfold proteins, a process that has been studied only for a few proteins. Here we address the interaction of β2-microglobulin, a paradigmatic protein for the study of amyloidogenesis, with hydrophobic surfaces. A system with 27 copies of the protein surrounded by a model cubic hydrophobic box is studied by implicit solvent molecular dynamics simulations. Most proteins adsorb on the walls of the box without major distortions in local geometry, whereas free molecules maintain proper structures and fluctuations as observed in explicit solvent molecular dynamics simulations. The major conclusions from the simulations are as follows: (i) the adopted implicit solvent model is adequate to describe protein dynamics and thermodynamics; (ii) adsorption occurs readily and is irreversible on the simulated timescale; (iii) the regions most involved in molecular encounters and stable interactions with the walls are the same as those that are important in protein-protein and protein-nanoparticle interactions; (iv) unfolding following adsorption occurs at regions found to be flexible by both experiments and simulations; (v) thermodynamic analysis suggests a very large contribution from van der Waals interactions, whereas unfavorable electrostatic interactions are not found to contribute much to adsorption energy. Surfaces with different degrees of hydrophobicity may occur in vivo. Our simulations show that adsorption is a fast and irreversible process which is accompanied by partial unfolding. The results and the thermodynamic analysis presented here are consistent with and rationalize previous experimental work.

Allosteric effects of gold nanoparticles on human serum albumin.

The ability of nanoparticles to alter protein structure and dynamics plays an important role in their medical and biological applications. We investigate allosteric effects of gold nanoparticles on human serum albumin protein using molecular simulations. The extent to which bound nanoparticles influence the structure and dynamics of residues distant from the binding site is analyzed. The root mean square deviation, root mean square fluctuation and variation in the secondary structure of individual residues on a human serum albumin protein are calculated for four protein-gold nanoparticle binding complexes. The complexes are identified in a brute-force search process using an implicit-solvent coarse-grained model for proteins and nanoparticles. They are then converted to atomic resolution and their structural and dynamic properties are investigated using explicit-solvent atomistic molecular dynamics simulations. The results show that even though the albumin protein remains in a folded structure, the presence of a gold nanoparticle can cause more than 50% of the residues to decrease their flexibility significantly, and approximately 10% of the residues to change their secondary structure. These affected residues are distributed on the whole protein, even on regions that are distant from the nanoparticle. We analyze the changes in structure and flexibility of amino acid residues on a variety of binding sites on albumin and confirm that nanoparticles could allosterically affect the ability of albumin to bind fatty acids, thyroxin and metals. Our simulations suggest that allosteric effects must be considered when designing and deploying nanoparticles in medical and biological applications that depend on protein-nanoparticle interactions.

Comparative Assessment of Nonlocal Continuum Solvent Models Exhibiting Overscreening

AbstractNonlocal continua have been proposed to offer a more realistic model for the electrostatic response of solutions such as the electrolyte solvents prominent in biology and electrochemistry. In this work, we review three nonlocal models based on the Landau-Ginzburg framework which have been proposed but not directly compared previously, due to different expressions of the nonlocal constitutive relationship. To understand the relationships between these models and the underlying physical insights from which they are derive, we situate these models into a single, unified Landau-Ginzburg framework. One of the models offers the capacity to interpret how temperature changes affect dielectric response, and we note that the variations with temperature are qualitatively reasonable even though predictions at ambient temperatures are not quantitatively in agreement with experiment. Two of these models correctly reproduce overscreening (oscillations between positive and negative polarization charge densities), and we observe small differences between them when we simulate the potential between parallel plates held at constant potential. These computations require reformulating the two models as coupled systems of local partial differential equations (PDEs), and we use spectral methods to discretize both problems. We propose further assessments to discriminate between the models, particularly in regards to establishing boundary conditions and comparing to explicit-solvent molecular dynamics simulations.

Interstrand cross-linking implies contrasting structural consequences for DNA: insights from molecular dynamics.

Oxidatively-generated interstrand cross-links rank among the most deleterious DNA lesions. They originate from abasic sites, whose aldehyde group can form a covalent adduct after condensation with the exocyclic amino group of purines, sometimes with remarkably high yields. We use explicit solvent molecular dynamics simulations to unravel the structures and mechanical properties of two DNA sequences containing an interstrand cross-link. Our simulations palliate the absence of experimental structural and stiffness information for such DNA lesions and provide an unprecedented insight into the DNA embedding of lesions that represent a major challenge for DNA replication, transcription and gene regulation by preventing strand separation. Our results based on quantum chemical calculations also suggest that the embedding of the ICL within the duplex can tune the reaction profile, and hence can be responsible for the high difference in yields of formation.

Microscopic mechanisms that govern the titration response and pKa values of buried residues in staphylococcal nuclease mutants.

To probe the microscopic mechanisms that govern the titration behavior of buried ionizable groups, microsecond explicit solvent molecular dynamics simulations are carried out for several mutants of Staphylococcal nuclease using both fixed charge and polarizable force fields. While the ionization of Asp 66, Glu 66, and Lys 125 lead to enhanced structural fluctuations and partial unfolding of adjacent α-helical regions, the ionization of Lys 25 causes local unfolding of adjacent β sheets. Using the sampled conformational ensembles, good agreement with experimental pKa values is obtained with Poisson-Boltzmann calculations using a protein dielectric constant of 2-4 for V66D/E; slightly larger dielectric constants are needed for Lys mutants especially L25K, suggesting that structural responses beyond microseconds are involved in ionization of Lys 25. Overall, the set of unbiased simulations provides insights into the spatial and temporal scales of protein and solvent motions that dictate the diverse titration behaviors of buried protein residues. Proteins 2017; 85:268-281. © 2016 Wiley Periodicals, Inc.