Thermal unfolding simulations of apo-calmodulin using leap-dynamics

Molecular Simulations Guidelines for Biological Nanomaterials: From Peptides to Membranes

In studying biological processes and focusing on the molecular mechanisms at the basis of these, molecular dynamics (MD) simulations have demonstrated to be a very useful tool for the past 50 years. This suite of computational methods calculates the time-dependent evolution of a molecular system using physics-based first principles. In this chapter, we give a brief introduction to the theory and practical use of molecular dynamics simulations, highlighting the different models and algorithms that have been developed to tackle specific problems, with a special focus on classical force fields. Some examples of how simulations have been used in the past will help the reader in discerning their power, limitations, and significance.

Uncovering the Early Stages of Domain Melting in Calmodulin with Ultrafast Temperature-Jump Infrared Spectroscopy

The signaling protein calmodulin (CaM) undergoes a well-known change in secondary structure upon binding Ca²⁺, but the structural plasticity of the Ca²⁺-free apo state is linked to CaM functionality. Variable temperature studies of apo-CaM indicate two structural transitions at 46 and 58 °C that are assigned to melting of the C- and N-terminal domains, respectively, but the molecular mechanism of domain unfolding is unknown. We report temperature-jump time-resolved infrared (IR) spectroscopy experiments designed to target the first steps in the C-terminal domain melting transition of human apo-CaM. A comparison of the nonequilibrium relaxation of apo-CaM with the more thermodynamically stable holo-CaM, with 4 equiv of Ca²⁺ bound, shows that domain melting of apo-CaM begins on microsecond time scales with α-helix destabilization. These observations enable the assignment of previously reported dynamics of CaM on hundreds of microsecond time scales to thermally activated melting, producing a complete mechanism for thermal unfolding of CaM.

Cotranslocational processing of the protein substrate calmodulin by an AAA+ unfoldase occurs via unfolding and refolding intermediates

Protein remodeling by AAA+ enzymes is central for maintaining proteostasis in a living cell. However, a detailed structural description of how this is accomplished at the level of the substrate molecules that are acted upon is lacking. Here, we combine chemical cross-linking and methyl transverse relaxation-optimized NMR spectroscopy to study, at atomic resolution, the stepwise unfolding and subsequent refolding of the two-domain substrate calmodulin by the VAT AAA+ unfoldase from Thermoplasma acidophilum By engineering intermolecular disulphide bridges between the substrate and VAT we trap the substrate at different stages of translocation, allowing structural studies throughout the translocation process. Our results show that VAT initiates substrate translocation by pulling on intrinsically unstructured N or C termini of substrate molecules without showing specificity for a particular amino acid sequence. Although the B1 domain of protein G is shown to unfold cooperatively, translocation of calmodulin leads to the formation of intermediates, and these differ on an individual domain level in a manner that depends on whether pulling is from the N or C terminus. The approach presented generates an atomic resolution picture of substrate unfolding and subsequent refolding by unfoldases that can be quite different from results obtained via in vitro denaturation experiments.

The Ca(2+) influence on calmodulin unfolding pathway: a steered molecular dynamics simulation study

The force-induced unfolding of calmodulin (CaM) was investigated at atomistic details with steered molecular dynamics. The two isolated CaM domains as well as the full-length CaM were simulated in N-C-terminal pulling scheme, and the isolated N-lobe of CaM was studied specially in two other pulling schemes to test the effect of pulling direction and compare with relevant experiments. Both Ca²⁺-loaded CaM and Ca²⁺-free CaM were considered in order to define the Ca²⁺ influence to the CaM unfolding. The results reveal that the Ca²⁺ significantly affects the stability and unfolding behaviors of both the isolated CaM domains and the full-length CaM. In Ca²⁺-loaded CaM, N-terminal domain unfolds in priori to the C-terminal domain. But in Ca²⁺-free CaM, the unfolding order changes, and C-terminal domain unfolds first. The force-extension curves of CaM unfolding indicate that the major unfolding barrier comes from conquering the interaction of two EF-hand motifs in both N- and C- terminal domains. Our results provide the atomistic-level insights in the force-induced CaM unfolding and explain the observation in recent AFM experiments.

