The Role of Conformational Dynamics and Allostery in the Disease Development of Human Ferritin

Determining the three-dimensional structure of myoglobin, the first solved structure of a protein, fundamentally changed the way protein function was understood. Even more revolutionary was the information that came afterward: protein dynamics play a critical role in biological functions. Therefore, understanding conformational dynamics is crucial to obtaining a more complete picture of protein evolution. We recently analyzed the evolution of different protein families including green fluorescent proteins (GFPs), β-lactamase inhibitors, and nuclear receptors, and we observed that the alteration of conformational dynamics through allosteric regulation leads to functional changes. Moreover, proteome-wide conformational dynamics analysis of more than 100 human proteins showed that mutations occurring at rigid residue positions are more susceptible to disease than flexible residue positions. These studies suggest that disease-associated mutations may impair dynamic allosteric regulations, leading to loss of function. Thus, in this study, we analyzed the conformational dynamics of the wild-type light chain subunit of human ferritin protein along with the neutral and disease forms. We first performed replica exchange molecular dynamics simulations of wild-type and mutants to obtain equilibrated dynamics and then used perturbation response scanning (PRS), where we introduced a random Brownian kick to a position and computed the fluctuation response of the chain using linear response theory. Using this approach, we computed the dynamic flexibility index (DFI) for each position in the chain for the wild-type and the mutants. DFI quantifies the resilience of a position to a perturbation and provides a flexibility/rigidity measurement for a given position in the chain. The DFI analysis reveals that neutral variants and the wild-type exhibit similar flexibility profiles in which experimentally determined functionally critical sites act as hinges in controlling the overall motion. However, disease mutations alter the conformational dynamic profile, making hinges more loose (i.e., softening the hinges), thus impairing the allosterically regulated dynamics.     

The Role of Conformational Dynamics and Allostery in the Disease Development of Human Ferritin

Determining the three-dimensional structure of myoglobin, the first solved structure of a protein, fundamentally changed the way protein function was understood. Even more revolutionary was the information that came afterward: protein dynamics play a critical role in biological functions. Therefore, understanding conformational dynamics is crucial to obtaining a more complete picture of protein evolution. We recently analyzed the evolution of different protein families including green fluorescent proteins (GFPs), β-lactamase inhibitors, and nuclear receptors, and we observed that the alteration of conformational dynamics through allosteric regulation leads to functional changes. Moreover, proteome-wide conformational dynamics analysis of more than 100 human proteins showed that mutations occurring at rigid residue positions are more susceptible to disease than flexible residue positions. These studies suggest that disease-associated mutations may impair dynamic allosteric regulations, leading to loss of function. Thus, in this study, we analyzed the conformational dynamics of the wild-type light chain subunit of human ferritin protein along with the neutral and disease forms. We first performed replica exchange molecular dynamics simulations of wild-type and mutants to obtain equilibrated dynamics and then used perturbation response scanning (PRS), where we introduced a random Brownian kick to a position and computed the fluctuation response of the chain using linear response theory. Using this approach, we computed the dynamic flexibility index (DFI) for each position in the chain for the wild-type and the mutants. DFI quantifies the resilience of a position to a perturbation and provides a flexibility/rigidity measurement for a given position in the chain. The DFI analysis reveals that neutral variants and the wild-type exhibit similar flexibility profiles in which experimentally determined functionally critical sites act as hinges in controlling the overall motion. However, disease mutations alter the conformational dynamic profile, making hinges more loose (i.e., softening the hinges), thus impairing the allosterically regulated dynamics.     

Dynamic features of homodimer interfaces calculated by normal‐mode analysis

