Molecular origin of Gerstmann-Sträussler-Scheinker syndrome: insight from computer simulation of an amyloidogenic prion peptide

Prion proteins become pathogenic through misfolding. Here, we characterize the folding of a peptide consisting of residues 109-122 of the Syrian hamster prion protein (the H1 peptide) and of a more amyloidogenic A117V point mutant that leads in humans to an inheritable form of the Gerstmann-Sträussler-Scheinker syndrome. Atomistic molecular dynamics simulations are performed for 2.5 μs. Both peptides lose their α-helical starting conformations and assume a β-hairpin that is structurally similar in both systems. In each simulation several unfolding/refolding events occur, leading to convergence of the thermodynamics of the conformational states to within 1 kJ/mol. The similar stability of the β-hairpin relative to the unfolded state is observed in the two peptides. However, substantial differences are found between the two unfolded states. A local minimum is found within the free energy unfolded basin of the A117V mutant populated by misfolded collapsed conformations of comparable stability to the β-hairpin state, consistent with increased amyloidogenicity. This population, in which V117 stabilizes a hydrophobic core, is absent in the wild-type peptide. These results are supported by simulations of oligomers showing a slightly higher stability of the associated structures and a lower barrier to association for the mutated peptide. Hence, a single point mutation carrying only two additional methyl groups is here shown to be responsible for rather dramatic differences of structuring within the unfolded (misfolded) state.

Optimal use of data in parallel tempering simulations for the construction of discrete-state Markov models of biomolecular dynamics

Parallel tempering (PT) molecular dynamics simulations have been extensively investigated as a means of efficient sampling of the configurations of biomolecular systems. Recent work has demonstrated how the short physical trajectories generated in PT simulations of biomolecules can be used to construct the Markov models describing biomolecular dynamics at each simulated temperature. While this approach describes the temperature-dependent kinetics, it does not make optimal use of all available PT data, instead estimating the rates at a given temperature using only data from that temperature. This can be problematic, as some relevant transitions or states may not be sufficiently sampled at the temperature of interest, but might be readily sampled at nearby temperatures. Further, the comparison of temperature-dependent properties can suffer from the false assumption that data collected from different temperatures are uncorrelated. We propose here a strategy in which, by a simple modification of the PT protocol, the harvested trajectories can be reweighted, permitting data from all temperatures to contribute to the estimated kinetic model. The method reduces the statistical uncertainty in the kinetic model relative to the single temperature approach and provides estimates of transition probabilities even for transitions not observed at the temperature of interest. Further, the method allows the kinetics to be estimated at temperatures other than those at which simulations were run. We illustrate this method by applying it to the generation of a Markov model of the conformational dynamics of the solvated terminally blocked alanine peptide.

In silico partitioning and transmembrane insertion of hydrophobic peptides under equilibrium conditions

Nascent transmembrane (TM) polypeptide segments are recognized and inserted into the lipid bilayer by the cellular translocon machinery. The recognition rules, described by a biological hydrophobicity scale, correlate strongly with physical hydrophobicity scales that describe the free energy of insertion of TM helices from water. However, the exact relationship between the physical and biological scales is unknown, because solubility problems limit our ability to measure experimentally the direct partitioning of hydrophobic peptides across lipid membranes. Here we use microsecond molecular dynamics (MD) simulations in which monomeric polyleucine segments of different lengths are allowed to partition spontaneously into and out of lipid bilayers. This approach directly reveals all states populated at equilibrium. For the hydrophobic peptides studied here, only surface-bound and transmembrane-inserted helices are found. The free energy of insertion is directly obtained from the relative occupancy of these states. A water-soluble state was not observed, consistent with the general insolubility of hydrophobic peptides. The approach further allows determination of the partitioning pathways and kinetics. Surprisingly, the transfer free energy appears to be independent of temperature, which implies that surface-to-bilayer peptide insertion is a zero-entropy process. We find that the partitioning free energy of the polyleucine segments correlates strongly with values from translocon experiments but reveals a systematic shift favoring shorter peptides, suggesting that translocon-to-bilayer partitioning is not equivalent but related to spontaneous surface-to-bilayer partitioning.

