Solvent quality and nonbiological oligomer folding: revisiting conventional paradigms.

We report on extensive molecular dynamics atomistic simulations of a meta-substituted poly-phenylacetylene (pPA) foldamer dispersed in three solvents, water H2O, cyclohexane cC6H12, and n-hexane nC6H14, and for three oligomer lengths 12mer, 16mer and 20mer. At room temperature, we find a tendency of the pPA foldamer to collapse into a helical structure in all three solvents but with rather different stability character, stable in water, marginally stable in n-hexane, unstable in cyclohexane. In the case of water, the initial and final number of hydrogen bonds of the foldamer with water molecules is found to be unchanged, with no formation of intrachain hydrogen bonding, thus indicating that hydrogen bonding plays no role in the folding process. In all three solvents, the folding is found to be mainly driven by electrostatics, nearly identical in the three cases, and largely dominant compared to van der Waals interactions that are different in the three cases. This scenario is also supported by the analysis of distribution of the bond and dihedral angles and by a direct calculation of the solvation and transfer free energies via thermodynamic integration. The different stability in the case of cyclohexane and n-hexane notwithstanding their rather similar chemical composition can be traced back to the different entropy-enthalpy compensation that is found similar for water and n-hexane, and very different for cyclohexane. A comparison with the same properties for poly-phenylalanine oligomers underscores the crucial differences between pPA and peptides. To highlight how these findings can hardly be interpreted in terms of a simple "good" and "poor" solvent picture, a molecular dynamics study of a bead-spring polymer chain in a Lennard-Jones fluid is also included.

Water will Find Its Way: Transport through Narrow Tunnels in Hydrolases.

An aqueous environment is vital for life as we know it, and water is essential for nearly all biochemical processes at the molecular level. Proteins utilize water molecules in various ways. Consequently, proteins must transport water molecules across their internal network of tunnels to reach the desired action sites, either within them or by functioning as molecular pipes to control cellular osmotic pressure. Despite water playing a crucial role in enzymatic activity and stability, its transport has been largely overlooked, with studies primarily focusing on water transport across membrane proteins. The transport of molecules through a protein's tunnel network is challenging to study experimentally, making molecular dynamics simulations the most popular approach for investigating such events. In this study, we focused on the transport of water molecules across three different α/ß-hydrolases: haloalkane dehalogenase, epoxide hydrolase, and lipase. Using a 5 µs adaptive simulation per system, we observed that only a few tunnels were responsible for the majority of water transport in dehalogenase, in contrast to a higher diversity of tunnels in other enzymes. Interestingly, water molecules could traverse narrow tunnels with subangstrom bottlenecks, which is surprising given the commonly accepted water molecule radius of 1.4 Å. Our analysis of the transport events in such narrow tunnels revealed a markedly increased number of hydrogen bonds formed between the water molecules and protein, likely compensating for the steric penalty of the process. Overall, these commonly disregarded narrow tunnels accounted for ∼20% of the total water transport observed, emphasizing the need to surpass the standard geometrical limits on the functional tunnels to properly account for the relevant transport processes. Finally, we demonstrated how the obtained insights could be applied to explain the differences in a mutant of the human soluble epoxide hydrolase associated with a higher incidence of ischemic stroke.

Rinmaker: a fast, versatile and reliable tool to determine residue interaction networks in proteins.

BACKGROUND: Residue Interaction Networks (RINs) map the crystallographic description of a protein into a graph, where amino acids are represented as nodes and non-covalent bonds as edges. Determination and visualization of a protein as a RIN provides insights on the topological properties (and hence their related biological functions) of large proteins without dealing with the full complexity of the three-dimensional description, and hence it represents an invaluable tool of modern bioinformatics. RESULTS: We present RINmaker, a fast, flexible, and powerful tool for determining and visualizing RINs that include all standard non-covalent interactions. RINmaker is offered as a cross-platform and open source software that can be used either as a command-line tool or through a web application or a web API service. We benchmark its efficiency against the main alternatives and provide explicit tests to show its performance and its correctness. CONCLUSIONS: RINmaker is designed to be fully customizable, from a simple and handy support for experimental research to a sophisticated computational tool that can be embedded into a large computational pipeline. Hence, it paves the way to bridge the gap between data-driven/machine learning approaches and numerical simulations of simple, physically motivated, models.

