Conformational Fluctuations Research Articles

Heterogeneous and Allosteric Role of Surface Hydration for Protein-Ligand Binding.

Atomistic-level understanding of surface hydration mediating protein-protein interactions and ligand binding has been a challenge due to the dynamic nature of water molecules near the surface. We develop a computational method to evaluate the solvation free energy based on the density map of the first hydration shell constructed from all-atom molecular dynamics simulation and use it to examine the binding of two intrinsically disordered ligands to their target protein domain. One ligand is from the human protein, and the other is from the 1918 Spanish flu virus. We find that the viral ligand incurs a 6.9 kcal/mol lower desolvation penalty upon binding to the target, which is consistent with its stronger binding affinity. The difference arises from the spatially fragmented and nonuniform water density profiles of the first hydration shell. In particular, residues that are distal from the ligand-binding site contribute to a varying extent to the desolvation penalty, among which the "entropy hotspot" residues contribute significantly. Thus, ligand binding alters hydration on remote sites in addition to affecting the binding interface. The nonlocal effect disappears when the conformational motion of the protein is suppressed. The present results elucidate the interplay between protein conformational dynamics and surface hydration. Our approach of measuring the solvation free energy based on the water density of the first hydration shell is tolerant of the conformational fluctuation of protein, and we expect it to be applicable to investigating a broad range of biomolecular interfaces.

Fast Motions Dominate Dynamics of Intrinsically Disordered Tau Protein at High Temperatures.

Reorientational dynamics of intrinsically disordered proteins (IDPs) contain multiple motions often clustered around three motional modes: ultrafast librational motions of amide groups, fast local backbone conformational fluctuations and slow chain segmental motions. This dynamic picture is mainly based on 15 N NMR relaxation studies of IDPs at relatively low temperatures where the amide-water proton exchange rates are sufficiently small. Less is known, however, about the dynamics of IDPs at more physiological temperatures. Here, we investigate protein dynamics in a 441-residue long IDP, tau protein, in the temperature range from 0-25 °C, using 15 N NMR relaxation rates and spectral density analysis. While at these temperatures relaxation rates are still better described in terms of amide group librational motions, local backbone dynamics and chain segmental motions, the temperature-dependent trend of spectral densities suggests that the timescales of fast backbone conformational fluctuations and slower chain segmental motions might become inseparable at higher temperatures. Our data demonstrate the remarkable dynamic plasticity of this prototypical IDP and highlight the need for dynamic studies of IDPs at multiple temperatures.

Biophysical Correlates of Enhanced Immunogenicity of a Stabilized Variant of the Receptor Binding Domain of SARS-CoV-2.

The receptor binding domain (RBD) of SARS-CoV-2 is the primary target of neutralizing antibodies. We have previously reported the design and characterization of a mammalian cell expressed RBD derivative, mRBD1-3.2, that has higher thermal stability and greatly enhanced immunogenicity relative to the wild type mRBD. The protein is highly thermotolerant and immunogenic and is being explored for use in room temperature stable Covid-19 vaccine formulations. In the current study, we have investigated the folding pathway of both WT and stabilized RBD. It was found that chemical denaturation of RBD proceeds through a stable equilibrium intermediate. Thermal and chemical denaturation is reversible, as assayed by binding to the receptor ACE2. Unusually, in its native state, RBD binds to the hydrophobic probe ANS, and enhanced ANS binding is observed for the equilibrium intermediate state. Further characterization of the folding of mRBD1-3.2, both in solution and after reconstitution of lyophilized protein stored for a month at 37 °C, revealed a higher stability represented by higher Cm, faster refolding, slower unfolding, and enhanced resistance to proteolytic cleavage relative to WT. In contrast to WT RBD, the mutant showed decreased interaction with the hydrophobic moiety linoleic acid. Collectively, these data suggest that the enhanced immunogenicity results from reduced conformational fluctuations that likely enhance in vivo half-life as well as reduce the exposure of irrelevant non-neutralizing epitopes to the immune system.

Re-Balancing Replica Exchange with Solute Tempering for Sampling Dynamic Protein Conformations.

