All-atom Explicit-solvent Molecular Dynamics Simulations Research Articles

Probing the Nucleic Acid Flexibility to Disarm the Viral Counter-Defense Machinery: Design and Characterization of Potent p19 Inhibitors.

Plant viruses are highly destructive and significant contributors to several global pandemics and epidemics in plants. A viral disease outbreak in plants can cause a scarcity of food supply and is a severe concern to humanity. The siRNA (small interfering RNA)-mediated RNA-induced silencing complex (RISC) formation is a primary defense mechanism in plants against viruses, where the RISC binds and degrades viral mRNAs. As a counter-defense, many viruses encode RNA-silencing suppressor proteins (e.g., the p19 protein from the Tombusviridae family) for viral proliferation in plants. The functional form of p19 (homodimer) binds to plant siRNA with high affinities, thereby interrupting the RISC formation and thus preventing the viral mRNA silencing in plants. By altering the RISC formation, the p19 protein helps the virus invasion in the plant and ultimately stunts host growth. In this study, we designed several modified siRNA-based molecules for p19 inhibition. The viral p19 protein is known to interact predominantly through H-bonds with 2'-OH and phosphates of the plant siRNA. We utilized this information and in silico-designed flexible substituents of siRNA, where we removed the C2'-C3' bond in each nucleotide unit. We performed all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicates each) for control/modified siRNA─p19 complexes (8 in total) followed by energetic estimations. Strikingly, in a few modified complexes, the siRNA not only retained the double-helical structural integrity but also displayed remarkably enhanced p19 binding compared to the control siRNA; hence, we consider it important to perform biological and chemical in vitro and in vivo studies on proposed flexible nucleic acids as p19 inhibitors for crop protection.

Elucidating the dynamics of individual IDPs in a biomolecular condensate using single-molecule FRET and molecular simulation.

A wide range of biomolecules in solution can phase-separate and form membraneless organelles in the cell. These assemblies often have liquid-like properties, and the corresponding dynamics and exchange of molecules with the environment are important for biological function. The dynamics and materials properties of these systems are commonly assessed based on translational diffusion or rheological properties, typically covering timescales of milliseconds and longer. However, information on the structure and dynamics at the molecular level is lacking. We have studied coacervates of two intrinsically disordered highly oppositely charged nuclear proteins as a prototypical example of heterotypic liquid-liquid phase separation of components abundant in the cell nucleus. Using single-molecule fluorescence spectroscopy, we show that the proteins in the condensates remain disordered, and that their chain dynamics occur on timescale surprisingly close to the ∼100-nanosecond dynamics of the 1:1 complex in dilute solution, even though the condensate is almost three orders of magnitude more viscous than the dilute solution and translational diffusion of the same proteins is ∼50 times slower in the dense than in the dilute phase. The experimental results are in good agreement with large-scale all-atom explicit-solvent molecular dynamics simulations consisting of ∼200 proteins and ∼4 million atoms in total, which can thus provide detailed insights into the protein interactions and dynamics within the condensate. The simulations indicate that the rapid protein chain dynamics in the dense phase are linked to the similarity of the local environment experienced by the proteins in the dilute and dense phases and the short lifetime of intermolecular contacts in the dense phase.

Harmonizing Interstrand Electrostatic Repulsion by Conformational Rigidity in Counterion-Deprived Z-DNA: A Molecular Dynamics Study.

Deoxyribonucleic acid (DNA) is a vital biomacromolecule. Although the right-handed B-DNA type helical structure is the most abundant and extensively studied form of DNA, several noncanonical forms, such as triplex, quadruplex, Z-DNA, A-DNA, and ss-DNA, have been probed from time to time to gain insights into the DNA's function. Z-DNA was recently found to be involved in cancer and several autoimmune diseases. In the present Article, we evaluated the conformational stability of locked-sugar-based Z-DNA via all-atom explicit-solvent molecular dynamics simulations and found that the modified DNA maintained the left-handed conformation even in the absence of counterions, wherein the structural rigidity dominates over the electrostatic repulsion between the complementary strands. The control Z-DNA without counterions, as expected, instantaneously resulted in unfolded states. The remarkable stability of the conformationally locked model system was thoroughly investigated via structural and energetic perspectives and was probably the result of the backbone widening in tandem with enhanced electrostatics between complementary strands. We believe that the design of the proposed modified Z-DNA construct could help understand the otherwise delicate Z-DNA conformation even in salt-deprived conditions. The design could also motivate the medicinal use of short segments of such modified nucleotides and could be utilized in more advanced modeling techniques, such as DNA origami which has gained popularity in recent years.