Three classes of motion in the dynamic neutron-scattering susceptibility of a globular protein

A simplified description of the 295 K dynamics of a globular protein over a wide frequency range (1-1000 GHz) is obtained by combining neutron scattering of lysozyme with molecular dynamics simulation. The molecular dynamics simulation agrees quantitatively with experiment for both the protein and the hydration water and shows that, whereas the hydration water molecules subdiffuse, the protein atoms undergo confined motion decomposable into three distinct classes: localized diffusion, methyl group rotations, and jumps. Each of the three classes gives rise to a characteristic neutron susceptibility signal.

Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum

Clostridium thermocellum is a thermophilic, obligately anaerobic, gram-positive bacterium that is a candidate microorganism for converting cellulosic biomass into ethanol through consolidated bioprocessing. Ethanol intolerance is an important metric in terms of process economics, and tolerance has often been described as a complex and likely multigenic trait for which complex gene interactions come into play. Here, we resequence the genome of an ethanol-tolerant mutant, show that the tolerant phenotype is primarily due to a mutated bifunctional acetaldehyde-CoA/alcohol dehydrogenase gene (adhE), hypothesize based on structural analysis that cofactor specificity may be affected, and confirm this hypothesis using enzyme assays. Biochemical assays confirm a complete loss of NADH-dependent activity with concomitant acquisition of NADPH-dependent activity, which likely affects electron flow in the mutant. The simplicity of the genetic basis for the ethanol-tolerant phenotype observed here informs rational engineering of mutant microbial strains for cellulosic ethanol production.

Engaging fathers in behavioral parent training: listening to fathers' voices

This is a qualitative study documenting the experiences of fathers who participated in the Helping Our Toddlers, Developing Our Children's Skills (HOT DOCS) behavioral parent training (BPT) series and later agreed to participate in a focus group. Focus groups methodology was used to capture the voices and perspectives of fathers regarding the benefits and barriers to their participation in BPT. The focus group interviews were conducted in both English and in Spanish, with three cohorts of male caregivers who were participants in HOT DOCS from 2006 to 2008. An analysis of their responses coded from transcripts of the focus groups identified five major themes, which are discussed as (a) motivational factors for joining BPT, (b) fathers' experiences with BPT, (c) barriers regarding fathers' participation, (d) changes in parenting as a result of BPT, and (e) perceived changes in children's behavior as a result of BPT. In addition, recommendations for improvement of BPT are presented. This research may be helpful in improving efforts to engage male caregivers in BPT and thereby reduce children's challenging behavior problems and improve program outcomes.

Small angle neutron scattering reveals pH-dependent conformational changes in Trichoderma reesei cellobiohydrolase I: implications for enzymatic activity

Cellobiohydrolase I (Cel7A) of the fungus Trichoderma reesei (now classified as an anamorph of Hypocrea jecorina) hydrolyzes crystalline cellulose to soluble sugars, making it of key interest for producing fermentable sugars from biomass for biofuel production. The activity of the enzyme is pH-dependent, with its highest activity occurring at pH 4-5. To probe the response of the solution structure of Cel7A to changes in pH, we measured small angle neutron scattering of it in a series of solutions having pH values of 7.0, 6.0, 5.3, and 4.2. As the pH decreases from 7.0 to 5.3, the enzyme structure remains well defined, possessing a spatial differentiation between the cellulose binding domain and the catalytic core that only changes subtly. At pH 4.2, the solution conformation of the enzyme changes to a structure that is intermediate between a properly folded enzyme and a denatured, unfolded state, yet the secondary structure of the enzyme is essentially unaltered. The results indicate that at the pH of optimal activity, the catalytic core of the enzyme adopts a structure in which the compact packing typical of a fully folded polypeptide chain is disrupted and suggest that the increased range of structures afforded by this disordered state plays an important role in the increased activity of Cel7A through conformational selection.

