Nonspecific Attractive Interactions Research Articles

Fly casting with ligand sliding and orientational selection supporting complex formation of a GPCR and a middle sized flexible molecule

A GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was conducted to elucidate binding mechanisms of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB). GA-mD-VcMD is a generalized ensemble method that produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stable maintenance of the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model. Then sampling was conducted from conformations where bosentan was distant from the binding site in the hETB binding pocket. The deepest basin in the resultant free-energy landscape was assigned to native-like complex conformation. The following binding mechanism was inferred. First, bosentan fluctuating randomly in solution is captured using a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly casting). Bosentan then slides occasionally from the tip to the root of the N-terminal tail (ligand–sliding). During this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with an accompanying rapid reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand–receptor attractive native contacts are formed. Eventually, the native-like complex is completed. The bosentan-captured conformations by the tip-region and root-region of the N-terminal tail correspond to two basins in the free-energy landscape. The ligand-sliding corresponds to overcoming of a free-energy barrier between the basins.

Seeding the Self-Assembly of DNA Origamis at Surfaces.

Unlike supramolecular self-assembly methods that can organize many distinct components into designer shapes in a homogeneous solution (e.g., DNA origami), only relatively simple, symmetric structures consisting of a few distinct components have been self-assembled at solid surfaces. As the self-assembly process is confined to the surface/interface by mostly nonspecific attractive interactions, an open question is how these interfacial interactions affect multicomponent self-assembly. To gain a mechanistic understanding of the roles of the surface environment in DNA origami self-assembly, here we studied the oligonucleotide-assisted folding of a long single-stranded DNA (ssDNA scaffold) that was end-tethered to a dynamic surface, which could actively regulate the DNA-surface interactions. The results showed that even weak surface attractions can lead to defective structures by inhibiting the merging of multiple domains into complete structures. A combination of surface anchoring and deliberate regulation of DNA-surface interactions allowed us to depart from the existing paradigm of surface confinement via nonspecific interactions and enabled DNA origami folding to proceed in a solution-like environment. Importantly, our strategy retains the key advantages of surface-mediated self-assembly. For example, surface-anchored oligonucleotides could sequence-specifically initiate the growth of DNA origamis of specific sizes and shapes. Our work enables information to be encoded into a surface and expressed into complex DNA surface architectures for potential nanoelectronic and nanophotonic applications. In addition, our approach to surface confinement may facilitate the 2D self-assembly of other molecular components, such as proteins, as maintaining conformational freedom may be a general challenge in the self-assembly of complex structures at surfaces.

Gluonic and Regulatory Solvents: A Paradigm for Tunable Phase Segregation in Polymers

A theoretical concept for phase segregation of polymers in the presence of multicomponent selective solvents is presented. Phase separation is caused by nonspecific attractive interactions between the polymers and a smaller component in the solution instead of repulsion between monomers and solvent molecules. We call the component that adsorbs on the polymers and thus causes condensation “gluonic”. It is shown that a discontinuous phase transition from a diluted or semidiluted state to the condensed state takes place if the fraction of the gluonic component is increased. The location of the phase coexistence can be shifted by influencing the binding efficiency of the gluonic component. This can be achieved by adding a regulator component to the solution. The latter is assumed to bind specifically to the polymer and interferes with the gluonic component. In this way a switching of the state of the polymer from dissolved to condensed can be achieved without changing the chemical properties of the system. Ap...

Nonexponential Kinetics of Loop Formation in Proteins and Peptides: A Signature of Rugged Free Energy Landscapes?

The kinetics of loop formation, i.e., the occurrence of contact between two atoms of a polypeptide, remains the focus of continuing interest. One of the reasons is that contact formation is the elementary event underlying processes such as folding and binding. More importantly, it is experimentally measurable and can be predicted theoretically for ideal polymers. Deviations from single exponential kinetics have sometimes been interpreted as a signature of rugged, protein-like, free energy landscapes. Here we present simulations, with different atomistic models, of short peptides with varied structural propensity, and of a structured protein. Results show exponential contact formation kinetics (or relaxation) at long times, and a power law relaxation at very short times. At intermediate times, a deviation from either power law or simple exponential kinetics is observed that appears to be characteristic of polypeptides with either specific or nonspecific attractive interactions but disappears if attractive interactions are absent. Our results agree with recent experimental measurements on peptides and proteins and offer a comprehensive interpretation for them.

