Glycan flexibility: insights into nanosecond dynamics from a microsecond molecular dynamics simulation explaining an unusual nuclear Overhauser effect

An atomistic all-atom molecular dynamics simulation of the trisaccharide β-d-ManpNAc-(1→4)[α-d-Glcp-(1→3)]-α-l-Rhap-OMe with explicit solvent molecules has been carried out. The trisaccharide represents a model for the branching region of the O-chain polysaccharide of a strain from Aeromonas salmonicida. The extensive MD simulations having a 1-μs duration revealed a conformational dynamics process on the nanosecond time scale, that is, a ‘time window’ not extensively investigated for carbohydrates to date. The results obtained from the MD simulation underscore the predictive power of molecular simulations in studies of biomolecular systems and also explain an unusual nuclear Overhauser effect originating from conformational exchange.     

Striking HIV-1 Entry by Targeting HIV-1 gp41. But,Where Should We Target?

HIV-1 gp41 facilitates the viral fusion through a conformational switch involving the association of three C-terminal helices along the conserved hydrophobic grooves of three N-terminal helices coiled-coil. The control of these structural rearrangements is thought to be central to HIV-1 entry and, therefore, different strategies of intervention are being developed. Herewith, we describe a procedure to simulate the folding of an HIV-1 gp41 simplified model. This procedure is based on the construction of plausible conformational pathways, which describe protein transition between non-fusogenic and fusogenic conformations. The calculation of the paths started with 100 molecular dynamics simulations of the non-fusogenic conformation, which were found to converge to different intermediate states. Those presenting defined criteria were selected for separate targeted molecular dynamics simulations, subjected to a force constant imposing a movement towards the gp41 fusogenic conformation. Despite significant diversity, a preferred sequence of events emerged when the simulations were analyzed in terms of the formation, breakage and evolution of the contacts. We pointed out 29 residues as the most relevant for the movement of gp41; also, 2696 possible interactions were reduced to only 48 major interactions, which reveals the efficiency of the method. The analysis of the evolution of the main interactions lead to the detection of four main behaviors for those contacts: stable, increasing, decreasing and repulsive interactions. Altogether, these results suggest a specific small cavity of the HIV-1 gp41 hydrophobic groove as the preferred target to small molecules.     

High performance computing in biology: multimillion atom simulations of nanoscale systems

Computational methods have been used in biology for sequence analysis (bioinformatics), all-atom simulation (molecular dynamics and quantum calculations), and more recently for modeling biological networks (systems biology). Of these three techniques, all-atom simulation is currently the most computationally demanding, in terms of compute load, communication speed, and memory load. Breakthroughs in electrostatic force calculation and dynamic load balancing have enabled molecular dynamics simulations of large biomolecular complexes. Here, we report simulation results for the ribosome, using approximately 2.64 million atoms, the largest all-atom biomolecular simulation published to date. Several other nano-scale systems with different numbers of atoms were studied to measure the performance of the NAMD molecular dynamics simulation program on the Los Alamos National Laboratory Q Machine. We demonstrate that multimillion atom systems represent a 'sweet spot' for the NAMD code on large supercomputers. NAMD displays an unprecedented 85% parallel scaling efficiency for the ribosome system on 1024 CPUs. We also review recent targeted molecular dynamics simulations of the ribosome that prove useful for studying conformational changes of this large biomolecular complex in atomic detail.     

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.     

Excluded volume approximation to protein-solvent interaction. The solvent contact model.

Important properties of globular proteins, such as the stability of its folded state, depend sensitively on interactions with solvent molecules. Existing methods for estimating these interactions, such as the geometrical surface model, are either physically misleading or too time consuming to be applied routinely in energy calculations. As an alternative, we derive here a simple model for the interactions between protein atoms and solvent atoms in the first hydration layer, the solvent contact model, based on the conservation of the total number of atomic contacts, a consequence of the excluded-volume effect. The model has the conceptual advantage that protein-protein contacts and protein-solvent contacts are treated in the same language and the technical advantage that the solvent term becomes a particularly simple function of interatomic distances. The model allows rapid calculation of any physical property that depends only on the number and type of protein-solvent nearest-neighbor contacts. We propose use of the method in the calculation of protein solvation energies, conformational energy calculations, and molecular dynamics simulations.     

