Thermal Fluctuations in Amphipol A8-35 Particles: A Neutron Scattering and Molecular Dynamics Study

Amphipols are a class of polymeric surfactants that can stabilize membrane proteins in aqueous solutions as compared to detergents. A8-35, the best-characterized amphipol to date, is composed of a polyacrylate backbone with ~35 % of the carboxylates free, ~25 % grafted with octyl side-chains, and ~40 % with isopropyl ones. In aqueous solutions, A8-35 self-organizes into globular particles with a molecular mass of ~40 kDa. The thermal dynamics of A8-35 particles was measured by neutron scattering in the 10-picosecond, 18-picosecond, and 1-nanosecond time-scales on natural abundance and deuterium-labeled molecules, which permitted to separate backbone and side-chain motions. A parallel analysis was performed on molecular dynamics trajectories (Perlmutter et al., Langmuir 27:10523–10537, 2011). Experimental results and simulations converge, from their respective time-scales, to show that A8-35 particles feature a more fluid hydrophobic core, predominantly containing the octyl chains, and a more rigid solvent-exposed surface, made up predominantly of the hydrophilic polymer backbone. The fluidity of the core is comparable to that of the lipid environment around proteins in the center of biological membranes, as also measured by neutron scattering. The biological activity of proteins depends sensitively on molecular dynamics, which itself is strongly dependent on the immediate macromolecular environment. In this context, the characterization of A8-35 particle dynamics constitutes a step toward understanding the effect of amphipols on membrane protein stability and function.     

Coupling Molecular Dynamics and Deep Learning to Mine Protein Conformational Space

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

     

Rotation-Vibration Coupling in Liquid Water Investigated by Molecular Dynamics Simulation

Abstract

The standard rotation-vibration decoupling approximation for liquid water is tested using a method based on the split of the molecular dipole fluctuation correlation function into pure rotational, pure vibrational and crossed terms. Our results corroborate the total decoupling between vibrational and rotational motions in liquid water at ambient conditions and provide a tool able to test the reliability of this approximation for molecular liquids at different states.     

Coupling Molecular Dynamics Simulations with Experiments for the Rational Design of Indolicidin-Analogous Antimicrobial Peptides

Antimicrobial peptides (AMPs) have attracted much interest in recent years because of their potential use as new-generation antibiotics. Indolicidin (IL) is a 13-residue cationic AMP that is effective against a broad spectrum of bacteria, fungi, and even viruses. Unfortunately, its high hemolytic activity retards its clinical applications. In this study, we adopted molecular dynamics (MD) simulations as an aid toward the rational design of IL analogues exhibiting high antimicrobial activity but low hemolysis. We employed long-timescale, multi-trajectory all-atom MD simulations to investigate the interactions of the peptide IL with model membranes. The lipid bilayer formed by the zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was chosen as the model erythrocyte membrane; lipid bilayers formed from a mixture of POPC and the negatively charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol were chosen to model bacterial membranes. MD simulations with a total simulation time of up to 4 μs revealed the mechanisms of the processes of IL adsorption onto and insertion into the membranes. The packing order of these lipid bilayers presumably correlated to the membrane stability upon IL adsorption and insertion. We used the degree of local membrane thinning and the reduction in the order parameter of the acyl chains of the lipids to characterize the membrane stability. The order of the mixed 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol/POPC lipid bilayer reduced significantly upon the adsorption of IL. On the other hand, although the order of the pure-POPC lipid bilayer was perturbed slightly during the adsorption stage, the value was reduced more dramatically upon the insertion of IL into the membrane's hydrophobic region. The results imply that enhancing IL adsorption on the microbial membrane may amplify its antimicrobial activity, while the degree of hemolysis may be reduced through inhibition of IL insertion into the hydrophobic region of the erythrocyte membrane. In addition, through simulations, we identified the amino acids that are most responsible for the adsorption onto or insertion into the two model membranes. Positive charges are critical to the peptide's adsorption, whereas the presence of hydrophobic Trp8 and Trp9 leads to its deeper insertion. Combining the hypothetical relationships between the membrane disordering and the antimicrobial and hemolytical activities with the simulated results, we designed three new IL-analogous peptides: IL-K7 (Pro7 → Lys), IL-F89 (Trp8 and Trp9 → Phe), and IL-K7F89 (Pro7 → Lys; Trp8 and Trp9 → Phe). The hemolytic activity of IL-F89 is considerably lower than that of IL, whereas the antimicrobial activity of IL-K7 is greatly enhanced. In particular, the de novo peptide IL-K7F89 exhibits higher antimicrobial activity against Escherichia coli; its hemolytic activity decreased to only 10% of that of IL. Our simulated and experimental results correlated well. This approach—coupling MD simulations with experimental design—is a useful strategy toward the rational design of AMPs for potential therapeutic use.     

