The effect of end constraints on protein loop kinematics

Morphing and docking visualisation of biomolecular structures using Multi-Dimensional Scaling

Protein structures are often solved at atomic resolution in two states defining a functional movement but intervening conformations are usually unknown. Morphing methods generate intervening conformations between two known structures. When viewed as an animation using molecular graphics, a smooth, direct morph enables the eye to track changes in structure that might be otherwise missed. We present a morphing method that aims to linearly interpolate interatomic distances and which uses SMACOF (Scaling by MAjorisation of COmplicated Function) and multigrid techniques with a cut-off distance based weighting that optimizes the MolProbity score of intervening structures. The all-atom morphs are smooth, move directly between the two structures, and are shown, in general, to pass closer to a set of known intermediates than those generated using other methods. The techniques are also used for docking by putting the unbound structures in a "near-approach pose" and then morphing to the bound complex. The resulting GPU-accelerated tools are available on a webserver, Morphit_Pro, at http://morphit-pro.cmp.uea.ac.uk/ and more than 5000 domains movements available at the DynDom website can now be viewed as morphs http://morphit-pro.cmp.uea.ac.uk/dyndom/.

Loop A is critical for the functional interaction of two Beta vulgaris PIP aquaporins

Research done in the last years strongly support the hypothesis that PIP aquaporin can form heterooligomeric assemblies, specially combining PIP2 monomers with PIP1 monomers. Nevertheless, the structural elements involved in the ruling of homo versus heterooligomeric organization are not completely elucidated. In this work we unveil some features of monomer-monomer interaction in Beta vulgaris PIP aquaporins. Our results show that while BvPIP2;2 is able to interact with BvPIP1;1, BvPIP2;1 shows no functional interaction. The lack of functional interaction between BvPIP2;1 and BvPIP1;1 was further corroborated by dose-response curves of water permeability due to aquaporin activity exposed to different acidic conditions. We also found that BvPIP2;1 is unable to translocate BvPIP1;1-ECFP from an intracellular position to the plasma membrane when co-expressed, as BvPIP2;2 does. Moreover we postulate that the first extracellular loop (loop A) of BvPIP2;1, could be relevant for the functional interaction with BvPIP1;1. Thus, we investigate BvPIP2;1 loop A at an atomic level by Molecular Dynamics Simulation (MDS) and by direct mutagenesis. We found that, within the tetramer, each loop A presents a dissimilar behavior. Besides, BvPIP2;1 loop A mutants restore functional interaction with BvPIP1;1. This work is a contribution to unravel how PIP2 and PIP1 interact to form functional heterooligomeric assemblies. We postulate that BvPIP2;1 loop A is relevant for the lack of functional interaction with BvPIP1;1 and that the monomer composition of PIP assemblies determines their functional properties.

Assessing protein loop flexibility by hierarchical Monte Carlo sampling

Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.