Communication: Role of explicit water models in the helix folding/unfolding processes

Kinetics from Metadynamics: Principles, Applications, and Outlook

Metadynamics is a popular enhanced sampling algorithm for computing the free energy landscape of rare events by using molecular dynamics simulation. Ten years ago, Tiwary and Parrinello introduced the infrequent metadynamics approach for calculating the kinetics of transitions across free energy barriers. Since then, metadynamics-based methods for obtaining rate constants have attracted significant attention in computational molecular science. Such methods have been applied to study a wide range of problems, including protein-ligand binding, protein folding, conformational transitions, chemical reactions, catalysis, and nucleation. Here, we review the principles of elucidating kinetics from metadynamics-like approaches, subsequent methodological developments in this area, and successful applications on chemical, biological, and material systems. We also highlight the challenges of reconstructing accurate kinetics from enhanced sampling simulations and the scope of future developments.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

Dynamics and Thermodynamics of Ibuprofen Conformational Isomerism at the Crystal/Solution Interface

Conformational flexibility of molecules involved in crystal growth and dissolution is rarely investigated in detail and usually considered to be negligible in the formulation of mesoscopic models of crystal growth. In this work, we set out to investigate the conformational isomerism of ibuprofen as it approaches and is incorporated in the morphologically dominant {100} crystal face, in a range of different solvents: water, 1-butanol, toluene, cyclohexanone, cyclohexane, acetonitrile, and trichloromethane. To this end, we combine extensive molecular dynamics and well-tempered metadynamics simulations to estimate the equilibrium distribution of conformers, compute conformer-conformer transition rates, and extract the characteristic relaxation time of the conformer population in solution, adsorbed at the solid/liquid interface, incorporated in the crystal in contact with the mother solution, and in the crystal bulk. We find that, while the conformational equilibrium distribution is weakly dependent on the solvent, relaxation times are instead significantly affected by it. Furthermore, differences in the relaxation dynamics are enhanced on the crystal surface, where conformational transitions become slower and specific conformational transition pathways are hindered. This leads us to observe that the dominant mechanisms of conformational transition can also change significantly moving from the bulk solution to the crystal interface, even for a small molecule with limited conformational flexibility such as ibuprofen. Our findings suggest that understanding conformational flexibility is key to provide an accurate description of the solid/liquid interface during crystal dissolution and growth, and therefore, its relevance should be systematically assessed in the formulation of mesoscopic growth models.

How Does Solvation Layer Mobility Affect Protein Structural Dynamics?

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.