Replica exchange with solute tempering (REST) is a highly effective variant of replica exchange for enhanced sampling in explicit solvent simulations of biomolecules. By scaling the Hamiltonian for a selected "solute" region of the system, REST effectively applies tempering only to the degrees of freedom of interest but not the rest of the system ("solvent"), allowing fewer replicas for covering the same temperature range. A key consideration of REST is how the solute-solvent interactions are scaled together with the solute-solute interactions. Here, we critically evaluate the performance of the latest REST2 protocol for sampling large-scale conformation fluctuations of intrinsically disordered proteins (IDPs). The results show that REST2 promotes artificial protein conformational collapse at high effective temperatures, which seems to be a designed feature originally to promote the sampling of reversible folding of small proteins. The collapse is particularly severe with larger IDPs, leading to replica segregation in the effective temperature space and hindering effective sampling of large-scale conformational changes. We propose that the scaling of the solute-solvent interactions can be treated as free parameters in REST, which can be tuned to control the solute conformational properties (e.g., chain expansion) at different effective temperatures and achieve more effective sampling. To this end, we derive a new REST3 protocol, where the strengths of the solute-solvent van der Waals interactions are recalibrated to reproduce the levels of protein chain expansion at high effective temperatures. The efficiency of REST3 is examined using two IDPs with nontrivial local and long-range structural features, including the p53 N-terminal domain and the kinase inducible transactivation domain of transcription factor CREB. The results suggest that REST3 leads to a much more efficient temperature random walk and improved sampling efficiency, which also further reduces the number of replicas required. Nonetheless, our analysis also reveals significant challenges of relying on tempering alone for sampling large-scale conformational fluctuations of disordered proteins. It is likely that more efficient sampling protocols will require incorporating more sophisticated Hamiltonian replica exchange schemes in addition to tempering.

Ligand-dependent conformational dynamics of MPER-binding HIV-1 antibody Fab235.

The glycoprotein gp160 of HIV-1 is targeted by many antibodies that differ in neutralizing capacity. Due to the high mutation rate of gp160, its relatively conserved membrane-proximal external region (MPER) is an important epitope for eliciting broadly neutralizing antibodies. However, the MPER is located underneath the gp160 ectodomain, making it difficult to access. Thus, the conformational property of the Fab domain of the antibody is likely a key factor in the binding to MPER and HIV-1 neutralization. As a test system, we use vaccine-elicited MPER antibody Fab235, to investigate its ligand-dependent conformational behavior by performing all-atom molecular dynamics simulations. We perform simulations with or without a bound MPER, 1-µs each. Compared to the APO state, the MPER-bound state had increased inter-domain contacts within Fab235 and attendant reduction in conformational fluctuation which propagated down to the constant domains, in particular to CH1. On the other hand, removing MPER from the MPER-bound structure led to reduction of inter-domain contacts and an increase in fluctuation of the CH1 domain, similar to the original APO structure. The compliance of Fab in the APO state may be beneficial for approaching and binding to the MPER, while stiffening upon binding to MPER would support interaction between the antibody and receptors on innate cells or complement proteins interacting via the Fc domain. Our results also indicate that ligand-dependent contacts between the variable and constant domains enable allosteric control of ligand binding.

Using single molecule polarization-sweep spectroscopy to monitor site-specific fluctuations of internal (iCy3)2 dimer-labeled DNA constructs at ss-ds DNA junctions.

Local fluctuations of the sugar-phosphate backbones of DNA (a form of DNA ‘breathing’) play key roles in protein-DNA assembly and enzymatic function. In this work, we monitored spectroscopic signals from single-molecules of DNA fork constructs labeled with two cyanine optical probes [(iCy3)2], which were rigidly inserted into the sugar-phosphate backbones at opposite positions of complementary single-strands. We thus studied the local conformational fluctuations of the DNA sugar-phosphate backbones at specific (iCy3)2 dimer-labeled positions relative to a single-stranded (ss) double-stranded (ds) DNA junction. These experiments are based on a newly-developed single-molecule spectroscopic method that uses a linearly polarized continuous wave laser in which the polarization direction is continuously rotated at radio frequencies. The emitted fluorescence from the single-molecule sample contains information about conformational changes of the exciton-coupled (iCy3)2 dimer probe, which monitor conformational fluctuations of the vicinal DNA sugar-phosphate backbones. Our results indicate that the local conformation of the DNA sugar-phosphate backbones at positions near ss-dsDNA junctions adopts four topologically-relevant macrostates. The probability that a particular macrostate is occupied depends on backbone position. We apply a kinetic network approach to interpret our observations of DNA backbone fluctuations at and near ss-dsDNA junctions, and we propose conformational assignments based on exciton-coupling models for the macrostates that we observe. This polarization-sweep single-molecule fluorescence method and analysis can be applied to study protein binding and macromolecular function at and near ss-dsDNA junctions, thus providing site-specific information about the free energy landscapes at various functionally relevant positions within such macromolecular complexes.

