Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

Arginine kinase: joint crystallographic and NMR RDC analyses link substrate-associated motions to intrinsic flexibility

The phosphagen kinase family, including creatine and arginine kinases (AKs), catalyzes the reversible transfer of a “high-energy” phosphate between ATP and a phosphoguanidino substrate. They have become a model for the study of both substrate-induced conformational change and intrinsic protein dynamics. Prior crystallographic studies indicated large substrate-induced domain rotations, but differences among a recent set of AK structures were interpreted as a plastic deformation. Here, the structure of Limulus substrate-free AK is refined against high-resolution crystallographic data and compared quantitatively with NMR chemical shifts and residual dipolar couplings (RDCs). This demonstrates the feasibility of this type of RDC analysis of proteins that are large by NMR standards (42 kDa) and illuminates the solution structure, free from crystal-packing constraints. Detailed comparison of the 1.7 Å resolution substrate-free crystal structure against the 1.7 Å transition-state analog complex shows large substrate-induced domain motions that can be broken down into movements of smaller quasi-rigid bodies. The solution-state structure of substrate-free AK is most consistent with an equilibrium of substrate-free and substrate-bound structures, with the substrate-free form dominating, but with varying displacements of the quasi-rigid groups. Rigid-group rotations evident from the crystal structures are about axes previously associated with intrinsic millisecond dynamics using NMR relaxation dispersion. Thus, “substrate-induced” motions are along modes that are intrinsically flexible in the substrate-free enzyme and likely involve some degree of conformational selection.     

Characterization of protein dynamics from residual dipolar couplings using the three dimensional Gaussian axial fluctuation model

Residual dipolar couplings are potentially very powerful probes of slower protein motions, providing access to dynamic events occurring on functionally important timescales up to the millisecond. One recent approach uses the three dimensional Gaussian Axial Fluctuation model (3D GAF) to determine the major directional modes and associated amplitudes of motions along the peptide chain. In this study we have used standard and accelerated molecular dynamics simulations to determine the accuracy of 3D GAF-based approaches in characterizing the nature and extent of local molecular motions. We compare modes determined directly from the trajectories with motional parameterization derived from RDCs simulated from the same trajectories. Three approaches are tested, that either suppose a known three-dimensional structure, simultaneously determine backbone structure and dynamics, or determine dynamic modes in the absence of a structural model. The results demonstrate the robustness of the 3D GAF analysis even in the presence of large-scale motions, and illustrate the remarkably quantitative nature of the extracted amplitudes. These observations suggest that the approach can be generally used for the study of functionally interesting biomolecular motions.     

Effects of cyclization on conformational dynamics and binding properties of Lys48-linked di-ubiquitin

In solution, Lys48-linked di-ubiquitin exists in dynamic equilibrium between closed and open conformations. To understand the effect of interdomain motion in polyubiquitin chains on their ability to bind ligands, we cyclized di-ubiquitin by cross-linking the free C terminus of the proximal ubiquitin with the side chain of residue 48 in the distal ubiquitin, using a chemical cross-linker, 1,6-Hexane-bis-vinylsulfone. Our NMR studies confirm that the cyclization affects conformational dynamics in di-ubiquitin by restricting opening of the interface and shifting the conformational equilibrium toward closed conformations. The cyclization, however, did not rigidly lock di-ubiquitin in a single closed conformation: The chain undergoes slow exchange between at least two closed conformations, characterized by interdomain contacts involving the same hydrophobic patch residues (Leu8-Ile44-Val70) as in the uncyclized di-ubiquitin. Lowering the pH changes the relative populations of these conformations, but in contrast with the uncyclized di-ubiquitin, does not lead to opening of the interface. This restriction of domain motions inhibits direct access of protein molecules to the hydrophobic patch residues located at the very center of the interdomain interface in di-ubiquitin, although the residual motions are sufficient to allow access of small molecules to the interface. This renders di-ubiquitin unable to bind protein molecules (e.g., UBA2 domain) in the normal manner, and thus could interfere with Ub(2) recognition by various downstream effectors. These results emphasize the importance of the opening/closing domain motions for the recognition and function of di-ubiquitin and possibly longer polyubiquitin chains.     

A new approach for applying residual dipolar couplings as restraints in structure elucidation

Residual dipolar couplings are useful global structural restraints. The dipolar couplings define the orientation of a vector with respect to the alignment tensor. Although the size of the alignment tensor can be derived from the distribution of the experimental dipolar couplings, its orientation with respect to the coordinate system of the molecule is unknown at the beginning of structure determination. This causes convergence problems in the simulated annealing process. We therefore propose a protocol that translates dipolar couplings into intervector projection angles, which are independent of the orientation of the alignment tensor with respect to the molecule. These restraints can be used during the whole simulated annealing protocol.     

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.     