Knowledge of the dynamic features of protein interfaces is necessary for a deeper understanding of protein–protein interactions. We performed normal‐mode analysis (NMA) of 517 nonredundant homodimers and their protomers to characterize dimer interfaces from a dynamic perspective. The motion vector calculated by NMA for each atom of a dimer was decomposed into internal and external motion vectors in individual component subunits, followed by the averaging of time‐averaged correlations between these vectors over atom pairs in the interface. This averaged correlation coefficient (ACC) was defined for various combinations of vectors and investigated in detail. ACCs decrease exponentially with an increasing interface area and r‐value, that is, interface area divided by the entire subunit surface area. As the r‐value reflects the nature of dimer formation, the result suggests that both the interface area and the nature of dimer formation are responsible for the dynamic properties of dimer interfaces. For interfaces with small or medium r‐values and without intersubunit entanglements, ACCs are found to increase on dimer formation when compared with those in the protomer state. In contrast, ACCs do not increase on dimer formation for interfaces with large r‐values and intersubunit entanglements such as in interwinding dimers. Furthermore, relationships between ACCs for intrasubunit atom pairs and for intersubunit atom pairs are found to significantly differ between interwinding and noninterwinding dimers for external motions. External motions are considered as an important factor for characterizing dimer interfaces.     

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site‐directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson‐Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ～5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures. Proteins 2014; 82:1128–1141. © 2013 Wiley Periodicals, Inc.     

Prediction of phenotypes of missense mutations in human proteins from biological assemblies

Single nucleotide polymorphisms (SNPs) are the most frequent variation in the human genome. Nonsynonymous SNPs that lead to missense mutations can be neutral or deleterious, and several computational methods have been presented that predict the phenotype of human missense mutations. These methods use sequence‐based and structure‐based features in various combinations, relying on different statistical distributions of these features for deleterious and neutral mutations. One structure‐based feature that has not been studied significantly is the accessible surface area within biologically relevant oligomeric assemblies. These assemblies are different from the crystallographic asymmetric unit for more than half of X‐ray crystal structures. We find that mutations in the core of proteins or in the interfaces in biological assemblies are significantly more likely to be disease‐associated than those on the surface of the biological assemblies. For structures with more than one protein in the biological assembly (whether the same sequence or different), we find the accessible surface area from biological assemblies provides a statistically significant improvement in prediction over the accessible surface area of monomers from protein crystal structures (P = 6e‐5). When adding this information to sequence‐based features such as the difference between wildtype and mutant position‐specific profile scores, the improvement from biological assemblies is statistically significant but much smaller (P = 0.018). Combining this information with sequence‐based features in a support vector machine leads to 82% accuracy on a balanced dataset of 50% disease‐associated mutations from SwissVar and 50% neutral mutations from human/primate sequence differences in orthologous proteins. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A superposition free method for protein conformational ensemble analyses and local clustering based on a differential geometry representation of backbone

Here a differential geometry (DG) representation of protein backbone is explored on the analyses of protein conformational ensembles. The protein backbone is described by curvature, κ, and torsion, τ, values per residue and we propose 1) a new dissimilarity and protein flexibility measurement and 2) a local conformational clustering method. The methods were applied to Ubiquitin and c-Myb-KIX protein conformational ensembles and results show that κ\τ metric space allows to properly judge protein flexibility by avoiding the superposition problem. The d_max measurement presents equally good or superior results when compared to RMSF, especially for the intrinsically unstructured protein. The clustering method is unique as it relates protein global to local dynamics by providing a global clustering solutions per residue. The methods proposed can be especially useful to the analyses of highly flexible proteins. The software written for the analyses presented here is available at https://github.com/AMarinhoSN/FleXgeo for academic usage only.     

Conformational Dynamics of Light-Harvesting Complex II in a Native Membrane Environment

