Binding Mechanism between Platelet Glycoprotein and Cyclic Peptide Elucidated by McMD-Based Dynamic Docking.

The cyclic peptide OS1 (amino acid sequence: CTERMALHNLC), which has a disulfide bond between both termini cysteine residues, inhibits complex formation between the platelet glycoprotein Ibα (GPIbα) and the von Willebrand factor (vWF) by forming a complex with GPIbα. To study the binding mechanism between GPIbα and OS1 and, therefore, the inhibition mechanism of the protein-protein GPIbα-vWF complex, we have applied our multicanonical molecular dynamics (McMD)-based dynamic docking protocol starting from the unbound state of the peptide. Our simulations have reproduced the experimental complex structure, although the top-ranking structure was an intermediary one, where the peptide was bound in the same location as in the experimental structure; however, the ß-switch of GPIbα attained a different conformation. Our analysis showed that subsequent refolding of the ß-switch results in a more stable binding configuration, although the transition to the native configuration appears to take some time, during which OS1 could dissociate. Our results show that conformational changes in the ß-switch are crucial for successful binding of OS1. Furthermore, we identified several allosteric binding sites of GPIbα that might also interfere with vWF binding, and optimization of the peptide to target these allosteric sites might lead to a more effective inhibitor, as these are not dependent on the ß-switch conformation.

Mutual induced-fit mechanism drives binding between intrinsically disordered Bim and cryptic binding site of Bcl-xL.

The intrinsically disordered region (IDR) of Bim binds to the flexible cryptic site of Bcl-xL, a pro-survival protein involved in cancer progression that plays an important role in initiating apoptosis. However, their binding mechanism has not yet been elucidated. We have applied our dynamic docking protocol, which correctly reproduced both the IDR properties of Bim and the native bound configuration, as well as suggesting other stable/meta-stable binding configurations and revealed the binding pathway. Although the cryptic site of Bcl-xL is predominantly in a closed conformation, initial binding of Bim in an encounter configuration leads to mutual induced-fit binding, where both molecules adapt to each other; Bcl-xL transitions to an open state as Bim folds from a disordered to an α-helical conformation while the two molecules bind each other. Finally, our data provides new avenues to develop novel drugs by targeting newly discovered stable conformations of Bcl-xL.

Accurate Binding Configuration Prediction of a G-Protein-Coupled Receptor to Its Antagonist Using Multicanonical Molecular Dynamics-Based Dynamic Docking.

We have performed dynamic docking between a prototypic G-protein-coupled receptor (GPCR) system, the ß2-adrenergic receptor, and its antagonist, alprenolol, using one of the enhanced conformation sampling methods, multicanonical molecular dynamics (McMD), which does not rely on any prior knowledge for the definition of the reaction coordinate. Although we have previously applied our McMD-based dynamic docking protocol to various globular protein systems, its application to GPCR systems would be difficult because of their complicated design, which include a lipid bilayer, and because of the difficulty in sampling the configurational space of a binding site that exists deep inside the GPCR. Our simulations sampled a wide array of ligand-bound and ligand-unbound structures, and we measured 427 binding events during our 48 µs production run. Analysis of the ensemble revealed several stable and meta-stable structures, where the most stable structure at the global free energy minimum matches the experimental one. Additional canonical MD simulations were used for refinement and validation of the structures, revealing that most of the intermediates are sufficiently stable to trap the ligand in these intermediary states and furthermore validated our prediction results. Given the difficulty in reaching the orthosteric binding site, chemical optimization of the compound for the second ranking configuration, which binds near the pocket's entrance, might lead to a high-affinity allosteric inhibitor. Accordingly, we show that the application of our methodology can be used to provide crucial insights for the rational design of drugs that target GPCRs.

Engineering of Aeromonas caviae Polyhydroxyalkanoate Synthase Through Site-Directed Mutagenesis for Enhanced Polymerization of the 3-Hydroxyhexanoate Unit.

