Mechanical response and deformation mechanics of Type IV pili investigated using steered coarse-grained molecular dynamics simulation

Type IV pili are long filamentous structures on the surface of bacteria, which can be rapidly assembled or disassembled with pilin subunits by molecular motors. They can generate force during retraction and are involved in many bacterial functions. Steered molecular dynamics simulations with coarse-grained MARTINI models are carried out to investigate the mechanical behaviors of pili under tension. Our study is the first to report a Young's modulus of 0.80 ± 0.07 GPa and a spring constant of 1294.6 ± 116.5 kJ mol⁻¹ nm⁻² for pilus. Our results show the mechanical responses of pili are different from those described by the worm-like chain model and the van der Waal's interactions play a critical role in the mechanical responses. Moreover, the effects of pulling rates and virtual spring constants of pilus on Young's modulus are studied and two distinct morphological stages with the conformational changes appear during the extension of pilus are observed. This work provide insight into the mechanics and the deformation mechanism of pilus assembly.     

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.     

Elucidation of the glucose transport pathway in glucose transporter 4 via steered molecular dynamics simulations

Background

GLUT4 is a predominant insulin regulated glucose transporter expressed in major glucose disposal tissues such as adipocytes and muscles. Under the unstimulated state, GLUT4 resides within intracellular vesicles. Various stimuli such as insulin translocate this protein to the plasma membrane for glucose transport. In the absence of a crystal structure for GLUT4, very little is known about the mechanism of glucose transport by this protein. Earlier we proposed a homology model for GLUT4 and performed a conventional molecular dynamics study revealing the conformational rearrangements during glucose and ATP binding. However, this study could not explain the transport of glucose through the permeation tunnel.

Methodology/Principal Findings

To elucidate the molecular mechanism of glucose transport and its energetic, a steered molecular dynamics study (SMD) was used. Glucose was pulled from the extracellular end of GLUT4 to the cytoplasm along the pathway using constant velocity pulling method. We identified several key residues within the tunnel that interact directly with either the backbone ring or the hydroxyl groups of glucose. A rotation of glucose molecule was seen near the sugar binding site facilitating the sugar recognition process at the QLS binding site.

Conclusions/Significance

This study proposes a possible glucose transport pathway and aids the identification of several residues that make direct interactions with glucose during glucose transport. Mutational studies are required to further validate the observation made in this study.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

Mechanism of CDK5 activation revealed by steered molecular dynamics simulations and energy calculations

In the current work, CDK5/p25 complexes were pulled apart by applying external forces with steered molecular dynamics (SMD) simulations. The crucial interactions between the kinase and the activation protein were investigated and the SMD simulations showed that several activation-relevant motifs of CDK5 leave p25 in sequence during the pulling and lead to an apo-CDK2 like CDK5 structure after separation. Based on systematic examination of hydrogen bond breaking and classical MD/molecular mechanics-generalized Born/surface area) (MM-GBSA) calculations, a CDK5 activation mechanism by p25 is suggested. This is the first step towards the systemic development of CDK inhibitors and the mechanism proposed could lead to a better understanding of the protein–protein recognition characteristics between the kinase and its activator.     

Conformational properties of glucose-based disaccharides investigated using molecular dynamics simulations with local elevation umbrella sampling

Explicit-solvent molecular dynamics (MD) simulations of the 11 glucose-based disaccharides in water at 300 K and 1 bar are reported. The simulations were carried out with the GROMOS 45A4 force-field and the sampling along the glycosidic dihedral angles ? and ψ was artificially enhanced using the local elevation umbrella sampling (LEUS) method. The trajectories are analyzed in terms of free-energy maps, stable and metastable conformational states (relative free energies and estimated transition timescales), intramolecular H-bonds, single molecule configurational entropies, and agreement with experimental data. All disaccharides considered are found to be characterized either by a single stable (overwhelmingly populated) state ((1→n)-linked disaccharides with n = 1, 2, 3, or 4) or by two stable (comparably populated and differing in the third glycosidic dihedral angle ; gg or gt) states with a low interconversion barrier ((1→6)-linked disaccharides). Metastable (anti-? or anti-ψ) states are also identified with relative free energies in the range of 8-22 kJ mol⁻¹. The 11 compounds can be classified into four families: (i) the α(1→1)α-linked disaccharide trehalose (axial-axial linkage) presents no metastable state, the lowest configurational entropy, and no intramolecular H-bonds; (ii) the four α(1→n)-linked disaccharides (n = 1, 2, 3, or 4; axial-equatorial linkage) present one metastable (anti-ψ) state, an intermediate configurational entropy, and two alternative intramolecular H-bonds; (iii) the four β(1→n)-linked disaccharides (n = 1, 2, 3, or 4; equatorial-equatorial linkage) present two metastable (anti-? and anti-ψ) states, an intermediate configurational entropy, and one intramolecular H-bond; (iv) the two (1→6)-linked disaccharides (additional glycosidic dihedral angle) present no (isomaltose) or a pair of (gentiobiose) metastable (anti-?) states, the highest configurational entropy, and no intramolecular H-bonds. The observed conformational preferences appear to be dictated by four main driving forces (ring conformational preferences, exo-anomeric effect, steric constraints, and possible presence of a third glycosidic dihedral angle), leaving a secondary role to intramolecular H-bonding and specific solvation effects. In spite of the weak conformational driving force attributed to solvent-exposed H-bonds in water (highly polar protic solvent), intramolecular H-bonds may still have a significant influence on the physico-chemical properties of the disaccharide by decreasing its hydrophilicity. Along with previous work, the results also complete the suggestion of a spectrum of approximate transition timescales for carbohydrates up to the disaccharide level, namely: ∼30 ps (hydroxyl groups), ∼1 ns (free lactol group, free hydroxymethyl groups, glycosidic dihedral angle in (1→6)-linked disaccharides), ∼10 ns to 2 μs (ring conformation, glycosidic dihedral angles ? and ψ). The calculated average values of the glycosidic torsional angles agree well with the available experimental data, providing validation for the force-field and simulation methodology employed.     

