Protein cavities and tunnels are critical for function. Ligand recognition and binding, transport, and enzyme catalysis require cavities rearrangements. Therefore, the flexibility of cavities should be guaranteed by protein vibrational dynamics. Molecular dynamics simulations provide a framework to explore conformational plasticity of protein cavities. Herein, we present a novel procedure to characterize the dynamics of protein cavities in terms of their volume gradient vector. For this purpose, we make use of algorithms for calculation of the cavity volume that result robust for numerical differentiations. Volume gradient vector is expressed in terms of principal component analysis obtained from equilibrated molecular dynamics simulations. We analyze contributions of principal component modes to the volume gradient vector according to their frequency and degree of delocalization. In all our test cases, we find that low frequency modes play a critical role together with minor contributions of high frequency modes. These modes involve concerted motions of significant fractions of the total residues lining the cavities. We make use of variations of the potential energy of a protein in the direction of the volume gradient vector as a measure of flexibility of the cavity. We show that proteins whose collective low frequency fluctuations contribute the most to changes of cavity volume exhibit more flexible cavities.

Epidermal Growth Factor Receptor (EGFR), a tyrosine kinase receptor, is one of the main tumor markers in different types of cancers. The kinase native state is mainly composed of two populations of conformers: active and inactive. Several sequence variations in EGFR kinase region promote the differential enrichment of conformers with higher activity. Some structural characteristics have been proposed to differentiate kinase conformations, but these considerations could lead to ambiguous classifications. We present a structural characterisation of EGFR kinase conformers, focused on active site pocket comparisons, and the mapping of known pathological sequence variations. A structural based clustering of this pocket accurately discriminates active from inactive, well-characterised conformations. Furthermore, this main pocket contains, or is in close contact with, ≈65% of cancer-related variation positions. Although the relevance of protein dynamics to explain biological function has been extensively recognised, the usage of the ensemble of conformations in dynamic equilibrium to represent the functional state of proteins and the importance of pockets, cavities and/or tunnels was often neglected in previous studies. These functional structures and the equilibrium between them could be structurally analysed in wild type as well as in sequence variants. Our results indicate that biologically important pockets, as well as their shape and dynamics, are central to understanding protein function in wild-type, polymorphic or disease-related variations.

Antibodies have become the Swiss Army tool for molecular biology and nanotechnology. Their outstanding ability to specifically recognise molecular antigens allows their use in many different applications from medicine to the industry. Moreover, the improvement of conventional structural biology techniques (e.g., X-ray, NMR) as well as the emergence of new ones (e.g., Cryo-EM), have permitted in the last years a notable increase of resolved antibody-antigen structures. This offers a unique opportunity to perform an exhaustive structural analysis of antibody-antigen interfaces by employing the large amount of data available nowadays. To leverage this factor, different geometric as well as chemical descriptors were evaluated to perform a comprehensive characterization.

Lipid-binding proteins (LBPs) are soluble proteins responsible for the uptake, transport, and storage of a large variety of hydrophobic lipophilic molecules including fatty acids, steroids, and other lipids in the cellular environment. Among the LBPs, fatty acid binding proteins (FABPs) present preferential binding affinities for long-chain fatty acids. While most of FABPs in vertebrates and invertebrates present similar β-barrel structures with ligands accommodated in their central cavity, parasitic nematode worms exhibit additional unusual α-helix rich fatty acid-and retinol-binding proteins (FAR). Herein, we report the comparison of extended molecular dynamics (MD) simulations performed on the ligand-free and palmitic acid-bond states of the Necator americanus FAR-1 (Na-FAR-1) with respect to other classical β-barrel FABPs. Principal component analysis (PCA) has been used to identify the different conformations adopted by each system during MD simulations. The αhelix fold encompasses a complex internal ligand-binding cavity with a remarkable conformational plasticity that allows reversible switching between distinct states in the holo-Na-FAR-1. The cavity can change up to one-third of its size affected by conformational changes of the protein−ligand complex. Besides, the ligand inside the cavity is not fixed but experiences large conformational changes between bent and stretched conformations. These changes in the ligand conformation follow changes in the cavity size dictated by the transient protein conformation. On the contrary, protein−ligand complex in β-barrel FABPs fluctuates around a unique conformation. The significantly more flexible holo-Na-FAR-1 ligand-cavity explains its larger ligand multiplicity respect to β-barrel FABPs. 38 perform very important functions for parasite growth and 39 development. 40 Fatty acid and retinol binding proteins (FARs) are LBPs that 41 have been described as components of E/S fluids from 42 parasitic nematodes 4−7 and they are hypothesized to play 43 essential roles in lipid acquisition and distribution of nutrients 44 as well as potential dampening of host's immune response. 8,9 45 FARs together with nematode polyprotein/allergens (NPAs) 10 46 are small (14−20 kDa), helix rich proteins that bind retinol 47 and fatty acids and have no recognizable counterparts in other 48 animal groups. 2 Given these characteristics, FARs have been 49 proven to be useful for serodiagnosis and experimental