Pressure dependence of vibrational optical activity of model biomolecules. A computational study

Change of molecular properties with pressure is an attracting means to regulate molecular reactivity or biological activity. However, the effect is usually small and so far explored rather scarcely. To obtain a deeper insight and estimate the sensitivity of vibrational optical activity spectra to pressure-induced conformational changes, we investigate small model molecules. The Ala-Ala dipeptide, isomaltose disaccharide and adenine-uracil dinucleotide were chosen to represent three different biomolecular classes. The pressure effects were modeled by molecular dynamics and density functional theory simulations. The dinucleotide was found to be the most sensitive to the pressure, whereas for the disaccharide the smallest changes are predicted. Pressure-induced relative intensity changes in vibrational circular dichroism and Raman optical activity spectra are predicted to be 2–3-times larger than for non-polarized IR and Raman techniques.     

Combined scanning probe and total internal reflection fluorescence microscopy

Combining scanning probe and optical microscopy represents a powerful approach for investigating structure-function relationships and dynamics of biomolecules and biomolecular assemblies, often in situ and in real-time. This platform technology allows us to obtain three-dimensional images of individual molecules with nanometer resolution, while simultaneously characterizing their structure and interactions though complementary techniques such as optical microscopy and spectroscopy. We describe herein the practical strategies for the coupling of scanning probe and total internal reflection fluorescence microscopy along with challenges and the potential applications of such platforms, with a particular focus on their application to the study of biomolecular interactions at membrane surfaces.     

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA''s secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.     

Multiscale methods for macromolecular simulations

In this article we review the key modeling tools available for simulating biomolecular systems. We consider recent developments and representative applications of mixed quantum mechanics/molecular mechanics (QM/MM), elastic network models (ENMs), coarse-grained molecular dynamics, and grid-based tools for calculating interactions between essentially rigid protein assemblies. We consider how the different length scales can be coupled, both in a sequential fashion (e.g. a coarse-grained or grid model using parameterization from MD simulations), and via concurrent approaches, where the calculations are performed together and together control the progression of the simulation. We suggest how the concurrent coupling approach familiar in the context of QM/MM calculations can be generalized, and describe how this has been done in the CHARMM macromolecular simulation package.     

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.     

Transportsomes and channelsomes: are they functional units for physiological responses?

Channels and transporters play essential biological roles primarily through the transportation of ions and small molecules that are required to maintain cellular activities across the biomembrane. Secondary to transportation, channels and transporters also integrate and coordinate biological functions at different levels, ranging from the subcellular (nm) to multicellular (μm) scales. This is underpinned by efficient functional coupling within molecular assemblies of channels, transporters, proteins, small molecules, and lipids. Molecular interactions create local microenvironments that, in some cases, uniquely modify the functional properties of the channels and transporters. These molecular assemblies built around a transporter or channel (“transportsomes” and “channelsomes”) can be considered as physiological functional units. In this special issue, we provide an overview of recent progress in our understanding of protein-protein and molecular interactions in transportsomes and channelsomes, which occur through both direct molecular contacts and more distal functional coupling, and examine the validity of these “somes”.     

Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange

Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.     

MC-PHS: A Monte Carlo Implementation of the Primary Hydration Shell for Protein Folding and Design

