Myosin dynamics on the millisecond time scale

Neural/Bayes network predictor for inheritable cardiac disease pathogenicity and phenotype

The cardiac muscle sarcomere contains multiple proteins contributing to contraction energy transduction and its regulation during a heartbeat. Inheritable heart disease mutants affect most of them but none more frequently than the ventricular myosin motor and cardiac myosin binding protein c (mybpc3). These co-localizing proteins have mybpc3 playing a regulatory role to the energy transducing motor. Residue substitution and functional domain assignment of each mutation in the protein sequence decides, under the direction of a sensible disease model, phenotype and pathogenicity. The unknown model mechanism is decided here using a method combing neural and Bayes networks. Missense single nucleotide polymorphisms (SNPs) are clues for the disease mechanism summarized in an extensive database collecting mutant sequence location and residue substitution as independent variables that imply the dependent disease phenotype and pathogenicity characteristics in 4 dimensional data points (4ddps). The SNP database contains entries with the majority having one or both dependent data entries unfulfilled. A neural network relating causes (mutant residue location and substitution) and effects (phenotype and pathogenicity) is trained, validated, and optimized using fulfilled 4ddps. It then predicts unfulfilled 4ddps providing the implicit disease model. A discrete Bayes network interprets fulfilled and predicted 4ddps with conditional probabilities for phenotype and pathogenicity given mutation location and residue substitution thus relating the neural network implicit model to explicit features of the motor and mybpc3 sequence and structural domains. Neural/Bayes network forecasting automates disease mechanism modeling by leveraging the world wide human missense SNP database that is in place and expanding.

Early events in helix unfolding under external forces: a milestoning analysis

