Designing Molecular Dynamics Simulations to Shift Populations of the Conformational States of Calmodulin

We elucidate the mechanisms that lead to population shifts in the conformational states of calcium-loaded calmodulin (Ca²⁺-CaM). We design extensive molecular dynamics simulations to classify the effects that are responsible for adopting occupied conformations available in the ensemble of NMR structures. Electrostatic interactions amongst the different regions of the protein and with its vicinal water are herein mediated by lowering the ionic strength or the pH. Amino acid E31, which is one of the few charged residues whose ionization state is highly sensitive to pH differences in the physiological range, proves to be distinctive in its control of population shifts. E31A mutation at low ionic strength results in a distinct change from an extended to a compact Ca²⁺-CaM conformation within tens of nanoseconds, that otherwise occur on the time scales of microseconds. The kinked linker found in this particular compact form is observed in many of the target-bound forms of Ca²⁺-CaM, increasing the binding affinity. This mutation is unique in controlling C-lobe dynamics by affecting the fluctuations between the EF-hand motif helices. We also monitor the effect of the ionic strength on the conformational multiplicity of Ca²⁺-CaM. By lowering the ionic strength, the tendency of nonspecific anions in water to accumulate near the protein surface increases, especially in the vicinity of the linker. The change in the distribution of ions in the vicinal layer of water allows N- and C- lobes to span a wide variety of relative orientations that are otherwise not observed at physiological ionic strength. E31 protonation restores the conformations associated with physiological environmental conditions even at low ionic strength.     

Conformational substates of calmodulin revealed by single-pair fluorescence resonance energy transfer: influence of solution conditions and oxidative modification

A calmodulin (CaM) mutant (T34,110C-CaM) doubly labeled with fluorescence probes AlexaFluor 488 and Texas Red in opposing domains (CaM-DA) has been used to examine conformational heterogeneity in CaM by single-pair fluorescence resonance energy transfer (spFRET). Burst-integrated FRET efficiencies of freely diffusing CaM-DA single molecules yielded distributions of distance between domains of CaM-DA. We recently reported distinct conformational substates of Ca(2+)-CaM-DA and apoCaM-DA, with peaks in the distance distributions centered at approximately 28 A, 34-38 A, and 55 A [Slaughter et al. (2004) J. Phys. Chem. B 108, 10388-10397]. In the present study, shifts in the amplitudes and center distances of the conformational substates were detected with variation in solution conditions. The amplitude of an extended conformation was observed to change as a function of Ca(2+) over a free Ca(2+) range that is consistent with binding to the high affinity, C-terminal Ca(2+) binding sites, suggesting the existence of communication between lobes of CaM. Lowering pH shifted the relative amplitudes of the conformations, with a marked increase in the presence of the compact conformations and an almost complete absence of the extended conformation. In addition, the single-molecule distance distribution of apoCaM-DA at reduced ionic strength was shifted to longer distance and showed evidence of an increase in conformational heterogeneity relative to apoCaM-DA at physiological ionic strength. Oxidation of methionine residues in CaM-DA produced a substantial increase in the amplitude of the extended conformation relative to the more compact conformation. The results are considered in light of a hypothesis that suggests that electrostatic interactions between charged amino acid side chains play an important role in determining the most stable CaM conformation under varying solution conditions.     

Calmodulin conformational changes in the activation of protein kinases

The conformation of Ca2+/calmodulin changes from extended when free in solution to compact when bound in peptide complexes. The extent and kinetics of calmodulin compaction in association with Ca2+/calmodulin-dependent protein kinases (CaMKs), as well as target peptides, were investigated by fluorescence, resonance energy transfer and stopped-flow kinetics. Compaction of Ca2+/calmodulin labelled with resonance energy-transfer probes in association with target peptides was rapid (>350 s(-1)). With the target enzymes smooth-muscle myosin light-chain kinase, CaMKIV and CaMKII, the rates of calmodulin compaction were one-two orders of magnitude lower compared with those of the peptides and in the case of alphaCaMKII, ATP binding and Thr(286) auto-phosphorylation were required for calmodulin compaction. In the absence of nucleotides, Ca2+/calmodulin bound to alphaCaMKII in extended conformations, initially probably attached by one lobe only. Kinetic data suggest that in the activation process of Ca2+/calmodulin-dependent protein kinases, productive as well as unproductive complexes are formed. The formation of productive complexes with Ca2+/calmodulin thus may determine the rate of activation.     

An activation switch in the rhodopsin family of G protein-coupled receptors: the thyrotropin receptor

