Local feature analysis: a statistical theory for reproducible essential dynamics of large macromolecules

Multivariate statistical methods are widely used to extract functional collective motions from macromolecular molecular dynamics (MD) simulations. In principal component analysis (PCA), a covariance matrix of positional fluctuations is diagonalized to obtain orthogonal eigenvectors and corresponding eigenvalues. The first few eigenvectors usually correspond to collective modes that approximate the functional motions in the protein. However, PCA representations are globally coherent by definition and, for a large biomolecular system, do not converge on the time scales accessible to MD. Also, the forced orthogonalization of modes leads to complex dependencies that are not necessarily consistent with the symmetry of biological macromolecules and assemblies. Here, we describe for the first time the application of local feature analysis (LFA) to construct a topographic representation of functional dynamics in terms of local features. The LFA representations are low dimensional, and like PCA provide a reduced basis set for collective motions, but they are sparsely distributed and spatially localized. This yields a more reliable assignment of essential dynamics modes across different MD time windows. Also, the intrinsic dynamics of local domains is more extensively sampled than that of globally coherent PCA modes.     

An NMR database for simulations of membrane dynamics

Computational methods are powerful in capturing the results of experimental studies in terms of force fields that both explain and predict biological structures. Validation of molecular simulations requires comparison with experimental data to test and confirm computational predictions. Here we report a comprehensive database of NMR results for membrane phospholipids with interpretations intended to be accessible by non-NMR specialists. Experimental 13C-1H and 2H NMR segmental order parameters (S(CH) or S(CD)) and spin-lattice (Zeeman) relaxation times (T(1Z)) are summarized in convenient tabular form for various saturated, unsaturated, and biological membrane phospholipids. Segmental order parameters give direct information about bilayer structural properties, including the area per lipid and volumetric hydrocarbon thickness. In addition, relaxation rates provide complementary information about molecular dynamics. Particular attention is paid to the magnetic field dependence (frequency dispersion) of the NMR relaxation rates in terms of various simplified power laws. Model-free reduction of the T(1Z) studies in terms of a power-law formalism shows that the relaxation rates for saturated phosphatidylcholines follow a single frequency-dispersive trend within the MHz regime. We show how analytical models can guide the continued development of atomistic and coarse-grained force fields. Our interpretation suggests that lipid diffusion and collective order fluctuations are implicitly governed by the viscoelastic nature of the liquid-crystalline ensemble. Collective bilayer excitations are emergent over mesoscopic length scales that fall between the molecular and bilayer dimensions, and are important for lipid organization and lipid-protein interactions. Future conceptual advances and theoretical reductions will foster understanding of biomembrane structural dynamics through a synergy of NMR measurements and molecular simulations.     

Conformational dynamics of cytochrome c: correlation to hydrogen exchange.

We study the dynamical fluctuations of horse heart cytochrome c by molecular dynamics (MD) simulations in aqueous solution, at four temperatures: 300 K, 360 K, 430 K, and 550 K. Each simulation covers a production time of at least 1.5 nanoseconds (ns). The conformational dynamics of the system is analyzed in terms of collective motions that involve the whole protein, and local motions that involve the formation and breaking of intramolecular hydrogen bonds. The character of the MD trajectories can be described within the framework of rugged energy landscape dynamics. The MD trajectories sample multiple conformational minima, with basins in protein conformational space being sampled for a few hundred picoseconds. The trajectories of the system in configurational space can be described in terms of diffusion of a particle in real space with a waiting time distribution due to partial trapping in shallow minima. As a consequence of the hierarchical nature of the dynamics, the mean square displacement autocorrelation function, <|x(t) - x(0)|2>, exhibits a power law dependence on time, with an exponent of around 0.5 for times shorter than 100 ps, and an exponent of 1.75 for longer times. This power law behavior indicates that the system exhibits suppressed diffusion (sub-diffusion) in sampling of configurational space at time scales shorter than 100 ps, and enhanced (super-diffusion) at longer time scales. The multi-basin feature of the trajectories is present at all temperatures simulated. Structural changes associated with inter-basin displacements correspond to collective motions of the Omega loops and coiled regions and relative motions of the alpha-helices as rigid bodies. Similar motions may be involved in experimentally observed amide hydrogen exchange. However, some groups showing large correlated motions do not expose the amino hydrogens to the solvent. We show that large fluctuations are not necessarily correlated to hydrogen exchange. For example, regions of the proteins forming alpha helices and turns show significant fluctuations, but as rigid bodies, and the hydrogen bonds involved in the formation of these structures do not break in proportion to these fluctuations. Proteins 1999;36:175-191. Published 1999 Wiley-Liss, Inc.     

Effect of membrane tension on the physical properties of DOPC lipid bilayer membrane