A Solvent-Free Coarse Grain Model for Crystalline and Amorphous Cellulose Fibrils

Understanding biomass structure and dynamics on a range of time and length scales is important for the development of cellulosic biofuels. Here, to enable length and time scale extension, we develop a coarse grain (CG) model for molecular dynamics (MD) simulations of cellulose. For this purpose, we use distribution functions from fully atomistic MD simulations as target observables. A single bead per monomer level coarse graining is found to be sufficient to successfully reproduce structural features of crystalline cellulose. Without the use of constraints the CG crystalline fibril is found to remain stable over the maximum simulation length explored in this study (>1 μs). We also extend the CG representation to model fully amorphous cellulose fibrils. This is done by using an atomistic MD simulation of fully solvated individual cellulose chains as a target for developing the corresponding fully amorphous CG force field. Fibril structures with different degrees of crystallinity are obtained using force fields derived using a parameter coupling the crystalline and amorphous potentials. The method provides an accurate and constraint-free approach to derive CG models for cellulose with a wide range of crystallinity, suitable for incorporation into large-scale models of lignocellulosic biomass.

Self-similar multiscale structure of lignin revealed by neutron scattering and molecular dynamics simulation

Lignin, a major polymeric component of plant cell walls, forms aggregates in vivo and poses a barrier to cellulosic ethanol production. Here, neutron scattering experiments and molecular dynamics simulations reveal that lignin aggregates are characterized by a surface fractal dimension that is invariant under change of scale from ~1-1000 Å. The simulations also reveal extensive water penetration of the aggregates and heterogeneous chain dynamics corresponding to a rigid core with a fluid surface.

Vibrational softening of a protein on ligand binding

Neutron scattering experiments have demonstrated that binding of the cancer drug methotrexate softens the low-frequency vibrations of its target protein, dihydrofolate reductase (DHFR). Here, this softening is fully reproduced using atomic detail normal-mode analysis. Decomposition of the vibrational density of states demonstrates that the largest contributions arise from structural elements of DHFR critical to stability and function. Mode-projection analysis reveals an increase of the breathing-like character of the affected vibrational modes consistent with the experimentally observed increased adiabatic compressibility of the protein on complexation.

Structural modeling and molecular dynamics simulation of the actin filament

Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility. Here, we present a model of the actin filament (F-actin) that not only incorporates the global structure of the recently published model by Oda et al. but also conserves internal stereochemistry. A comparison is made using molecular dynamics simulation of the model with other recent F-actin models. A number of structural determents such as the protomer propeller angle, the number of hydrogen bonds, and the structural variation among the protomers are analyzed. The MD comparison is found to reflect the evolution in quality of actin models over the last 6 years. In addition, simulations of the model are carried out in states with both ADP or ATP bound and local hydrogen-bonding differences characterized.

Task-parallel message passing interface implementation of Autodock4 for docking of very large databases of compounds using high-performance super-computers

A message passing interface (MPI)-based implementation (Autodock4.lga.MPI) of the grid-based docking program Autodock4 has been developed to allow simultaneous and independent docking of multiple compounds on up to thousands of central processing units (CPUs) using the Lamarkian genetic algorithm. The MPI version reads a single binary file containing precalculated grids that represent the protein-ligand interactions, i.e., van der Waals, electrostatic, and desolvation potentials, and needs only two input parameter files for the entire docking run. In comparison, the serial version of Autodock4 reads ASCII grid files and requires one parameter file per compound. The modifications performed result in significantly reduced input/output activity compared with the serial version. Autodock4.lga.MPI scales up to 8192 CPUs with a maximal overhead of 16.3%, of which two thirds is due to input/output operations and one third originates from MPI operations. The optimal docking strategy, which minimizes docking CPU time without lowering the quality of the database enrichments, comprises the docking of ligands preordered from the most to the least flexible and the assignment of the number of energy evaluations as a function of the number of rotatable bounds. In 24 h, on 8192 high-performance computing CPUs, the present MPI version would allow docking to a rigid protein of about 300K small flexible compounds or 11 million rigid compounds.