Perspective: On the importance of hydrodynamic interactions in the subcellular dynamics of macromolecules.

An outstanding challenge in computational biophysics is the simulation of a living cell at molecular detail. Over the past several years, using Stokesian dynamics, progress has been made in simulating coarse grained molecular models of the cytoplasm. Since macromolecules comprise 20%-40% of the volume of a cell, one would expect that steric interactions dominate macromolecular diffusion. However, the reduction in cellular diffusion rates relative to infinite dilution is due, roughly equally, to steric and hydrodynamic interactions, HI, with nonspecific attractive interactions likely playing rather a minor role. HI not only serve to slow down long time diffusion rates but also cause a considerable reduction in the magnitude of the short time diffusion coefficient relative to that at infinite dilution. More importantly, the long range contribution of the Rotne-Prager-Yamakawa diffusion tensor results in temporal and spatial correlations that persist up to microseconds and for intermolecular distances on the order of protein radii. While HI slow down the bimolecular association rate in the early stages of lipid bilayer formation, they accelerate the rate of large scale assembly of lipid aggregates. This is suggestive of an important role for HI in the self-assembly kinetics of large macromolecular complexes such as tubulin. Since HI are important, questions as to whether continuum models of HI are adequate as well as improved simulation methodologies that will make simulations of more complex cellular processes practical need to be addressed. Nevertheless, the stage is set for the molecular simulations of ever more complex subcellular processes.

Equilibrium simulation of trp-cage in the presence of protein crowders.

While steric crowders tend to stabilize globular proteins, it has been found that protein crowders can have an either stabilizing or destabilizing effect, where a destabilization may arise from nonspecific attractive interactions between the test protein and the crowders. Here, we use Monte Carlo replica-exchange methods to explore the equilibrium behavior of the miniprotein trp-cage in the presence of protein crowders. Our results suggest that the surrounding crowders prevent trp-cage from adopting its global native fold, while giving rise to a stabilization of its main secondary-structure element, an α-helix. With the crowding agent used (bovine pancreatic trypsin inhibitor), the trp-cage-crowder interactions are found to be specific, involving a few key residues, most of which are prolines. The effects of these crowders are contrasted with those of hard-sphere crowders.

Kinetic of Loop Formation in Polypeptides and Free Energy Landscapes

The kinetics of contact formation between the two ends of a polypeptide has been the focus of extended interest in the recent and less recent past. One of the reasons is that contact formation is the elementary event underlying processes such as folding and binding. More importantly, it is experimentally measurable with increasing time and space resolution; it can be predicted theoretically for ideal polymers. Deviations from single exponential kinetics has been sometimes interpreted as a signature of a rugged, protein-like, free energy landscape. Here we present simulations of short peptides with different structural propensity and of a structured protein with different atomistic models. Results show a remarkable deviation from exponentiality for times shorter than a characteristic time that depends on the sequence and length of polypeptides. Deviation from exponentiality cannot be simply explained with the non-Markovian polymer dynamics at very short times, but appears to be characteristic of polypeptides with either specific or non-specific attractive interactions. Our results agree with recent experimental measurements and offer a comprehensive interpretation for them.

Molecular Crowding Drives Active Pin1 into Nonspecific Complexes with Endogenous Proteins Prior to Substrate Recognition