Use of SANS and biophysical techniques to reveal subtle conformational differences between native apo-calmodulin and its unfolded states

Apo-calmodulin, a small, mainly α, soluble protein is a calcium-dependent protein activator. It is made of two N- and C-terminal domains having a sequence homology of 70%, an identical folding but different stabilities, and is thus an interesting system for unfolding studies. The use of small angle neutron scattering (SANS) and other biophysical techniques has permitted to reveal conformational difference between native and thermal denatured states of apo-calmodulin. The results show that secondary and tertiary structures of apo-calmodulin evolve in a synchronous way, indicating the absence in the unfolding pathway of molten-globule state sufficiently stable to affect transition curves. From SANS experiments, at 85 °C, apo-calmodulin adopts a polymer chain conformation with some residual local structures. After cooling down, apo-calmodulin recovers a compact state, with a secondary structure close to the native one but with a higher radius of gyration and a different tyrosine environment. In fact on a timescale of few minutes, heat denaturation of apo-calmodulin is partially reversible, but on a time scale of hours (for SANS experiments), the long exposure to heat may lead to a non-reversibility due to some chemical perturbation of the protein. In fact, from Mass Spectrometry measurements, we got evidence of dehydration and deamidation of heated apo-calmodulin.

Implicit Solvation Parameters Derived from Explicit Water Forces in Large-Scale Molecular Dynamics Simulations

Implicit solvation is a mean force approach to model solvent forces acting on a solute molecule. It is frequently used in molecular simulations to reduce the computational cost of solvent treatment. In the first instance, the free energy of solvation and the associated solvent–solute forces can be approximated by a function of the solvent-accessible surface area (SASA) of the solute and differentiated by an atom–specific solvation parameter σ_i^SASA. A procedure for the determination of values for the σ_i^SASA parameters through matching of explicit and implicit solvation forces is proposed. Using the results of Molecular Dynamics simulations of 188 topologically diverse protein structures in water and in implicit solvent, values for the σ_i^SASA parameters for atom types i of the standard amino acids in the GROMOS force field have been determined. A simplified representation based on groups of atom types σ_g^SASA was obtained via partitioning of the atom–type σ_i^SASA distributions by dynamic programming. Three groups of atom types with well separated parameter ranges were obtained, and their performance in implicit versus explicit simulations was assessed. The solvent forces are available at http://mathbio.nimr.mrc.ac.uk/wiki/Solvent_Forces.

Protein-water interactions in MD simulations: POPS/POPSCOMP solvent accessibility analysis, solvation forces and hydration sites

The effects of solvation on molecular recognition are investigated from different perspectives, ranging from methods to analyse explicit solvent dynamical behaviour at the protein surface to methods for the implicit treatment of solvent effects associated with the conformational behaviour of biomolecules. The here presented implicit solvation method is based on an analytical approximation of the Solvent Accessible Surface Area (SASA) of solute molecules, which is computationally efficient and easy to parametrise. The parametrised SASA solvation method is discussed in the light of protein design and ligand binding studies. The POPS program for the SASA computation on single molecules and complex interfaces is described in detail. Explicit solvent behaviour is described here in the form of solvent density maps at the protein surface. We highlight the usefulness of that approach in defining the organisation of specific water molecules at functional sites and in determining hydrophobicity scores for the identification of potential interaction patches.

The effect of macromolecular crowding, ionic strength and calcium binding on calmodulin dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM's target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.

Proteins are workhorses for driving biological functions inside cells. Calmodulin (CaM) is a protein that can carry cellular signals by triggered conformational changes due to calcium binding that alters target binding. Interestingly, CaM is able to bind over 300 targets. One of the challenges in characterizing CaM's ability to bind multiple targets lies in that CaM is a flexible protein and its structure is easily modulated by the physicochemical changes in its surroundings, particularly inside a complex cellular milieu. In order to determine structure-function relationships of CaM, we employed a combined approach of experiments, computer simulations and statistical physics in the investigation of the effect of calcium-binding, salt concentration, and macromolecular crowding on CaM. The results revealed unique folding energy landscapes of CaM in the absence and presence of calcium ions and the structural implications of CaM are interpreted under cell-like conditions. Further, a large conformational change in CaM in response to environmental impacts, dictates the packing of local helices that may be critical to its function of target binding and recognition among vast target selections.