Polyhydroxyalkanoate (PHA) synthase is an enzyme that polymerizes the acyl group of hydroxyacyl-coenzyme A (CoA) substrates. Aeromonas caviae PHA synthase (PhaCAc) is an important biocatalyst for the synthesis of a useful PHA copolymer, poly[(R)-3-hydroxybutyrate-co-(R)-3-hydroxyhexanoate] [P(3HB-co-3HHx)]. Previously, a PhaCAc mutant with double mutations in asparagine 149 (replaced by serine [N149S]) and aspartate 171 (replaced by glycine [D171G]) was generated to synthesize a 3HHx-rich P(3HB-co-3HHx) and was named PhaCAc NSDG. In this study, to further increase the 3HHx fraction in biosynthesized PHA, PhaCAc was engineered based on the three-dimensional structural information of PHA synthases. First, a homology model of PhaCAc was built to target the residues for site-directed mutagenesis. Three residues, namely tyrosine 318 (Y318), serine 389 (S389), and leucine 436 (L436), were predicted to be involved in substrate recognition by PhaCAc. These PhaCAc NSDG residues were replaced with other amino acids, and the resulting triple mutants were expressed in the engineered strain of Ralstonia eutropha for application in PHA biosynthesis from palm kernel oil. The S389T mutation allowed the synthesis of P(3HB-co-3HHx) with an increased 3HHx fraction without a significant reduction in PHA yield. Thus, a new workhorse enzyme was successfully engineered for the biosynthesis of a higher 3HHx-fraction polymer.

High-Precision Atomic Charge Prediction for Protein Systems Using Fragment Molecular Orbital Calculation and Machine Learning.

Here, we have constructed neural network-based models that predict atomic partial charges with high accuracy at low computational cost. The models were trained using high-quality data acquired from quantum mechanics calculations using the fragment molecular orbital method. We have succeeded in obtaining highly accurate atomic partial charges for three representative molecular systems of proteins, including one large biomolecule (approx. 2000 atoms). The novelty of our approach is the ability to take into account the electronic polarization in the system, which is a system-dependent phenomenon, being important in the field of drug design. Our high-precision models are useful for the prediction of atomic partial charges and expected to be widely applicable in structure-based drug designs such as structural optimization, high-speed and high-precision docking, and molecular dynamics calculations.

Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking.

Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.

Dynamic Docking of a Medium-Sized Molecule to Its Receptor by Multicanonical MD Simulations.

A medium-sized and highly flexible inhibitor to the enzyme ß-secretase 1 (BACE), which produces the amyloid ß-peptide by cleavage of its precursor protein, was dynamically docked into the large and wide catalytic cleft of BACE that binds to the amyloid-precursor by employing multicanonical molecular dynamics (McMD) simulations. We applied our method to predict the native binding configuration and sample the intermediary structures connecting this natively bound state to the unbound one. Representative structures located at free energy minima obtained from McMD were taken and subjected to canonical simulations to refine and validate them, reproducing the native complex structure in agreement with the experimental data in the most stable structure, i.e., the one at the global minimum. In addition, the binding free energy was estimated by umbrella sampling (US) simulations along representative pathways obtained from the McMD ensemble, followed by weighted histogram analysis to estimate the affinity, which also reproduced the experimental inhibitory affinity. Interestingly, the loss of interactions between the two molecules along the pathway was clearly shown in the free energy landscape, reiterating the fundamental importance of atomistic interactions to the binding affinity between receptor and drug compound. The sampled ensemble by the US simulations smoothly connected the bound and unbound states, refining the binding pathway while staying true to the McMD ensemble.

Head-head interactions of resting myosin crossbridges in intact frog skeletal muscles, revealed by synchrotron x-ray fiber diffraction.