We aimed at understanding molecular events involved in the activation of a member of the G protein-coupled receptor family, the thyrotropin receptor. We have focused on the transmembrane region and in particular on a network of polar interactions between highly conserved residues. Using molecular dynamics simulations and site-directed mutagenesis techniques we have identified residue Asn-7.49, of the NPxxY motif of TM 7, as a molecular switch in the mechanism of thyrotropin receptor (TSHr) activation. Asn-7.49 appears to adopt two different conformations in the inactive and active states. These two states are characterized by specific interactions between this Asn and polar residues in the transmembrane domain. The inactive gauche+ conformation is maintained by interactions with residues Thr-6.43 and Asp-6.44. Mutation of these residues into Ala increases the constitutive activity of the receptor by factors of approximately 14 and approximately 10 relative to wild type TSHr, respectively. Upon receptor activation Asn-7.49 adopts the trans conformation to interact with Asp-2.50 and a putatively charged residue that remains to be identified. In addition, the conserved Leu-2.46 of the (N/S)LxxxD motif also plays a significant role in restraining the receptor in the inactive state because the L2.46A mutation increases constitutive activity by a factor of approximately 13 relative to wild type TSHr. As residues Leu-2.46, Asp-2.50, and Asn-7.49 are strongly conserved, this molecular mechanism of TSHr activation can be extended to other members of the rhodopsin-like family of G protein-coupled receptors.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Structural determinants of the conformations of medium-sized loops in proteins

Loops are integral components of protein structures, providing links between elements of secondary structure, and in many cases contributing to catalytic and binding sites. The conformations of short loops are now understood to depend primarily on their amino acid sequences. In contrast, the structural determinants of longer loops involve hydrogen-bonding and packing interactions within the loop and with other parts of the protein. By searching solved protein structures for regions similar in main chain conformation to the antigen-binding loops in immunoglobulins, we identified medium-sized loops of similar structure in unrelated proteins, and compared the determinants of their conformations. For loops that form compact substructures the major determinant of the conformation is the formation of hydrogen bonds to inward-pointing main chain atoms. For loops that have more extended conformations, the major determinant of their structure is the packing of a particular residue or residues against the rest of the protein. The following picture emerges: Medium-sized loops of similar conformation are stabilized by similar interactions. The groups that interact with the loop have very similar spatial dispositions with respect to the loop. However, the residues that provide these interactions may arise from dissimilar parts of the protein: The conformation of the loop requires certain interactions that the protein may provide in a variety of ways.     

Dependence of protein stability on the structure of the denatured state: free energy calculations of I56V mutation in human lysozyme.

Free energy calculations were carried out to understand the effect of the I56V mutation of human lysozyme on its thermal stability. In the simulation of the denatured state, a short peptide including the mutation site in the middle is employed. To study the dependence of the stability on the denatured-state structure, five different initial conformations, native-like, extended, and three random-coil-like conformations, were examined. We found that the calculated free energy difference, DeltaDeltaGcal, depends significantly on the structure of the denatured state. When native-like structure is employed, DeltaDeltaGcal is in good agreement with the experimental free energy difference, DeltaDeltaGexp, whereas in the other four models, DeltaDeltaGcal differs sharply from DeltaDeltaGexp. It is therefore strongly suggested that the structure around the mutation site takes a native-like conformation rather than an extended or random-coil conformation. From the free energy component analysis, it has been shown that free energy components originating from Lennard-Jones and covalent interactions dominantly determine DeltaDeltaGcal. The contribution of protein-protein interactions to the nonbonded component of DeltaDeltaGcal is about the same as that from protein-water interactions. The residues that are located in a hydrophobic core (F3, L8, Y38, N39, T40, and I89) contribute significantly to the nonbonded free energy component of DeltaDeltaGcal. We also propose a general computational strategy for the study of protein stability that is equally conscious of the denatured and native states.     

Development of an extended simulated annealing method: application to the modeling of complementary determining regions of immunoglobulins.

An extended simulated annealing process (ESAP) has been developed in order to obtain an ensemble of conformations of a peptide segment from a protein fluctuating at a given temperature. The annealing process was performed with a fast Monte Carlo method using the scaled collective variables developed by Noguti and Go. The system was divided into two parts: one consists of one or more peptide segments and is flexible around the main-chain and side-chain torsional angles; the other represents the rest of the molecule and was maintained fixed at the atomic positions determined by x-ray experiments. The target function included the nonbonding atomic interactions and a distance function to anchor the N and C terminal ends of each segment to the fixed part. Three systems of complementary determining regions (CDR) of antibodies were tested and compared to x-ray data: L2 loop (7 residues) of the light chain of lambda-type Bence-Jones protein, H1 and the H2 loops (14 residues) of McPC603, and H1 and H2 loops (12 residues) of HyHEL-5. Each state of CDR conformations was characterized at room temperature by the average of their coordinates (average conformation) and the internal energy. With a limited number of annealing processes (10), starting from the extended conformation, we have obtained states with conformations close to the observed x-ray structures, from 1.1 to 1.7 A root mean square deviation (rmsd) of main-chain atoms depending on the system. These states were identical or within 0.25 A rmsd of those of lowest internal energy. For unknown CDR structures the criteria of lowest internal energies from ESAP can be used to predict hypervariable loop structures in antibodies with an accuracy comparable to other methods.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations

