Local Conformational Dynamics Research Articles

Resolving Atomic-Level Dynamics and Interactions of High-Molecular-Weight Hyaluronic Acid by Multidimensional Solid-State NMR.

High-molecular-weight (HMW) hyaluronic acid (HA) is a highly abundant natural polysaccharide and a fundamental component of the extracellular matrix (ECM). Its size and concentration regulate tissues' macro- and microenvironments, and its upregulation is a hallmark feature of certain tumors. Yet, the conformational dynamics of HMW-HA and how it engages with the components of the ECM microenvironment remain poorly understood at the molecular level. Probing the molecular structure and dynamics of HMW polysaccharides in a hydrated, physiological-like environment is crucial and also technically challenging. Here, we deploy advanced magic-angle spinning (MAS) solid-state NMR spectroscopy in combination with isotopic enrichment to enable an in-depth study of HMW-HA to address this challenge. This approach resolves multiple coexisting HA conformations and dynamics as a function of environmental conditions. By combining 13C-labeled HA with unlabeled ECM components, we detect by MAS NMR HA-specific changes in global and local conformational dynamics as a consequence of hydration and ECM interactions. These measurements reveal atom-specific variations in the dynamics and structure of the N-acetylglucosamine moiety of HA. We discuss possible implications for interactions that stabilize the structure of HMW-HA and facilitate its recognition by HA-binding proteins. The described methods apply similarly to the studies of the molecular structure and dynamics of HA in tumor contexts and in other biological tissues as well as HMW-HA hydrogels and nanoparticles used for biomedical and/or pharmaceutical applications.

Studies of DNA breathing in exciton-coupled (iCy3)2 dimer-labeled DNA constructs by polarization-sweep single-molecule fluorescence (PS-SMF) microscopy.

Local fluctuations of the sugar-phosphate backbones and bases of DNA (a form of DNA 'breathing') play a central role in the assembly of protein-DNA complexes. We present a single-molecule fluorescence method to sensitively measure the local conformational fluctuations of exciton-coupled cyanine [(iCy3)2] dimer-labeled DNA fork constructs in which the dimer probes are placed at varying positions relative to the DNA fork junction. These systems exhibit spectroscopic signals that are sensitive to the local conformations adopted by the sugar-phosphate backbones and bases immediately surrounding the dimer probe label positions. The (iCy3)2 dimer has one symmetric (+) and one anti-symmetric (-) exciton with respective transition dipole moments oriented perpendicular to one another. We excite single molecule samples using a continuous-wave, linearly polarized laser with its polarization direction rotated at a frequency of 1 MHz. The ensuing fluorescence signal is modulated as the laser polarization alternately excites the symmetric and anti-symmetric excitons of the (iCy3)2 dimer probe. Phase-sensitive detection of the signal at the photon-counting level provides information about the distribution of local conformations and conformational dynamics. We analyze our data using a kinetic network model, which we use to parametrize the free energy surface of the system. In addition to observing DNA breathing at and near ss-dsDNA junctions, the approach can be used to study the effects of proteins that bind and function at these sites.

Isolation and conformational analysis of the Gα α-helical domain

Heterotrimeric G proteins (G proteins), composed of Gα, Gβ, and Gγ subunits, are the major downstream signaling molecules of the G protein-coupled receptors. Upon activation, Gα undergoes conformational changes both in the Ras-like domain (RD) and the α-helical domain (AHD), leading to the dissociation of Gα from Gβγ and subsequent regulation of downstream effector proteins. Gα RD mediate the most of classical functions of Gα. However, the role of Gα AHD is relatively not well elucidated despite its much higher sequence differences between Gα subtypes than those between Gα RD. Here, we isolated AHD from Gαs, Gαi1, and Gαq to provide tools for examining Gα AHD. We investigated the conformational dynamics of the isolated Gα AHD compared to those of the GDP-bound Gα. The results showed higher local conformational dynamics of Gα AHD not only at the domain interfaces but also in regions further away from the domain interfaces. This finding is consistent with the conformation of Gα AHD in the receptor-bound nucleotide-free state. Therefore, the isolated Gα AHD could provide a platform for studying the functions of Gα AHD, such as identification of the Gα AHD-binding proteins.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding.

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75±0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. II. Effect of Temperature and Internal Friction.