Mechanical features of various silkworm crystalline considering hydration effect via molecular dynamics simulations

Silk materials are receiving significant attention as base materials for various functional nanomaterials and nanodevices, due to its exceptionally high mechanical properties, biocompatibility, and degradable characteristics. Although crystalline silk regions are composed of various repetitive motifs with differing amino acid sequences, how the effect of humidity works differently on each of the motifs and their structural characteristics remains unclear. We report molecular dynamics (MD) simulations on various silkworm fibroins composed of major motifs (i.e. (GAGAGS)_n, (GAGAGA)_n, and (GAGAGY)_n) at varying degrees of hydration, and reveal how each major motifs of silk fibroins change at each degrees of hydration using MD simulations and their structural properties in mechanical perspective via steered molecular dynamics simulations. Our results explain what effects humidity can have on nanoscale materials and devices consisting of crystalline silk materials.     

Allostery in G protein-coupled receptors investigated by molecular dynamics simulations

1. Download : Download high-res image (101KB)
2. Download : Download full-size image

     

Refinement of protein structure homology models via long, all-atom molecular dynamics simulations

Accurate computational prediction of protein structure represents a longstanding challenge in molecular biology and structure-based drug design. Although homology modeling techniques are widely used to produce low-resolution models, refining these models to high resolution has proven difficult. With long enough simulations and sufficiently accurate force fields, molecular dynamics (MD) simulations should in principle allow such refinement, but efforts to refine homology models using MD have for the most part yielded disappointing results. It has thus far been unclear whether MD-based refinement is limited primarily by accessible simulation timescales, force field accuracy, or both. Here, we examine MD as a technique for homology model refinement using all-atom simulations, each at least 100 μs long-more than 100 times longer than previous refinement simulations-and a physics-based force field that was recently shown to successfully fold a structurally diverse set of fast-folding proteins. In MD simulations of 24 proteins chosen from the refinement category of recent Critical Assessment of Structure Prediction (CASP) experiments, we find that in most cases, simulations initiated from homology models drift away from the native structure. Comparison with simulations initiated from the native structure suggests that force field accuracy is the primary factor limiting MD-based refinement. This problem can be mitigated to some extent by restricting sampling to the neighborhood of the initial model, leading to structural improvement that, while limited, is roughly comparable to the leading alternative methods.     

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer’s disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.     

Aggregation of small peptides studied by molecular dynamics simulations

Peptides and proteins tend to aggregate under appropriate conditions. The amyloid fibrils that are ubiquitously found among these structures are associated with major human diseases like Alzheimer's disease, type II diabetes, and various prion diseases. Lately, it has been observed that even very short peptides like tetra and pentapeptides can form ordered amyloid structures. Here, we present aggregation studies of three such small polypeptide systems, namely, the two amyloidogenic peptides DFNKF and FF, and a control (nonamyloidogenic) one, the AGAIL. The respective aggregation process is studied by all-atom Molecular Dynamics simulations, which allow to shed light on the fine details of the association and aggregation process. Our analysis suggests that naturally aggregating systems exhibit significantly diverse overall cluster shape properties and specific intermolecular interactions. Additional analysis was also performed on the previously studied NFGAIL system.     

The pathway of oligomeric DNA melting investigated by molecular dynamics simulations