The intensities of the myosin-based layer lines in the x-ray diffraction patterns from live resting frog skeletal muscles with full thick-thin filament overlap from which partial lattice sampling effects had been removed were analyzed to elucidate the configurations of myosin crossbridges around the thick filament backbone to nanometer resolution. The repeat of myosin binding protein C (C-protein) molecules on the thick filaments was determined to be 45.33 nm, slightly longer than that of myosin crossbridges. With the inclusion of structural information for C-proteins and a pre-powerstroke head shape, modeling in terms of a mixed population of regular and perturbed regions of myosin crown repeats along the filament revealed that the myosin filament had azimuthal perturbations of crossbridges in addition to axial perturbations in the perturbed region, producing pseudo-six-fold rotational symmetry in the structure projected down the filament axis. Myosin crossbridges had a different organization about the filament axis in each of the regular and perturbed regions. In the regular region that lacks C-proteins, there were inter-molecular interactions between the myosin heads in axially adjacent crown levels. In the perturbed region that contains C-proteins, in addition to inter-molecular interactions between the myosin heads in the closest adjacent crown levels, there were also intra-molecular interactions between the paired heads on the same crown level. Common features of the interactions in both regions were interactions between a portion of the 50-kDa-domain and part of the converter domain of the myosin heads, similar to those found in the phosphorylation-regulated invertebrate myosin. These interactions are primarily electrostatic and the converter domain is responsible for the head-head interactions. Thus multiple head-head interactions of myosin crossbridges also characterize the switched-off state and have an important role in the regulation or other functions of myosin in thin filament-regulated muscles as well as in the thick filament-regulated muscles.

Structural alterations of thin actin filaments in muscle contraction by synchrotron X-ray fiber diffraction.

Strong evidence has been accumulated that the conformational changes of the thin actin filaments are occurring and playing an important role in the entire process of muscle contraction. The conformational changes and the mechanical properties of the thin actin filaments we have found by X-ray fiber diffraction on skeletal muscle contraction are explored. Recent studies on the conformational changes of regulatory proteins bound to actin filaments upon activation and in the force generation process are also described. Finally, the roles of structural alterations and dynamics of the actin filaments are discussed in conjunction with the regulation mechanism and the force generation mechanism.

Axial dispositions and conformations of myosin crossbridges along thick filaments in relaxed and contracting states of vertebrate striated muscles by X-ray fiber diffraction.

X-ray diffraction patterns from live vertebrate striated muscles were analyzed to elucidate the detailed structural models of the myosin crown arrangement and the axial disposition of two-headed myosin crossbridges along the thick filaments in the relaxed and contracting states. The modeling studies were based upon the previous notion that individual myosin filaments had a mixed structure with two regions, a "regular" and a "perturbed". In the relaxed state the distributions and sizes of the regular and perturbed regions on myosin filaments, each having its own axial periodicity for the arrangement of crossbridge crowns within the basic period, were similar to those reported previously. A new finding was that in the contracting state, this mixed structure was maintained but the length of each region, the periodicities of the crowns and the axial disposition of two heads of a crossbridge were altered. The perturbed regions of the crossbridge repeat shifted towards the Z-bands in the sarcomere without changing the lengths found in the relaxed state, but in which the intervals between three successive crowns within the basic period became closer to the regular 14.5-nm repeat in the contracting state. In high resolution modeling for a myosin head, the two heads of a crossbridge were axially tilted in opposite directions along the three-fold helical tracks of myosin filaments and their axial orientations were different from each other in perturbed and regular regions in both states. Under relaxing conditions, one head of a double-headed crossbridge pair appeared to be in close proximity to another head in a pair at the adjacent crown level in the axial direction in the regular region. In the perturbed region this contact between heads occurred only on the narrower inter-crown levels. During contraction, one head of a crossbridge oriented more perpendicular to the fiber axis and the partner head flared axially. Several factors that significantly influence the intensities of the myosin based-meridional reflections and their relative contributions are discussed.

Modeling analysis of myosin-based meridional X-ray reflections from frog skeletal muscles in relaxed and contracting states.

Analysis of the myosin-based meridional intensity data in the X-ray diffraction patterns of live frog skeletal muscles was performed to propose a more precise model for a myosin crown periodicity and an axial disposition of two-headed crossbridges along the thick filament in a sarcomere. Modeling studies revealed that the thick filament has a mixed structure of two different periodicities of the myosin crossbridge crown arrangement and that the crown periodicity and the axial disposition of crossbridges are altered when muscle goes from the relaxed state to the contracting state. Factors that primarily affect the meridional intensities were examined.