Internal dynamics of proteins are essential for protein folding and function. Dynamics in unfolded proteins are of particular interest since they are the basis for many cellular processes like folding, misfolding, aggregation, and amyloid formation and also determine the properties of intrinsically disordered proteins (IDPs). It is still an open question of what governs motions in unfolded proteins and whether they encounter major energy barriers. Here we use triplet-triplet energy transfer (TTET) in unfolded homopolypeptide chains and IDPs to characterize the barriers for local and long-range loop formation. The results show that the formation of short loops encounters major energy barriers with activation energies (Ea) up to 18 kJ/mol (corrected for effects of temperature on water viscosity) with very little dependence on amino acid sequence. For poly(Gly-Ser) and polySer chains the barrier decreases with increasing loop size and reaches a limiting value of 4.6 ± 0.4 kJ/mol for long and flexible chains. This observation is in accordance with the concept of internal friction encountered by chain motions due to steric effects, which is high for local motions and decreases with increasing loop size. Comparison with the results from the viscosity dependence of loop formation shows a negative correlation between Ea and the sensitivity of the reaction to solvent viscosity (α) in accordance with the Grote-Hynes theory of memory friction. The Arrhenius pre-exponential factor (A) also decreases with increasing loop size, indicating increased entropic costs for loop formation. Long-range loop formation in the investigated sequences derived from IDPs shows increased Ea and A compared with poly(Gly-Ser) and polySer chains. This increase is exclusively due to steric effects that cause additional internal friction, whereas intramolecular hydrogen bonds, dispersion forces, and charge interactions do not affect the activation parameters.

Investigating local DNA conformational trapping dynamics during DNA polymerase holoenzyme assembly and exonuclease proof-reading activity.

The sugar-phosphate backbones of DNA continuously undergo thermally-activated conformational fluctuations at equilibrium, which directly affect base-pair unstacking and hydrogen bond breaking. The influence of these motions on biomolecular function are particularly important at single-stranded (ss)—double-stranded (ds) DNA forks and primer-template (p/t) junctions. At these junctions, non-canonical structures that deviate from the Watson-Crick B-form DNA duplex may serve as recognition sites for protein binding, and the kinetics of these interactions likely play a central regulatory role for fundamental biological processes such as DNA replication and repair. The DNA polymerase/sliding clamp holoenzyme is the central protein complex that carries out the process of complementary nucleotide addition to nascent ssDNA template strands initiated at p/t junctions. In addition to playing a role in protein assembly, local DNA conformations resulting from nucleotide misincorporation regulate the exonuclease proof-reading activity of DNA polymerase. To study the details of these processive local DNA conformational dynamics, our group has developed a molecular probing method that relies on the rigid insertion of two fluorescent cyanine (Cy3) probes at opposite positions within the sugar-phosphate backbones of complementary DNA strands. Electrostatic coupling between probes results in delocalized optical properties that sensitively depend on relative probe orientations. Utilizing a series of complementary ensemble-based (i.e., fluorescence, absorbance, circular dichroism) and single-molecule spectroscopic methods, we study the local DNA backbone conformational dynamics as they relate to protein binding and activity. In this presentation, I will discuss recent studies of DNA polymerase holoenzyme assembly and exonuclease proofreading activity at p/t junctions. Specifically, results from microsecond-resolved polarization-sweep single-molecule spectroscopy (PS-SMS) demonstrate how protein-dependent conformational trapping of otherwise sparsely populated transient DNA local conformations is involved in processive holoenzyme formation.

Probing RNA Structures and Interactions Using Fluorescence Lifetime Analyses.

Structural analyses of large, complex noncoding RNAs continue to lag behind their rapid discovery and functional descriptions. Site-specifically incorporated, minimally invasive fluorescent probes such as 2-aminopurine (2AP) and pyrrolo-cytosine (PyC) have provided essential complementary information about local RNA structure, conformational dynamics, and interactions. Here I describe a protocol that benchmarks and correlates local RNA conformations with their respective fluorescence lifetimes, as a general technique that confers key advantages over fluorescence intensity-based methods. The observation that fluorescence lifetimes are more sensitive to local structures than sequence contexts suggests broad utility across diverse RNA and ribonucleoprotein systems.