Molecular Dynamics Simulations of Peptides from the Central Domain of Smooth Muscle Caldesmon

Abstract

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable α-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of α-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from α-helical structure in all the simulations arise mainly from the formation of long stretches of π-helix. Interconversion between α-helical and π-helical conformations occurs through insertion of water molecules into α-helical hydrogen bonds and subsequent formation of reverse turns. The α-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the π-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of α-helical structure.     

Thermal Stability of Apolipoprotein A-I in High-Density Lipoproteins by Molecular Dynamics

Apolipoprotein (apo) A-I is an unusually flexible protein whose lipid-associated structure is poorly understood. Thermal denaturation, which is used to measure the global helix stability of high-density lipoprotein (HDL)-associated apoA-I, provides no information about local helix stability. Here we report the use of temperature jump molecular dynamics (MD) simulations to scan the per-residue helix stability of apoA-I in phospholipid-rich HDL. When three 20 ns MD simulations were performed at 500 K on each of two particles created by MD simulations at 310 K, bilayers remained intact but expanded by 40%, and total apoA-I helicity decreased from 95% to 72%. Of significance, the conformations of the overlapping N- and C-terminal domains of apoA-I in the particles were unusually mobile, exposing hydrocarbon regions of the phospholipid to solvent; a lack of buried interhelical salt bridges in the terminal domains correlated with increased mobility. Nondenaturing gradient gels show that 40% expansion of the phospholipid surface of 100:2 particles by addition of palmitoyloleoylphosphatidylcholine exceeds the threshold of particle stability. As a unifying hypothesis, we propose that the terminal domains of apoA-I are phospholipid concentration-sensitive molecular triggers for fusion/remodeling of HDL particles. Since HDL remodeling is necessary for cholesterol transport, our model for remodeling has substantial biomedical implications.     

Molecular Dynamics of Thermoenzymes at High Temperature and Pressure: A Review

Lipases are known for their versatility in addition to their ability to digest fat. They can be used for the formulation of detergents, as food ingredients and as biocatalysts in many industrial processes. Because conventional enzymes are frangible at high temperatures, the replacement of conventional chemical routes with biochemical processes that utilize thermostable lipases is vital in the industrial setting. Recent theoretical studies on enzymes have provided numerous fundamental insights into the structures, folding mechanisms and stabilities of these proteins. The studies corroborate the experimental results and provide additional information regarding the structures that were determined experimentally. In this paper, we review the computational studies that have described how temperature affects the structure and dynamics of thermoenzymes, including the thermoalkalophilic L1 lipase derived from Bacillus stearothermophilus. We will also discuss the potential of using pressure for the analysis of the stability of thermoenzymes because high pressure is also important for the processing and preservation of foods.     

Gibbs-Ensemble Molecular Dynamics: A New Method for Simulations Involving Particle Exchange

We discuss a novel simulation method suitable for simulating phenomena involving particle exchange. The method is a molecular dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles or molecules in a thermodynamically consistent fashion. We discuss the general idea of the Gibbs-Ensemble Molecular Dynamics technique and present examples for different simple atomic and molecular fluids. Specifically we will discuss Gibbs-Ensemble Molecular Dynamics simulations of gas-liquid and liquid-solid equilibria in Lennard-Jones systems and in hexane as well as an application of the method to adsorption.     

分子动力学模拟与分子生物力学

生物大分子的微观结构动力学决定其生物学功能,其力学-化学耦合规律是分子生物力学的重点关注方向。分子动力学模拟是耦合生物大分子力学-化学性质微观结构动力学基础的有效手段,其结果可用于预测结构-功能关系、指导实验设计和诠释实验结果。本文简要介绍了分子动力学模拟的方法学特点、基本工作原理及其在分子生物力学中的应用,并展望了未来可能的发展方向和应用前景。     

Brownian Dynamics Simulations of Colloidal Liquids: Hydrodynamics and Stress Relaxation

Abstract