Atomistic Insights into A315E Mutation-Enhanced Pathogenicity of TDP-43 Core Fibrils.

The aggregation of TAR DNA-binding protein of 43 kDa (TDP-43) into fibrillary deposits is implicated in amyotrophic lateral sclerosis (ALS), and some hereditary mutations localized in the low complexity domain (LCD) facilitate the formation of pathogenic TDP-43 fibrils. A recent cryo-EM study reported the atomic-level structures of the A315E TDP-43 LCD (residues 288-319, TDP-43288-319) core fibril in which the protofilaments have R-shaped structures and hypothesized that A315E U-shaped protofilaments can readily convert to R-shaped protofilaments compared to the wild-type (WT) ones. There are no atomic structures of WT protofilaments available yet. Herein, we performed extensive all-atom explicit-solvent molecular dynamics simulations on A315E and WT protofilaments starting from both the cryo-EM-determined R-shaped and our constructed U-shaped structures. Our simulations show that WT protofilaments also adopt the R-shaped structures but are less stable than their A315E counterparts. Except for R293-E315 salt bridges, N312-F316 hydrophobic interactions and F316-F316 π-π stacking interactions are also crucial for the stabilization of the neck region of the R-shaped A315E protofilaments. The loss of R293-E315 salt bridges and the weakened interactions of N312-F316 and F316-F316 result in the reduced stability of the R-shaped WT protofilaments. Simulations starting from U-shaped folds reveal that A315E protofilaments can spontaneously convert to the cryo-EM-derived R-shaped protofilaments, whereas WT protofilaments convert to R-shape-like structures with remodeled neck regions. The R-shape-like WT protofilaments could act as intermediate states slowing down the U-to-R transition. This study reveals that A315E mutation can not only enhance the structural stability of the R-shaped TDP-43288-319 protofilaments but also promote the U-to-R transition, which provides atomistic insights into the A315E mutation-enhanced TDP-43 pathogenicity in ALS.

Assessing the DNA structural integrity via selective annihilation of Watson-Crick hydrogen bonds: Insights from molecular dynamics simulations.

Understanding the role of base pairing and stacking displayed by polynucleotide chains interwind together resulting in a double-helical B-DNA type structure is crucial to gaining access to the sophisticated structural arrangement of DNA. Several computational and experimental studies hinted towards the dominance of base pairing over stacking for duplex stability. To find out how significant the individual Watson-Crick hydrogen bonds are in maintaining the double-helical integrity of the DNA, in the present article, we selectively switched off the hydrogen bonds (one specific bond or their combinations in all the base pairs at a time) via manipulating the force fields for A-T and G-C base pairs. We studied 12 systems in total via all-atom explicit-solvent molecular dynamics simulations (200ns each). The MD output structures were compared with the control system by means of structural, dynamic, and energetic properties to monitor the overall consequences of removing H-bond(s) on the B-DNA characteristics of the model systems. Our findings suggest that all the individual hydrogen bonds involved in base pairing are vital for maintaining the DNA structural integrity as any possible alteration in Watson-Crick hydrogen bond(s) leads to the disintegration/collapse of DNA strands resulting in unfolded states.

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins.

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

Combinatorial Virtual Library Screening Study of Transforming Growth Factor-β2-Chondroitin Sulfate System.

Transforming growth factor-beta (TGF-β), a member of the TGF-β cytokine superfamily, is known to bind to sulfated glycosaminoglycans (GAGs), but the nature of this interaction remains unclear. In a recent study, we found that preterm human milk TGF-β2 is sequestered by chondroitin sulfate (CS) in its proteoglycan form. To understand the molecular basis of the TGF-β2–CS interaction, we utilized the computational combinatorial virtual library screening (CVLS) approach in tandem with molecular dynamics (MD) simulations. All possible CS oligosaccharides were generated in a combinatorial manner to give 24 di- (CS02), 192 tetra- (CS04), and 1536 hexa- (CS06) saccharides. This library of 1752 CS oligosaccharides was first screened against TGF-β2 using the dual filter CVLS algorithm in which the GOLDScore and root-mean-square-difference (RMSD) between the best bound poses were used as surrogate markers for in silico affinity and in silico specificity. CVLS predicted that both the chain length and level of sulfation are critical for the high affinity and high specificity recognition of TGF-β2. Interestingly, CVLS led to identification of two distinct sites of GAG binding on TGF-β2. CVLS also deduced the preferred composition of the high specificity hexasaccharides, which were further assessed in all-atom explicit solvent MD simulations. The MD results confirmed that both sites of binding form stable GAG–protein complexes. More specifically, the highly selective CS chains were found to engage the TGF-β2 monomer with high affinity. Overall, this work present key principles of recognition with regard to the TGF-β2–CS system. In the process, it led to the generation of the in silico library of all possible CS oligosaccharides, which can be used for advanced studies on other protein–CS systems. Finally, the study led to the identification of unique CS sequences that are predicted to selectively recognize TGF-β2 and may out-compete common natural CS biopolymers.