Dynamical fingerprints for probing individual relaxation processes in biomolecular dynamics with simulations and kinetic experiments

There is a gap between kinetic experiment and simulation in their views of the dynamics of complex biomolecular systems. Whereas experiments typically reveal only a few readily discernible exponential relaxations, simulations often indicate complex multistate behavior. Here, a theoretical framework is presented that reconciles these two approaches. The central concept is "dynamical fingerprints" which contain peaks at the time scales of the dynamical processes involved with amplitudes determined by the experimental observable. Fingerprints can be generated from both experimental and simulation data, and their comparison by matching peaks permits assignment of structural changes present in the simulation to experimentally observed relaxation processes. The approach is applied here to a test case interpreting single molecule fluorescence correlation spectroscopy experiments on a set of fluorescent peptides with molecular dynamics simulations. The peptides exhibit complex kinetics shown to be consistent with the apparent simplicity of the experimental data. Moreover, the fingerprint approach can be used to design new experiments with site-specific labels that optimally probe specific dynamical processes in the molecule under investigation.

Dirac point degenerate with massive bands at a topological quantum critical point

The quasilinear bands in the topologically trivial skutterudite insulator CoSb(3) are studied under adiabatic, symmetry-conserving displacement of the Sb sublattice. In this cubic, time-reversal and inversion symmetric system, a transition from trivial insulator to topological point Fermi surface system occurs through a critical point in which massless (Dirac) bands appear, and moreover are degenerate with massive bands. Spin-orbit coupling, while small due to the type of band character, coupled with tetragonal strain opens the gap required to give the topological insulator. The mineral skutterudite (CoSb(3)) is very near the critical point in its natural state.

Configurational subdiffusion of peptides: a network study

Molecular dynamics (MD) simulation of linear peptides reveals configurational subdiffusion at equilibrium extending from 10⁻¹² to 10⁻⁸ s. Rouse chain and continuous-time random walk models of the subdiffusion are critically discussed. Network approaches to analyzing MD simulations are shown to reproduce the time dependence of the subdiffusive mean squared displacement, which is found to arise from the fractal-like geometry of the accessible volume in the configuration space. Convergence properties of the simulation pertaining to the subdiffusive dynamics are characterized and the effect on the subdiffusive properties of representing the solvent explicitly or implicitly is compared. Non-Markovianity and other factors limiting the range of applicability of the network models are examined.

Catalytic mechanism of cellulose degradation by a cellobiohydrolase, CelS

The hydrolysis of cellulose is the bottleneck in cellulosic ethanol production. The cellobiohydrolase CelS from Clostridium thermocellum catalyzes the hydrolysis of cello-oligosaccharides via inversion of the anomeric carbon. Here, to examine key features of the CelS-catalyzed reaction, QM/MM (SCCDFTB/MM) simulations are performed. The calculated free energy profile for the reaction possesses a 19 kcal/mol barrier. The results confirm the role of active site residue Glu87 as the general acid catalyst in the cleavage reaction and show that Asp255 may act as the general base. A feasible position in the reactant state of the water molecule responsible for nucleophilic attack is identified. Sugar ring distortion as the reaction progresses is quantified. The results provide a computational approach that may complement the experimental design of more efficient enzymes for biofuel production.

Probing the mechanism of cellulosome attachment to the Clostridium thermocellum cell surface: computer simulation of the Type II cohesin-dockerin complex and its variants

