Establishment of a deep-learning system to diagnose BI-RADS4a or higher using breast ultrasound for clinical application.

Although the categorization of ultrasound using the Breast Imaging Reporting and Data System (BI-RADS) has become widespread worldwide, the problem of inter-observer variability remains. To maintain uniformity in diagnostic accuracy, we have developed a system in which artificial intelligence (AI) can distinguish whether a static image obtained using a breast ultrasound represents BI-RADS3 or lower or BI-RADS4a or higher to determine the medical management that should be performed on a patient whose breast ultrasound shows abnormalities. To establish and validate the AI system, a training dataset consisting of 4028 images containing 5014 lesions and a test dataset consisting of 3166 images containing 3656 lesions were collected and annotated. We selected a setting that maximized the area under the curve (AUC) and minimized the difference in sensitivity and specificity by adjusting the internal parameters of the AI system, achieving an AUC, sensitivity, and specificity of 0.95, 91.2%, and 90.7%, respectively. Furthermore, based on 30 images extracted from the test data, the diagnostic accuracy of 20 clinicians and the AI system was compared, and the AI system was found to be significantly superior to the clinicians (McNemar test, p < 0.001). Although deep-learning methods to categorize benign and malignant tumors using breast ultrasound have been extensively reported, our work represents the first attempt to establish an AI system to classify BI-RADS3 or lower and BI-RADS4a or higher successfully, providing important implications for clinical actions. These results suggest that the AI diagnostic system is sufficient to proceed to the next stage of clinical application.

High-performance drug discovery: computational screening by combining docking and molecular dynamics simulations.

Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, this method is not completely reliable and therefore unsatisfactory. In this study, we used massive molecular dynamics simulations of protein-ligand conformations obtained by molecular docking in order to improve the enrichment performance of molecular docking. Our screening approach employed the molecular mechanics/Poisson-Boltzmann and surface area method to estimate the binding free energies. For the top-ranking 1,000 compounds obtained by docking to a target protein, approximately 6,000 molecular dynamics simulations were performed using multiple docking poses in about a week. As a result, the enrichment performance of the top 100 compounds by our approach was improved by 1.6-4.0 times that of the enrichment performance of molecular dockings. This result indicates that the application of molecular dynamics simulations to virtual screening for lead discovery is both effective and practical. However, further optimization of the computational protocols is required for screening various target proteins.

FGF9 monomer-dimer equilibrium regulates extracellular matrix affinity and tissue diffusion.

The spontaneous dominant mouse mutant, Elbow knee synostosis (Eks), shows elbow and knee joint synosotsis, and premature fusion of cranial sutures. Here we identify a missense mutation in the Fgf9 gene that is responsible for the Eks mutation. Through investigation of the pathogenic mechanisms of joint and suture synostosis in Eks mice, we identify a key molecular mechanism that regulates FGF9 signaling in developing tissues. We show that the Eks mutation prevents homodimerization of the FGF9 protein and that monomeric FGF9 binds to heparin with a lower affinity than dimeric FGF9. These biochemical defects result in increased diffusion of the altered FGF9 protein (FGF9(Eks)) through developing tissues, leading to ectopic FGF9 signaling and repression of joint and suture development. We propose a mechanism in which the range of FGF9 signaling in developing tissues is limited by its ability to homodimerize and its affinity for extracellular matrix heparan sulfate proteoglycans.

Folding dynamics of 10-residue beta-hairpin peptide chignolin.

Short peptides that fold into beta-hairpins are ideal model systems for investigating the mechanism of protein folding because their folding process shows dynamics typical of proteins. We performed folding, unfolding, and refolding molecular dynamics simulations (total of 2.7 micros) of the 10-residue beta-hairpin peptide chignolin, which is the smallest beta-hairpin structure known to be stable in solution. Our results revealed the folding mechanism of chignolin, which comprises three steps. First, the folding begins with hydrophobic assembly. It brings the main chain together; subsequently, a nascent turn structure is formed. The second step is the conversion of the nascent turn into a tight turn structure along with interconversion of the hydrophobic packing and interstrand hydrogen bonds. Finally, the formation of the hydrogen-bond network and the complete hydrophobic core as well as the arrangement of side-chain-side-chain interactions occur at approximately the same time. This three-step mechanism appropriately interprets the folding process as involving a combination of previous inconsistent explanations of the folding mechanism of the beta-hairpin, that the first event of the folding is formation of hydrogen bonds and the second is that of the hydrophobic core, or vice versa.