Details of the reaction coordinate for DNA melting are fundamental to much of biology and biotechnology. Recently, it has been shown experimentally that there are at least three states involved. To clarify the reaction mechanism of the melting transition of DNA, we perform 100-ns molecular dynamics simulations of a homo-oligomeric, 12-basepair DNA duplex, d(A₁₂)·d(T₁₂), with explicit salt water at 400 K. Analysis of the trajectory reveals the various biochemically important processes that occur on different timescales. Peeling (including fraying from the ends), searching for Watson-Crick complements, and dissociation are recognizable processes. However, we find that basepair searching for Watson-Crick complements along a strand is not mechanistically tied to or directly accessible from the dissociation steps of strand melting. A three-step melting mechanism is proposed where the untwisting of the duplex is determined to be the major component of the reaction coordinate at the barrier. Though the observations are limited to the characteristics of the system being studied, they provide important insight into the mechanism of melting of other more biologically relevant forms of DNA, which will certainly differ in details from those here.     

Free-energy landscape of glycerol permeation through aquaglyceroporin GlpF determined from steered molecular dynamics simulations

The free-energy landscape of glycerol permeation through the aquaglyceroporin GlpF has been estimated in the literature by the nonequilibrium method of steered molecular dynamics (SMD) simulations and by the equilibrium method of adaptive biasing force (ABF) simulations. However, the ABF results qualitatively disagree with the SMD results that were based on the Jarzynski equality (JE) relating the equilibrium free-energy difference to the nonequilibrium work of the irreversible pulling experiments. In this paper, I present a new SMD study of the glycerol permeation through GlpF to explore the free-energy profile of glycerol along the permeation channel. Instead of the JE in terms of thermodynamic work, I use the fluctuation-dissipation theorem (FDT) of Brownian dynamics (BD), in terms of mechanical work, for extracting the free-energy difference from the nonequilibrium work of irreversible pulling experiments. The results of this new SMD-BD-FDT study are in agreement with the experimental data and with the ABF results.     

Transmembrane peptides from tyrosine kinase receptor. Mutation-related behavior in a lipid bilayer investigated by molecular dynamics simulations

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.     

Exploring the energetics of water permeation in photosystem II by multiple steered molecular dynamics simulations

The Mn(4)Ca cluster of the oxygen-evolving complex (OEC) of photosynthesis catalyzes the light-driven splitting of water into molecular oxygen, protons, and electrons. The OEC is buried within photosystem II (PSII), a multisubunit integral membrane protein complex, and water must find its way to the Mn(4)Ca cluster by moving through protein. Molecular dynamics simulations were used to determine the energetic barriers for water permeation though PSII extrinsic proteins. Potentials of mean force (PMFs) for water were derived by using the technique of multiple steered molecular dynamics (MSMD). Calculation of free energy profiles for water permeation allowed us to characterize previously identified water channels, and discover new pathways for water movement toward the Mn(4)Ca cluster. Our results identify the main constriction sites in these pathways which may serve as selectivity filters that restrict both the access of solutes detrimental to the water oxidation reaction and loss of Ca(2+) and Cl(-) from the active site.     

Nanocontraction flows of short-chain polyethylene via molecular dynamics simulations

Nanocontraction flows of liquid short-chain polyethylene ([CH₂]₅₀) that were uniformly extruded by a constant-speed piston into a surrounding vacuum from a reservoir through an abrupt contraction nozzle were performed by employing molecular dynamics simulations. The extrudate exhibits a similar die swell phenomenon around the exit of the nozzle. In addition, numerous molecular chains are strongly adsorbed on the external surface of the nozzle. At high extrusion speeds, the velocity and temperature profiles in the nozzle show convex and concave parabolic curves, respectively, whereas the profiles are relatively flat at lower speeds. Near the internal boundary of the nozzle, the wall slip is inspected. Significantly, during the flow, the molecular chains undergo structural deformation, including compressed, stretched and shrunk motions. Comparisons with related experimental observations show that the simulated probability distributions of the bending and dihedral angles, and variations of the squared radius of gyration and orientations, are in reasonable agreement.     

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics‐based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge‐based or implicit membrane‐based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid‐facing residues. Scoring with knowledge‐based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane‐based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.     

Investigation of formate transport through the substrate channel of formate dehydrogenase by steered molecular dynamics simulations

Steered molecular dynamics simulation has revealed the mechanism of formate transport via the substrate channel of formate dehydrogenase. It is shown that the structural organization of the channel promotes the transport of formate anion in spite of the fact that the channel is too narrow even for such a small molecule. The conformational mobility of Arg284 residue, one of the residues forming the wall of the substrate channel, provides for the binding and delivery of formate to the active site.