Molecular dynamics simulations of a dioleoylphosphocholine (DOPC) lipid bilayer were performed to explore its mechanosensitivity. Variations in the bilayer properties, such as area per lipid, volume, thickness, hydration depth (HD), hydration thickness (HT), lateral diffusion coefficient, and changes in lipid structural order were computed in the membrane tension range 0 to 15dyn/cm. We determined that an increase in membrane tension results in a decrease in the bilayer thickness and HD of ~5% and ~5.7% respectively, whereas area per lipid, volume, and HT/HD increased by 6.8%, 2.4%, and 5% respectively. The changes in lipid conformation and orientation were characterized using orientational (S(2)) and deuterium (S(CD)) order parameters. Upon increase of membrane tension both order parameters indicated an increase in lipid disorder by 10-20%, mostly in the tail end region of the hydrophobic chains. The effect of membrane tension on lipid lateral diffusion in the DOPC bilayer was analyzed on three different time scales corresponding to inertial motion, anomalous diffusion and normal diffusion. The results showed that lateral diffusion of lipid molecules is anomalous in nature due to the non-exponential distribution of waiting times. The anomalous and normal diffusion coefficients increased by 20% and 52% when the membrane tension changed from 0 to 15dyn/cm, respectively. In conclusion, our studies showed that membrane tension causes relatively significant changes in the area per lipid, volume, polarity, membrane thickness, and fluidity of the membrane suggesting multiple mechanisms by which mechanical perturbation of the membrane could trigger mechanosensitive response in cells.     

Solvent effects in the slow dynamics of proteins

The influence of solvent on the slow internal dynamics of proteins is studied by comparing molecular dynamics simulations of solvated and unsolvated lysozyme. The dynamical trajectories are projected onto the protein's normal modes in order to obtain a separate analysis for each of the associated time scales. The results show that solvent effects are important for the slowest motions (below approximately 1 ps(-1)) but negligible for faster motions. The damping effects seen in the latter show that the principal source of friction in protein dynamics is not the solvent, but the protein itself.     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Protein dynamics in a family of laboratory evolved thermophilic enzymes

Molecular dynamics simulations were employed to study how protein solution structure and dynamics are affected by adaptation to high temperature. Simulations were carried out on a para-nitrobenzyl esterase (484 residues) and two thermostable variants that were generated by laboratory evolution. Although these variants display much higher melting temperatures than wild-type (up to 18 degrees C higher) they are both >97% identical in sequence to the wild-type. In simulations at 300 K the thermostable variants remain closer to their crystal structures than wild-type. However, they also display increased fluctuations about their time-averaged structures. Additionally, both variants show a small but significant increase in radius of gyration relative to wild-type. The vibrational density of states was calculated for each of the esterases. While the density of states profiles are similar overall, both thermostable mutants show increased populations of the very lowest frequency modes (<10 cm(-1)), with the more stable mutant showing the larger increase. This indicates that the thermally stable variants experience increased concerted motions relative to wild-type. Taken together, these data suggest that adaptation for high temperature stability has resulted in a restriction of large deviations from the native state and a corresponding increase in smaller scale fluctuations about the native state. These fluctuations contribute to entropy and hence to the stability of the native state. The largest changes in localized dynamics occur in surface loops, while other regions, particularly the active site residues, remain essentially unchanged. Several mutations, most notably L313F and H322Y in variant 8G8, are in the region showing the largest increase in fluctuations, suggesting that these mutations confer more flexibility to the loops. As a validation of our simulations, the fluctuations of Trp102 were examined in detail, and compared with Trp102 phosphorescence lifetimes that were previously measured. Consistent with expectations from the theory of phosphorescence, an inverse correlation between out-of-plane fluctuations on the picosecond time scale and phosphorescence lifetime was observed.     

Computer simulations of the flexibility of a series of synthetic cyclic peptide analogues

The flexibility of a series of cyclic peptides derived from the epitope of a snake toxin is investigated using computer simulation techniques. Molecular dynamics (MD) simulations and vibrational analyses are performed on chemically constrained peptides modeled on the parent protein loop. In the 50 ps MD simulations, large variations in the atomic fluctuations are seen between the peptides, and can be related to the nature of the chemical constraints present in the molecules. Normal mode analyses are performed on energy-minimized configurations derived from the dynamics trajectories. The atomic fluctuations calculated from the normal modes are about 30% of those of the molecular dynamics for the more flexible peptides and 70% for the more constrained molecules. The calculated differences in flexibility between the molecules are much less significant in the harmonic approximation. © 1993 John Wiley & Sons, Inc.     

Molecular dynamics simulations of phospholipid bilayers with cholesterol