We describe the statistical mechanics background and additional algorithmic features of a recently proposed simple mean-field Brownian Dynamics algorithm formulated to include many-body hydrodynamics, using a local density approximation for the friction coefficient. We show that the equations of motion satisfy the incompressibility of phase space. We make further developments to the model, computing the hydrodynamic effects on the shear stress relaxation function. We show that stress relaxation takes place over two well-defined regimes, in both cases with and without mean field hydrodynamics, MFH. At short times ta ²/D ₀ < 10^?3, where a is the radius of the colloidal particle and D ₀ is the self-diffusion coefficient at infinite dilution, decay of the stress autocorrelation function, C_s(t) is essentially independent of volume fraction and does not fit to a simple analytic form. At longer times than ta ²/D ₀ < 10^?2 the decay has the fractional exponential form ～exp(-t ^β) with β ? 1. The transition between these two regimes coincides with a rapid fall in the time-dependent diffusion coefficient from the so-called short-time to long-time values. We do not find any evidence for power law decay in the C_s(t) as predicted by recent mode-coupling based analytical expansions.     

Capturing Intracellular pH Dynamics by Coupling Its Molecular Mechanisms within a Fully Tractable Mathematical Model

We describe the construction of a fully tractable mathematical model for intracellular pH. This work is based on coupling the kinetic equations depicting the molecular mechanisms for pumps, transporters and chemical reactions, which determine this parameter in eukaryotic cells. Thus, our system also calculates the membrane potential and the cytosolic ionic composition. Such a model required the development of a novel algebraic method that couples differential equations for slow relaxation processes to steady-state equations for fast chemical reactions. Compared to classical heuristic approaches based on fitted curves and ad hoc constants, this yields significant improvements. This model is mathematically self-consistent and allows for the first time to establish analytical solutions for steady-state pH and a reduced differential equation for pH regulation. Because of its modular structure, it can integrate any additional mechanism that will directly or indirectly affect pH. In addition, it provides mathematical clarifications for widely observed biological phenomena such as overshooting in regulatory loops. Finally, instead of including a limited set of experimental results to fit our model, we show examples of numerical calculations that are extremely consistent with the wide body of intracellular pH experimental measurements gathered by different groups in many different cellular systems.     

Collective Fluctuations in the Dynamics of Adaptation and Other Traveling Waves

The dynamics of adaptation are difficult to predict because it is highly stochastic even in large populations. The uncertainty emerges from random genetic drift arising in a vanguard of particularly fit individuals of the population. Several approaches have been developed to analyze the crucial role of genetic drift on the expected dynamics of adaptation, including the mean fitness of the entire population, or the fate of newly arising beneficial deleterious mutations. However, little is known about how genetic drift causes fluctuations to emerge on the population level, where it becomes palpable as variations in the adaptation speed and the fitness distribution. Yet these phenomena control the decay of genetic diversity and variability in evolution experiments and are key to a truly predictive understanding of evolutionary processes. Here, we show that correlations induced by these emergent fluctuations can be computed at any arbitrary order by a suitable choice of a dynamical constraint. The resulting linear equations exhibit fluctuation-induced terms that amplify short-distance correlations and suppress long-distance ones. These terms, which are in general not small, control the decay of genetic diversity and, for wave-tip dominated (“pulled”) waves, lead to anticorrelations between the tip of the wave and the lagging bulk of the population. While it is natural to consider the process of adaptation as a branching random walk in fitness space subject to a constraint (due to finite resources), we show that other traveling wave phenomena in ecology and evolution likewise fall into this class of constrained branching random walks. Our methods, therefore, provide a systematic approach toward analyzing fluctuations in a wide range of population biological processes, such as adaptation, genetic meltdown, species invasions, or epidemics.     

Direct Measurements of Local Coupling between Myosin Molecules Are Consistent with a Model of Muscle Activation

Muscle contracts due to ATP-dependent interactions of myosin motors with thin filaments composed of the proteins actin, troponin, and tropomyosin. Contraction is initiated when calcium binds to troponin, which changes conformation and displaces tropomyosin, a filamentous protein that wraps around the actin filament, thereby exposing myosin binding sites on actin. Myosin motors interact with each other indirectly via tropomyosin, since myosin binding to actin locally displaces tropomyosin and thereby facilitates binding of nearby myosin. Defining and modeling this local coupling between myosin motors is an open problem in muscle modeling and, more broadly, a requirement to understanding the connection between muscle contraction at the molecular and macro scale. It is challenging to directly observe this coupling, and such measurements have only recently been made. Analysis of these data suggests that two myosin heads are required to activate the thin filament. This result contrasts with a theoretical model, which reproduces several indirect measurements of coupling between myosin, that assumes a single myosin head can activate the thin filament. To understand this apparent discrepancy, we incorporated the model into stochastic simulations of the experiments, which generated simulated data that were then analyzed identically to the experimental measurements. By varying a single parameter, good agreement between simulation and experiment was established. The conclusion that two myosin molecules are required to activate the thin filament arises from an assumption, made during data analysis, that the intensity of the fluorescent tags attached to myosin varies depending on experimental condition. We provide an alternative explanation that reconciles theory and experiment without assuming that the intensity of the fluorescent tags varies.     