Biophysical study of thermal denaturation of apo-calmodulin: dynamics of native and unfolded states

Apo-calmodulin, a small, mainly alpha, soluble protein is a calcium-dependent protein activator. This article presents a study of internal dynamics of native and thermal unfolded apo-calmodulin, using quasi-elastic neutron scattering. This technique can probe protein internal dynamics in the picosecond timescale and in the nanometer length-scale. It appears that a dynamical transition is associated with thermal denaturation of apo-calmodulin. This dynamical transition goes together with a decrease of the confinement of hydrogen atoms, a decrease of immobile protons proportion and an increase of dynamical heterogeneity. The comparison of native and unfolded states dynamics suggests that the dynamics of protein atoms is more influenced by their distance to the backbone than by their solvent exposure.

On searching in, sampling of, and dynamically moving through conformational space of biomolecular systems: A review

Methods to search for low-energy conformations, to generate a Boltzmann-weighted ensemble of configurations, or to generate classical-dynamical trajectories for molecular systems in the condensed liquid phase are briefly reviewed with an eye to application to biomolecular systems. After having chosen the degrees of freedom and method to generate molecular configurations, the efficiency of the search or sampling can be enhanced in various ways: (i) efficient calculation of the energy function and forces, (ii) application of a plethora of search enhancement techniques, (iii) use of a biasing potential energy term, and (iv) guiding the sampling using a reaction or transition pathway. The overview of the available methods should help the reader to choose the combination that is most suitable for the biomolecular system, degrees of freedom, interaction function, and molecular or thermodynamic properties of interest.

The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.

Plasticity of influenza haemagglutinin fusion peptides and their interaction with lipid bilayers

A detailed molecular dynamics study of the haemagglutinin fusion peptide (N-terminal 20 residues of the HA2 subunits) in a model bilayer has yielded useful information about the molecular interactions leading to insertion into the lipids. Simulations were performed on the native sequence, as well as a number of mutant sequences, which are either fusogenic or nonfusogenic. For the native sequence and fusogenic mutants, the N-terminal 11 residues of the fusion peptides are helical and insert with a tilt angle of approximately 30 degrees with respect to the membrane normal, in very good agreement with experimental data. The tilted insertion of the native sequence peptide leads to membrane bilayer thinning and the calculated order parameters show larger disorder of the alkyl chains. These results indicate that the lipid packing is perturbed by the fusion peptide and could be used to explain membrane fusion. For the nonfusogenic sequences investigated, it was found that most of them equilibrate parallel to the interface plane and do not adopt a tilted conformation. The presence of a charged residue at the beginning of the sequence (G1E mutant) resulted in a more difficult case, and the outcomes do not fall straightforwardly into the general picture. Sequence searches have revealed similarities of the fusion peptide of influenza haemagglutinin with peptide sequences such as segments of porin, amyloid alpha eta peptide, and a peptide from the prion sequence. These results confirm that the sequence can adopt different folds in different environments. The plasticity and the conformational dependence on the local environment could be used to better understand the function of fusion peptides.

Crystal structures of apocalmodulin and an apocalmodulin/SK potassium channel gating domain complex

Small conductance Ca2+-activated K+ channels (SK channels) are composed of the pore-forming alpha subunit and calmodulin (CaM). CaM binds to a region of the alpha subunit called the CaM binding domain (CaMBD), located intracellular and immediately C-terminal to the inner helix gate, in either the presence or absence of Ca2+. SK gating occurs when Ca2+ binds the N lobe of CaM thereby transmitting the signal to the attached inner helix gate to open. Here we present crystal structures of apoCaM and apoCaM/SK2 CaMBD complex. Several apoCaM crystal forms with multiple (12) packing environments reveal the same EF hand domain-swapped dimer providing potentially new insight into CaM regulation. The apoCaM/SK2 CaMBD structure, combined with our Ca2+/CaM/CaMBD structure suggests that Ca2+ binding induces folding and dimerization of the CaMBD, which causes large CaMBD-CaM C lobe conformational changes, including a >90 degrees rotation of the region of the CaMBD directly connected to the gate.