To investigate the microscopic interactions between cholesterol and lipids in biological membranes, we have performed a series of molecular dynamics simulations of large membranes with different levels of cholesterol content. The simulations extend to 10 ns, and were performed with hydrated dipalmitoylphosphatidylcholine (DPPC) bilayers. The bilayers contain 1024 lipids of which 0-40% were cholesterol and the rest DPPC. The effects of cholesterol on the structure and mesoscopic dynamics of the bilayer were monitored as a function of cholesterol concentration. The main effects observed are a significant ordering of the DPPC chains (as monitored by NMR type order parameters), a reduced fraction of gauche bonds, a reduced surface area per lipid, less undulations--corresponding to an increased bending modulus for the membrane, smaller area fluctuations, and a reduced lateral diffusion of DPPC-lipids as well as cholesterols.     

Thermal fluctuations of large cylindrical phospholipid vesicles

The time correlation function of the shape fluctuations of large (greater than 10 micron), cylindrical, hydrated, phospholipid-membrane vesicles consisting of one bimolecular layer was measured. The restoring force of the membrane was due to the excess curvature of a membrane element. A value for the curvature elastic modulus, Kc, was obtained from the mean-square amplitude of the normal modes of the fluctuations using the equipartition theorem. An expression for the correlation time was found by solving the dynamics of the membrane's relaxation against the low Reynolds number viscous drag of the surrounding fluid. The amplitudes and correlation times of the fundamental bending mode of the cylindrical vesicles both yield Kc = 1-2 X 10(-12) ergs.     

Molecular dynamics of myoglobin at 298 degrees K. Results from a 300-ps computer simulation.

We have carried out a very long (300 ps) molecular dynamics simulation of the protein myoglobin. This trajectory is approximately three times longer than the longest previous molecular dynamics simulation of a protein, and ten times longer than protein simulations of comparable size (1,423 atoms in our model). Here we report results from this long simulation concerning the average structure, the mean square fluctuations of atoms about the average structure, and the nuclear magnetic resonance order parameters for various groups in myoglobin. The results demonstrate that the average coordinates change very slowly during the simulation. The relative atomic mobilities are well described by the simulation. For both the mean square atomic fluctuations and the order parameters, however, there are significant quantitative differences when values calculated using shorter portions of the trajectory are compared with results obtained for the entire 300-ps simulation. The implications of this result for obtaining converged properties from protein molecular dynamics simulations for comparison with experiment are discussed.     

The Temporal Spectrum of Adult Mosquito Population Fluctuations: Conceptual and Modeling Implications

An improved understanding of mosquito population dynamics under natural environmental forcing requires adequate field observations spanning the full range of temporal scales over which mosquito abundance fluctuates in natural conditions. Here we analyze a 9-year daily time series of uninterrupted observations of adult mosquito abundance for multiple mosquito species in North Carolina to identify characteristic scales of temporal variability, the processes generating them, and the representativeness of observations at different sampling resolutions. We focus in particular on Aedes vexans and Culiseta melanura and, using a combination of spectral analysis and modeling, we find significant population fluctuations with characteristic periodicity between 2 days and several years. Population dynamical modelling suggests that the observed fast fluctuations scales (2 days-weeks) are importantly affected by a varying mosquito activity in response to rapid changes in meteorological conditions, a process neglected in most representations of mosquito population dynamics. We further suggest that the range of time scales over which adult mosquito population variability takes place can be divided into three main parts. At small time scales (indicatively 2 days-1 month) observed population fluctuations are mainly driven by behavioral responses to rapid changes in weather conditions. At intermediate scales (1 to several month) environmentally-forced fluctuations in generation times, mortality rates, and density dependence determine the population characteristic response times. At longer scales (annual to multi-annual) mosquito populations follow seasonal and inter-annual environmental changes. We conclude that observations of adult mosquito populations should be based on a sub-weekly sampling frequency and that predictive models of mosquito abundance must include behavioral dynamics to separate the effects of a varying mosquito activity from actual changes in the abundance of the underlying population.     

Small- and large-scale conformational changes of adenylate kinase: a molecular dynamics study of the subdomain motion and mechanics

Adenylate kinase, an enzyme that catalyzes the phosphoryl transfer between ATP and AMP, can interconvert between the open and catalytically potent (closed) forms even without binding ligands. Several aspects of the enzyme elasticity and internal dynamics are analyzed here by atomistic molecular dynamics simulations covering a total time span of 100 ns. This duration is sufficiently long to reveal a partial conversion of the enzyme that proceeds through jumps between structurally different substates. The intra- and intersubstates contributions to the enzyme's structural fluctuations are analyzed and compared both in magnitude and directionality. It is found that, despite the structural heterogeneity of the visited conformers, the generalized directions accounting for conformational fluctuations within and across the substates are mutually consistent and can be described by a limited set of collective modes. The functional-oriented nature of the consensus modes is suggested by their good overlap with the deformation vector bridging the open and closed crystal structures. The consistency of adenylate kinase's internal dynamics over timescales wide enough to capture intra- and intersubstates fluctuations adds elements in favor of the recent proposal that the free (apo) enzyme possesses an innate ability to sustain the open/close conformational changes.     