Structure and dynamics of RNA polymerase II elongation complex.

RNA polymerase (Pol) II is a fundamental and important enzyme in the transcription process. However, two mysterious questions have remained unsolved: how an unwound bubble of DNA is established and maintained, and how the enzyme moves along the DNA. To answer these questions, we constructed a model structure of the Pol II elongation complex with the 50 base pairs of DNA-24 bases of RNA including the unwound bubble of DNA and performed a molecular dynamics simulation. We obtained a reliable model structure of the Pol II elongation complex in the pre-translocation state which has not yet been determined by the X-ray crystallographic study. The model structure revealed that multiple protein loops work concertedly to form and maintain the bubble structure. We also found that the conformational change of a loop in the Pol II, fork loop 1, couples with the unidirectional movement of the Pol II along the DNA.

Novel mechanism of interaction of p85 subunit of phosphatidylinositol 3-kinase and ErbB3 receptor-derived phosphotyrosyl peptides.

Ligand-activated and tyrosine-phosphorylated ErbB3 receptor binds to the SH2 domain of the p85 subunit of phosphatidylinositol 3-kinase and initiates intracellular signaling. Here, we studied the interactions between the N- (N-SH2) and C- (C-SH2) terminal SH2 domains of the p85 subunit of the phosphatidylinositol 3-kinase and eight ErbB3 receptor-derived phosphotyrosyl peptides (P-peptides) by using molecular dynamics, free energy, and surface plasmon resonance (SPR) analyses. In SPR analysis, these P-peptides showed no binding to the C-SH2 domain, but P-peptides containing a phospho-YXXM or a non-phospho-YXXM motif did bind to the N-SH2 domain. The N-SH2 domain has two phosphotyrosine binding sites in its N- (N1) and C- (N2) terminal regions. Interestingly, we found that P-peptides of pY1180 and pY1241 favored to bind to the N2 site, although all other P-peptides showed favorable binding to the N1 site. Remarkably, two phosphotyrosines, pY1178 and pY1243, which are just 63 amino acids apart from the pY1241 and pY1180, respectively, showed favorable binding to the N1 site. These findings indicate a possibility that the pair of phosphotyrosines, pY1178-pY1241 or pY1243-pY1180, will fold into an appropriate configuration for binding to the N1 and N2 sites simultaneously. Our model structures of the cytoplasmic C-terminal domain of ErbB3 receptor also strongly supported the speculation. The calculated binding free energies between the N-SH2 domain and P-peptides showed excellent qualitative agreement with SPR data with a correlation coefficient of 0.91. The total electrostatic solvation energy between the N-SH2 domain and P-peptide was the dominant factor for its binding affinity.

Cooperative motions of protein and hydration water molecules: molecular dynamics study of scytalone dehydratase.

Two molecular dynamics (MD) simulations totaling 25 ns of simulation time of monomeric scytalone dehydratase (SD) were performed. The enzyme has a ligand-binding pocket containing a cone-shaped alpha+beta barrel, and the C-terminal region covers the binding pocket. Our simulations clarified the difference in protein dynamics and conformation between the liganded protein and the unliganded protein. The liganded protein held the ligand molecule tightly and the initial structure was maintained during the simulation. The unliganded protein, on the other hand, fluctuated dynamically and its structure changed largely from the initial structure. In the equilibrium state, the binding pocket was fully solvated by opening of the C-terminal region, and the protein dynamics was connected with hydration water molecules entry into and release from the binding pocket. In addition, the cooperative motions of the unliganded protein and the hydration water molecules produced the path through the protein interior for ligand binding.

