Ryuhei Harada

Protein crowding affects hydration structure and dynamics.

The effect of protein crowding on the structure and dynamics of water was examined from explicit solvent molecular dynamics simulations of a series of protein G and protein G/villin systems at different protein concentrations. Hydration structure was analyzed in terms of radial distribution functions, three-dimensional hydration sites, and preservation of tetrahedral coordination. Analysis of hydration dynamics focused on self-diffusion rates and dielectric constants as a function of crowding. The results show significant changes in both structure and dynamics of water under highly crowded conditions. The structure of water is altered mostly beyond the first solvation shell. Diffusion rates and dielectric constants are significantly reduced following linear trends as a function of crowding reflecting highly constrained water in crowded environments. The reduced dynamics of diffusion is expected to be strongly related to hydrodynamic properties of crowded cellular environments while the reduced dielectric constant under crowded conditions has implications for the stability of biomolecules in crowded environments. The results from this study suggest a prescription for modeling solvation in simulations of cellular environments.

Reduced native state stability in crowded cellular environment due to protein-protein interactions.

The effect of cellular crowding environments on protein structure and stability is a key issue in molecular and cellular biology. The classical view of crowding emphasizes the volume exclusion effect that generally favors compact, native states. Here, results from molecular dynamics simulations and NMR experiments show that protein crowders may destabilize native states via protein-protein interactions. In the model system considered here, mixtures of villin head piece and protein G at high concentrations, villin structures become increasingly destabilized upon increasing crowder concentrations. The denatured states observed in the simulation involve partial unfolding as well as more subtle conformational shifts. The unfolded states remain overall compact and only partially overlap with unfolded ensembles at high temperature and in the presence of urea. NMR measurements on the same systems confirm structural changes upon crowding based on changes of chemical shifts relative to dilute conditions. An analysis of protein-protein interactions and energetic aspects suggests the importance of enthalpic and solvation contributions to the crowding free energies that challenge an entropic-centered view of crowding effects.

Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm

Biological macromolecules function in highly crowded cellular environments. The structure and dynamics of proteins and nucleic acids are well characterized in vitro, but in vivo crowding effects remain unclear. Using molecular dynamics simulations of a comprehensive atomistic model cytoplasm we found that protein-protein interactions may destabilize native protein structures, whereas metabolite interactions may induce more compact states due to electrostatic screening. Protein-protein interactions also resulted in significant variations in reduced macromolecular diffusion under crowded conditions, while metabolites exhibited significant two-dimensional surface diffusion and altered protein-ligand binding that may reduce the effective concentration of metabolites and ligands in vivo. Metabolic enzymes showed weak non-specific association in cellular environments attributed to solvation and entropic effects. These effects are expected to have broad implications for the in vivo functioning of biomolecules. This work is a first step towards physically realistic in silico whole-cell models that connect molecular with cellular biology. DOI: http://dx.doi.org/10.7554/eLife.19274.001

Complete atomistic model of a bacterial cytoplasm for integrating physics, biochemistry, and systems biology.

A model for the cytoplasm of Mycoplasma genitalium is presented that integrates data from a variety of sources into a physically and biochemically consistent model. Based on gene annotations, core genes expected to be present in the cytoplasm were determined and a metabolic reaction network was reconstructed. The set of cytoplasmic genes and metabolites from the predicted reactions were assembled into a comprehensive atomistic model consisting of proteins with predicted structures, RNA, protein/RNA complexes, metabolites, ions, and solvent. The resulting model bridges between atomistic and cellular scales, between physical and biochemical aspects, and between structural and systems views of cellular systems and is meant as a starting point for a variety of simulation studies.

Exploring the Folding Free Energy Landscape of a β-Hairpin Miniprotein, Chignolin, Using Multiscale Free Energy Landscape Calculation Method

The folding process for a β-hairpin miniprotein, chignolin, was investigated by free energy landscape (FEL) calculations using the recently proposed multiscale free energy landscape calculation method (MSFEL). First, coarse-grained molecular dynamics simulations searched a broad conformational space, then multiple independent, all-atom molecular dynamics simulations with explicit solvent determined the detailed local FEL using massively distributed computing. The combination of the two models enabled efficient calculation of the free energy landscapes. The MSFEL analysis showed that chignolin has an intermediate state as well as a misfolded state. The folding process is initiated by the formation of a β-hairpin turn, followed by the formation of contacts in the hydrophobic core between Tyr2 and Trp9. Furthermore, mutation of Tyr2 shifts the population to the misfolded conformation. The results indicate that the hydrophobic core plays an important role in stabilizing the native state of chignolin.

On the induced‐fit mechanism of substrate‐enzyme binding structures of nylon‐oligomer hydrolase