Photosynthetic light-harvesting complexes (LHCs) of higher plants, moss, and green algae can undergo dynamic conformational transitions, which have been correlated to their ability to adapt to fluctuations in the light environment. Herein, we demonstrate the application of solid-state NMR spectroscopy on native, heterogeneous thylakoid membranes of Chlamydomonas reinhardtii (Cr) and on Cr light-harvesting complex II (LHCII) in thylakoid lipid bilayers to detect LHCII conformational dynamics in its native membrane environment. We show that membrane-reconstituted LHCII contains selective sites that undergo fast, large-amplitude motions, including the phytol tails of two chlorophylls. Protein plasticity is also observed in the N-terminal stromal loop and in protein fragments facing the lumen, involving sites that stabilize the xanthophyll-cycle carotenoid violaxanthin and the two luteins. The results report on the intrinsic flexibility of LHCII pigment-protein complexes in a membrane environment, revealing putative sites for conformational switching. In thylakoid membranes, fast dynamics of protein and pigment sites is significantly reduced, which suggests that in their native organelle membranes, LHCII complexes are locked in specific conformational states.     

Exploration of human serum albumin binding sites by docking and molecular dynamics flexible ligand–protein interactions

Five‐nanosecond molecular dynamics (MD) simulations were performed on human serum albumin (HSA) to study the conformational features of its primary ligand binding sites (I and II). Additionally, 11 HSA snapshots were extracted every 0.5 ns to explore the binding affinity (K_d) of 94 known HSA binding drugs using a blind docking procedure. MD simulations indicate that there is considerable flexibility for the protein, including the known sites I and II. Movements at HSA sites I and II were evidenced by structural analyses and docking simulations. The latter enabled the study and analysis of the HSA–ligand interactions of warfarin and ketoprofen (ligands binding to sites I and II, respectively) in greater detail. Our results indicate that the free energy values by docking (K_d observed) depend upon the conformations of both HSA and the ligand. The 94 HSA–ligand binding K_d values, obtained by the docking procedure, were subjected to a quantitative structure‐activity relationship (QSAR) study by multiple regression analysis. The best correlation between the observed and QSAR theoretical (K_d predicted) data was displayed at 2.5 ns. This study provides evidence that HSA binding sites I and II interact specifically with a variety of compounds through conformational adjustments of the protein structure in conjunction with ligand conformational adaptation to these sites. These results serve to explain the high ligand‐promiscuity of HSA. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 161–170, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Constrained geometric dynamics of the Fenna–Matthews–Olson complex: the role of correlated motion in reducing uncertainty in excitation energy transfer

The trimeric Fenna–Mathews–Olson (FMO) complex of green sulphur bacteria is a well-studied example of a photosynthetic pigment–protein complex, in which the electronic properties of the pigments are modified by the protein environment to promote efficient excitonic energy transfer from antenna complexes to the reaction centres. By a range of simulation methods, many of the electronic properties of the FMO complex can be extracted from knowledge of the static crystal structure. However, the recent observation and analysis of long-lasting quantum dynamics in the FMO complex point to protein dynamics as a key factor in protecting and generating quantum coherence under laboratory conditions. While fast inter- and intra-molecular vibrations have been investigated extensively, the slow, conformational dynamics which effectively determine the optical inhomogeneous broadening of experimental ensembles has received less attention. The following study employs constrained geometric dynamics to study the flexibility in the protein network by efficiently generating the accessible conformational states from the published crystal structure. Statistical and principle component analyses reveal highly correlated low frequency motions between functionally relevant elements, including strong correlations between pigments that are excitonically coupled. Our analysis reveals a hierarchy of structural interactions which enforce these correlated motions, from the level of monomer-monomer interfaces right down to the α-helices, β-sheets and pigments. In addition to inducing strong spatial correlations across the conformational ensemble, we find that the overall rigidity of the FMO complex is exceptionally high. We suggest that these observations support the idea of highly correlated inhomogeneous disorder of the electronic excited states, which is further supported by the remarkably low variance (typically <5 %) of the excitonic couplings of the conformational ensemble.     

A multiscale approach to predicting affinity changes in protein–protein interfaces