TransportTools: a library for high-throughput analyses of internal voids in biomolecules and ligand transport through them.

SUMMARY: Information regarding pathways through voids in biomolecules and their roles in ligand transport is critical to our understanding of the function of many biomolecules. Recently, the advent of high-throughput molecular dynamics simulations has enabled the study of these pathways, and of rare transport events. However, the scale and intricacy of the data produced requires dedicated tools in order to conduct analyses efficiently and without excessive demand on users. To fill this gap, we developed the TransportTools, which allows the investigation of pathways and their utilization across large, simulated datasets. TransportTools also facilitates the development of custom-made analyses. AVAILABILITY AND IMPLEMENTATION: TransportTools is implemented in Python3 and distributed as pip and conda packages. The source code is available at https://github.com/labbit-eu/transport_tools. Data are available in a repository and can be accessed via a link: https://doi.org/10.5281/zenodo.5642954. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Solvent quality and solvent polarity in polypeptides.

Using molecular dynamics and thermodynamic integration, we report on the solvation process of seven polypeptides (GLY, ALA, ILE, ASN, LYS, ARG, GLU) in water and in cyclohexane. The polypeptides are selected to cover the full hydrophobic scale while varying their chain length from tri- to undeca-homopeptides, providing indications on possible non-additivity effects as well as the role of the peptide backbone in the overall stability of the polypeptides. The use of different solvents and different polypeptides allows us to investigate the relation between solvent quality - the capacity of a given solvent to fold/unfold a given biopolymer often described on a scale ranging from "good" to "poor"; and solvent polarity - related to the specific interactions of any solvent with respect to a reference solvent. Undeca-glycine is found to be the only polypeptide to have a stable collapse in water (polar solvent), with the other hydrophobic polypeptides displaying repeated folding and unfolding events in water, with polar polypeptides presenting even more complex behavior. By contrast, all polypeptides are found to keep an extended conformation in cyclohexane, irrespective of their polarity. All considered polypeptides are also found to have favorable solvation free energy independent of the solvent polarity and their intrinsic hydrophobicity, clearly highlighting the prominent stabilizing role of the peptide backbone - with the solvation process largely enthalpically dominated in polar polypeptides and partially entropically driven for hydrophobic polypeptides. Our study thus reveals the complexity of the solvation process of polypeptides defying the common view "like dissolves like", with the solute polarity playing the most prominent role. The absence of mirror symmetry upon the inversion of polarities of both the solvent and the polypeptides is confirmed.

Insights on peptide topology in the computational design of protein ligands: the example of lysozyme binding peptides.

Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.

Can the roles of polar and non-polar moieties be reversed in non-polar solvents?

Using thermodynamic integration, we study the solvation free energy of 18 amino acid side chain equivalents in solvents with different polarities, ranging from the most polar water to the most non-polar cyclohexane. The amino acid side chain equivalents are obtained from the 20 natural amino acids by replacing the backbone part with a hydrogen atom, and discarding proline and glycine that have special properties. A detailed analysis of the relative solvation free energies suggests how it is possible to achieve a robust and unambiguous hydrophobic scale for the amino acids. By discriminating the relative contributions of the entropic and enthalpic terms, we find strong negative correlations in water and ethanol, associated with the well-known entropy-enthalpy compensation, and a much reduced correlation in cyclohexane. This shows that in general the role of the polar and non-polar moieties cannot be reversed in a non-polar solvent. Our findings are compared with past experimental as well as numerical results, and may shed additional light on the unique role of water as a biological solvent.

Computational design of cyclic peptides for the customized oriented immobilization of globular proteins.