Tyr-317 phosphorylation increases Shc structural rigidity and reduces coupling of domain motions remote from the phosphorylation site as revealed by molecular dynamics simulations.

Activated receptor tyrosine kinases bind the Shc adaptor protein through its N-terminal phosphotyrosine-binding (PTB) and C-terminal Src homology 2 (SH2) domains. After binding, Shc is phosphorylated within the central collagen-homology (CH) linker region on Tyr-317, a residue remote to both the PTB and SH2 domains. Shc phosphorylation plays a pivotal role in the initiation of mitogenic signaling through the Ras/Raf/MEK/ERK pathway, but it is unclear if Tyr-317 phosphorylation affects Shc-receptor interactions through the PTB and SH2 domains. To investigate the structural impact of Shc phosphorylation, molecular dynamics simulations were carried out using special-purpose Molecular Dynamics Machine-Grape computers. After a 1-nanosecond equilibration, atomic motions in the structures of unphosphorylated Shc and Shc phosphorylated on Tyr-317 were calculated during a 2-nanosecond period. The results reveal larger phosphotyrosine-binding domain fluctuations and more structural flexibility of unphosphorylated Shc compared with phosphorylated Shc. Collective motions between the PTB-SH2, PTB-CH, and CH-SH2 domains were highly correlated only in unphosphorylated Shc. Dramatic changes in domain coupling and structural rigidity, induced by Tyr-317 phosphorylation, may alter Shc function, bringing about marked differences in the association of unphosphorylated and phosphorylated Shc with its numerous partners, including activated membrane receptors.

Increased rigidity of domain structures enhances the stability of a mutant enzyme created by directed evolution.

A mutant of kanamycin nucleotidyltransferase (KNT) was previously created by directed evolution. This mutant, HTK, has 19 amino acid substitutions, which increase the thermostability by 20 degrees C. In this study, we have examined to what extent each mutation contributes to the increased stability and analyzed how the mutations affect the structure of KNT at 72 degrees C using molecular dynamics simulations. The effects of some mutations on the stability are simply additive, but those of others are cooperative. Mutations with large effects on the stability are introduced into the N-terminal domain, which appears to be less stable than the C-terminal domain. Results of the molecular dynamics simulations have indicated that the rigidity of the domain structures is increased by the mutations: at 72 degrees C, the intradomain fluctuations of HTK are decreased, and in turn, its interdomain motions are pronounced, whereas the structure of the preevolved KNT fluctuates randomly. Chemical modification experiments of cysteine residues have shown that the cysteine residues of HTK are less accessible to an SH reagent than those of the preevolved KNT. The present results suggest that the 19 mutations of HTK stabilize KNT by affecting the dynamic behavior of the structure of this enzyme without significantly changing its static overall structure.

Molecular dynamics, free energy, and SPR analyses of the interactions between the SH2 domain of Grb2 and ErbB phosphotyrosyl peptides.

We studied the interactions between the SH2 domain of growth factor receptor binding protein 2 (Grb2) and ErbB receptor-derived phosphotyrosyl peptides using molecular dynamics, free energy calculations, and surface plasmon resonance (SPR) analysis. Binding free energies for nine phosphotyrosyl peptides were calculated using the MM-PBSA continuum solvent method, and excellent qualitative agreement with the SPR experimental data, with a correlation coefficient of 0.92, was obtained. Consistent with previous experimental findings, phosphotyrosyl peptides with the consensus sequence pYXNX showed favorable binding affinity for the Grb2. Unexpectedly, phosphotyrosyl peptides with the consensus sequence pYQQD, which had not shown any specific binding affinity for the Grb2 in earlier studies, also showed favorable binding affinity for the Grb2 in our experimental and computational analyses. Component analysis of the calculated binding free energies revealed that van der Waals interaction between the Grb2 and the phosphotyrosyl peptide was the dominant factor for specificity and binding affinity. These results indicate that current methods of estimating binding free energies are efficient for obtaining important information about protein-protein interactions, which are essential for the transmission of signals in cellular signaling pathways.