Conformational dynamics of glucagon receptor domains revealed by single-molecule fluorescence microscopy.

G protein-coupled receptors (GPCRs) are integral membrane proteins involved in many cellular signaling processes. Representing the largest family of proteins with approximately 800 members in humans, GPCRs are the target for more than 40% of current pharmaceutical drugs. The structural and dynamic information is vital to understanding the activation mechanism of these receptors. Despite the availability of high-resolution structures of GPCRs at different conformational states, the dynamics of those conformational states at the molecular level are missing. We used total internal reflection fluorescence (TIRF) imaging to investigate the conformational dynamics of the glucagon receptor (GCGR), a class B family GPCR involved in glucose homeostasis. We expressed and purified GCGR from mammalian cells and analyzed them at a single molecule level. We used single-molecule Förster resonance energy transfer (smFRET) to observe the conformational dynamics of the donor and acceptor fluorophore-labeled GCGR molecules. We observed ligand-dependent slow and reversible conformational fluctuation of N-terminal extra cellular domain (ECD) without ligand binding. The addition of a saturated concentration of peptide ligand (Glucagon) altered the dynamics towards the lower FRET states. These results provide new insight into the structural dynamics of glucagon receptor and receptor activation.

Investigating the binding affinity of andrographolide against human SARS-CoV-2 spike receptor-binding domain through docking and molecular dynamics simulations

SARS-CoV-2 is a positive-sense single-stranded RNA virus that causes a deadly coronavirus disease (COVID-19) in humans. The infection of SARS-CoV-2 in humans involves a viral surface spike glycoprotein containing the receptor-binding domain (RBD). The interactions of SARS-CoV-2 with the host angiotensin-converting enzyme 2 (ACE2) receptor are mediated by RBD. It binds to the host ACE2 and influences viral replication and disease pathogenesis. Therefore, targeting the RBD to prevent SARS-CoV-2 infections is of utmost importance. In this study, we used docking and molecular dynamics simulations to understand the binding effect of andrographolide on the SARS-CoV-2 spike protein. During docking, a strong binding affinity was observed between the ligand and the target receptor protein. MD results demonstrated higher conformational fluctuations in the ligand-free protein compared to the bound form. Several residues in the active sites make conformational rearrangements for the S protein to interact with the ligand. While RBD experiences conformational transition to gain more stability upon binding with the ligand. This binding is strengthened via several non-covalent interactions that make the complex structure more stable with higher binding affinity. Overall findings of the study may shed some valuable insights concerning the development of potential therapeutics in the strategies for COVID-19 prevention.

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.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins.

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

De novo design of anti-variant COVID-19 vaccine.

Recent studies highlight the effectiveness of hybrid Severe Acute RespiratorySyndrome-Coronavirus-2 (SARS-CoV-2) vaccines combining wild-type nucleocapsid and Spike proteins. We have further enhanced this strategy by incorporating delta and omicron variants' spike protein mutations. Both delta and omicron mark the shifts in viral transmissibility and severity in unvaccinated and vaccinated patients. So their mutations are highly crucial for future viral variants also. Omicron is particularly adept at immune evasion by mutating spike epitopes. The rapid adaptations of Omicron and sub-variants to spike-based vaccines and simultaneous transmissibility underline the urgency for new vaccines in the continuous battle against SARS-CoV-2. Therefore, we have added three persistent T-cell-stimulating nucleocapsid peptides similar to homologous sequences from seasonal Human Coronaviruses (HuCoV) and an envelope peptide that elicits a strong T-cell immune response. These peptides are clustered in the hybrid spike's cytoplasmic region with non-immunogenic linkers, enabling systematic arrangement. AlphaFold (Artificial intelligence-based model building) analysis suggests omitting the transmembrane domain enhances these cytoplasmic epitopes' folding efficiency which can ensure persistent immunity for CD4+ structural epitopes. Further molecular dynamics simulations validate the compact conformation of the modeled structures and a flexible C-terminus region. Overall, the structures show stability and less conformational fluctuation throughout the simulation. Also, the AlphaFold predicted structural epitopes maintained their folds during simulation to ensure the specificity of CD4+ T-cell response after vaccination. Our proposed approach may provide options for incorporating diverse anti-viral T-cell peptides, similar to HuCoV, into linker regions. This versatility can be promising to address outbreaks and challenges posed by various viruses for effective management in this era of innovative vaccines.