Virus-Bacterium Coupling Driven by both Turbidity and Hydrodynamics in an Amazonian Floodplain Lake

The importance of viruses in aquatic ecosystem functioning has been widely described. However, few studies have examined tropical aquatic ecosystems. Here, we evaluated for the first time viruses and their relationship with other planktonic communities in an Amazonian freshwater ecosystem. Coupling between viruses and bacteria was studied, focusing both on hydrologic dynamics and anthropogenic forced turbidity in the system (Lake Batata). Samples were taken during four hydrologic seasons at both natural and impacted sites to count virus-like particles (VLP) and bacteria. In parallel, virus-infected bacteria were identified and quantified by transmission electron microscopy (TEM). Viral abundance ranged from 0.5 × 10⁷ ± 0.2 × 10⁷ VLP ml⁻¹ (high-water season, impacted site) to 1.7 × 10⁷ ± 0.4 × 10⁷ VLP ml⁻¹ (low-water season, natural site). These data were strongly correlated with the bacterial abundance (r² = 0.84; P < 0.05), which ranged from 1.0 × 10⁶ ± 0.5 × 10⁶ cells ml⁻¹ (high water, impacted site) to 3.4 × 10⁶ ± 0.7 × 10⁶ cells ml⁻¹ (low water, natural site). Moreover, the viral abundance was weakly correlated with chlorophyll a, suggesting that most viruses were bacteriophages. TEM quantitative analyses revealed that the frequency of visibly infected cells was 20%, with 10 ± 3 phages per cell section. In general, we found a low virus-bacterium ratio (<7). Both the close coupling between the viral and bacterial abundances and the low virus-bacterium ratio suggest that viral abundance tends to be driven by the reduction of hosts for viral infection. Our results demonstrate that viruses are controlled by biological substrates, whereas in addition to grazing, bacteria are regulated by physical processes caused by turbidity, which affect underwater light distribution and dissolved organic carbon availability.Viruses are the most abundant and dynamic components of the aquatic microbial community (6, 31, 32). Viruses influence many biogeochemical and ecological processes, including nutrient cycling, system respiration, particle-size distribution, bacterial and algal biodiversity, species distribution, algal blooms, and genetic transfer between microorganisms (21, 49). In addition, viruses play a major role in aquatic microbial food webs by releasing carbon trapped in host cells to the dissolved organic carbon (DOC) pool and ultimately back to the bacterial community (11, 21). The action of viruses is an important mechanism of bacterial regulation in aquatic ecosystems, acting directly on bacterial populations and indirectly on bacterial diversity by decreasing the density of dominant bacterial species (31). Studies based on viral decay rates and electron microscopy analyses have shown that viruses can cause up to 40% of bacterial mortality and more than 10% of phytoplankton mortality in aquatic systems (11, 22, 48, 50, 54, 55). It also has been suggested that viral lysis and protistan grazing cause similar bacterial mortality in aquatic ecosystems (22, 40).Several environmental factors, including solar radiation and temperature, can influence viral abundance. Exposure to solar radiation decreases viral abundance in aquatic ecosystems, while low temperatures decrease their virulence (33, 56). However, the majority of the studies on virus ecology have been performed in temperate or polar regions, where seasonal changes in solar radiation and water temperature are more pronounced (26, 30, 32). Viral abundances have been little investigated in tropical aquatic ecosystems (6, 39) and particularly in the Amazonian region, where the abundance and activity of aquatic viruses have not been studied.The greatest watershed in the world is located in the Amazonian region. It is composed of clear-water, black-water, and turbid freshwater ecosystems, which are seasonally influenced by the flood pulse. The hydrologic pulse is characterized by a pronounced change in water level, defining the flood seasons. Nutrient sources and stocks and species dynamics vary according to the water level (25). During the high-water season (flood season), the tight connection between terrestrial and aquatic environments results in an increase in allochthonous DOC input and the dilution of inorganic nutrients and organisms. During the low-water season, there is an increase in nutrient concentrations, organism abundances, and the importance of autochthonous DOC. These seasonal changes differently impact ecosystem functions and aquatic community dynamics (5, 9, 18). For instance, bacterioplankton abundance is less changeable than phytoplankton abundance throughout the hydrological cycle due to the alternative sources of DOC (allochthonous in the high-water period and autochthonous in the low-water period) for bacterial communities (5, 24).Lake Batata is a clear-water Amazonian floodplain lake located in the watershed of the Trombetas River, a tributary of the Amazon River. As a clear-water Amazonian ecosystem, it contains low concentrations of suspended particles and inorganic nutrients (47). Lake Batata is distinct because it was impacted by bauxite tailings for 10 years (1979 to 1989), affecting 30% of the lake''s area. The tailings caused a huge increase in turbidity; large amounts of tailings settled on the sediment surface and often are resuspended by physical mixing or biotic movements (28). The presence of tailings resulted in a clear spatial variation in the lake, forming impacted and natural sites. Furthermore, tailing particles can directly act as a substrate for attaching bacteria and also can adsorb DOC (5).Previous studies on bacterio-, phyto-, and zooplankton communities have shown that flood pulse acts as the primary driver of plankton community structure in Lake Batata (5, 9, 24). Bauxite tailings also affect microbial processes in impacted sites, such as bacterial growth and production (5), photosynthesis rates and primary production (43), or the availability of food for zooplankton (8). However, there still is no evidence indicating a complementary (flood pulse and forced turbidity) effect among these factors in any microbial community. Based on published data, we assume that (i) bacterioplankton abundance is less changeable through the hydrologic cycle than phytoplankton abundance (5, 24), and (ii) tailing particles can act as a substrate for attaching bacteria and also can adsorb organic matter, which is controlled by the flood pulse (5). Therefore, we hypothesized that the relationship between viruses and bacteria in Lake Batata is modulated by a synergistic effect between the hydrological cycle and turbidity.     