Energy transfer fluctuation observed by single-molecule spectroscopy of red-shifted bacteriochlorophyll in the homodimeric photosynthetic reaction center.

The photosynthetic reaction center of heliobacteria (hRC) is a homodimeric chromoprotein responsible for light harvesting and photoelectric conversion. The fluorescence of the hRC is radiated from a bacteriochlorophyll (Bchl) g having the lowest energy level, called red-Bchl g. The homodimeric architecture of the hRC indicates that it includes two red-Bchls g arranged symmetrically in pairs. Red-Bchl g is a fluorescent probe useful for monitoring the energy transfer network in the RC. Here, we show the fluorescence polarization dependences of two red-Bchls g, individually measured with selective excitation of chlorophyll a serving as the primary electron acceptor. The two red-Bchls g exhibit almost the same polarization dependences. Based on the polarization dependence and structural data of the hRC, we propose a candidate molecule for red-Bchl g. The fluorescence spectra of single hRCs represent the spectral heterogeneity reflecting the local conformational inhomogeneity. A time series of the fluorescence spectra indicates occasional peak shifts between blue- and red-shifted states without significant changes in the fluorescence intensity. The spectral fluctuation is interpreted to be due to the local conformational dynamics around a Bchl g mediating the energy transfer, switching the terminal energy acceptor between two red-Bchls g. In conclusion, while the energy transfer network in the RC can be perturbed by microscopic dynamics, the total energy transfer efficiency, i.e., the light-harvesting function, is rather robust. The functional robustness may be due to multiple energy transfer pathways composed of many antenna pigments in the RC.

RNA modifications stabilize the tertiary structure of tRNAfMet by locally increasing conformational dynamics.

A plethora of modified nucleotides extends the chemical and conformational space for natural occurring RNAs. tRNAs constitute the class of RNAs with the highest modification rate. The extensive modification modulates their overall stability, the fidelity and efficiency of translation. However, the impact of nucleotide modifications on the local structural dynamics is not well characterized. Here we show that the incorporation of the modified nucleotides in tRNAfMet from Escherichia coli leads to an increase in the local conformational dynamics, ultimately resulting in the stabilization of the overall tertiary structure. Through analysis of the local dynamics by NMR spectroscopic methods we find that, although the overall thermal stability of the tRNA is higher for the modified molecule, the conformational fluctuations on the local level are increased in comparison to an unmodified tRNA. In consequence, the melting of individual base pairs in the unmodified tRNA is determined by high entropic penalties compared to the modified. Further, we find that the modifications lead to a stabilization of long-range interactions harmonizing the stability of the tRNA’s secondary and tertiary structure. Our results demonstrate that the increase in chemical space through introduction of modifications enables the population of otherwise inaccessible conformational substates.

Investigating the role of local DNA conformational dynamics in processive DNA polymerase holoenzyme-primer-template DNA complex formation

Single-strand-double-strand (ss-ds) DNA junctions play important roles in protein recognition and binding interactions in DNA replication and repair. The functionality of ss-dsDNA junctions is derived, at least in part, from intrinsic disorder at these sites that results in non-canonical local DNA conformations and conformational dynamics. Primer/template (P/T) DNA junctions serve as sites for the multistep assembly of the DNA polymerase/sliding clamp holoenzyme that drives DNA replication and proof-reading.

Ligand modulation of the conformational dynamics of the A2A adenosine receptor revealed by single-molecule fluorescence