Elucidating the Sluggish Water Dynamics at the Ice-Binding Surface of the Hyperactive Tenebrio molitor Antifreeze Protein.

Quasi-ice-like hydration waters on the ice-binding surface (IBS) of an antifreeze protein (AFP) commonly exhibit sluggish dynamics especially at low temperatures. In this work, we have analyzed molecular dynamics (MD) simulation trajectories at two different temperatures for Tenebrio molitor antifreeze protein (TmAFP) to explore whether the unique quasi-ice-like structuring of hydration water has any impact on making their dynamics slower on the IBS of the protein. Our calculation reveals that, as translational dynamics is coupled with the conformational fluctuations, hydration water on the IBS exhibits sluggish translational motion due to reduced flexibility of the IBS compared to that on the non-ice-binding surface (NIBS) of the protein. Interestingly, it is noticed that rotational motion of hydration water is not coupled with the conformational fluctuations of the surfaces. In that case, structural relaxations of the protein-water (PW) and water-water (WW) hydrogen bonds compete with each other to make the rotational dynamics of hydration water around the IBS either faster or slower with respect to those around the NIBS. At low temperature, the slower structural relaxation of water-water hydrogen bonds dominates and imparts sluggish rotational motion of the hydration water on the IBS of the protein. The slower structural relaxation of water-water hydrogen bonds and hence the retarded rotational dynamics, despite the weak short-lived PW hydrogen bonds on the IBS, is clearly a manifestation of the rigid quasi-ice-like structure of the hydration shell on that surface.

5-Formylcytosine weakens the G-C pair and imparts local conformational fluctuations to DNA duplexes.

DNA epigenetic modifications such as 5-methyl (5mC), 5-hydroxymethyl (5hmC), 5-formyl (5fC) and 5-carboxyl (5caC) cytosine have unique and specific biological roles. Crystallographic studies of 5mC containing duplexes were conducted in the A-, B- or the intermediate E-DNA polymorphic forms. 5fC-modified duplexes initially observed in the disputed F-DNA architecture were subsequently crystallized in the A-form, suggesting that epigenetic modifications enable DNA sequences to adopt diverse conformational states that plausibly contribute to their function. Solution-state studies of these modifications were found in the B-DNA form, with marked differences in the conformational flexibility of 5fC containing duplexes in comparison to C/5mC containing duplexes, compromising the DNA duplex's stability. Herein, we systematically evaluate sensitive and commonly inaccessible NMR parameters to map the subtle differences between C, 5mC, and their oxidized (5hmC/5fC) counterparts. We observe that 15N/1H chemical shifts effectively report on the weakening of 5fC-G Watson-Crick base-pair H-bonding, extending the instability beyond any achievable within the sequence-specific changes in DNA. Triple 5fC containing sequences propagate the destabilization farther from the site of modifications, explaining reduced duplex stability upon multiple modifications. Additionally, scalar and residual dipolar coupling measurements unravel local sugar pucker fluctuations. One-bond 13C-1H scalar coupling measurements point towards a significant deviation away from the anticipated C2'-endo pucker for the 5fC modified nucleotide. Structural models obtained employing 13C-1H residual dipolar couplings and inter-proton distances corroborate the sugar pucker's deviation for 5fC modified DNA duplexes. The changes in the sugar pucker equilibria remain local to the 5fC modified nucleotide sans additive/long-range effects arising from multiple contiguous modifications. These observations highlight the impact of a major groove modification that alters the physical properties of DNA duplex without disturbing the Watson-Crick face. The changes observed in our studies for the 5fC containing DNA contrast with the perturbations induced by damage/lesion highlight the varied conformational preferences that modified nucleobases impart to the DNA duplex. As sequence-specific DNA transactions are rooted in the base-pair stability and pucker deviations, the observed structural perturbations for 5fC-modified DNA potentially play critical functional roles, such as protein-DNA recognition and interactions.