The oriented immobilization of proteins, key for the development of novel responsive biomaterials, relies on the availability of effective probes. These are generally provided by standard approaches based on in vivo maturation and in vitro selection of antibodies and/or aptamers. These techniques can suffer technical problems when a non-immunogenic epitope needs to be targeted. Here we propose a strategy to circumvent this issue by in silico design. In our method molecular binders, in the form of cyclic peptides, are computationally evolved by stochastically exploring their sequence and structure space to identify high-affinity peptides for a chosen epitope of a target globular protein: here a solvent-exposed site of ß2-microglobulin (ß2m). Designed sequences were screened by explicit solvent molecular dynamics simulations (MD) followed by experimental validation. Five candidates gave dose-response surface plasmon resonance signals with dissociation constants in the micromolar range. One of them was further analyzed by means of isothermal titration calorimetry, nuclear magnetic resonance, and 250 ns of MD. Atomic-force microscopy imaging showed that this peptide is able to immobilize ß2m on a gold surface. In short, we have shown by a variety of experimental techniques that it is possible to capture a protein through an epitope of choice by computational design.

How perfluoroalkyl substances modify fluorinated self-assembled monolayer architectures: An electrochemical and computational study.

There is an urgent need for sensing strategies to screen perfluoroalkyl substances (PFAS) in aqueous matrices. These strategies must be applicable in large-scale monitoring plans to face the ubiquitous use of PFAS, their wide global spread, and their fast evolution towards short-chain, branched molecules. To this aim, the changes in fluorinated self-assembled monolayers (SAM) with different architectures (pinholes/defects-free and with randomized pinholes/defects) were studied upon exposure to both long and short-chain PFAS. The applicability of fluorinated SAM in PFAS sensing was evaluated. Changes in the SAM structures were characterised combining electrochemical impedance spectroscopy and voltammetric techniques. The experimental data interpretation was supported by molecular dynamics simulations to gain a more in-depth understanding of the interaction mechanisms involved. Pinhole/defect-free fluorinated SAM were found to be applicable to long-chain PFAS screening within switch-on sensing strategy, while a switch-off sensing strategy was reported for screening of both short/long-chain PFAS. These strategies confirmed the possibility to play on fluorophilic interactions when designing PFAS screening methods.

Dynamics and Thermodynamics of Transthyretin Association from Molecular Dynamics Simulations.

Molecular dynamics simulations are used in this work to probe the structural stability and the dynamics of engineered mutants of transthyretin (TTR), i.e., the double mutant F87M/L110M (MT-TTR) and the triple mutant F87M/L110M/S117E (3M-TTR), in relation to wild-type. Free energy analysis from end-point simulations and statistical effective energy functions are used to analyze trajectories, revealing that mutations do not have major impact on protein structure but rather on protein association, shifting the equilibria towards dissociated species. The result is confirmed by the analysis of 3M-TTR which shows dissociation within the first 10 ns of the simulation, indicating that contacts are lost at the dimer-dimer interface, whereas dimers (formed by monomers which pair to form two extended ß-sheets) appear fairly stable. Overall the simulations provide a detailed view of the dynamics and thermodynamics of wild-type and mutant transthyretins and a rationale of the observed effects.

Molecular dynamics simulations of ß2-microglobulin interaction with hydrophobic surfaces.

Hydrophobic surfaces are known to adsorb and unfold proteins, a process that has been studied only for a few proteins. Here we address the interaction of ß2-microglobulin, a paradigmatic protein for the study of amyloidogenesis, with hydrophobic surfaces. A system with 27 copies of the protein surrounded by a model cubic hydrophobic box is studied by implicit solvent molecular dynamics simulations. Most proteins adsorb on the walls of the box without major distortions in local geometry, whereas free molecules maintain proper structures and fluctuations as observed in explicit solvent molecular dynamics simulations. The major conclusions from the simulations are as follows: (i) the adopted implicit solvent model is adequate to describe protein dynamics and thermodynamics; (ii) adsorption occurs readily and is irreversible on the simulated timescale; (iii) the regions most involved in molecular encounters and stable interactions with the walls are the same as those that are important in protein-protein and protein-nanoparticle interactions; (iv) unfolding following adsorption occurs at regions found to be flexible by both experiments and simulations; (v) thermodynamic analysis suggests a very large contribution from van der Waals interactions, whereas unfavorable electrostatic interactions are not found to contribute much to adsorption energy. Surfaces with different degrees of hydrophobicity may occur in vivo. Our simulations show that adsorption is a fast and irreversible process which is accompanied by partial unfolding. The results and the thermodynamic analysis presented here are consistent with and rationalize previous experimental work.