G protein-coupled receptors (GPCRs) are the largest class of transmembrane proteins, making them an important target for therapeutics. Activation of these receptors is modulated by orthosteric ligands, which stabilize one or several states within a complex conformational ensemble. The intra- and inter-state dynamics, however, is not well documented. Here, we used single-molecule fluorescence to measure ligand-modulated conformational dynamics of the adenosine A2A receptor (A2AR) on nanosecond to millisecond timescales. Experiments were performed on detergent-purified A2R in either the ligand-free (apo) state, or when bound to an inverse, partial or full agonist ligand. Single-molecule Förster resonance energy transfer (smFRET) was performed on detergent-solubilized A2AR to resolve active and inactive states via the separation between transmembrane (TM) helices 4 and 6. The ligand-dependent changes of the smFRET distributions are consistent with conformational selection and with inter-state exchange lifetimes ≥ 3 ms. Local conformational dynamics around residue 2296.31 on TM6 was measured using fluorescence correlation spectroscopy (FCS), which captures dynamic quenching due to photoinduced electron transfer (PET) between a covalently-attached dye and proximal aromatic residues. Global analysis of PET-FCS data revealed fast (150–350 ns), intermediate (50–60 μs) and slow (200–300 μs) conformational dynamics in A2AR, with lifetimes and amplitudes modulated by ligands and a G-protein mimetic (mini-Gs). Most notably, the agonist binding and the coupling to mini-Gs accelerates and increases the relative contribution of the sub-microsecond phase. Molecular dynamics simulations identified three tyrosine residues (Y112, Y2887.53, and Y2907.55) as being responsible for the dynamic quenching observed by PET-FCS and revealed associated helical motions around residue 2296.31 on TM6. This study provides a quantitative description of conformational dynamics in A2AR and supports the idea that ligands bias not only GPCR conformations but also the dynamics within and between distinct conformational states of the receptor.

Using Exciton-Coupled Probe Spectroscopy to Investigate the Roles of Local DNA Conformational Dynamics in Processive DNA Polymerase Holoenzyme Formation

Using Exciton-Coupled Probe Spectroscopy to Investigate the Roles of Local DNA Conformational Dynamics in Processive DNA Polymerase Holoenzyme Formation

Effects of electrostatic interactions on global folding and local conformational dynamics of a multidomain Y-family DNA polymerase.

The Y-family DNA polymerases specialize in translesion DNA synthesis, which is essential for replicating damaged DNA. The Y-family polymerases, which are made up of four stable domains, exhibit extensive distributions of charged residues, and are responsible for the tight formation of the protein-DNA complex. However, it is still unclear how the electrostatic interactions influence the conformational dynamics of the polymerases. Here, we focus on the case of a prototype Y-family DNA polymerase, Dpo4. Using coarse-grained models including a salt-dependent electrostatic potential, we investigate the effects of the electrostatic interactions on the folding process of Dpo4. Our simulations show that strong electrostatic interactions result in a three-state folding of Dpo4, consistent with the experimental observations. This folding process exhibits low cooperativity led by low salt concentration, where the individual domains fold one by one through one single pathway. Since the refined folding order of domains in multidomain proteins can shrink the configurational space, we suggest that the electrostatic interactions facilitate the Dpo4 folding. In addition, we study the local conformational dynamics of Dpo4 in terms of fluctuation and frustration analyses. We show that the electrostatic interactions can exaggerate the local conformational properties, which are in favor of the large-scale conformational transition of Dpo4 during the functional DNA binding. Our results underline the importance of electrostatic interactions in the conformational dynamics of Dpo4 at both the global and local scale, providing useful guidance in protein engineering at the multidomain level.

Requirements for efficient ligand-gated co-transcriptional switching in designed variants of the B. subtilis pbuE adenine-responsive riboswitch in E. coli.

Riboswitches, generally located in the 5’-leader of bacterial mRNAs, direct expression via a small molecule-dependent structural switch that informs the transcriptional or translational machinery. While the structure and function of riboswitch effector-binding (aptamer) domains have been intensely studied, only recently have the requirements for efficient linkage between small molecule binding and the structural switch in the cellular and co-transcriptional context begun to be actively explored. To address this aspect of riboswitch function, we have performed a structure-guided mutagenic analysis of the B. subtilis pbuE adenine-responsive riboswitch, one of the simplest riboswitches that employs a strand displacement switching mechanism to regulate transcription. Using a cell-based fluorescent protein reporter assay to assess ligand-dependent regulatory activity in E. coli, these studies revealed previously unrecognized features of the riboswitch. Within the aptamer domain, local and long-range conformational dynamics influenced by sequences within helices have a significant effect upon efficient regulatory switching. Sequence features of the expression platform including the pre-aptamer leader sequence, a toehold helix and an RNA polymerase pause site all serve to promote strong ligand-dependent regulation. By optimizing these features, we were able to improve the performance of the B. subtilis pbuE riboswitch in E. coli from 5.6-fold induction of reporter gene expression by the wild type riboswitch to over 120-fold in the top performing designed variant. Together, these data point to sequence and structural features distributed throughout the riboswitch required to strike a balance between rates of ligand binding, transcription and secondary structural switching via a strand exchange mechanism and yield new insights into the design of artificial riboswitches.