We present a detailed computational investigation of the induced‐fit motion in a nylon‐oligomer hydrolase (NylB) upon substrate binding. To this aim, we resort on the recently introduced parallel cascade selection molecular dynamics approach, allowing for an accelerated access to the set of conformational changes from an open‐ to a closed‐state structure to form the enzyme‐substrate complex in a specific induce‐fit mechanism. The structural investigation is quantitatively complemented by free energy analyses within the umbrella sampling algorithm accompanied by weighted histogram analysis. We find that the stabilization free energy is about 1.4 kcal/mol, whereas the highest free energy barrier to be overcome is about 2.3 kcal/mol. Conversely, the energetic contribution for the substrate binding is about 20 kcal/mol, as estimated from Generalized Born/Surface Area. This means that the open‐close induced‐fit motion could occur frequently once the substrate binds to the open state of NylB.

Protein folding pathways extracted by OFLOOD: Outlier FLOODing method

The Outlier FLOODing method (OFLOOD) is proposed as an efficient conformational sampling method to extract biologically rare events such as protein folding. In OFLOOD, sparse distributions (outliers in the conformational space) were regarded as relevant states for the transitions. Then, the transitions were enhanced through conformational resampling from the outliers. This evidence indicates that the conformational resampling of the sparse distributions might increase chances for promoting the transitions from the outliers to other meta‐stable states, which resembles a conformational flooding from the outliers to the neighboring clusters. OFLOOD consists of (i) detections of outliers from conformational distributions and (ii) conformational resampling from the outliers by molecular dynamics (MD) simulations. Cycles of (i) and (ii) are simply repeated. As demonstrations, OFLOOD was applied to folding of Chignolin and HP35. In both cases, OFLOOD automatically extracted folding pathways from unfolded structures with ns‐order computational costs, although µs‐order canonical MD failed to extract them.

The Fast-Folding Mechanism of Villin Headpiece Subdomain Studied by Multiscale Distributed Computing.

The fast-folding mechanism of a 35-residue mini-protein, villin headpiece subdomain (HP35), was investigated using folding free energy landscape analysis with the multiscale free energy landscape calculation method (MSFEL). A major and a minor folding pathway were deduced from the folding free energy landscape. In the major folding pathway, the formation of helices II and III was the rate-limiting step in the transition to an intermediate state, triggered by the folding of the PLWK motif. HP35 then folds into the native structure through the formation of the hydrophobic core located at the center of the three-helix bundle. Mutations in the motif and hydrophobic core that suppressed folding into the native state drastically changed the folding free energy landscape compared to the wild type protein. In the minor folding pathway, nucleation of the hydrophobic core preceded formation of the motif.

Enhanced conformational sampling method for proteins based on the TaBoo SeArch algorithm: Application to the folding of a mini‐protein, chignolin

The conformational samplings are indispensible for obtaining reliable canonical ensembles, which provide statistical averages of physical quantities such as free energies. However, the samplings of vast conformational space of biomacromolecules by conventional molecular dynamics (MD) simulations might be insufficient, due to their inadequate accessible time‐scales for investigating biological functions. Therefore, the development of methodologies for enhancing the conformational sampling of biomacromolecules still remains as a challenging issue in computational biology. To tackle this problem, we newly propose an efficient conformational search method, which is referred as TaBoo SeArch (TBSA) algorithm. In TBSA, an inverse energy histogram is used to select seeds for the conformational resampling so that states with high frequencies are inhibited, while states with low frequencies are efficiently sampled to explore the unvisited conformational space. As a demonstration, TBSA was applied to the folding of a mini‐protein, chignolin, and automatically sampled the native structure (Cα root mean square deviation < 1.0 Å) with nanosecond order computational costs started from a completely extended structure, although a long‐time 1‐µs normal MD simulation failed to sample the native structure. Furthermore, a multiscale free energy landscape method based on the conformational sampling of TBSA were quantitatively evaluated through free energy calculations with both implicit and explicit solvent models, which enable us to find several metastable states on the folding landscape.

Nontargeted Parallel Cascade Selection Molecular Dynamics for Enhancing the Conformational Sampling of Proteins.

Nontargeted parallel cascade selection molecular dynamics (nt-PaCS-MD) is proposed as an efficient conformational sampling method to enhance the conformational transitions of proteins, which is an extension of the original targeted PaCS-MD (t-PaCS-MD). The original PaCS-MD comprises cycles of (i) selection of initial structures for multiple independent MD simulations toward a predetermined target and (ii) conformational sampling by the independent MDs. In nt-PaCS-MD, structures that significantly deviate from an average are regarded as candidates that have high potential to address other metastable states and are chosen as the initial structures in the selection. To select significantly deviated structures, we examine the root-mean-square deviation (RMSD) of snapshots generated from the average structure based on Gram-Schmidt orthogonalization. nt-PaCS-MD was applied to the folding of the mini-protein chignolin in implicit solvent and to the open-closed conformational transitions of T4 lysozyme (T4L) and glutamine binding protein (QBP) in explicit solvent. We show that nt-PaCS-MD can reach chignolins native state and can also cause the open-closed transition of T4L and QBP on a nanosecond time scale, which are very efficient in terms of conformational sampling and comparable to that with t-PaCS-MD.