Sequence Tendency for the Interaction between Low-Complexity Intrinsically Disordered Proteins.

Reversible interaction between intrinsically disordered proteins (IDPs) is considered as the driving force for liquid-liquid phase separation (LLPS), while the detailed description of such a transient interaction process still remains a challenge. And the mechanisms underlying the behavior of IDP interaction, for example, the possible relationship with its inherent conformational fluctuations and sequence features, remain elusive. Here, we use atomistic molecular dynamics (MD) simulation to investigate the reversible association of the LAF-1 RGG domain, the IDP with ultra-low LLPS concentration (0.06 mM). We find that the duration of the association between two RGG domains is highly heterogeneous, and the sustained associations essentially dominate the IDP interaction. More interestingly, such sustained associations are mediated by a finite region, that is, the C-terminal region 138-168 (denoted as a contact-prone region). We noticed that such sequence tendency is attributed to the extended conformation of the RGG domain during its inherent conformational fluctuations. Hence, our results suggest that there is a certain region in this low-complexity IDP which can essentially dominate their interaction and should be also important to the LLPS. And the inherent conformational fluctuations are actually essential for the emergence of such a hot region of IDP interaction. The importance of this hot region to LLPS is verified by experiment.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

Strand 1A variant in neuroserpin shows increased aggregation and no loss of inhibition: implication in ameliorating polymerization to retain activity.

Neuroserpin (NS) is predominantly expressed in the brain and is the primary inhibitor of tissue plasminogen activator (tPA). NS variants are associated with the neurogenerative disease termed familial encephalopathy with neuroserpin inclusion bodies (FENIB). The disease is characterized by variable age of onset and severity. The reactive center loop (RCL) insertion-based inhibitory mechanism of NS requires a coordinated conformational change leading to a shift in the strands of the β-sheet A and movement of helix F. Strand 1A is connected to the helix F at its C terminal end and with the strand 2A at its N terminal, both these domain move for accommodating the inserting loop; therefore, a variant that influences their movement may alter the inhibition rates. A molecular dynamic simulation analysis of a H138C NS variant from strand 1A showed a large decrease in conformational fluctuations as compared with wild-type NS. H138 was mutated, expressed, purified and a native-PAGE and transmission electron microscopy (TEM) analysis showed that this variant forms large molecular weight aggregates on a slight increase in temperature. However, a circular dichroism analysis showed its secondary structure to be largely conserved. Surprisingly, its tPA inhibition activity and complex formation remain unhindered even after the site-specific labeling of H138C with Alexa fluor C5 maleimide. Further, a helix F-strand 1A (W154C-H138C) double variant still shows appreciable inhibitory activity. Increasingly, it appears that aggregation and not loss of inhibition is the more likely cause of shutter region-based variants phenotypes, indicating that hindering polymer formation using small molecules may retain inhibitory activity in pathological variants of NS.

Protein as evolvable functionally constrained amorphous matter

We explore current ideas around the representation of a protein as an amorphous material, in turn represented by an abstract graph G with edges weighted by elastic stiffnesses. By embedding this graph in physical space, we can map every graph to a spectrum of conformational fluctuations and responses (as a result of, say, ligandbinding). This sets up a 'genotype-phenotype' map, which we use to evolve the amorphous material to select for fitness. Using this, we study the emergence of allosteric interaction, hinge joint, crack formation and a slide bolt in functional proteins such as adenylate kinase, HSP90, calmodulin and GPCR proteins. We find that these emergent features are associated with specific geometries and mode spectra of floppy or liquid-like regions. Our analysis provides insight into understanding the architectural demands on a protein that enable a prescribed function and its stability to mutations.

Explicit analytical form for memory kernel in the generalized Langevin equation for end-to-end vector of Rouse chains.