Quantum dynamical simulations of intra-chain exciton diffusion in an oligo (para-phenylene vinylene) chain at finite temperature.

We report on quantum dynamical simulations of exciton diffusion in an oligo(para-phenylene vinylene) chain segment with 20 repeat units (OPV-20) at finite temperature, complementary to our recent study of the same system at T = 0 K [R. Binder and I. Burghardt, J. Chem. Phys. 152, 204120 (2020)]. Accurate quantum dynamical simulations are performed using the multi-layer multi-configuration time-dependent Hartree method as applied to a site-based Hamiltonian comprising 20 electronic states of Frenkel type and 460 vibrational modes, including site-local quinoid-distortion modes along with site-correlated bond-length alternation (BLA) modes, ring torsional modes, and an explicit harmonic-oscillator bath. A first-principles parameterized Frenkel-Holstein type Hamiltonian is employed, which accounts for correlations between the ring torsional modes and the anharmonically coupled BLA coordinates located at the same junction. Thermally induced fluctuations of the torsional modes are described by a stochastic mean-field approach, and their impact on the excitonic motion is characterized in terms of the exciton mean-squared displacement. A normal diffusion regime is observed under periodic boundary conditions, apart from transient localization features. Even though the polaronic exciton species are comparatively weakly bound, exciton diffusion is found to be a coherent-rather than hopping type-process, driven by the fluctuations of the soft torsional modes. Similar to the previous observations for oligothiophenes, the evolution for the most part exhibits a near-adiabatic dynamics of local exciton ground states (LEGSs) that adjust to the local conformational dynamics. However, a second mechanism, involving resonant transitions between neighboring LEGSs, gains importance at higher temperatures.

A Hydrogen-Deuterium Exchange Mass Spectrometry (HDX-MS) Platform for Investigating Peptide Biosynthetic Enzymes.

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is a powerful method for the biophysical characterization of enzyme conformational changes and enzyme-substrate interactions. Among its many benefits, HDX-MS consumes only small amounts of material, can be performed under near native conditions without the need for enzyme/substrate labeling, and can provide spatially resolved information on enzyme conformational dynamics-even for large enzymes and multiprotein complexes. The method is initiated by the dilution of the enzyme of interest into buffer prepared in D2O. This triggers the exchange of protium in peptide bond amides (N-H) with deuterium (N-D). At the desired exchange time points, reaction aliquots are quenched, the enzyme is proteolyzed into peptides, the peptides are separated by ultra-performance liquid chromatography (UPLC), and the change in mass of each peptide (due to the exchange of hydrogen for deuterium) is recorded by MS. The amount of deuterium uptake by each peptide is strongly dependent on the local hydrogen bonding environment of that peptide. Peptides present in very dynamic regions of the enzyme exchange deuterium very rapidly, whereas peptides derived from well-ordered regions undergo exchange much more slowly. In this manner, the HDX rate reports on local enzyme conformational dynamics. Perturbations to deuterium uptake levels in the presence of different ligands can then be used to map ligand binding sites, identify allosteric networks, and to understand the role of conformational dynamics in enzyme function. Here, we illustrate how we have used HDX-MS to better understand the biosynthesis of a type of peptide natural products called lanthipeptides. Lanthipeptides are genetically encoded peptides that are post-translationally modified by large, multifunctional, conformationally dynamic enzymes that are difficult to study with traditional structural biology approaches. HDX-MS provides an ideal and adaptable platform for investigating the mechanistic properties of these types of enzymes.

Electrostatic Switching Controls Channel Dynamics of the Sensor Protein VirB10 in A. tumefaciens Type IV Secretion System.