Coarse-Grained Modeling of Coronavirus Spike Proteins and ACE2 Receptors

We developed coarse-grained models of spike proteins in SARS-CoV-2 coronavirus and angiotensin-converting enzyme 2 (ACE2) receptor proteins to study the endocytosis of a whole coronavirus under physiologically relevant spatial and temporal scales. We first conducted all-atom explicit-solvent molecular dynamics simulations of the recently characterized structures of spike and ACE2 proteins. We then established coarse-grained models using the shape-based coarse-graining approach based on the protein crystal structures and extracted the force field parameters from the all-atom simulation trajectories. To further analyze the coarse-grained models, we carried out normal mode analysis of the coarse-grained models to refine the force field parameters by matching the fluctuations of the internal coordinates with the original all-atom simulations. Finally, we demonstrated the capability of these coarse-grained models by simulating the endocytosis of a whole coronavirus through the host cell membrane. We embedded the coarse-grained models of spikes on the surface of the virus envelope and anchored ACE2 receptors on the host cell membrane, which is modeled using a one-particle-thick lipid bilayer model. The coarse-grained simulations show the spike proteins adopt bent configurations due to their unique flexibility during their interaction with the ACE2 receptors, which makes it easier for them to attach to the host cell membrane than rigid spikes.

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions.

The effects of sulfation and calcium cations (Ca2+) on the atomic-resolution conformational properties of chondroitin carbohydrate polymers in aqueous solutions are not well studied owing to experimental challenges. Here, we compare all-atom explicit-solvent molecular dynamics simulations results for pairs of O-type (nonsulfated) and A-type (GlcNAc 4-O-sulfated) chondroitin 20-mers in 140 mM NaCl with and without Ca2+ and find that both sulfation and Ca2+ favor more compact polymer conformations. We also show that subtle differences in force-field parametrization can have dramatic effects on Ca2+ binding to chondroitin carboxylate and sulfate functional groups and thereby determine Ca2+-mediated intra- and interstrand association. In addition to providing an atomic-resolution picture of the interaction of Ca2+ with sulfated and nonsulfated chondroitin polymers, the molecular dynamics data emphasize the importance of careful force-field parametrization to balance ion–water and ion–chondroitin interactions and suggest additional parametrization efforts to tune interactions involving sulfate.

Efficient Construction of Atomic-Resolution Models of Non-Sulfated Chondroitin Glycosaminoglycan Using Molecular Dynamics Data.

Glycosaminoglycans (GAGs) are linear, structurally diverse, conformationally complex carbohydrate polymers that may contain up to 200 monosaccharides. These characteristics present a challenge for studying GAG conformational thermodynamics at atomic resolution using existing experimental methods. Molecular dynamics (MD) simulations can overcome this challenge but are only feasible for short GAG polymers. To address this problem, we developed an algorithm that applies all conformational parameters contributing to GAG backbone flexibility (i.e., bond lengths, bond angles, and dihedral angles) from unbiased all-atom explicit-solvent MD simulations of short GAG polymers to rapidly construct models of GAGs of arbitrary length. The algorithm was used to generate non-sulfated chondroitin 10- and 20-mer ensembles which were compared to MD-generated ensembles for internal validation. End-to-end distance distributions in constructed and MD-generated ensembles have minimal differences, suggesting that our algorithm produces conformational ensembles that mimic the backbone flexibility seen in simulation. Non-sulfated chondroitin 100- and 200-mer ensembles were constructed within a day, demonstrating the efficiency of the algorithm and reduction in time and computational cost compared to simulation.