The recalcitrance of lignocellulosic biomass to hydrolysis is the bottleneck in cellulosic ethanol production. Efficient degradation of biomass by the anaerobic bacterium Clostridium thermocellum is carried out by the multicomponent cellulosome complex. The bacterial cell-surface attachment of the cellulosome is mediated by high-affinity protein-protein interactions between the Type II cohesin domain borne by the cell envelope protein and the Type II dockerin domain, together with neighboring X-module present at the C-terminus of the scaffolding protein (Type II coh-Xdoc). Here, the Type II coh-Xdoc interaction is probed using molecular dynamics simulations, free-energy calculations and essential dynamics analyses on both the wild type and various mutants of the C. thermocellum Type II coh-Xdoc in aqueous solution. The simulations identify the hot spots, i.e. the amino acid residues that may lead to a dramatic decrease in binding affinity upon mutation and also probe the effects of mutations on the mode of binding. The results suggest that bulky and hydrophobic residues at the protein interface, which make specific contacts with their counterparts, may play essential roles in retaining a rigid cohesin-dockerin interface. Moreover, dynamical cross-correlation analysis indicates that the X-module has a dramatic effect on the cohesin-dockerin interaction and is required for the dynamical integrity of the interface.

Intersegmental vessel formation in zebrafish: requirement for VEGF but not BMP signalling revealed by selective and non-selective BMP antagonists

BACKGROUND AND PURPOSE

Bone morphogenetic proteins (BMPs) were first identified through their role in inducing bone and cartilage formation, but many other important functions have since been ascribed to BMPs, including dorsoventral patterning, angiogenesis and tissue homeostasis. Using dorsomorphin and LDN193189, selective small molecule inhibitors of BMP signalling, we investigated the role of BMP signalling in early vascular patterning in zebrafish.

EXPERIMENTAL APPROACH

The effects of dorsomorphin and LDN193189 on vascular endothelial growth factor-a (VEGF) and BMP signalling in developing zebrafish and in human pulmonary artery endothelial cells were determined using confocal microscopy, Western blotting and quantitative PCR.

KEY RESULTS

We showed that dorsomorphin, similar to the VEGF inhibitor SU5416, strongly inhibits intersegmental vessel formation in zebrafish and that this is due to inhibition of VEGF activation of VEGF receptor 2 (VEGFR2), leading to reduced VEGF-induced phospho-ERK (extracellular regulated kinase) 1/2 and VEGF target gene transcription. These effects occurred at concentrations of dorsomorphin that block BMP signalling. We also showed that LDN193189, an analogue of dorsomorphin, more potently blocks BMP signalling but has no effect on VEGF signalling in zebrafish and does not disrupt early vascular patterning.

CONCLUSIONS AND IMPLICATIONS

Dorsomorphin inhibits both BMP and VEGF signalling, whereas LDN193189 is a more selective BMP antagonist. Results obtained in cardiovascular studies using dorsomorphin need to be interpreted with caution, and use of LDN193189 would be preferable due to its selectivity. Our data also suggest that BMP signalling is dispensable for early patterning of intersegmental vessels in zebrafish.

Autotrophic picoplankton in the tropical ocean

In phytoplankton of the eastern tropical Pacific Ocean from 25 to 90 percent of the biomass (measured as chlorophyll a) and 20 to 80 percent of the inorganic carbon fixation were attributable to particles that could pass a screen with a 1-micrometer pore diameter. Evidence is presented that these are indeed autotrophic cells and not cell fragments.

Common folding mechanism of a beta-hairpin peptide via non-native turn formation revealed by unbiased molecular dynamics simulations

The folding of a 15-residue beta-hairpin peptide (Peptide 1) is characterized using multiple unbiased, atomistic molecular dynamics (MD) simulations. Fifteen independent MD trajectories, each 2.5 micros-long for a total of 37.5 micros, are performed of the peptide in explicit solvent, at room temperature, and without the use of enhanced sampling techniques. The computed folding time of 1-1.5 micros obtained from the simulations is in good agreement with experiment [Xu, Y.; et al. J. Am. Chem. Soc. 2003, 125, 15388-15394]. A common folding mechanism is observed, in which the turn is always found to be the major determinant in initiating the folding process, followed by cooperative formation of the interstrand hydrogen bonds and the side-chain packing. Furthermore, direct transition to the folded state from fully unstructured conformations does not take place. Instead, the peptide is always observed to form partially structured conformations involving a non-native (ESYI) turn from which the native (NPDG) turn forms, triggering the folding to the beta-hairpin.