The generalized Langevin equation (GLE) provides an attractive theoretical framework for investigating the dynamics of conformational fluctuations of polymeric systems. While the memory kernel is a central function in the GLE, explicit analytical forms for this function have been challenging to obtain, even for the simple models of polymer dynamics. Here, we achieve an explicit analytical expression for the memory kernel in the GLE for the end-to-end vector of Rouse chains in the overdamped limit. Our derivation takes advantage of the finding that the dynamics of the end-to-end vector of Rouse chains with both free ends are equivalent to those of Rouse chains with one free end and the other fixed. For the latter model, we first show that the equations of motion of the Rouse modes as well as their statistical properties can be obtained under the boundary conditions where the free end is held fixed temporarily. We then analytically solve the terms associated with intrachain interactions in the GLE. By formally comparing these terms with the GLE based on the Rouse modes, we obtain an explicit expression for the memory kernel, along with analytical forms for the potential field and the random colored noise force. Our analytical memory kernel is confirmed by numerical calculations in the Laplace space and is shown to yield asymptotic behaviors that are consistent with previous studies. Finally, we utilize our analytical result to simulate the cyclization dynamics of Rouse chains and discuss the scaling of the cyclization time with chain length.

Tenuous transcriptional threshold of human sex determination. II. SRY exploits water-mediated clamp at the edge of ambiguity.

Y-encoded transcription factor SRY initiates male differentiation in therian mammals. This factor contains a high-mobility-group (HMG) box, which mediates sequence-specific DNA binding with sharp DNA bending. A companion article in this issue described sex-reversal mutations at box position 72 (residue 127 in human SRY), invariant as Tyr among mammalian orthologs. Although not contacting DNA, the aromatic ring seals the domain's minor wing at a solvent-exposed junction with a basic tail. A seeming paradox was posed by the native-like biochemical properties of inherited Swyer variant Y72F: its near-native gene-regulatory activity is consistent with the father's male development, but at odds with the daughter's XY female somatic phenotype. Surprisingly, aromatic rings (Y72, F72 or W72) confer higher transcriptional activity than do basic or polar side chains generally observed at solvated DNA interfaces (Arg, Lys, His or Gln). Whereas biophysical studies (time-resolved fluorescence resonance energy transfer and heteronuclear NMR spectroscopy) uncovered only subtle perturbations, dissociation of the Y72F complex was markedly accelerated relative to wild-type. Studies of protein-DNA solvation by molecular-dynamics (MD) simulations of an homologous high-resolution crystal structure (SOX18) suggest that Y72 para-OH anchors a network of water molecules at the tail-DNA interface, perturbed in the variant in association with nonlocal conformational fluctuations. Loss of the Y72 anchor among SRY variants presumably "unclamps" its basic tail, leading to (a) rapid DNA dissociation despite native affinity and (b) attenuated transcriptional activity at the edge of sexual ambiguity. Conservation of Y72 suggests that this water-mediated clamp operates generally among SRY and metazoan SOX domains.

Computational investigation of conformational fluctuations of Aβ42 monomers in aqueous ionic liquid mixtures

Aggregation of intrinsically disordered proteins (IDPs) is often triggered by protein–protein interaction, leading to several neurodegenerative diseases.In this study, we aim to understand the role of 1-butyl-3-methylimidazolium tetrafluoroborate [BMIM][BF4] ionic liquid (IL) on altering secondary structure of the peptide and study the interactions responsible for their conformational modification. Here we report structural features of an ensemble of amyloid beta (Aβ42) peptide monomers employing all-atom molecular dynamics simulations in aqueous solution of IL at room temperature. Our calculations reveal, Aβ42 peptide exhibit rapidly inter-converting conformations adopting extended chain as well as collapsed coil-like structures with a propensity to form few transient α-helices and β-sheets at intermediate times. We have correlated this higher flexibility of the peptide with formation of less stable salt-bridge in presence of IL in aqueous solution. It is found that water molecules around the polar regions are more structured imparting structural rigidity to the peptide backbone than around the hydrophobic domains. Thus identifying the effect of electrostatic interaction of IL on expulsion of water molecules from the polar domains on the peptide surface can shed light on preventing the aggregation process. We have also determined the most energetically stable peptide conformations from various accessible intermediate states using metadynamics simulations that are characterized by the rugged free energy surface having multiple energy basins.