Go reconfigure: how fish change shape as they swim and evolve

The bodies of fish change shape over propulsive, behavioral, developmental, and evolutionary time scales, a general phenomenon that we call "reconfiguration". Undulatory, postural, and form-reconfiguration can be distinguished, studied independently, and examined in terms of mechanical interactions and evolutionary importance. Using a combination of live, swimming fishes and digital robotic fish that are autonomous and self-propelled, we examined the functional relation between undulatory and postural reconfiguration in forward swimming, backward swimming, and yaw turning. To probe how postural and form reconfiguration interact, the yaw turning of leopard sharks was examined using morphometric and kinematic analyses. To test how undulatory reconfiguration might evolve, the digital robotic fish were subjected to selection for enhanced performance in a simulated ecology in which each individual had to detect and move towards a food source. In addition to the general issue of reconfiguration, these investigations are united by the fact that the dynamics of undulatory and postural reconfigurations are predicted to be determined, in part, by the structural stiffness of the fish's body. Our method defines undulatory reconfiguration as the combined, point-by-point periodic motion of the body, leaving postural reconfiguration as the combined deviations from undulatory reconfiguration. While undulatory reconfiguration appears to be the sole or primary propulsive driver, postural reconfiguration may contribute to propulsion in hagfish and it is correlated with differences in forward, and backward, swimming in lamprey. Form reconfigures over developmental time in leopard sharks in a manner that is consistent with an allometric scaling theory in which structural stiffness of the body is held constant. However, correlation of a form proxy for structural stiffness of the body suggests that body stiffness may scale in order to limit maximum postural reconfiguration during routine yaw turns. When structural stiffness and undulatory frequency are modeled as determining the tail's undulatory wave speed, both factors evolve under selection for enhanced foraging behavior in the digital fish-like robots. The methods used in making these distinctions between kinds of reconfiguration have broad applicability in fish biology, especially for quantifying complex motor behaviors in the wild and for simulating selection on behavior that leads to directional evolution of functional phenotypes.     

Influence of Structural Symmetry on Protein Dynamics

Structural symmetry in homooligomeric proteins has intrigued many researchers over the past several decades. However, the implication of protein symmetry is still not well understood. In this study, we performed molecular dynamics (MD) simulations of two forms of trp RNA binding attenuation protein (TRAP), the wild-type 11-mer and an engineered 12-mer, having two different levels of circular symmetry. The results of the simulations showed that the inter-subunit fluctuations in the 11-mer TRAP were significantly smaller than the fluctuations in the 12-mer TRAP while the internal fluctuations were larger in the 11-mer than in the 12-mer. These differences in thermal fluctuations were interpreted by normal mode analysis and group theory. For the 12-mer TRAP, the wave nodes of the normal modes existed at the flexible interface between the subunits, while the 11-mer TRAP had its nodes within the subunits. The principal components derived from the MD simulations showed similar mode structures. These results demonstrated that the structural symmetry was an important determinant of protein dynamics in circularly symmetric homooligomeric proteins.     

A Forced Damped Oscillation Framework for Undulatory Swimming Provides New Insights into How Propulsion Arises in Active and Passive Swimming

A fundamental issue in locomotion is to understand how muscle forcing produces apparently complex deformation kinematics leading to movement of animals like undulatory swimmers. The question of whether complicated muscle forcing is required to create the observed deformation kinematics is central to the understanding of how animals control movement. In this work, a forced damped oscillation framework is applied to a chain-link model for undulatory swimming to understand how forcing leads to deformation and movement. A unified understanding of swimming, caused by muscle contractions (“active” swimming) or by forces imparted by the surrounding fluid (“passive” swimming), is obtained. We show that the forcing triggers the first few deformation modes of the body, which in turn cause the translational motion. We show that relatively simple forcing patterns can trigger seemingly complex deformation kinematics that lead to movement. For given muscle activation, the forcing frequency relative to the natural frequency of the damped oscillator is important for the emergent deformation characteristics of the body. The proposed approach also leads to a qualitative understanding of optimal deformation kinematics for fast swimming. These results, based on a chain-link model of swimming, are confirmed by fully resolved computational fluid dynamics (CFD) simulations. Prior results from the literature on the optimal value of stiffness for maximum speed are explained.     

Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor

The dynamics of alpha-amylase inhibitors has been investigated using molecular dynamics (MD) simulations and two analytical approaches, the Gaussian network model (GNM) and anisotropic network model (ANM). MD simulations use a full atomic approach with empirical force fields, while the analytical approaches are based on a coarse-grained single-site-per-residue model with a single-parameter harmonic potential between sufficiently close (r 相似文献