Dynamic Interreceptor Coupling Contributes to the Consistent Open Duration of Ryanodine Receptors

Ca²⁺ spark is the elementary Ca²⁺ signaling event in muscle excitation-contraction coupling. The rise time of Ca²⁺ spark is rather stable under different conditions, suggesting consistent open duration of ryanodine receptors (RyRs) in vivo. It has been proposed that the array-based behavior of RyRs plays an important role in shaping Ca²⁺ spark dynamics, particularly in controlling the open duration of RyR clusters. Therefore, we investigated the possible role of inter-RyR coupling in stabilization of the open time of arrayed RyRs under several potential perturbations, for instance, array size, inter-RyR coupling noise, and up-regulation or down-regulation of the activity of partial RyRs in the array. We found that RyR arrays with dynamic coupling showed consistent open duration against the perturbations, whereas the RyR array with constant coupling did not. On the other hand, the open probability and amplitude of RyR arrays with dynamic interreceptor coupling were sensitive to the perturbations. These two points were consistent with experimental observations, indicating that the RyR array with dynamic coupling could recapitulate in vivo open properties of RyRs. Our findings support the idea that dynamic coupling is a feasible in vivo working mechanism of RyR arrays.     

Molecular Dynamics Simulation Studies of the Density and Temperature Dependence of Self-diffusion in a Cylindrical Micropore

We report Molecular Dynamics calculations of radial density profiles and self-diffusion coefficients of Lennard-Jones fluids in a cylindrical pore of radius 2σ, for a wide range of temperatures and densities. At n _p σ³ = 0.825 the self-diffusion coefficient parallel to the pore walls D _∥*. follows a monotonic (nearly linear) increase with kT/ε and is very similar to that of the bulk self-diffusion coefficient D _b *. At n _p σ³ = 0.4 and kT/ε ≤ 1.0 the curve of D _∥* vs. kT/ε shows a distinct inflection in the region 0.7 ≤ kT/ε ≤ 0.9 and values of D _∥* are much less than D _b * decreasing to near solid state values at very low temperatures. At the highest temperature studied, kT/ε = 2.98, D _∥* is almost inversely proportional to density and in a fairly close agreement with that of D _b *. At KT/ε = 0.49, D _∥* is much smaller than D _b *. The motion of adsorbate particles normal to the walls is also discussed.