A Physical Picture of Protein Dynamics and Conformational Changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an α-helix. Molecules performing these types of motion are in the “flexible state”. This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.     

Function of interdomain α-helix in human brain hexokinase: covalent linkage and catalytic regulation between N- and C-terminal halves

Human brain relies on a steady supply of glucose as the source of fuel, and type I hexokinase is the major isozyme governing the introduction of glucose to glycolysis in the brain. One unique regulatory property associated with type I isozyme is the alleviation of product inhibition by inorganic phosphate which binds to the N-terminal half, and the conformational change induced by inorganic phosphate must be propagated to the active site in the C-terminal half. With a single interdomain α-helix as the only covalent connection between the N- and C-terminal halves, the question arises as what role the interdomain α-helix plays at the interdomain signal transduction. Two mutants were constructed in an attempt to answer this question. The first mutant, A464P/E465G, with a helix breaker embedded in the interdomain α-helix had a smaller magnitude of phosphate alleviation than the wild type. The second mutant, with an insertion of seven additional residues between Gln 466 and His 467, had this phosphate relief property further diminished. Neither mutant showed dramatic changes nor the other kinetic properties. It is speculated that the interdomain α-helix is important for keeping the proper non-covalent contact so that transmission of the conformational changes across the N- and C-terminal half boundary can be achieved.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

Conformational dynamics of supramolecular protein assemblies

Supramolecular protein assemblies including molecular motors, cytoskeletal filaments, chaperones, and ribosomes play a central role in a broad array of cellular functions ranging from cell division and motility to RNA and protein synthesis and folding. Single-particle reconstructions of such assemblies have been growing rapidly in recent years, providing increasingly high resolution structural information under native conditions. While the static structure of these assemblies provides essential insight into their mechanism of biological function, their dynamical motions provide additional important information that cannot be inferred from structure alone. Here we present an unsupervised computational framework for the analysis of high molecular weight protein assemblies and use it to analyze the conformational dynamics of structures deposited in the Electron Microscopy Data Bank. Protein assemblies are modeled using a recently introduced coarse-grained modeling framework based on the finite element method, which is used to compute equilibrium thermal fluctuations, elastic strain energy distributions associated with specific conformational transitions, and dynamical correlations in distant molecular domains. Results are presented in detail for the ribosome-bound termination factor RF2 from Escherichia coli, the nuclear pore complex from Dictyostelium discoideum, and the chaperonin GroEL from E. coli. Elastic strain energy distributions reveal hinge-regions associated with specific conformational change pathways, and correlations in collective molecular motions reveal dynamical coupling between distant molecular domains that suggest new, as well as confirm existing, allosteric mechanisms. Results are publically available for use in further investigation and interpretation of biological function including cooperative transitions, allosteric communication, and molecular mechanics, as well as in further classification and refinement of electron microscopy based structures.     

Detection of initiation sites in protein folding of the four helix bundle ACBP by chemical shift analysis

A simple alternative method for obtaining "random coil" chemical shifts by intrinsic referencing using the protein's own peptide sequence is presented. These intrinsic random coil backbone shifts were then used to calculate secondary chemical shifts, that provide important information on the residual secondary structure elements in the acid-denatured state of an acyl-coenzyme A binding protein. This method reveals a clear correlation between the carbon secondary chemical shifts and the amide secondary chemical shifts 3-5 residues away in the primary sequence. These findings strongly suggest transient formation of short helix-like segments, and identify unique sequence segments important for protein folding.     

Some Surprising Implications of NMR-directed Simulations of Substrate Recognition and Binding by Cytochrome P450cam (CYP101A1)

Cytochrome P450_cam (CYP101A1) catalyzes the stereospecific 5-exo hydroxylation of d-camphor by molecular oxygen. Previously, residual dipolar couplings measured for backbone amide ¹H–¹⁵N correlations in both substrate-free and bound forms of CYP101A1 were used as restraints in soft annealing molecular dynamic simulations in order to identify average conformations of the enzyme with and without substrate bound. Multiple substrate-dependent conformational changes remote from the enzyme active site were identified, and site-directed mutagenesis and activity assays confirmed the importance of these changes in substrate recognition. The current work makes use of perturbation response scanning (PRS) and umbrella sampling molecular dynamic of the residual dipolar coupling-derived CYP101A1 structures to probe the roles of remote structural features in enforcing the regio- and stereospecific nature of the hydroxylation reaction catalyzed by CYP101A1. An improper dihedral angle Ψ was defined and used to maintain substrate orientation in the CYP101A1 active site, and it was observed that different values of Ψ result in different PRS response maps. Umbrella sampling methods show that the free energy of the system is sensitive to Ψ, and bound substrate forms an important mechanical link in the transmission of mechanical coupling through the enzyme structure. Finally, a qualitative approach to interpreting PRS maps in terms of the roles of secondary structural features is proposed.     