Implicit Micelle Model for Membrane Proteins Using Superellipsoid Approximation.

Surfactant micelles are often utilized as membrane mimetics for structure determination and functional analysis of membrane proteins. The curved-surface effects of the micelle can perturb membrane protein structure. However, it is difficult to assess such effects and membrane mimetic artifacts by experimental and theoretical methods. Here, we propose an implicit micelle model (IMIC) to be used in molecular dynamics (MD) simulations of membrane proteins. IMIC is an extension of the IMM1 implicit membrane model and additionally introduces a superellipsoid approximation to represent the curved-surface effects. Most of the IMIC parameters are obtained from all-atom explicit solvent MD simulations of 12 membrane proteins in various micelles. The HIV envelope protein gp41, M13 major coat protein gp8, and amyloid precursor protein (APP) dimer are simulated via MD simulations with IMIC. These simulations clearly show how the micelle influences membrane protein structures compared to the bilayer environments. The MD simulations with IMIC provide reliable membrane protein structures in various micelle environments quickly with smaller computational cost than that for an explicit solvent/micelle model.

TRNA Dissociation from EF-Tu after GTP Hydrolysis: Primary Steps and Antibiotic Inhibition

In each round of ribosomal translation, the translational GTPase elongation factor Tu (EF-Tu) delivers a transfer RNA (tRNA) to the ribosome. After successful decoding, EF-Tu hydrolyzes GTP, which triggers a conformational change that ultimately results in the release of the tRNA from EF-Tu. To identify the primary steps of these conformational changes and how they are prevented by the antibiotic kirromycin, we employed all-atom explicit-solvent molecular dynamics simulations of the full ribosome-EF-Tu complex. Our results suggest that after GTP hydrolysis and Pi release, the loss of interactions between the nucleotide and the switch 1 loop of EF-Tu allows domain D1 of EF-Tu to rotate relative to domains D2 and D3 and leads to an increased flexibility of the switch 1 loop. This rotation induces a closing of the D1-D3 interface and an opening of the D1-D2 interface. We propose that the opening of the D1-D2 interface, which binds the CCA tail of the tRNA, weakens the crucial EF-Tu-tRNA interactions, which lowers tRNA binding affinity, representing the first step of tRNA release. Kirromycin binds within the D1-D3 interface, sterically blocking its closure, but does not prevent hydrolysis. The resulting increased flexibility of switch 1 explains why it is not resolved in kirromycin-bound structures.

Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides.

Chondroitin sulfate (CS) is one of several glycosaminoglycans that are major components of proteoglycans. A linear polymer consisting of repeats of the disaccharide -4GlcAβ1-3GalNAcβ1-, CS undergoes differential sulfation resulting in five unique sulfation patterns. Because of the dimer repeat, the CS glycosidic "backbone" has two distinct sets of conformational degrees of freedom defined by pairs of dihedral angles: (ϕ1, ψ1) about the β1-3 glycosidic linkage and (ϕ2, ψ2) about the β1-4 glycosidic linkage. Differential sulfation and the possibility of cation binding, combined with the conformational flexibility and biological diversity of CS, complicate experimental efforts to understand CS three-dimensional structures at atomic resolution. Therefore, all-atom explicit-solvent molecular dynamics simulations with Adaptive Biasing Force sampling of the CS backbone were applied to obtain high-resolution, high-precision free energies of CS disaccharides as a function of all possible backbone geometries. All 10 disaccharides (β1-3 vs β1-4 linkage × five different sulfation patterns) were studied; additionally, ion effects were investigated by considering each disaccharide in the presence of either neutralizing sodium or calcium cations. GlcAβ1-3GalNAc disaccharides have a single, broad, thermodynamically important free-energy minimum, whereas GalNAcβ1-4GlcA disaccharides have two such minima. Calcium cations but not sodium cations bind to the disaccharides, and binding is primarily to the GlcA -COO(-) moiety as opposed to sulfate groups. This binding alters the glycan backbone thermodynamics in instances where a calcium cation bound to -COO(-) can act to bridge and stabilize an interaction with an adjacent sulfate group, whereas, in the absence of this cation, the proximity of a sulfate group to -COO(-) results in two like charges being both desolvated and placed adjacent to each other and is found to be destabilizing. In addition to providing information on sulfation and cation effects, the present results can be applied to building models of CS polymers and as a point of comparison in studies of CS polymer backbone dynamics and thermodynamics.