A primary hydration shell (PHS) approach is developed for Monte Carlo simulations of conformationally rich macromolecular systems in an environment that efficiently captures principal solvation effects. It has been previously demonstrated that molecular dynamics using PHS is an efficient method to study peptide structure and dynamics in aqueous solution. Here, we extend the PHS approach to Monte Carlo simulations, whereby a stable shell of water molecules is maintained with a flexible, nonspherical, half-harmonic potential, tuned to maintain a constant restraining energy, with the difference between the restraint and shell energies used to dynamically adjust the shell radius. Examination of the shell and system size dependence of the restraining potential reveals its robustness. Moreover, its suitability for biomolecular simulations is evaluated using small spheres of water, hydration properties of small biological molecules, and configurational sampling of β-hairpin pentapeptide YPGDV. This method, termed MC-PHS, appears to provide efficient representation of dominant solvation effects and should prove useful in the study of protein folding and design.     

Measurement of Mesoscale Conformational Dynamics of Freely Diffusing Molecules with Tracking FCS

Few techniques are suited to probe the structure and dynamics of molecular complexes at the mesoscale level ($～$100–1000 nm). We have developed a single-molecule technique that uses tracking fluorescence correlation spectroscopy (tFCS) to probe the conformation and dynamics of mesoscale molecular assemblies. tFCS measures the distance fluctuations between two fluorescently labeled sites within an untethered, freely diffusing biomolecule. To achieve subdiffraction spatial resolution, we developed a feedback scheme that allows us to maintain the molecule at an optimal position within the laser intensity gradient for fluorescence correlation spectroscopy. We characterized tFCS spatial sensitivity by measuring the Brownian end-to-end dynamics of DNA molecules as short as 1000 bp. We demonstrate that tFCS detects changes in the compaction of reconstituted nucleosome arrays and can assay transient protein-mediated interactions between distant sites in an individual DNA molecule. Our measurements highlight the applicability of tFCS to a wide variety of biochemical processes involving mesoscale conformational dynamics.     

Phosphorylation-induced conformational changes in a mitogen-activated protein kinase substrate. Implications for tyrosine hydroxylase activation

Mitogen-activated protein (MAP) kinase-mediated phosphorylation of specific residues in tyrosine hydroxylase leads to an increase in enzyme activity. However, the mechanism whereby phosphorylation affects enzyme turnover is not well understood. We used a combination of fluorescence resonance energy transfer (FRET) measurements and molecular dynamics simulations to explore the conformational free energy landscape of a 10-residue MAP kinase substrate found near the N terminus of the enzyme. This region is believed to be part of an autoregulatory sequence that overlies the active site of the enzyme. FRET was used to measure the effect of phosphorylation on the ensemble of peptide conformations, and molecular dynamics simulations generated free energy profiles for both the unphosphorylated and phosphorylated peptides. We demonstrate how FRET transfer efficiencies can be calculated from molecular dynamics simulations. For both the unphosphorylated and phosphorylated peptides, the calculated FRET efficiencies are in excellent agreement with the experimentally determined values. Moreover, the FRET measurements and molecular simulations suggest that phosphorylation causes the peptide backbone to change direction and fold into a compact structure relative to the unphosphorylated state. These results are consistent with a model of enzyme activation where phosphorylation of the MAP kinase substrate causes the N-terminal region to adopt a compact structure away from the active site. The methods we employ provide a general framework for analyzing the accessible conformational states of peptides and small molecules. Therefore, they are expected to be applicable to a variety of different systems.     

Non-Ewald methods: theory and applications to molecular systems

Several non-Ewald methods for calculating electrostatic interactions have recently been developed, such as the Wolf method, the reaction field method, the pre-averaging method, and the zero-dipole summation method, for molecular dynamics simulations of various physical systems, including biomolecular systems. We review the theories of these approaches and their potential applications to molecular simulations, and discuss their relationships.     

Fuzziness in Protein Interactions—A Historical Perspective