Proteins and nucleic acids maintain the crowded interior of a living cell and can reach concentrations in the order of 200-400 g/L which affects the physicochemical parameters of the environment, such as viscosity and hydrodynamic as well as nonspecific strong repulsive and weak attractive interactions. Dynamics, structure, and activity of macromolecules were demonstrated to be affected by these parameters. However, it remains controversially debated, which of these factors are the dominant cause for the observed alterations in vivo. In this study we investigated the globular folded peptidyl-prolyl isomerase Pin1 in Xenopus laevis oocytes and in native-like crowded oocyte extract by in-cell NMR spectroscopy. We show that active Pin1 is driven into nonspecific weak attractive interactions with intracellular proteins prior to substrate recognition. The substrate recognition site of Pin1 performs specific and nonspecific attractive interactions. Phosphorylation of the WW domain at Ser16 by PKA abrogates both substrate recognition and the nonspecific interactions with the endogenous proteins. Our results validate the hypothesis formulated by McConkey that the majority of globular folded proteins with surface charge properties close to neutral under physiological conditions reside in macromolecular complexes with other sticky proteins due to molecular crowding. In addition, we demonstrate that commonly used synthetic crowding agents like Ficoll 70 are not suitable to mimic the intracellular environment due to their incapability to simulate biologically important weak attractive interactions.

Quantitative assessment of the relative contributions of steric repulsion and chemical interactions to macromolecular crowding

The term "macromolecular crowding" denotes the combined effects of high volume fractions of nominally unrelated macromolecules upon the equilibrium and transport properties of all macrosolutes, dilute as well as concentrated, in the crowded medium. We present a formal partitioning of the total crowding effect into contributions from steric exclusion (excluded volume) and weak, nonspecific attractive interactions between a concentrated "crowding agent" and reactant and product species present at trace concentration. A numerical example of the combined effect of both steric and chemical interactions between crowder and tracer upon the reversible dimerization of tracer is presented, based upon reasonable estimates of the magnitude of both repulsive and attractive interactions between tracer and crowder species.

Viscosity Behavior of High-Concentration Monoclonal Antibody Solutions: Correlation with Interaction Parameter and Electroviscous Effects

The purpose of this work was to understand the viscosity behavior of high-concentration monoclonal antibody (mAb) solutions in a wide range of solution conditions and generate guidelines helpful to formulate products with manageable viscosity. The zeta potential and effective isoelectric point (pI) were determined from electrophoretic mobility measurements. High-frequency rheology studies characterized viscoelasticity at high concentrations. The interaction parameter (k(D) ) obtained from dynamic light scattering quantified intermolecular interactions. Circular dichroism characterized conformational stability upon change in solution pH. Except for mAb-1, all other mAb solutions were found to be more viscous at solution pHs closer to the molecular pI. For mAb-2, mAb-3, and mAb-10,the k(D) indicated intermolecular attractions at the pI, wherein the net molecular charge (Z) was zero, whereas repulsions dominated at pHs away from the pI. At the pI, Z and, hence, the charge-induced repulsions are minimal, whereas the charge distribution becomes most conspicuous. The resulting dominance of nonspecific attractive interactions at the pI increases the self-association or aggregation behavior of protein molecules, leading to a higher viscosity at the pI. mAb-1 is an exception to this general behavior. The k(D) could serve as a qualitative screening tool to predict high-concentration viscosity behavior, whereas the correlation with net charge was inconsistent. A higher negative k(D) generally resulted in a more viscous solution at high concentrations; however, direct quantitative assessment was not possible.

Crowding and hydrodynamic interactions likely dominate in vivo macromolecular motion

To begin to elucidate the principles of intermolecular dynamics in the crowded environment of cells, employing brownian dynamics (BD) simulations, we examined possible mechanism(s) responsible for the great reduction in diffusion constants of macromolecules in vivo from that at infinite dilution. In an Escherichia coli cytoplasm model comprised of 15 different macromolecule types at physiological concentrations, BD simulations of molecular-shaped and equivalent sphere representations were performed with a soft repulsive potential. At cellular concentrations, the calculated diffusion constant of GFP is much larger than experiment, with no significant shape dependence. Next, using the equivalent sphere system, hydrodynamic interactions (HI) were considered. Without adjustable parameters, the in vivo experimental GFP diffusion constant was reproduced. Finally, the effects of nonspecific attractive interactions were examined. The reduction in diffusivity is very sensitive to macromolecular radius with the motion of the largest macromolecules dramatically slowed down; this is not seen if HI dominate. In addition, long-lived clusters involving the largest macromolecules form if attractions dominate, whereas HI give rise to significant, size independent intermolecular dynamic correlations. These qualitative differences provide a testable means of differentiating the importance of HI vs. nonspecific attractive interactions on macromolecular motion in cells.