Parametrization of Backbone Flexibility in a Coarse-Grained Force Field for Proteins (COFFDROP) Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of All Possible Two-Residue Peptides.

Recently, we reported the parametrization of a set of coarse-grained (CG) nonbonded potential functions, derived from all-atom explicit-solvent molecular dynamics (MD) simulations of amino acid pairs and designed for use in (implicit-solvent) Brownian dynamics (BD) simulations of proteins; this force field was named COFFDROP (COarse-grained Force Field for Dynamic Representations Of Proteins). Here, we describe the extension of COFFDROP to include bonded backbone terms derived from fitting to results of explicit-solvent MD simulations of all possible two-residue peptides containing the 20 standard amino acids, with histidine modeled in both its protonated and neutral forms. The iterative Boltzmann inversion (IBI) method was used to optimize new CG potential functions for backbone-related terms by attempting to reproduce angle, dihedral, and distance probability distributions generated by the MD simulations. In a simple test of the transferability of the extended force field, the angle, dihedral, and distance probability distributions obtained from BD simulations of 56 three-residue peptides were compared to results from corresponding explicit-solvent MD simulations. In a more challenging test of the COFFDROP force field, it was used to simulate eight intrinsically disordered proteins and was shown to quite accurately reproduce the experimental hydrodynamic radii (Rhydro), provided that the favorable nonbonded interactions of the force field were uniformly scaled downward in magnitude. Overall, the results indicate that the COFFDROP force field is likely to find use in modeling the conformational behavior of intrinsically disordered proteins and multidomain proteins connected by flexible linkers.

COFFDROP: A Coarse-Grained Nonbonded Force Field for Proteins Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of Amino Acids.

We describe the derivation of a set of bonded and nonbonded coarse-grained (CG) potential functions for use in implicit-solvent Brownian dynamics (BD) simulations of proteins derived from all-atom explicit-solvent molecular dynamics (MD) simulations of amino acids. Bonded potential functions were derived from 1 μs MD simulations of each of the 20 canonical amino acids, with histidine modeled in both its protonated and neutral forms; nonbonded potential functions were derived from 1 μs MD simulations of every possible pairing of the amino acids (231 different systems). The angle and dihedral probability distributions and radial distribution functions sampled during MD were used to optimize a set of CG potential functions through use of the iterative Boltzmann inversion (IBI) method. The optimized set of potential functions—which we term COFFDROP (COarse-grained Force Field for Dynamic Representation Of Proteins)—quantitatively reproduced all of the “target” MD distributions. In a first test of the force field, it was used to predict the clustering behavior of concentrated amino acid solutions; the predictions were directly compared with the results of corresponding all-atom explicit-solvent MD simulations and found to be in excellent agreement. In a second test, BD simulations of the small protein villin headpiece were carried out at concentrations that have recently been studied in all-atom explicit-solvent MD simulations by Petrov and Zagrovic (PLoS Comput. Biol.2014, 5, e1003638). The anomalously strong intermolecular interactions seen in the MD study were reproduced in the COFFDROP simulations; a simple scaling of COFFDROP’s nonbonded parameters, however, produced results in better accordance with experiment. Overall, our results suggest that potential functions derived from simulations of pairwise amino acid interactions might be of quite broad applicability, with COFFDROP likely to be especially useful for modeling unfolded or intrinsically disordered proteins.

Atomic Resolution Mechanism of Cd44-Hyaluronan Protein-Carbohydrate Binding

All-atom explicit-solvent molecular dynamics simulations with advanced sampling methods were applied to completely determine, at atomic resolution, the mechanism of binding of the protein CD44 with the glycosaminoglycan carbohydrate hyaluronan. The results include descriptions of: formation of the protein-carbohydrate complex from separate protein and carbohydrate; conformational switching at the binding site to enhance binding affinity; additional binding affinity enhancement through an order-to-disorder transition in the protein; and inhibition of binding by glycosylation of CD44. These computational results are discussed in the context of and illuminate existing experimental data.

Conformational and Dynamic Properties of Extracellular Domains of Cell Adhesion Molecules