The goal of the current study is to utilize molecular dynamic (MD) simulations to investigate the dynamic behavior of 16S rRNA in the presence and absence of S15 and to identify the binding interactions between these two molecules. The simulations show that: (i) 16S rRNA remains in a highly folded structure when it is bound to S15; (ii) in the absence of S15, 16S rRNA significantly alters its conformation and transiently forms conformations that are similar to the bound structure that make it available for binding with S15; (iii) the unbound rRNA spends the majority of its time in extended conformations. The formation of the extended conformations is a result of the molecule reaching a lower electrostatic energy and the formation of the highly folded, crystal-like conformation is a result of achieving a lower solvation energy. In addition, our MD simulations show that 16S rRNA and S15 bind across the major groove of helix 22 (H22) via electrostatic interactions. The negatively charged phosphate groups of G658, U740, G741 and G742 bind to the positively charged S15 residues Lys7, Arg34 and Arg37. The current study provides a dynamic view of the binding of 16S rRNA with S15.     

Substrate-Specific Reorganization of the Conformational Ensemble of CSK Implicates Novel Modes of Kinase Function

Protein kinases use ATP as a phosphoryl donor for the posttranslational modification of signaling targets. It is generally thought that the binding of this nucleotide induces conformational changes leading to closed, more compact forms of the kinase domain that ideally orient active-site residues for efficient catalysis. The kinase domain is oftentimes flanked by additional ligand binding domains that up- or down-regulate catalytic function. C-terminal Src kinase (Csk) is a multidomain tyrosine kinase that is up-regulated by N-terminal SH2 and SH3 domains. Although the X-ray structure of Csk suggests the enzyme is compact, X-ray scattering studies indicate that the enzyme possesses both compact and open conformational forms in solution. Here, we investigated whether interactions with the ATP analog AMP-PNP and ADP can shift the conformational ensemble of Csk in solution using a combination of small angle x-ray scattering and molecular dynamics simulations. We find that binding of AMP-PNP shifts the ensemble towards more extended rather than more compact conformations. Binding of ADP further shifts the ensemble towards extended conformations, including highly extended conformations not adopted by the apo protein, nor by the AMP-PNP bound protein. These ensembles indicate that any compaction of the kinase domain induced by nucleotide binding does not extend to the overall multi-domain architecture. Instead, assembly of an ATP-bound kinase domain generates further extended forms of Csk that may have relevance for kinase scaffolding and Src regulation in the cell.     

Single molecule analyses of the conformational substates of calmodulin bound to the phosphorylase kinase complex

The four integral delta subunits of the phosphorylase kinase (PhK) complex are identical to calmodulin (CaM) and confer Ca(2+) sensitivity to the enzyme, but bind independently of Ca(2+). In addition to binding Ca(2+), an obligatory activator of PhK's phosphoryltransferase activity, the delta subunits transmit allosteric signals to PhK's remaining alpha, beta, and gamma subunits in activating the enzyme. Under mild conditions about 10% of the delta subunits can be exchanged for exogenous CaM. In this study, a CaM double-mutant derivatized with a fluorescent donor-acceptor pair (CaM-DA) was exchanged for delta to assess the conformational substates of PhKdelta by single molecule fluorescence resonance energy transfer (FRET) +/-Ca(2+). The exchanged subunits were determined to occupy distinct conformations, depending on the absence or presence of Ca(2+), as observed by alterations of the compact, mid-length, and extended populations of their FRET distance distributions. Specifically, the combined predominant mid-length and less common compact conformations of PhKdelta became less abundant in the presence of Ca(2+), with the delta subunits assuming more extended conformations. This behavior is in contrast to the compact forms commonly observed for many of CaM's Ca(2+)-dependent interactions with other proteins. In addition, the conformational distributions of the exchanged PhKdelta subunits were distinct from those of CaM-DA free in solution, +/-Ca(2+), as well as from exogenous CaM bound to the PhK complex as delta'. The distinction between delta and delta' is that the latter binds only in the presence of Ca(2+), but stoichiometrically and at a different location in the complex than delta.     

Effects of the Protonation State of the EEEE Motif of a Bacterial Na-channel on Conduction and Pore Structure

A distinctive feature of prokaryotic Na⁺-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na⁺ ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na⁺ ions. These higher-energy barriers preclude Na⁺ ions to permeate the selectivity filter of prokaryotic Na⁺-channels when one or more Glu177 residues are protonated.     