The proposal that coupled folding to binding is not an obligatory mechanism for intrinsically disordered (ID) proteins was put forward 10 years ago. The notion of fuzziness implies that conformational heterogeneity can be maintained upon interactions of ID proteins, which has a functional impact either on regulated assembly or activity of the corresponding complexes. Here I review how the concept has evolved in the past decade, via increasing experimental data providing insights into the mechanisms, pathways and regulatory modes. The effects of structural diversity and transient contacts on protein assemblies have been collected and systematically analyzed (Fuzzy Complexes Database, http://protdyn-database.org). Fuzziness has also been exploited as a framework to decipher molecular organization of higher-order protein structures. Quantification of conformational heterogeneity opens exciting future perspectives for drug discovery from small molecule–ID protein interactions to supramolecular assemblies.     

Searching for the molecular arrangement of transmembrane ceramide channels

Ceramides have been implicated in the initiation of apoptosis by permeabilizing the mitochondrial outer membrane to small proteins, including cytochrome c. In addition, ceramides were shown to form large metastable channels in planar membranes and liposomes, indicating that these lipids permeabilize membranes directly. Here we analyze molecular models of ceramide channels and test their stability in molecular dynamics simulations. The structural units are columns of four to six ceramides H-bonded via amide groups and arranged as staves in either a parallel or antiparallel manner. Two cylindrical assemblies of 14 columns (four or six molecules per column) were embedded in a fully hydrated palmitoyloleoyl-phosphatidylcholine phospholipid bilayer, and simulated for 24 ns in total. After equilibration, the water-filled pore adopted an hourglass-like shape as headgroups of ceramides and phospholipids formed a smooth continuous interface. The structure-stabilizing interactions were both hydrogen bonds between the headgroups (including water-mediated interactions) and packing of the hydrocarbon tails. Ceramide's essential double bond reduced the mobility of the hydrocarbon tails and stabilized their packing. The six-column assembly remained stable throughout a 10-ns simulation. During simulations of four-column assemblies, pairs of columns displayed the tendency of splitting out from the channels, consistent with the previously proposed mechanism of channel disassembly.     

CN-GELFrEE - Clear Native Gel-eluted Liquid Fraction Entrapment Electrophoresis

Protein complexes perform an array of crucial cellular functions. Elucidating their non-covalent interactions and dynamics is paramount for understanding the role of complexes in biological systems. While the direct characterization of biomolecular assemblies has become increasingly important in recent years, native fractionation techniques that are compatible with downstream analysis techniques, including mass spectrometry, are necessary to further expand these studies. Nevertheless, the field lacks a high-throughput, wide-range, high-recovery separation method for native protein assemblies. Here, we present clear native gel-eluted liquid fraction entrapment electrophoresis (CN-GELFrEE), which is a novel separation modality for non-covalent protein assemblies. CN-GELFrEE separation performance was demonstrated by fractionating complexes extracted from mouse heart. Fractions were collected over 2 hr and displayed discrete bands ranging from ~30 to 500 kDa. A consistent pattern of increasing molecular weight bandwidths was observed, each ranging ~100 kDa. Further, subsequent reanalysis of native fractions via SDS-PAGE showed molecular-weight shifts consistent with the denaturation of protein complexes. Therefore, CN-GELFrEE was proved to offer the ability to perform high-resolution and high-recovery native separations on protein complexes from a large molecular weight range, providing fractions that are compatible with downstream protein analyses.     

Atomistic View of the Conformational Activation of Src Kinase Using the String Method with Swarms-of-Trajectories

The inactive-to-active conformational transition of the catalytic domain of human c-Src tyrosine kinase is characterized using the string method with swarms-of-trajectories with all-atom explicit solvent molecular dynamics simulations. The activation process occurs in two main steps in which the activation loop (A-loop) opens first, followed by the rotation of the αC helix. The computed potential of mean force energy along the activation pathway displays a local minimum, which allows the identification of an intermediate state. These results show that the string method with swarms-of-trajectories is an effective technique to characterize complex and slow conformational transitions in large biomolecular systems.