Due to their specific binding properties, cell adhesion molecules (CAMs), such as integrin, cadherin and the immunoglobulin superfamily CAMs are of primary importance in cell-cell and cell-substrate interactions, signaling pathways and other crucial biological processes. We investigate the molecular structures and conformational dynamics of the extracellular domains of the Aplysia CAM (apCAM), by constructing and studying their atomically detailed structural model based on homology with human neural CAM (NCAM). The stability and dynamic properties of the apCAM domains, individually and as a connected in the molecule, are probed and analyzed using all-atom explicit-solvent molecular dynamics simulations. The refined structural model reveals a specific pattern of amino acid interactions that controls the stability of their beta-sheet reach structures, and which could affect their specific response to environmental changes.

Protein High-Force Pulling Simulations Yield Low-Force Results

All-atom explicit-solvent molecular dynamics simulations are used to pull with extremely large constant force (750–3000 pN) on three small proteins. The introduction of a nondimensional timescale permits direct comparison of unfolding across all forces. A crossover force of approximately 1100 pN divides unfolding dynamics into two regimes. At higher forces, residues sequentially unfold from the pulling end while maintaining the remainder of the protein force-free. Measurements of hydrodynamic viscous stresses are made easy by the high speeds of unfolding. Using an exact low-Reynolds-number scaling, these measurements can be extrapolated to provide, for the first time, an estimate of the hydrodynamic force on low-force unfolding. Below 1100 pN, but surprisingly still at extremely large applied force, intermediate states and cooperative unfoldings as seen at much lower forces are observed. The force-insensitive persistence of these structures indicates that decomposition into unfolded fragments requires a large fluctuation. This finding suggests how proteins are constructed to resist transient high force. The progression of helix and sheet unfolding is also found to be insensitive to force. The force-insensitivity of key aspects of unfolding opens the possibility that numerical simulations can be accelerated by high applied force while still maintaining critical features of unfolding.

Binding Site Flexibility in the CD44:Hyaluronan Protein:Carbohydrate Interaction

The CD44:hyaluronan protein:carbohydrate binding interaction is an important component of inflammation and cancer cell invasiveness and metastasis. Association of a key arginine sidechain in the beta1-alpha1 loop of the binding site with the carbohydrate ligand has been postulated to lead to a high-affinity state of the complex based on atomic-resolution x-ray crystal studies. NMR studies point to coupling between an order-to-disorder transition of the beta9-alpha3 region of the CD44 protein, which is distant from the binding site, and the formation of the high-affinity state; however, this region of the protein is well-ordered in all crystal structures. All-atom explicit-solvent molecular dynamics simulations of the CD44:hyaluronan complex revealed beta1-alpha1 loop backbone instability when arginine is in close association with hyaluronan, whereas the loop was stable in the absence of sidechain association. Additionally, a conformational transition in the backbone phi angle of a single residue in the beta1-alpha1 loop was demonstrated to be sufficient to induce reversible arginine sidechain:hyaluronan association, leading to the description of a molecular switching mechanism. The instability of the beta1-alpha1 loop in the context of a well-ordered beta9-alpha3 region suggests possible allosteric coupling between the binding site conformation and the order-to-disorder transition.

Mechanism of Binding Site Conformational Switching in the CD44–Hyaluronan Protein–Carbohydrate Binding Interaction

The transmembrane protein CD44, which has been implicated in cancer biology and inflammation, mediates cell adhesion through multimeric interactions with the linear extracellular glycosaminoglycan hyaluronan (HA; in megadaltons). Affinity switching of CD44 from a low-affinity state to a high-affinity state is required for normal CD44 physiological function; crystal structures of the CD44 hyaluronan binding domain complexed with HA oligomers point to a conformational rearrangement at a binding site loop, leading to the formation of direct contact between the oligomer and an arginine side chain as a molecular basis for affinity switching. Here, all-atom explicit-solvent molecular dynamics simulations were used to characterize the dynamics and thermodynamics of oligomeric hyaluronan (oHA) and its two crystallographic complexes with the CD44 hyaluronan binding domain: the “A-form,” which lacks arginine–HA close contact, and the “B-form,” which has direct arginine side-chain–HA contact. From the simulations, the conformational properties of oHA are essentially unaltered in going from the unbound state to either the A-form or the B-form bound state, with the oligomer retaining its flexibility when bound and with only two of the eight monosaccharides in the oligomer maintaining uninterrupted contact with the protein. Biased simulations revealed that altering the backbone conformation of a tyrosine residue in the arginine loop can induce the A-form → B-form conformational transition and that a large free-energy barrier prevents ready interconversion between the two forms, thereby suggesting that the tyrosine backbone forms a molecular switch.