Explicit Solvent Molecular Dynamics Research Articles

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.

Molecular dynamic simulations of glycine amino acid association with potassium and sodium ions in explicit solvent

Salt solutions are the natural environment in which biological molecules act, and dissolved ions are actively involved in biochemical processes. With metal ions, the membrane potentials are maintained. Ions are crucial for the activity of many enzymes, and their ability to coordinate with chemical groups modulates protein-protein interactions. Here we present a comparative study of sodium and potassium coordination with zwitterionic glycine, by means of explicit solvent molecular dynamics. We demonstrated that contact ion pair of cations and carboxylate group splits into two distinct coordination states. Sodium binding is significantly stronger than for potassium. These results can shed light on the different roles of sodium and potassium ions in abiogenic peptide synthesis.

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.

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.

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.

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.

Strengthening of the Coordination Shell by Counter Ions in Aqueous Th4+ Solutions.

The presence of counterions in solutions containing highly charged metal cations can trigger processes such as ion-pair formation, hydrogen bond breakages and subsequent re-formation, and ligand exchanges. In this work, it is shown how halide (Cl-, Br-) and perchlorate (ClO4-) anions affect the strength of the primary solvent coordination shells around Th4+ using explicit-solvent and finite-temperature ab initio molecular dynamics modeling methods. The 9-fold solvent geometry was found to be the most stable hydration structure in each aqueous solution. Relative to the dilute aqueous solution, the presence of the counterions did not significantly alter the geometry of the primary hydration shell. However, the free energy analyses indicated that the 10-fold hydrated states were thermodynamically accessible in dilute and bromide aqueous solutions within 1 kcal/mol. Analysis of the results showed that the hydrogen bond lifetimes were longer and solvent exchange energy barriers were larger in solutions with counterions in comparison with the solution with no counterions. This implies that the presence of the counterions induces a strengthening of the Th4+ hydration shell.

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.

Predicting the orientation of protein G B1 on hydrophobic surfaces using Monte Carlo simulations

A Monte Carlo algorithm was developed to predict the most likely orientations of protein G B1, an immunoglobulin G (IgG) antibody-binding domain of protein G, adsorbed onto a hydrophobic surface. At each Monte Carlo step, the protein was rotated and translated as a rigid body. The assumption about rigidity was supported by quartz crystal microbalance with dissipation monitoring experiments, which indicated that protein G B1 adsorbed on a polystyrene surface with its native structure conserved and showed that its IgG antibody-binding activity was retained. The Monte Carlo simulations predicted that protein G B1 is likely adsorbed onto a hydrophobic surface in two different orientations, characterized as two mutually exclusive sets of amino acids contacting the surface. This was consistent with sum frequency generation (SFG) vibrational spectroscopy results. In fact, theoretical SFG spectra calculated from an equal combination of the two predicted orientations exhibited reasonable agreement with measured spectra of protein G B1 on polystyrene surfaces. Also, in explicit solvent molecular dynamics simulations, protein G B1 maintained its predicted orientation in three out of four runs. This work shows that using a Monte Carlo approach can provide an accurate estimate of a protein orientation on a hydrophobic surface, which complements experimental surface analysis techniques and provides an initial system to study the interaction between a protein and a surface in molecular dynamics simulations.

Microsecond molecular dynamics simulations provide insight into the ATP-competitive inhibitor-induced allosteric protection of Akt kinase phosphorylation.

Akt is a serine/threonine protein kinase, a critical mediator of growth factor-induced survival in key cellular pathways. Allosteric signaling between protein intramolecular domains requires long-range communication mediated by hotspot residues, often triggered by ligand binding. Here, based on extensive 3μs explicit solvent molecular dynamics (MD) simulations of Akt1 kinase domain in the unbound (apo) and ATP-competitive inhibitor, GDC-0068-bound states, we propose a molecular mechanism for allosteric regulation of Akt1 kinase phosphorylation by GDC-0068 binding to the ATP-binding site. MD simulations revealed that the apo Akt1 is flexible with two disengaged N- and C-lobes, equilibrated between the open and closed conformations. GDC-0068 occupancy of the ATP-binding site shifts the conformational equilibrium of Akt1 from the open conformation toward the closed conformation and stabilizes the closed state. This effect enables allosteric signal propagation from the GDC-0068 to the phosphorylated T308 (pT308) in the activation loop and restrains phosphatase access to pT308, thereby protecting the pT308 in the GDC-0068-bound Akt1. Importantly, functional hotspots involved in the allosteric communication from the GDC-0068 to the pT308 are identified. Our analysis of GDC-0068-induced allosteric protection of Akt kinase phosphorylation yields important new insights into the molecular mechanism of allosteric regulation of Akt kinase activity.

