The Treatment of Solvent in Multiscale Biophysical Modeling

Elucidating the multi-targeted anti-amyloid activity and enhanced islet amyloid polypeptide binding of β-wrapins

β-wrapins are engineered binding proteins stabilizing the β-hairpin conformations of amyloidogenic proteins islet amyloid polypeptide (IAPP), amyloid-β, and α-synuclein, thus inhibiting their amyloid propensity. Here, we use computational and experimental methods to investigate the molecular recognition of IAPP by β-wrapins. We show that the multi-targeted, IAPP, amyloid-β, and α-synuclein, binding properties of β-wrapins originate mainly from optimized interactions between β-wrapin residues and sets of residues in the three amyloidogenic proteins with similar physicochemical properties. Our results suggest that IAPP is a comparatively promiscuous β-wrapin target, probably due to the low number of charged residues in the IAPP β-hairpin motif. The sub-micromolar affinity of β-wrapin HI18, specifically selected against IAPP, is achieved in part by salt-bridge formation between HI18 residue Glu10 and the IAPP N-terminal residue Lys1, both located in the flexible N-termini of the interacting proteins. Our findings provide insights towards developing novel protein-based single- or multi-targeted therapeutics.

Evaluating the dynamics and electrostatic interactions of folded proteins in implicit solvents

Three implicit solvent models, namely GBMVII, FACTS, and SCPISM, were evaluated for their abilities to emulate an explicit solvent environment by comparing the simulated conformational ensembles, dynamics, and electrostatic interactions of the Src SH2 domain and the Lyn kinase domain. This assessment in terms of structural features in folded proteins expands upon the use of hydration energy as a metric for comparison. All-against-all rms coordinate deviation, average positional fluctuations, and ion-pair distance distribution were used to compare the implicit solvent models with the TIP3P explicit solvent model. Our study shows that the Src SH2 domains solvated with TIP3P, GBMVII, and FACTS sample similar global conformations. Additionally, the Src SH2 ion-pair distance distributions of solvent-exposed side chains corresponding to TIP3P, GBMVII, and FACTS do not differ substantially, indicating that GBMVII and FACTS are capable of modeling these electrostatic interactions. The ion-pair distance distributions of SCPISM are distinct from others, demonstrating that these electrostatic interactions are not adequately reproduced with the SCPISM model. On the other hand, for the Lyn kinase domain, a non-globular protein with bilobal structure and a large concavity on the surface, implicit solvent does not accurately model solvation to faithfully reproduce partially buried electrostatic interactions and lobe-lobe conformations. Our work reveals that local structure and dynamics of small, globular proteins are modeled well using FACTS and GBMVII. Nonetheless, global conformations and electrostatic interactions in concavities of multi-lobal proteins resulting from simulations with implicit solvent models do not match those obtained from explicit water simulations.

Flexible Boundaries for Multiresolution Solvation: An Algorithm for Spatial Multiscaling in Molecular Dynamics Simulations

An algorithm is proposed for performing molecular dynamics (MD) simulations of a biomolecular solute represented at atomistic resolution surrounded by a surface layer of atomistic fine-grained (FG) solvent molecules within a bulk represented by coarse-grained (CG) solvent beads. The method, called flexible boundaries for multiresolution solvation (FBMS), is based on: (i) a three-region layering of the solvent around the solute, involving an FG layer surrounded by a mixed FG-CG buffer layer, itself surrounded by a bulk CG region; (ii) a definition of the layer boundary that relies on an effective distance to the solute surface and is thus adapted to the shape of the solute as well as adjusts to its conformational changes. The effective surface distance is defined by inverse-nth power averaging over the distances to all non-hydrogen solute atoms, and the layering is enforced by means of half-harmonic distance restraints, attractive for the FG molecules and repulsive for the CG beads. A restraint-free region at intermediate distances enables the formation of the buffer layer, where the FG and CG solvents can mix freely. The algorithm is tested and validated using the GROMOS force field and the associated FG (SPC) and CG (polarizable CGW) water models. The test systems include pure-water systems, where one FG molecule plays the role of a solute, and a deca-alanine peptide with two widely different solute shapes considered, α-helical and fully extended. In particular, as the peptide unfolds, the number of FG molecules required to fill its close-range solvation layer increases, with the additional molecules being provided by the buffer layer. Further validation involves simulations of four proteins in multiresolution FG/CG mixtures. The resulting structural, energetic, and solvation properties are found to be similar to those observed in corresponding pure FG simulations.

Modeling Protein-Micelle Systems in Implicit Water

Several membrane proteins and numerous membrane-active peptides have been studied in detergent micelles by solution NMR. However, the detailed structure of these complexes remains unknown. We propose a modeling approach that treats the protein and detergent in atomistic detail and the solvent implicitly. The model is based on previous work on dodecylphosphocholine micelles, adapted for use with the CHARMM36 force field and extended to sodium dodecyl sulfate micelles. Solvation parameters were slightly adjusted to reproduce experimental data on aggregation numbers and critical micelle concentrations. To test the approach, several membrane-active peptides and three β-barrel membrane proteins were subjected to molecular dynamics simulations in the presence of a large number of detergent molecules. Their experimentally determined secondary structure was maintained and the RMSD values were less than 2 Å. Deformations were commonly observed in the N or C termini. The atomistic view of the protein-micelle systems that this approach provides could be useful in interpreting biophysical experiments carried out in the presence of detergent.