Type IV secretion systems are large nanomachines assembled across the bacterial cell envelope for effector translocation and conjugation. VirB10 traverses the inner and outer membranes, sensing cellular signals for coordinating the conformational switch for pilus biogenesis and/or secretion. Mutations uncoupling secretion from pilus biogenesis were identified in Agrobacterium tumefaciens VirB10 including a gating defect mutation G272R that made VirB10 unresponsive to intracellular ATP, causing unregulated secretion of VirE2 in a contact-independent manner. Comparative long-timescale molecular dynamics of the wild type and G272R mutant of the A. tumefaciens VirB10CTD tetradecamer reveals how the G272R mutation locks the oligomer in a rigid conformation by swapping the ionic interactions between the loops from the β-barrel close to the inner leaflet of the outer membrane. This electrostatic switching changes the allosteric communication pathway from the extracellular loop to the base of the barrel, suggesting that the local conformational dynamics in the loops can gate information across VirB10.

An Investigation of the YidC-Mediated Membrane Insertion of a Pf3 Coat Protein Using MD Simulations

YidC is a bacterial member of the YidC/Alb3/Oxa1 insertase family which mediates the insertion of newly synthesized peptides and proteins in the membrane. Proteins can be inserted into lipid bilayers via the SecYEG complex or through Sec-independent mechanisms. In this molecular dynamics (MD) study, we have developed different docking models of YidC with a nascent Pf3 coat protein in a POPE membrane. Our equilibrium MD simulations show that sequential global conformational changes occur in the transmembrane (TM) helices, cytoplasmic loops and EH1 domains of the YidC which have a prominent role at different stages of insertion mechanism. We have also performed a steered molecular dynamics (SMD) simulation to characterize these conformational changes, using our docking models as a reference. The SMD simulations reveal that a positively charged arginine residue in the hydrophilic groove establishes strong electrostatic interactions with negatively charged Pf3 coat residues over the course of insertion process. Therefore, we propose that the formation of strong interactions between the Pf3 coat protein and YidC hydrophilic groove are essential for the insertion of Pf3 into the lipid bilayer. In addition, our results from equilibrium simulations, such as Pf3 coat interactions with lipids, global, and local conformational changes in YidC domains at different conformational stages are also reappeared in SMD simulation process very similarly. This study thus characterizes the local conformational dynamics of the YidC-Pf3 interaction in the context of a Sec-independent insertion mechanism.

Comparative effects of oncogenic mutations G12C, G12V, G13D, and Q61H on local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated protein in human cancers. However, until very recently, its oncogenic mutants were viewed as undruggable. To develop inhibitors that directly target oncogenic K-Ras mutants, we need to understand both their mutant-specific and pan-mutant dynamics and conformations. Recently, we have investigated how the most frequently observed K-Ras mutation in cancer patients, G12D, changes its local dynamics and conformations (Vatansever et al., 2019). Here, we extend our analysis to study and compare the local effects of other frequently observed oncogenic mutations, G12C, G12V, G13D and Q61H. For this purpose, we have performed Molecular Dynamics (MD) simulations of each mutant when active (GTP-bound) and inactive (GDP-bound), analyzed their trajectories, and compared how each mutant changes local residue conformations, inter-protein distance distributions, local flexibility and residue pair correlated motions. Our results reveal that in the four active oncogenic mutants we have studied, the α2 helix moves closer to the C-terminal of the α3 helix. However, P-loop mutations cause α3 helix to move away from Loop7, and only G12 mutations change the local conformational state populations of the protein. Furthermore, the motions of coupled residues are mutant-specific: G12 mutations lead to new negative correlations between residue motions, while Q61H destroys them. Overall, our findings on the local conformational states and protein dynamics of oncogenic K-Ras mutants can provide insights for both mutant-selective and pan-mutant targeted inhibition efforts.

Local conformational dynamics regulating transport properties of a Cl- /H+ antiporter.

ClC-ec1 is a Cl- /H+ antiporter that exchanges Cl- and H+ ions across the membrane. Experiments have demonstrated that several mutations, including I109F, decrease the Cl- and H+ transport rates by an order of magnitude. Using reactive molecular dynamics simulations of explicit proton transport across the central region in the I109F mutant, a two-dimensional free energy profile has been constructed that is consistent with the experimental transport rates. The importance of a phenylalanine gate formed by F109 and F357 and its influence on hydration connectivity through the central proton transport pathway is revealed. This work demonstrates how seemingly subtle changes in local conformational dynamics can dictate hydration changes and thus transport properties. © 2019 Wiley Periodicals, Inc.