Structure and conformational dynamics of the metalloregulator MerR upon binding of Hg(II)

The bacterial metalloregulator MerR is the index case of an eponymous family of regulatory proteins, which controls the transcription of a set of genes (the mer operon) conferring mercury resistance in many bacteria. Homodimeric MerR represses transcription in the absence of mercury and activates transcription upon Hg(II) binding. Here, the average structures of the apo and Hg(II)-bound forms of MerR in aqueous solution are examined using small-angle X-ray scattering, indicating an extended conformation of the metal-bound protein and revealing the existence of a novel compact conformation in the absence of Hg(II). Molecular dynamics (MD) simulations are performed to characterize the conformational dynamics of the Hg(II)-bound form. In both small-angle X-ray scattering and MD, the average torsional angle between DNA-binding domains is approximately 65 degrees. Furthermore, in MD, interdomain motions on a timescale of approximately 10 ns involving large-amplitude (approximately 20 A) domain opening-and-closing, coupled to approximately 40 degrees variations of interdomain torsional angle, are revealed. This correlated domain motion may propagate allosteric changes from the metal-binding site to the DNA-binding site while maintaining DNA contacts required to initiate DNA underwinding.

Nucleotide-dependence of G-actin conformation from multiple molecular dynamics simulations and observation of a putatively polymerization-competent superclosed state

The assembly of monomeric G-actin into filamentous F-actin is nucleotide dependent: ATP-G-actin is favored for filament growth at the "barbed end" of F-actin, whereas ADP-G-actin tends to dissociate from the "pointed end." Structural differences between ATP- and ADP-G-actin are examined here using multiple molecular dynamics simulations. The "open" and "closed" conformational states of G-actin in aqueous solution are characterized, with either ATP or ADP in the nucleotide binding pocket. With both ATP and ADP bound, the open state closes in the absence of actin-bound profilin. The position of the nucleotide in the protein is found to be correlated with the degree of opening of the active site cleft. Further, the simulations reveal the existence of a structurally well-defined, compact, "superclosed" state of ATP-G-actin, as yet unseen crystallographically and absent in the ADP-G-actin simulations. The superclosed state resembles structurally the actin monomer in filament models derived from fiber diffraction and is putatively the polymerization competent conformation of ATP-G-actin.

Hydrogen-bond driven loop-closure kinetics in unfolded polypeptide chains

Characterization of the length dependence of end-to-end loop-closure kinetics in unfolded polypeptide chains provides an understanding of early steps in protein folding. Here, loop-closure in poly-glycine-serine peptides is investigated by combining single-molecule fluorescence spectroscopy with molecular dynamics simulation. For chains containing more than 10 peptide bonds loop-closing rate constants on the 20–100 nanosecond time range exhibit a power-law length dependence. However, this scaling breaks down for shorter peptides, which exhibit slower kinetics arising from a perturbation induced by the dye reporter system used in the experimental setup. The loop-closure kinetics in the longer peptides is found to be determined by the formation of intra-peptide hydrogen bonds and transient β-sheet structure, that accelerate the search for contacts among residues distant in sequence relative to the case of a polypeptide chain in which hydrogen bonds cannot form. Hydrogen-bond-driven polypeptide-chain collapse in unfolded peptides under physiological conditions found here is not only consistent with hierarchical models of protein folding, that highlights the importance of secondary structure formation early in the folding process, but is also shown to speed up the search for productive folding events.