Substitution mutations in protein–protein interfaces can have a substantial effect on binding, which has consequences in basic and applied biomedical research. Experimental expression, purification, and affinity determination of protein complexes is an expensive and time‐consuming means of evaluating the effect of mutations, making a fast and accurate in silico method highly desirable. When the structure of the wild‐type complex is known, it is possible to economically evaluate the effect of point mutations with knowledge based potentials, which do not model backbone flexibility, but these have been validated only for single mutants. Substitution mutations tend to induce local conformational rearrangements only. Accordingly, ZEMu (Zone Equilibration of Mutants) flexibilizes only a small region around the site of mutation, then computes its dynamics under a physics‐based force field. We validate with 1254 experimental mutants (with 1–15 simultaneous substitutions) in a wide variety of different protein environments (65 protein complexes), and obtain a significant improvement in the accuracy of predicted ΔΔG. Proteins 2014; 82:2681–2690. © 2014 Wiley Periodicals, Inc.     

A hybrid NMR/SAXS‐based approach for discriminating oligomeric protein interfaces using Rosetta

Oligomeric proteins are important targets for structure determination in solution. While in most cases the fold of individual subunits can be determined experimentally, or predicted by homology‐based methods, protein–protein interfaces are challenging to determine de novo using conventional NMR structure determination protocols. Here we focus on a member of the bet‐V1 superfamily, Aha1 from Colwellia psychrerythraea. This family displays a broad range of crystallographic interfaces none of which can be reconciled with the NMR and SAXS data collected for Aha1. Unlike conventional methods relying on a dense network of experimental restraints, the sparse data are used to limit conformational search during optimization of a physically realistic energy function. This work highlights a new approach for studying minor conformational changes due to structural plasticity within a single dimeric interface in solution. Proteins 2015; 83:309–317. © 2014 Wiley Periodicals, Inc.     

Flap flexibility amongst plasmepsins I,II, III,IV, and V: Sequence,structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I – V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap‐structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap‐structure in plasmepsins is made up of a flip tip in the N‐terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C‐terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap‐structure in apo plasmepsins. Namely, the distance, d₁, between the flap tip and the flexible region; the dihedral angle, ?, to account for the twisting motion; and the TriCα angle, θ₁. Simulations have shown that as the flap‐structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP. Proteins 2015; 83:1693–1705. © 2015 Wiley Periodicals, Inc.     

Investigating the inter‐subunit/subdomain interactions and motions relevant to disease mutations in the N‐terminal domain of ryanodine receptors by molecular dynamics simulation

The ryanodine receptors (RyR) are essential to calcium signaling in striated muscles, and numerous disease mutations have been identified in two RyR isoforms, RyR1 in skeletal muscle and RyR2 in cardiac muscle. A deep understanding of the activation/regulation mechanisms of RyRs has been hampered by the shortage of high‐resolution structures and dynamic information for this giant tetrameric complex in different functional states. Toward elucidating the molecular mechanisms of disease mutations in RyRs, we performed molecular dynamics simulation of the N‐terminal domain (NTD) which is not only the best‐resolved structural component of RyRs, but also a hotspot of disease mutations. First, we simulated the tetrameric NTD of wild‐type RyR1 and three disease mutants (K155E, R157Q, and R164Q) that perturb the inter‐subunit interfaces. Our simulations identified a dynamic network of salt bridges involving charged residues at the inter‐subunit/subdomain interfaces and disease‐mutation sites. By perturbing this key network, the above three mutations result in greater flexibility with the highest inter‐subunit opening probability for R157Q. Next, we simulated the monomeric NTD of RyR2 in the presence or absence of a central Cl^– anion which is known to stabilize the interfaces between the three NTD subdomains (A, B, and C). We found that the loss of Cl^– restructures the salt‐bridge network near the Cl^–‐binding site, leading to rotations of subdomain A/B relative to subdomain C and enhanced mobility between the subdomains. This finding supports a mechanism for disease mutations in the NTD of RyR2 via perturbation of the Cl^– binding. The rich structural and dynamic information gained from this study will guide future mutational and functional studies of the NTD of RyRs. Proteins 2017; 85:1633–1644. © 2017 Wiley Periodicals, Inc.     