Insight into the redox partner interaction mechanism in cytochrome P450BM‐3 using molecular dynamics simulations

Flavocytochrome P450BM‐3 is a soluble bacterial reductase composed of two flavin (FAD/FMN) and one HEME domains. In this article, we have performed molecular dynamics simulations on both the isolated FMN and HEME domains and their crystallographic complex, with the aim to study their binding modes and to garner insight into the interdomain electron transfer (ET) mechanism. The results evidenced an interdomain conformational rearrangement that reduces the average distance between the FMN and HEME cofactors from 1.81 nm, in the crystal structure, to an average value of 1.41 ± 0.09 nm along the simulation. This modification is in agreement with previously proposed hypotheses suggesting that the crystallographic FMN/HEME complex is not in the optimal arrangement for favorable ET rate under physiological conditions. The calculation of the transfer rate along the simulation, using the Pathways Path method, demonstrated the occurrence of seven ET pathways between the two redox centers, with three of them providing ET rates (K_ET) comparable with the experimental one. The sampled ET pathways comprise the amino acids N319, L322, F390, K391, P392, F393, A399, C400, and Q403 of the HEME domain and M490 of the FMN domain. The values of K_ET closer to the experiment were found along the pathways FMN(C7) → F390 → K391 → P392 → HEME(Fe) and FMN(C8) → M490 → F393 → HEME(Fe). Finally, the analysis of the collective modes of the protein complex evidences a clear correlation of the first two essential modes with the activation of the most effective ET pathways along the trajectory. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 197–209, 2014.     

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein “breathing”, that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.     

Energy landscape of a native protein: Jumping-among-minima model

We have investigated energy landscape of human lysozyme in its native state by using principal component analysis and a model, jumping-among-minima (JAM) model. These analyses are applied to 1 nsec molecular dynamics trajectory of the protein in water. An assumption embodied in the JAM model allows us to divide protein motions into intra-substate and inter-substate motions. By examining intra-substate motions, it is shown that energy surfaces of individual conformational substates are nearly harmonic and mutually similar. As a result of principal component analysis and JAM model analysis, protein motions are shown to consist of three types of collective modes, multiply hierarchical modes, singly hierarchical modes, and harmonic modes. Multiply hierarchical modes, the number of which accounts only for 0.5% of all modes, dominate contributions to total mean-square atomic fluctuation. Inter-substate motions are observed only in a small-dimensional subspace spanned by the axes of multiplyhierarchical and singly hierarchical modes. Inter-substate motions have two notable time components: faster component seen within 200 psec and slower component. The former involves transitions among the conformational substates of the low-level hierarchy, whereas the latter involves transitions of the higher level substates observed along the first four multiply hierarchical modes. We also discuss dependence of the subspace, which contains conformational substates, on time duration of simulation. Proteins 33:496–517, 1998. © 1998 Wiley-Liss, Inc.     

Segment swapping aided the evolution of enzyme function: The case of uroporphyrinogen III synthase

In an earlier study, we showed that two‐domain segment‐swapped proteins can evolve by domain swapping and fusion, resulting in a protein with two linkers connecting its domains. We proposed that a potential evolutionary advantage of this topology may be the restriction of interdomain motions, which may facilitate domain closure by a hinge‐like movement, crucial for the function of many enzymes. Here, we test this hypothesis computationally on uroporphyrinogen III synthase, a two‐domain segment‐swapped enzyme essential in porphyrin metabolism. To compare the interdomain flexibility between the wild‐type, segment‐swapped enzyme (having two interdomain linkers) and circular permutants of the same enzyme having only one interdomain linker, we performed geometric and molecular dynamics simulations for these species in their ligand‐free and ligand‐bound forms. We find that in the ligand‐free form, interdomain motions in the wild‐type enzyme are significantly more restricted than they would be with only one interdomain linker, while the flexibility difference is negligible in the ligand‐bound form. We also estimated the entropy costs of ligand binding associated with the interdomain motions, and find that the change in domain connectivity due to segment swapping results in a reduction of this entropy cost, corresponding to ～20% of the total ligand binding free energy. In addition, the restriction of interdomain motions may also help the functional domain‐closure motion required for catalysis. This suggests that the evolution of the segment‐swapped topology facilitated the evolution of enzyme function for this protein by influencing its dynamic properties. Proteins 2016; 85:46–53. © 2016 Wiley Periodicals, Inc.