How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes?

We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.

Comparative Assessment of Different RNA Tetranucleotides from the DFT-D3 and Force Field Perspective.

Classical force field (FF) molecular dynamics (MD) simulations of RNA tetranucleotides have substantial problems in reproducing conformer populations indicated by NMR experiments. To provide more information about the possible sources of errors, we performed quantum mechanical (QM, TPSS-D3/def2-TZVP) and molecular mechanics (MM, AMBER parm99bsc0+χOL3) calculations of different r(CCCC), r(GACC), and r(UUUU) conformers obtained from explicit solvent MD simulations. Solvent effects in the static QM and MM calculations were mimicked using implicit solvent models (COSMO and Poisson-Boltzmann, respectively). The comparison of QM and MM geometries and energies revealed that the two methodologies provide qualitatively consistent results in most of the cases. Even though we found some differences, these were insufficient to indicate any systematic corrections of the RNA FF terms that could improve the performance of classical MD in simulating tetranucleotides. On the basis of these findings, we inferred that the overpopulation of intercalated conformers in the MD simulations of RNA tetramers, which were not observed experimentally, might be predominantly caused by imbalanced water-solvent and water-water interactions. Apart from the large-scale QM calculations performed to assess the performance of the AMBER FF, a representative spectrum of faster QM methods was tested.

Simulation study of the ability of a computationally-designed peptide to recognize target tRNALys3 and other decoy tRNAs.

In this paper, we investigate the ability of our computationally-designed peptide, Pept10 (PNWNGNRWLNNCLRG), to recognize the anticodon stem and loop (ASL) domain of the hypermodified tRNALys3 (mcm5 s2 U34 ,ms2 t6 A37 ), a reverse transcription primer of HIV replication. Five other ASLs, the singly modified ASLLys3 (ms2 t6 A37 ), ASLLys3 (s2 U34 ), ASLLys3 (Ψ39 ), ASLLys1,2 (t6 A37 ), and ASLGlu (s2 U34 ), were used as decoys. Explicit-solvent atomistic molecular dynamics simulations were performed to examine the process of binding of Pept10 with the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) and the decoy ASLs. Simulation results demonstrated that Pept10 is capable of recognizing the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) as well as one of the decoys, ASLLys3 (Ψ39 ), but screens out the other four decoy ASLs. The interchain van der Waals (VDW) and charge-charge (ELE + EGB) energies for the two best complexes were evaluated to shed light on the molecular recognition mechanism between Pept10 and ASLs. The results indicated that Pept10 recognizes and binds to the target ASLLys3 (mcm5 s2 U34 ,ms2 t6 A37 ) through residues W3 and R7 which interact with the nucleotides mcm5 s2 U34 , U35 , and ms2 t6 A37 via the interchain VDW energy. Pept10 also recognizes the decoy ASLLys3 (Ψ39 ) through residue R14 which contacts the nucleotide U36 via the interchain VDW energy. Regardless of the type of ASL, the positively charged arginines on Pept10 are attracted to the negatively charged phosphate linkages on the ASL via the interchain ELE + EGB energy, thereby enhancing the binding affinity.

An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase.

The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.

Correlation of bistranded clustered abasic DNA lesion processing with structural and dynamic DNA helix distortion

Clustered apurinic/apyrimidinic (AP; abasic) DNA lesions produced by ionizing radiation are by far more cytotoxic than isolated AP lesion entities. The structure and dynamics of a series of seven 23-bp oligonucleotides featuring simple bistranded clustered damage sites, comprising of two AP sites, zero, one, three or five bases 3′ or 5′ apart from each other, were investigated through 400 ns explicit solvent molecular dynamics simulations. They provide representative structures of synthetically engineered multiply damage sites-containing oligonucleotides whose repair was investigated experimentally (Nucl. Acids Res. 2004, 32:5609-5620; Nucl. Acids Res. 2002, 30: 2800–2808). The inspection of extrahelical positioning of the AP sites, bulge and non Watson–Crick hydrogen bonding corroborates the experimental measurements of repair efficiencies by bacterial or human AP endonucleases Nfo and APE1, respectively. This study provides unprecedented knowledge into the structure and dynamics of clustered abasic DNA lesions, notably rationalizing the non-symmetry with respect to 3′ to 5′ position. In addition, it provides strong mechanistic insights and basis for future studies on the effects of clustered DNA damage on the recognition and processing of these lesions by bacterial or human DNA repair enzymes specialized in the processing of such lesions.