Initial events of helix breakage as a function of load are considered using molecular dynamics simulations and milestoning analysis. A helix length of ∼100 amino acids is considered as a model for typical helices found in molecular machines and as a model that minimizes end effects for early events of unfolding. Transitions of individual amino acids (averaged over the helix's interior residues) are examined and its surrounding hydrogen bonds are considered. Dense kinetic networks are constructed that, with milestoning analysis, provide the overall kinetics of early breakage events. Network analysis and selection of MaxFlux pathways illustrate that load impacts unfolding mechanisms in addition to time scales. At relatively high (100 pN) load levels, the principal intermediate is the 3(10)-helix, while at relatively low (10 pN) levels the π-helix is significantly populated, albeit not as an unfolding intermediate. Coarse variables are examined at different levels of resolution; the rate of unfolding illustrates remarkable stability under changes in the coarsening. Consistent prediction of about ∼5 ns for the time of a single amino-acid unfolding event are obtained. Hydrogen bonds are much faster coarse variables (by about 2 orders of magnitude) compared to backbone torsional transition, which gates unfolding and thereby provides the appropriate coarse variable for the initiation of unfolding. Results provide an atomic description of "catch-bond" behavior, based on alternative pathways, in which unfolding of a simple protein structural element occurs over longer timescales for intermediate (10 pN) loads than for zero (0 pN) or large (100 pN) loads.

Self-assembly of coiled coils in synthetic biology: inspiration and progress

Biological self-assembly is very complex and results in highly functional materials. In effect, it takes a bottom-up approach using biomolecular building blocks of precisely defined shape, size, hydrophobicity, and spatial distribution of functionality. Inspired by, and drawing lessons from self-assembly processes in nature, scientists are learning how to control the balance of many small forces to increase the complexity and functionality of self-assembled nanomaterials. The coiled-coil motif, a multipurpose building block commonly found in nature, has great potential in synthetic biology. In this review we examine the roles that the coiled-coil peptide motif plays in self-assembly in nature, and then summarize the advances that this has inspired in the creation of functional units, assemblies, and systems.

Myosin individualized: single nucleotide polymorphisms in energy transduction

BACKGROUND

Myosin performs ATP free energy transduction into mechanical work in the motor domain of the myosin heavy chain (MHC). Energy transduction is the definitive systemic feature of the myosin motor performed by coordinating in a time ordered sequence: ATP hydrolysis at the active site, actin affinity modulation at the actin binding site, and the lever-arm rotation of the power stroke. These functions are carried out by several conserved sub-domains within the motor domain. Single nucleotide polymorphisms (SNPs) affect the MHC sequence of many isoforms expressed in striated muscle, smooth muscle, and non-muscle tissue. The purpose of this work is to provide a rationale for using SNPs as a functional genomics tool to investigate structurefunction relationships in myosin. In particular, to discover SNP distribution over the conserved sub-domains and surmise what it implies about sub-domain stability and criticality in the energy transduction mechanism.

RESULTS

An automated routine identifying human nonsynonymous SNP amino acid missense substitutions for any MHC gene mined the NCBI SNP data base. The routine tested 22 MHC genes coding muscle and non-muscle isoforms and identified 89 missense mutation positions in the motor domain with 10 already implicated in heart disease and another 8 lacking sequence homology with a skeletal MHC isoform for which a crystallographic model is available. The remaining 71 SNP substitutions were found to be distributed over MHC with 22 falling outside identified functional sub-domains and 49 in or very near to myosin sub-domains assigned specific crucial functions in energy transduction. The latter includes the active site, the actin binding site, the rigid lever-arm, and regions facilitating their communication. Most MHC isoforms contained SNPs somewhere in the motor domain.

CONCLUSIONS

Several functional-crucial sub-domains are infiltrated by a large number of SNP substitution sites suggesting these domains are engineered by evolution to be too-robust to be disturbed by otherwise intrusive sequence changes. Two functional sub-domains are SNP-free or relatively SNP-deficient but contain many disease implicated mutants. These sub-domains are apparently highly sensitive to any missense substitution suggesting they have failed to evolve a robust sequence paradigm for performing their function.

The myosin C-loop is an allosteric actin contact sensor in actomyosin

Actin and myosin form the molecular motor in muscle. Myosin is the enzyme performing ATP hydrolysis under the allosteric control of actin such that actin binding initiates product release and force generation in the myosin power stroke. Non-equilibrium Monte Carlo molecular dynamics simulation of the power stroke suggested that a structured surface loop on myosin, the C-loop, is the actin contact sensor initiating actin activation of the myosin ATPase. Previous experimental work demonstrated C-loop binds actin and established the forward and reverse allosteric link between the C-loop and the myosin active site. Here, smooth muscle heavy meromyosin C-loop chimeras were constructed with skeletal (sCl) and cardiac (cCl) myosin C-loops substituted for the native sequence. In both cases, actin-activated ATPase inhibition is indicated mainly by the lower V(max). In vitro motility was also inhibited in the chimeras. Motility data were collected as a function of myosin surface density, with unregulated actin, and with skeletal and cardiac isoforms of tropomyosin-bound actin for the wild type, cCl, and sCl. Slow and fast subpopulations of myosin velocities in the wild-type species were discovered and represent geometrically unfavorable and favorable actomyosin interactions, respectively. Unfavorable interactions are detected at all surface densities tested. Favorable interactions are more probable at higher myosin surface densities. Cardiac tropomyosin-bound actin promotes the favorable actomyosin interactions by lowering the inhibiting geometrical constraint barriers with a structural effect on actin. Neither higher surface density nor cardiac tropomyosin-bound actin can accelerate motility velocity in cCl or sCl, suggesting the element initiating maximal myosin activation by actin resides in the C-loop.

Energetics of subdomain movements and fluorescence probe solvation environment change in ATP-bound myosin

Energetics of conformational changes experienced by an ATP-bound myosin head detached from actin was studied by all-atom explicit water umbrella sampling simulations. The statistics of coupling between large scale domain movements and smaller scale structural features were examined, including the closing of the ATP binding pocket, and a number of key hydrogen bond formations shown to play roles in structural and biochemical studies. The statistics for the ATP binding pocket open/close transition show an evolution of the relative stability from the open state in the early stages of the recovery stroke to the stable closed state after the stroke. The change in solvation environment of the fluorescence probe Trp507 (scallop numbering; 501 in Dictyostelium discoideum) indicates that the probe faithfully reflects the closing of the binding pocket as previously shown in experimental studies, while being directly coupled to roughly the early half of the overall large scale conformational change of the converter domain rotation. The free energy change of this solvation environment change, in particular, is -1.3 kcal/mol, in close agreement with experimental estimates. In addition, our results provide direct molecular level data allowing for interpretations of the fluorescence experiments of myosin conformational change in terms of the de-solvation of Trp side chain.