Mobility in the structure of E.coli recQ helicase upon substrate binding as seen from molecular dynamics simulations

RecQ helicases feature multiple domains in their structure, of which the helicase domain, the RecQ-Ct domain and the HRDC domains are well conserved among the SF2 helicases. The helicase domain and the RecQ-Ct domain constitute the catalytic core of the enzyme. The domain interfaces are the DNA binding sites which display significant conformational changes in our molecular dynamics simulation studies. The preferred conformational states of the DNA bound and unbound forms of RecQ appear to be quite different from each other. DNA binding induces inter-domain flexibility leading to hinge mobility between the domains. The divergence in the dynamics of the two structures is caused by changes in the interactions at the domain interface, which seems to propagate along the whole protein structure. This could be essential in ssDNA binding after strand separation, as well as aiding translocation of the RecQ protein like an inch-worm.     

Solution structure and interface‐driven self‐assembly of NC2, a new member of the Class II hydrophobin proteins

Hydrophobins are fungal proteins that self‐assemble spontaneously to form amphipathic monolayers at hydrophobic:hydrophilic interfaces. Hydrophobin assemblies facilitate fungal transitions between wet and dry environments and interactions with plant and animal hosts. NC2 is a previously uncharacterized hydrophobin from Neurospora crassa. It is a highly surface active protein and is able to form protein layers on a water:air interface that stabilize air bubbles. On a hydrophobic substrate, NC2 forms layers consisting of an ordered network of protein molecules, which dramatically decrease the water contact angle. The solution structure and dynamics of NC2 have been determined using nuclear magnetic resonance spectroscopy. The structure of this protein displays the same core fold as observed in other hydrophobin structures determined to date, including the Class II hydrophobins HFBI and HFBII from Trichoderma reesei, but certain features illuminate the structural differences between Classes I and II hydrophobins and also highlight the variations between structures of Class II hydrophobin family members. The unique properties of hydrophobins have attracted much attention for biotechnology applications. The insights obtained through determining the structure, biophysical properties and assembly characteristics of NC2 will facilitate the development of hydrophobin‐based applications. Proteins 2014; 82:990–1003. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulations of a K+ channel blocker: Tc1 toxin from Tityus cambridgei

Toxins that block voltage-gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel-blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.     

Specificity of a protein–protein interface: Local dynamics direct substrate recognition of effector caspases

Proteases are prototypes of multispecific protein–protein interfaces. Proteases recognize and cleave protein and peptide substrates at a well‐defined position in a substrate binding groove and a plethora of experimental techniques provide insights into their substrate recognition. We investigate the caspase family of cysteine proteases playing a key role in programmed cell death and inflammation, turning caspases into interesting drug targets. Specific ligand binding to one particular caspase is difficult to achieve, as substrate specificities of caspase isoforms are highly similar. In an effort to rationalize substrate specificity of two closely related caspases, we investigate the substrate promiscuity of the effector Caspases 3 and 7 by data mining (cleavage entropy) and by molecular dynamics simulations. We find a strong correlation between binding site rigidity and substrate readout for individual caspase subpockets explaining more stringent substrate readout of Caspase 7 via its narrower conformational space. Caspase 3 subpockets S3 and S4 show elevated local flexibility explaining the more unspecific substrate readout of that isoform in comparison to Caspase 7. We show by in silico exchange mutations in the S3 pocket of the proteases that a proline residue in Caspase 7 contributes to the narrowed conformational space of the binding site. These findings explain the substrate specificities of caspases via a mechanism of conformational selection and highlight the crucial importance of binding site local dynamics in substrate recognition of proteases. Proteins 2014; 82:546–555. © 2013 Wiley Periodicals, Inc.