Lowered pH Leads to Fusion Peptide Release and a Highly Dynamic Intermediate of Influenza Hemagglutinin.

Hemagglutinin (HA), the membrane-bound fusion protein of the influenza virus, enables the entry of virus into host cells via a structural rearrangement. There is strong evidence that the primary trigger for this rearrangement is the low pH environment of a late endosome. To understand the structural basis and the dynamic consequences of the pH trigger, we employed explicit-solvent molecular dynamics simulations to investigate the initial stages of the HA transition. Our results indicate that lowered pH destabilizes HA and speeds up the dissociation of the fusion peptides (FPs). A buried salt bridge between the N-terminus and Asp1122 of HA stem domain locks the FPs and may act as one of the pH sensors. In line with recent observations from simplified protein models, we find that, after the dissociation of FPs, a structural order-disorder transition in a loop connecting the central coiled-coil to the C-terminal domains produces a highly mobile HA. This motion suggests the existence of a long-lived asymmetric or "symmetry-broken" intermediate during the HA conformational change. This intermediate conformation is consistent with models of hemifusion, and its early formation during the conformational change has implications for the aggregation seen in HA activity.

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two-domain constructs: A computational characterization of novel chimeric proteins.

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one-domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two-domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all-atom explicit solvent molecular dynamics simulations. The well-characterized PDZ3 and SH3 domains of human zonula occludens (ZO-1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two-domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358-1374. © 2016 Wiley Periodicals, Inc.

Protein-protein structure prediction by scoring molecular dynamics trajectories of putative poses.

The prediction of protein-protein interactions and their structural configuration remains a largely unsolved problem. Most of the algorithms aimed at finding the native conformation of a protein complex starting from the structure of its monomers are based on searching the structure corresponding to the global minimum of a suitable scoring function. However, protein complexes are often highly flexible, with mobile side chains and transient contacts due to thermal fluctuations. Flexibility can be neglected if one aims at finding quickly the approximate structure of the native complex, but may play a role in structure refinement, and in discriminating solutions characterized by similar scores. We here benchmark the capability of some state-of-the-art scoring functions (BACH-SixthSense, PIE/PISA and Rosetta) in discriminating finite-temperature ensembles of structures corresponding to the native state and to non-native configurations. We produce the ensembles by running thousands of molecular dynamics simulations in explicit solvent starting from poses generated by rigid docking and optimized in vacuum. We find that while Rosetta outperformed the other two scoring functions in scoring the structures in vacuum, BACH-SixthSense and PIE/PISA perform better in distinguishing near-native ensembles of structures generated by molecular dynamics in explicit solvent. Proteins 2016; 84:1312-1320. © 2016 Wiley Periodicals, Inc.

Differential Water Thermodynamics Determine PI3K-Beta/Delta Selectivity for Solvent-Exposed Ligand Modifications

Phosphoinositide 3-kinases (PI3Ks) are involved in important cellular functions and represent desirable targets for drug discovery efforts, especially related to oncology; however, the four PI3K subtypes (α, β, γ, and δ) have highly similar binding sites, making the design of selective inhibitors challenging. A series of inhibitors with selectivity toward the β subtype over δ resulted in compound 3(S), which has entered a phase I/Ib clinical trial for patients with advanced PTEN-deficient cancer. Interestingly, X-ray crystallography revealed that the modifications making inhibitor 3(S) and related compounds selective toward the β-isoform do not interact directly with either PI3Kβ or PI3Kδ, thereby confounding rationalization of the SAR. Here, we apply explicit solvent molecular dynamics and solvent thermodynamic analysis using WaterMap in an effort to understand the unusual affinity and selectivity trends. We find that differences in solvent energetics and water networks, which are modulated upon binding of different ligands, explain the experimental affinity and selectivity trends. This study highlights the critical role of water molecules in molecular recognition and the importance of considering water networks in drug discovery efforts to rationalize and improve selectivity.