Effects of the Protonation State of the EEEE Motif of a Bacterial Na+-channel on Conduction and Pore Structure

A distinctive feature of prokaryotic Na⁺-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na⁺ ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na⁺ ions. These higher-energy barriers preclude Na⁺ ions to permeate the selectivity filter of prokaryotic Na⁺-channels when one or more Glu177 residues are protonated.     

Construction and molecular dynamics simulation of calmodulin in the extended and in a bent conformation.

Analysis of sequence similarity and comparison of the three-dimensional (3D) structures of troponin C and calmodulin have revealed a sequence in the central helix of calmodulin with a high probability for bending. The three amino acids known to form a bend in the N-terminal portion of troponin C are also found in the central helix of calmodulin. The modelling of a bent calmodulin structure, using the dihedral angles of the three residues in the bend of troponin C as a 3D template, results in a conformation of calmodulin where the N- and C-terminal domains are able to form contacts. Dynamics simulations starting from the X-ray structure of calmodulin and from the modelled bent calmodulin were carried out to compare flexibility and correlated movements of Ca2+ in the binding loops. Both conformations of calmodulin remained stable during the period of simulation. In the simulation of calmodulin in the extended form, the motions of the Ca2+ atoms in the two domains (Ca2+1 and Ca2+2 in one domain, and Ca2+3 and Ca2+4 in the other) are correlated. In the simulation of the bent form, an additional correlation between the Ca atoms in the two different domains is observed. The results are compatible with the occurrence of a bent conformation of calmodulin in the presence of targets, and with increased Ca2+ affinity and cooperativity of the Ca(2+)-binding loops in the calmodulin-peptide complexes.     

Indolicidin,a 13-residue basic antimicrobial peptide rich in tryptophan and proline,interacts with Ca(2+)-calmodulin

Indolicidin, ILPWKWPWWPWRR-NH(2), a short 13-residue antimicrobial and cytolytic peptide characterized from bovine neutrophils, has the calmodulin-recognition 1-5-10 hydrophobic pattern (indicated by amino acids in bold), is cationic, and thereby fulfills the requirements to interact with calmodulin. Hence, we have investigated the calmodulin-binding properties of indolicidin. Indolicidin interacted with calmodulin with fairly high affinity in a Ca(2+)-dependent manner. However, when bound, the peptide did not adopt helical conformation. Indolicidin also inhibited calmodulin-stimulated phosphodiesterase activity with IC(50) values in the nanomolar range. Replacement of either the proline residues of indolicidin with alanines or tryptophan residues with phenylalanines did not affect binding to calmodulin. However, these replacements had distinctive effects on the conformations of the bound peptides. While the alanine analog of indolicidin adopted predominantly alpha-helical conformation, the phenylalanine analog remained largely unordered. Differences in the ability of these analogs to inhibit the calmodulin-stimulated phosphodiesterase activity were observed. While the alanine analog was capable of inhibiting the activity with IC(50) values comparable to that of indolicidin, the phenylalanine analog did not inhibit the activity. Our results indicate that ability to adopt amphiphilic alpha-helical structure is not a prerequisite for binding to calmodulin and also binding does not necessarily result in inhibition of calmodulin-stimulated enzyme activities.     

Local propensities and statistical potentials of backbone dihedral angles in proteins

The following three issues concerning the backbone dihedral angles of protein structures are presented. (1) How do the dihedral angles of the 20 amino acids depend on the identity and conformation of their nearest residues? (2) To what extent are the native dihedral angles determined by local (dihedral) potentials? (3) How to build a knowledge-based potential for a residue's dihedral angles, considering the identity and conformation of its nearest residues? We find that the dihedral angle distribution for a residue can significantly depend on the identity and conformation of its adjacent residues. These correlations are in sharp contrast to the Flory isolated-pair hypothesis. Statistical potentials are built for all combinations of residue triplets and depend on the dihedral angles between consecutive residues. First, a low-resolution potential is obtained, which only differentiates between the main populated basins in the dihedral angle density plots. Minimization of the dihedral potential for 125 test proteins reveals that most native alpha-helical residues (89%) and a large fraction of native beta-sheet residues (47%) adopt conformations close to their native one. For native loop residues, the percentage is 48%. It is also found that this fraction is higher for residues away from the ends of alpha or beta secondary structure elements. In addition, a higher resolution potential is built as a function of dihedral angles by a smoothing procedure and continuous functions interpolations. Monte Carlo energy minimization with this potential results in a lower fraction for native beta-sheet residues. Nevertheless, because of the higher flexibility and entropy of beta structures, they could be preferred under the influence of non-local interactions. In general, most alpha-helices and many beta-sheets are strongly determined by the local potential, while the conformations in loops and near the end of beta-sheets are more influenced by non-local interactions.