In studies of protein folding evidence exists for early compaction in the unfolded state, although it is unclear whether these compact conformations contain specific secondary structures (through hydrophilic interactions) or whether compaction is a non-specific hydrophobic-driven effect. Here we combine single-molecule fluorescence spectroscopy and molecular dynamics simulation to demonstrate peptide hydrogen-bond-driven polypeptide-chain collapse involving secondary structure formation as the key process in the early stage of folding. Partial structuring in unfolded polypeptide chains is shown to lead to faster contact formation kinetics than would be expected if the unfolded state were populated by featureless random-coils.

Downstream of FGF during mesoderm formation in Xenopus: the roles of Elk-1 and Egr-1

Signalling by members of the FGF family is required for induction and maintenance of the mesoderm during amphibian development. One of the downstream effectors of FGF is the SRF-interacting Ets family member Elk-1, which, after phosphorylation by MAP kinase, activates the expression of immediate-early genes. Here, we show that Xenopus Elk-1 is phosphorylated in response to FGF signalling in a dynamic pattern throughout the embryo. Loss of XElk-1 function causes reduced expression of Xbra at neurula stages, followed by a failure to form notochord and muscle and then the partial loss of trunk structures. One of the genes regulated by XElk-1 is XEgr-1, which encodes a zinc finger transcription factor: we show that phosphorylated XElk-1 forms a complex with XSRF that binds to the XEgr-1 promoter. Superficially, Xenopus tropicalis embryos with reduced levels of XEgr-1 resemble those lacking XElk-1, but to our surprise, levels of Xbra are elevated at late gastrula stages in such embryos, and over-expression of XEgr-1 causes the down-regulation of Xbra both in whole embryos and in animal pole regions treated with activin or FGF. In contrast, the myogenic regulatory factor XMyoD is activated by XEgr-1 in a direct manner. We discuss these counterintuitive results in terms of the genetic regulatory network to which XEgr-1 contributes.

REACH coarse-grained normal mode analysis of protein dimer interaction dynamics

The REACH (realistic extension algorithm via covariance Hessian) coarse-grained biomolecular simulation method is a self-consistent multiscale approach directly mapping atomistic molecular dynamics simulation results onto a residue-scale model. Here, REACH is applied to calculate the dynamics of protein-protein interactions. The intra- and intermolecular fluctuations and the intermolecular vibrational densities of states derived from atomistic molecular dynamics are well reproduced by the REACH normal modes. The phonon dispersion relations derived from the REACH lattice dynamics model of crystalline ribonuclease A are also in satisfactory agreement with the corresponding all-atom results. The REACH model demonstrates that increasing dimer interaction strength decreases the translational and rotational intermolecular vibrational amplitudes, while their vibrational frequencies are relatively unaffected. A comparative study of functionally interacting biological dimers with crystal dimers, which are formed artificially via crystallization, reveals a relation between their static structures and the interprotein dynamics: i.e., the consequence of the extensive interfaces of biological dimers is reduction of the intermonomer translational and rotational amplitudes, but not the frequencies.

Response of small-scale, methyl rotors to protein-ligand association: a simulation analysis of calmodulin-peptide binding

Changes in the free energy barrier (DeltaE), entropy, and motional parameters associated with the rotation of methyl groups in a protein (calmodulin (CaM)) on binding a ligand (the calmodulin-binding domain of smooth-muscle myosin (smMLCKp)) are investigated using molecular dynamics simulation. In both the bound and uncomplexed forms of CaM, the methyl rotational free energy barriers follow skewed-Gaussian distributions that are not altered significantly upon ligand binding. However, site-specific perturbations are found. Around 11% of the methyl groups in CaM exhibit changes in DeltaE greater than 0.7 kcal/mol on binding. The rotational entropies of the methyl groups exhibit a nonlinear dependence on DeltaE. The relations are examined between motional parameters (the methyl rotational NMR order parameter and the relaxation time) and DeltaE. Low-barrier methyl group rotational order parameters deviate from ideal tetrahedrality by up to approximately 20%. There is a correlation between rotational barrier changes and proximity to the protein-peptide binding interface. Methyl groups that exhibit large changes in DeltaE are found to report on elements in the protein undergoing structural change on binding.