Calculation of hydrodynamic properties of globular proteins from their atomic-level structure

Quantitative analysis of the slow exchange process by 19F NMR in the presence of scalar and dipolar couplings: applications to the ribose 2'-19F probe in nucleic acids

Solution NMR spectroscopy is a particularly powerful technique for characterizing the functional dynamics of biomolecules, which is typically achieved through the quantitative characterization of chemical exchange processes via the measurement of spin relaxation rates. In addition to the conventional nuclei such as 15N and 13C, which are abundant in biomolecules, fluorine-19 (19F) has recently garnered attention and is being widely used as a site-specific spin probe. While 19F offers the advantages of high sensitivity and low background, it can be susceptible to artifacts in quantitative relaxation analyses due to a multitude of dipolar and scalar coupling interactions with nearby 1H spins. In this study, we focused on the ribose 2'-19F spin probe in nucleic acids and investigated the effects of 1H-19F spin interactions on the quantitative characterization of slow exchange processes on the millisecond time scale. We demonstrated that the 1H-19F dipolar coupling can significantly affect the interpretation of 19F chemical exchange saturation transfer (CEST) experiments when 1H decoupling is applied, while the 1H-19F interactions have a lesser impact on Carr-Purcell-Meiboom-Gill relaxation dispersion applications. We also proposed a modified CEST scheme to alleviate these artifacts along with experimental verifications on self-complementary RNA systems. The theoretical framework presented in this study can be widely applied to various 19F spin systems where 1H-19F interactions are operative, further expanding the utility of 19F relaxation-based NMR experiments.

NMR characterization of RNA binding property of the DEAD-box RNA helicase DDX3X and its implications for helicase activity

The DEAD-box RNA helicase (DDX) plays a central role in many aspects of RNA metabolism by remodeling the defined structure of RNA molecules. While a number of structural studies have revealed the atomistic details of the interaction between DDX and RNA ligands, the molecular mechanism of how this molecule unwinds a structured RNA into an unstructured single-stranded RNA (ssRNA) has largely remained elusive. This is due to challenges in structurally characterizing the unwinding intermediate state and the lack of thermodynamic details underlying this process. In this study, we use solution nuclear magnetic resonance (NMR) spectroscopy to characterize the interaction of human DDX3X, a member of the DDX family, with various RNA ligands. Our results show that the inherent binding affinity of DDX3X for ssRNA is significantly higher than that for structured RNA elements. This preferential binding, accompanied by the formation of a domain-closed conformation in complex with ssRNA, effectively stabilizes the denatured ssRNA state and thus underlies the unwinding activity of DDX3X. Our results provide a thermodynamic and structural basis for the DDX function, whereby DDX can recognize and remodel a distinct set of structured RNAs to participate in a wide range of physiological processes.

Effects of Aqueous Isotopic Substitution on the Adsorption Dynamics and Dilational Rheology of β-Lactoglobulin Layers at the Water/Air Interface

The effect of the degree of isotopic substitution of the aqueous medium on the adsorption kinetics and the surface dilational rheological behavior at the water/air interface of the globular protein β-lactoglobulin was investigated. Aqueous solutions with fixed concentrations of 1 μM protein and 10 mM hydrogenous buffer with controlled pH 7 were prepared in H2O, D2O, and an isotopic mixture of 8.1% v/v D2O in H2O (called air contrast matched water, ACMW). Using a bubble shape analysis tensiometer, we obtained various experimental dependencies of the dilational viscoelasticity modulus E as a function of the dynamic surface pressure and of the frequency and amplitude of bubble surface area oscillations, either in the course of adsorption or after having reached a steady state. In general, the results revealed virtually no effect from substituting H2O by ACMW but distinct albeit relatively weak effects for intermediate adsorption times for D2O as the aqueous phase. In the final stage of adsorption, established after around 10 h, the equilibrium adsorption and the dilational rheological behavior of all protein layers under investigation are only very weakly affected by the presence of D2O. The obtained results help to design experimental protocols for protein adsorption studies, for example, by neutron reflectivity.

Probing effects of site-specific aspartic acid isomerization on structure and stability of GB1 through chemical protein synthesis

Chemical modifications of long-lived proteins, such as isomerization and epimerization, have been evoked as prime triggers for protein-damage related diseases. Deamidation of Asn residues, which results in formation of a mixture of l- and d-Asp and isoAsp via an intermediate aspartyl succinimide, can result in the disruption of cellular proteostasis and toxic protein depositions. In contrast to extensive data on the biological prevalence and functional implications of aspartyl succinimide formation, much less is known about the impact of the resulting altered backbone composition on properties of individual proteins at a molecular level. Here, we report the total chemical synthesis, biophysical characterization, and NMR structural analysis of a series of variants of the B1 domain of protein G from Streptococcal bacteria (GB1) in which all possible Asp isomers as well as an aspartyl succinimide were individually incorporated at a defined position in a solvent-exposed loop. Subtle local structural effects were observed; however, these were accompanied by notable differences in thermodynamic folded stability. Surprisingly, the noncanonical backbone connectivity of d-isoAsp led to a variant that exhibited enhanced stability relative to the natural protein.

Enhanced TROSY Effect in [2-19 F, 2-13 C] Adenosine and ATP Analogs Facilitates NMR Spectroscopy of Very Large Biological RNAs in Solution

Large RNAs are central to cellular functions, but characterizing such RNAs remains challenging by solution NMR. We present two labeling technologies based on [2-19 F, 2-13 C]-adenosine, which allow the incorporation of aromatic 19 F-13 C spin pairs. The labels when coupled with the transverse relaxation optimized spectroscopy (TROSY) enable us to probe RNAs comprising up to 124 nucleotides. With our new [2-19 F, 2-13 C]-adenosine-phosphoramidite, all resonances of the human hepatitis B virus epsilon RNA could be readily assigned. With [2-19 F, 2-13 C]-adenosine triphosphate, the 124 nt pre-miR-17-NPSL1-RNA was produced via in vitro transcription and the TROSY spectrum of this 40 kDa [2-19 F, 2-13 C]-A-labeled RNA featured sharper resonances than the [2-1 H, 2-13 C]-A sample. The mutual cancelation of the chemical-shift-anisotropy and the dipole-dipole-components of TROSY-resonances leads to narrow linewidths over a wide range of molecular weights. With the synthesis of a non-hydrolysable [2-19 F, 2-13 C]-adenosine-triphosphate, we facilitate the probing of co-factor binding in kinase complexes and NMR-based inhibitor binding studies in such systems. Our labels allow a straightforward assignment for larger RNAs via a divide-and-conquer/mutational approach. The new [2-19 F, 2-13 C]-adenosine precursors are a valuable addition to the RNA NMR toolbox and will allow the study of large RNAs/RNA protein complexes in vitro and in cells.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Structure, Topology, and Stability of Multiple G-quadruplexes in Long Telomeric Overhangs

Telomeres and their single stranded overhangs gradually shorten with successive cell divisions, as part of the natural aging process, but can be elongated by telomerase, a nucleoprotein complex which is activated in the majority of cancers. This prominent implication in cancer and aging has made the repetitive telomeric sequences (TTAGGG repeats) and the G-quadruplex structures that form in their overhangs the focus of intense research in the past several decades. However, until recently most in vitro efforts to understand the structure, stability, dynamics, and interactions of telomeric overhangs had been focused on short sequences that are not representative of longer sequences encountered in a physiological setting. In this review, we will provide a broad perspective about telomeres and associated factors, and introduce the agents and structural characteristics involved in organizing, maintaining, and protecting telomeric DNA. We will also present a summary of recent research performed on long telomeric sequences, nominally defined as those that can form two or more tandem G-quadruplexes, i.e., which contain eight or more TTAGGG repeats. Results of experimental studies using a broad array of experimental tools, in addition to recent computational efforts will be discussed, particularly in terms of their implications for the stability, folding topology, and compactness of the tandem G-quadruplexes that form in long telomeric overhangs.

First coarse grain then scale: How to estimate diffusion coefficients of confined molecules

An approach for approximating position and orientation dependent translational and rotational diffusion coefficients of rigid molecules of any shape suspended in a viscous fluid under geometric confinement is proposed. It is an extension of the previously developed scheme for evaluating near-wall diffusion of macromolecules, now applied to any geometry of boundaries. The method relies on shape based coarse-graining combined with scaling of mobility matrix components by factors derived based on energy dissipation arguments for Stokes flows. Tests performed for a capsule shaped molecule and its coarse-grained model, a dumbbell, for three different types of boundaries (a sphere, an open cylinder, and two parallel planes) are described. An almost perfect agreement between mobility functions of the detailed and coarse-grained models, even close to boundary surfaces, is obtained. The proposed method can be used to simplify hydrodynamic calculations and reduce errors introduced due to coarse-graining of molecular shapes.

Protein translational diffusion as a way to detect intermolecular interactions

In this work, we analyze the information on the protein intermolecular interactions obtained from macromolecular diffusion. We have shown that the most hopeful results are given by our approach based on analysis of protein translational self-diffusion and collective diffusion obtained by dynamic light scattering and pulsed-field gradient NMR (PFG NMR) spectroscopy with the help of Vink's approach to analyze diffusion motion of particles by frictional formalism of non-equilibrium thermodynamics and the usage of the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloid particles interactions in electrolyte solutions. Early we have shown that integration of Vink's theory with DLVO provides a reliable basis for uniform interpreting of PFG NMR and DLS experiments on concentration dependence of diffusion coefficients. Basic details of theoretical and mathematical procedures and a broad analysis of experimental attestation of proposed conception on proteins of various structural form, size, and shape are presented. In the present review, the main capabilities of our approach obtain the details of intermolecular interactions of proteins with different shapes, internal structures, and mass. The universality of Vink's approach is experimentally shown, which gives the appropriate description of experimental results for proteins of complicated structure and shape.

Unveiling protein-protein interaction potential through Monte Carlo simulation combined with small-angle X-ray scattering

Protein interactions are investigated under different conditions of lysozyme concentration, temperature and ionic strength by means of in-solution small angle X-Ray scattering (SAXS) experiments and Monte Carlo (MC) simulations. Initially, experimental data were analysed through a Hard-Sphere Double Yukawa (HSDY) model combined with Random Phase Approximation (RPA), a closure relationship commonly used in the literature for monodisperse systems. We realized by means of MC that the HSDY/RPA modelling fails to describe the protein-protein pair potential for moderated and dense systems at low ionic strength, mainly due to inherent distortions of the RPA approximation. Our SAXS/MC results thus show that lysozyme concentrations between 2 (diluted) and 20 mg/mL (not crowded) present similar protein-protein pair potential preserving the values of surface net charge around 7 e, protein diameter of 28 Å, decay range of attractive well potential of 3 Å and a depth of the well potential varying from 1 to 5 kBT depending on temperature and salt addition. Noteworthy, we here propose a novel method to analyse the SAXS data from interacting proteins through MC simulations, which overcomes the deficiencies presented by the use of a closure relationship. Furthermore, this new methodology of combining SAXS with MC simulations gives a step forward to investigate more complex systems as those composed of a mixture of proteins of distinct species presenting different molecular weights (and hence sizes) and surface net charges at low, moderate and very dense systems.

Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides a Fast but Still Accurate Treatment of Hydrodynamic Interactions in Brownian Dynamics Simulations of Biological Macromolecules

The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HIs). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here, we explore the use of an alternative way to accelerate the calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA-RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA-RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast, approximate way of including HIs in BD simulations of intermediate-scale systems.

Modelling red blood cell optical trapping by machine learning improved geometrical optics calculations

Optically trapping red blood cells allows for the exploration of their biophysical properties, which are affected in many diseases. However, because of their nonspherical shape, the numerical calculation of the optical forces is slow, limiting the range of situations that can be explored. Here we train a neural network that improves both the accuracy and the speed of the calculation and we employ it to simulate the motion of a red blood cell under different beam configurations. We found that by fixing two beams and controlling the position of a third, it is possible to control the tilting of the cell. We anticipate this work to be a promising approach to study the trapping of complex shaped and inhomogeneous biological materials, where the possible photodamage imposes restrictions in the beam power.

Solution structure, dynamics and tetrahedral assembly of Anti-TRAP, a homo-trimeric triskelion-shaped regulator of tryptophan biosynthesis in Bacillus subtilis

Cellular production of tryptophan is metabolically expensive and tightly regulated. The small Bacillus subtilis zinc binding Anti-TRAP protein (AT), which is the product of the yczA/rtpA gene, is upregulated in response to accumulating levels of uncharged tRNATrp through a T-box antitermination mechanism. AT binds to the undecameric ring-shaped protein TRAP (trp RNA Binding Attenuation Protein), thereby preventing it from binding to the trp leader RNA. This reverses the inhibitory effect of TRAP on transcription and translation of the trp operon. AT principally adopts two symmetric oligomeric states, a trimer (AT3) featuring a three-helix bundle, or a dodecamer (AT12) comprising a tetrahedral assembly of trimers, whereas only the trimeric form has been shown to bind and inhibit TRAP. We demonstrate the utility of native mass spectrometry (nMS) and small-angle x-ray scattering (SAXS), together with analytical ultracentrifugation (AUC) for monitoring the pH and concentration-dependent equilibrium between the trimeric and dodecameric structural forms of AT. In addition, we report the use of solution nuclear magnetic resonance (NMR) spectroscopy to determine the solution structure of AT3, while heteronuclear 15N relaxation measurements on both oligomeric forms of AT provide insights into the dynamic properties of binding-active AT3 and binding-inactive AT12, with implications for TRAP inhibition.

An overview of descriptors to capture protein properties - Tools and perspectives in the context of QSAR modeling

Proteins are important ingredients in food and feed, they are the active components of many pharmaceutical products, and they are necessary, in the form of enzymes, for the success of many technical processes. However, production can be challenging, especially when using heterologous host cells such as bacteria to express and assemble recombinant mammalian proteins. The manufacturability of proteins can be hindered by low solubility, a tendency to aggregate, or inefficient purification. Tools such as in silico protein engineering and models that predict separation criteria can overcome these issues but usually require the complex shape and surface properties of proteins to be represented by a small number of quantitative numeric values known as descriptors, as similarly used to capture the features of small molecules. Here, we review the current status of protein descriptors, especially for application in quantitative structure activity relationship (QSAR) models. First, we describe the complexity of proteins and the properties that descriptors must accommodate. Then we introduce descriptors of shape and surface properties that quantify the global and local features of proteins. Finally, we highlight the current limitations of protein descriptors and propose strategies for the derivation of novel protein descriptors that are more informative.

Disulfide-Bond-Induced Structural Frustration and Dynamic Disorder in a Peroxiredoxin from MAS NMR

Disulfide bond formation is fundamentally important for protein structure and constitutes a key mechanism by which cells regulate the intracellular oxidation state. Peroxiredoxins (PRDXs) eliminate reactive oxygen species such as hydrogen peroxide through a catalytic cycle of Cys oxidation and reduction. Additionally, upon Cys oxidation PRDXs undergo extensive conformational rearrangements that may underlie their presently structurally poorly defined functions as molecular chaperones. Rearrangements include high molecular-weight oligomerization, the dynamics of which are, however, poorly understood, as is the impact of disulfide bond formation on these properties. Here we show that formation of disulfide bonds along the catalytic cycle induces extensive μs time scale dynamics, as monitored by magic-angle spinning NMR of the 216 kDa-large Tsa1 decameric assembly and solution-NMR of a designed dimeric mutant. We ascribe the conformational dynamics to structural frustration, resulting from conflicts between the disulfide-constrained reduction of mobility and the desire to fulfill other favorable contacts.

An Orientationally Averaged Version of the Rotne-Prager-Yamakawa Tensor Provides A Fast But Still Accurate Treatment Of Hydrodynamic Interactions In Brownian Dynamics Simulations Of Biological Macromolecules

The Brownian dynamics (BD) simulation technique is widely used to model the diffusive and conformational dynamics of complex systems comprising biological macromolecules. For the diffusive properties of macromolecules to be described correctly by BD simulations, it is necessary to include hydrodynamic interactions (HI). When modeled at the Rotne-Prager-Yamakawa (RPY) level of theory, for example, the translational and rotational diffusion coefficients of isolated macromolecules can be accurately reproduced; when HIs are neglected, however, diffusion coefficients can be underestimated by an order of magnitude or more. The principal drawback to the inclusion of HIs in BD simulations is their computational expense, and several previous studies have sought to accelerate their modeling by developing fast approximations for the calculation of the correlated random displacements. Here we explore the use of an alternative way to accelerate calculation of HIs, i.e., by replacing the full RPY tensor with an orientationally averaged (OA) version which retains the distance dependence of the HIs but averages out their orientational dependence. We seek here to determine whether such an approximation can be justified in application to the modeling of typical proteins and RNAs. We show that the use of an OA RPY tensor allows translational diffusion of macromolecules to be modeled with very high accuracy at the cost of rotational diffusion being underestimated by ∼25%. We show that this finding is independent of the type of macromolecule simulated and the level of structural resolution employed in the models. We also show, however, that these results are critically dependent on the inclusion of a non-zero term that describes the divergence of the diffusion tensor: when this term is omitted from simulations that use the OA RPY model, unfolded macromolecules undergo rapid collapse. Our results indicate that the orientationally averaged RPY tensor is likely to be a useful, fast approximate way of including HIs in BD simulations of intermediate-scale systems.

A conserved oligomerization domain in the disordered linker of coronavirus nucleocapsid proteins

The nucleocapsid (N-)protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a key role in viral assembly and scaffolding of the viral RNA. It promotes liquid-liquid phase separation (LLPS), forming dense droplets that support the assembly of ribonucleoprotein particles with as-of-yet unknown macromolecular architecture. Combining biophysical experiments, molecular dynamics simulations, and analysis of the mutational landscape, we describe a heretofore unknown oligomerization site that contributes to LLPS, is required for the assembly of higher-order protein-nucleic acid complexes, and is coupled to large-scale conformational changes of N-protein upon nucleic acid binding. The self-association interface is located in a leucine-rich sequence of the intrinsically disordered linker between N-protein folded domains and formed by transient helices assembling into trimeric coiled-coils. Critical residues stabilizing hydrophobic and electrostatic interactions between adjacent helices are highly protected against mutations in viable SARS-CoV-2 genomes, and the oligomerization motif is conserved across related coronaviruses, thus presenting a target for antiviral therapeutics.

Translational and rotational diffusion of short ribonucleic acids

Inevitable precondition for ribonucleic acids to regulate gene expression and to perform gene editing is diffusion. Free three-dimensional translational diffusion velocity of RNA of up to 200 nucleotides could be predicted with high accuracy by the empirical model D = 4.58 10^-10 N^-0.39 m²s^-1. Furthermore, the biological function of ribonucleic acids is determined by rotational diffusion. In the presented work, an empirical model is derived applying atom-level shell-modeling of electron density maps, D_r = 1.62 10⁹ N^-1.20 s^-1, to predict the rotational diffusion coefficient of short ribonucleic acids based on the polymer size.

Therapeutic phosphorodiamidate morpholino oligonucleotides: Physical properties, solution structures, and folding thermodynamics

Elucidating the structure-function relationships for therapeutic RNA mimicking phosphorodiamidate morpholino oligonucleotides (PMOs) is challenging due to the lack of information about their structures. While PMOs have been approved by the US Food and Drug Administration for treatment of Duchenne muscular dystrophy, no structural information on these unique, charge-neutral, and stable molecules is available. We performed circular dichroism and solution viscosity measurements combined with molecular dynamics simulations and machine learning to resolve solution structures of 22-mer, 25-mer, and 30-mer length PMOs. The PMO conformational dynamics are defined by the competition between non-polar nucleobases and uncharged phosphorodiamidate groups for shielding from solvent exposure. PMO molecules form non-canonical, partially helical, stable folded structures with a small 1.4- to 1.7-nm radius of gyration, low count of three to six base pairs and six to nine base stacks, characterized by -34 to -51 kcal/mol free energy, -57 to -103 kcal/mol enthalpy, and -23 to -53 kcal/mol entropy for folding. The 4.5- to 6.2-cm³/g intrinsic viscosity and Huggins constant of 4.5-9.9 are indicative of extended and aggregating systems. The results obtained highlight the importance of the conformational ensemble view of PMO solution structures, thermodynamic stability of their non-canonical structures, and concentration-dependent viscosity properties. These principles form a paradigm to understand the structure-properties-function relationship for therapeutic PMOs to advance the design of new RNA-mimic-based drugs.

Trimeric Architecture Ensures the Stability and Biological Activity of the Calf Purine Nucleoside Phosphorylase: In Silico and In Vitro Studies of Monomeric and Trimeric Forms of the Enzyme

Mammalian purine nucleoside phosphorylase (PNP) is biologically active as a homotrimer, in which each monomer catalyzes a reaction independently of the others. To answer the question of why the native PNP forms a trimeric structure, we constructed, in silico and in vitro, the monomeric form of the enzyme. Molecular dynamics simulations showed different geometries of the active site in the non-mutated trimeric and monomeric PNP forms, which suggested that the active site in the isolated monomer could be non-functional. To confirm this hypothesis, six amino acids located at the interface of the subunits were selected and mutated to alanines to disrupt the trimer and obtain a monomer (6Ala PNP). The effects of these mutations on the enzyme structure, stability, conformational dynamics, and activity were examined. The solution experiments confirmed that the 6Ala PNP mutant occurs mainly as a monomer, with a secondary structure almost identical to the wild type, WT PNP, and importantly, it shows no enzymatic activity. Simulations confirmed that, although the secondary structure of the 6Ala monomer is similar to the WT PNP, the positions of the amino acids building the 6Ala PNP active site significantly differ. These data suggest that a trimeric structure is necessary to stabilize the geometry of the active site of this enzyme.

Less is more: A simple methyl-TROSY based pulse scheme offers improved sensitivity in applications to high molecular weight complexes

The HMQC pulse sequence and variants thereof have been exploited in studies of high molecular weight protein complexes, taking advantage of the fact that fast and slow relaxing magnetization components are sequestered along two distinct magnetization transfer pathways. Despite the simplicity of the HMQC scheme an even shorter version can be designed, based on elimination of the terminal refocusing period, as a further means of increasing signal. Here we present such an experiment, and show that significant sensitivity gains, in some cases by factors of two or more, are realized in studies of proteins varying in molecular masses from 100 kDa to 1 MDa.

Effects of model membranes on lysozyme amyloid aggregation

The study of the interaction between lipid membranes and amyloidogenic peptides is a turning point for understanding the processes involving the cytotoxicity of peptides involved in neurodegenerative diseases. In this work, we perform an experimental study of model membrane-lysozyme interaction to understand how the formation of amyloid fibrils can be affected by the presence of polar and zwitterionic phospholipid molecules (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine [POPC] and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol [POPG]). The study was conducted above and below the critical micellar concentration (CMC) using dynamic light scattering (DLS), atomic force microscopy (AFM), UV-Vis spectrophotometry, and the quartz crystal microbalance (QCM). Our results show that the presence of phospholipids appears to be a factor favoring the formation of amyloid aggregates. Spectrophotometric and DLS data revealed that the quantity of β -structure increases in the presence of POPG and POPC at different concentrations. The presence of POPG and POPC increases the speed of the nucleation process, without altering the overall structures of the fibrillar final products.

Predicting molecular properties of α-synuclein using force fields for intrinsically disordered proteins

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 μs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 μs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and ³ J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.

RNA double strand hybridization measured at the single molecule level

RNA double strand hybridization is a hallmark for gene expression regulation. In this function, single stranded regulatory RNA forms Watson-Crick base pairs with complementary messenger RNA. In the presented work the dissociation constants of complementary equally sized RNA single strands were measured at the single molecule level applying fluorescence correlation spectroscopy (FCS). Dissociation constants of 3.2 nM, 1.4 nM and 1.0 nM were determined for 26 bp, 41 bp and 54 bp dsRNA, respectively. The translational diffusion coefficients of RNA, measured at infinite dilution, could be accurately predicted applying the model D = 4.58 × 10^-10 N^-0.39 m²s^-1.

Diffusion of small molecule drugs is affected by surface interactions and crowder proteins

Crowded environments are known to affect the diffusion of macromolecules, but their effects on the diffusion of small molecules are largely uncharacterized. We investigate how three protein crowders, bovine serum albumin (BSA), hen egg-white lysozyme, and myoglobin, influence the diffusion rates and interactions of four small molecules: fluorescein, and three drugs, doxorubicin, glycogen synthase kinase-3 inhibitor SB216763, and quinacrine. Using Line-FRAP measurements, Brownian dynamics simulations, and molecular docking, we find that the diffusion rates of the small molecules are highly affected by self-aggregation, interactions with the proteins, and surface adsorption. The diffusion of fluorescein is decreased because of its interactions with the protein crowders and their surface adsorption. Protein crowders increase the diffusion rates of doxorubicin and SB216763 by reducing surface interactions and self-aggregation, respectively. Quinacrine diffusion was not affected by protein crowders. The mechanistic insights gained here may assist in optimization of compounds for higher mobility in complex macromolecular environments.

Convergent views on disordered protein dynamics from NMR and computational approaches

Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.

Resolving molecular diffusion and aggregation of antibody proteins with megahertz X-ray free-electron laser pulses

X-ray free-electron lasers (XFELs) with megahertz repetition rate can provide novel insights into structural dynamics of biological macromolecule solutions. However, very high dose rates can lead to beam-induced dynamics and structural changes due to radiation damage. Here, we probe the dynamics of dense antibody protein (Ig-PEG) solutions using megahertz X-ray photon correlation spectroscopy (MHz-XPCS) at the European XFEL. By varying the total dose and dose rate, we identify a regime for measuring the motion of proteins in their first coordination shell, quantify XFEL-induced effects such as driven motion, and map out the extent of agglomeration dynamics. The results indicate that for average dose rates below 1.06 kGy μs^-1 in a time window up to 10 μs, it is possible to capture the protein dynamics before the onset of beam induced aggregation. We refer to this approach as correlation before aggregation and demonstrate that MHz-XPCS bridges an important spatio-temporal gap in measurement techniques for biological samples.

Poly(amidoamine) dendrimer immunosensor for ultrasensitive gravimetric and electrochemical detection of matrix metalloproteinase-9

Monitoring the level of matrix metalloproteinase-9 (MMP-9) and inhibiting its expression is important for the diagnosis and treatment of various diseases. However, the analysis of MMP-9 is challenging owing to its very low content in the blood, especially at the early stages of diseases. Therefore, we developed an ultrasensitive and easy-to-use immunosensor based on a three-dimensional (3D) bioplatform for the determination of the total MMP-9 concentration in plasma. The used 3D bioplatform (G2 poly(amidoamine) dendrimer; PAMAM) improved the sensitivity of the determination by significantly expanding the surface area of the receptor layer. The antigen-antibody recognition process was controlled by quartz crystal microbalance with dissipation (QCM-D) and electrochemical impedance spectroscopy (EIS). The effect of the orientation of antibody molecules in the sensing layer on the work parameters of the immunosensor was analyzed using unmodified PAMAM (PAMAM-NH₂) and PAMAM functionalized with -COOH groups (PAMAM-COOH). The developed immunosensor based on PAMAM-NH₂ was characterized by a lower detection limit (LOD = 2.0 pg⋅mL^-1) and wider analytical range (1·10^-4 - 5 μg⋅mL^-1 for EIS and QCM-D) compared to PAMAM-COOH immunosensor (EIS: 1·10^-4 - 0.5 μg⋅mL^-1; QCM-D: 5·10^-4 - 0.5 μg⋅mL^-1). The functionality of the proposed device was verified in spiked plasma. The recoveries determined in commercial human and rat plasma and noncommercial rat plasma were very close to the value of 100% and in the range of 96-120% for Au/PAMAM-NH₂/Ab and Au/PAMAM-COOH/Ab immunosensors, respectively. The designed analytical devices showed high selectivity and sensitivity without the use of any amplifiers such as metal nanoparticles or enzymes.

Effect of Protein–Protein Interactions on Translational Diffusion of Spheroidal Proteins

One of the commonly accepted approaches to estimate protein–protein interactions (PPI) in aqueous solutions is the analysis of their translational diffusion. The present review article observes a phenomenological approach to analyze PPI effects via concentration dependencies of self- and collective translational diffusion coefficient for several spheroidal proteins derived from the pulsed field gradient NMR (PFG NMR) and dynamic light scattering (DLS), respectively. These proteins are rigid globular α-chymotrypsin (ChTr) and human serum albumin (HSA), and partly disordered α-casein (α-CN) and β-lactoglobulin (β-Lg). The PPI analysis enabled us to reveal the dominance of intermolecular repulsion at low ionic strength of solution (0.003–0.01 M) for all studied proteins. The increase in the ionic strength to 0.1–1.0 M leads to the screening of protein charges, resulting in the decrease of the protein electrostatic potential. The increase of the van der Waals potential for ChTr and α-CN characterizes their propensity towards unstable weak attractive interactions. The decrease of van der Waals interactions for β-Lg is probably associated with the formation of stable oligomers by this protein. The PPI, estimated with the help of interaction potential and idealized spherical molecular geometry, are in good agreement with experimental data.

Conformational buffering underlies functional selection in intrinsically disordered protein regions

Many disordered proteins conserve essential functions in the face of extensive sequence variation, making it challenging to identify the mechanisms responsible for functional selection. Here we identify the molecular mechanism of functional selection for the disordered adenovirus early gene 1A (E1A) protein. E1A competes with host factors to bind the retinoblastoma (Rb) protein, subverting cell cycle regulation. We show that two binding motifs tethered by a hypervariable disordered linker drive picomolar affinity Rb binding and host factor displacement. Compensatory changes in amino acid sequence composition and sequence length lead to conservation of optimal tethering across a large family of E1A linkers. We refer to this compensatory mechanism as conformational buffering. We also detect coevolution of the motifs and linker, which can preserve or eliminate the tethering mechanism. Conformational buffering and motif-linker coevolution explain robust functional encoding within hypervariable disordered linkers and could underlie functional selection of many disordered protein regions.

Estimating near-wall diffusion coefficients of arbitrarily shaped rigid macromolecules

We developed a computationally efficient approach to approximate near-wall diffusion coefficients of arbitrarily shaped rigid macromolecules. The proposed method relies on extremum principles for Stokes flows produced by the motion of rigid bodies. In the presence of the wall, the rate of energy dissipation is decreased relative to the unbounded fluid. In our approach, the position- and orientation-dependent mobility matrix of a body suspended near a no-slip plane is calculated numerically using a coarse-grained molecular model and the Rotne-Prager-Yamakawa description of hydrodynamics. Effects of the boundary are accounted for via Blake's image construction. The matrix components are scaled using ratios of the corresponding bulk values evaluated for the detailed representation of the molecule and its coarse-grained model, leading to accurate values of the near-wall diffusion coefficients. We assess the performance of the approach for two biomolecules at different levels of coarse-graining.

Variation of Structural and Dynamical Flexibility of Myelin Basic Protein in Response to Guanidinium Chloride

Myelin basic protein (MBP) is intrinsically disordered in solution and is considered as a conformationally flexible biomacromolecule. Here, we present a study on perturbation of MBP structure and dynamics by the denaturant guanidinium chloride (GndCl) using small-angle scattering and neutron spin–echo spectroscopy (NSE). A concentration of 0.2 M GndCl causes charge screening in MBP resulting in a compact, but still disordered protein conformation, while GndCl concentrations above 1 M lead to structural expansion and swelling of MBP. NSE data of MBP were analyzed using the Zimm model with internal friction (ZIF) and normal mode (NM) analysis. A significant contribution of internal friction was found in compact states of MBP that approaches a non-vanishing internal friction relaxation time of approximately 40 ns at high GndCl concentrations. NM analysis demonstrates that the relaxation rates of internal modes of MBP remain unaffected by GndCl, while structural expansion due to GndCl results in increased amplitudes of internal motions. Within the model of the Brownian oscillator our observations can be rationalized by a loss of friction within the protein due to structural expansion. Our study highlights the intimate coupling of structural and dynamical plasticity of MBP, and its fundamental difference to the behavior of ideal polymers in solution.

Biophysical characterisation of human LincRNA-p21 sense and antisense Alu inverted repeats

Human Long Intergenic Noncoding RNA-p21 (LincRNA-p21) is a regulatory noncoding RNA that plays an important role in promoting apoptosis. LincRNA-p21 is also critical in down-regulating many p53 target genes through its interaction with a p53 repressive complex. The interaction between LincRNA-p21 and the repressive complex is likely dependent on the RNA tertiary structure. Previous studies have determined the two-dimensional secondary structures of the sense and antisense human LincRNA-p21 AluSx1 IRs using SHAPE. However, there were no insights into its three-dimensional structure. Therefore, we in vitro transcribed the sense and antisense regions of LincRNA-p21 AluSx1 Inverted Repeats (IRs) and performed analytical ultracentrifugation, size exclusion chromatography, light scattering, and small angle X-ray scattering (SAXS) studies. Based on these studies, we determined low-resolution, three-dimensional structures of sense and antisense LincRNA-p21. By adapting previously known two-dimensional information, we calculated their sense and antisense high-resolution models and determined that they agree with the low-resolution structures determined using SAXS. Thus, our integrated approach provides insights into the structure of LincRNA-p21 Alu IRs. Our study also offers a viable pipeline for combining the secondary structure information with biophysical and computational studies to obtain high-resolution atomistic models for long noncoding RNAs.

Curvature dependence of BAR protein membrane association and dissociation kinetics

BAR (Bin/Amphiphysin/Rvs) domain containing proteins function as lipid bilayer benders and curvature sensors, and they contribute to membrane shaping involved in cell signaling and metabolism. The mechanism for their membrane shape sensing has been investigated by both equilibrium binding and kinetic studies. In prior research, stopped-flow spectroscopy has been used to deduce a positive dependence on membrane curvature for the binding rate constant, k_on, of a BAR protein called endophilin. However, the impact of bulk diffusion of endophilin, on the kinetic binding parameters has not been thoroughly considered. Employing similar methods, and using lipid vesicles of multiple sizes, we obtained a linear dependence of k_on on vesicle curvature. However, we found that the observed relation can be explained without considering the local curvature sensing ability of endophilin in the membrane association process. In contrast, the diffusion-independent unbinding rate constant (k_off) obtained from stopped-flow measurements shows a negative dependence on membrane curvature, which is controlled/mediated by endophilin-membrane interactions. This latter dependency, in addition to protein-protein interactions on the membrane, explains the selective binding of BAR proteins to highly curved membranes in equilibrium binding experiments.

Complementary approaches to obtaining thermodynamic parameters from protein ligand systems-challenges and opportunities

Protein ligand interactions play an important role in biology. Increasingly the aim is to understand and influence protein ligand binding. The binding process is heavily influenced by its thermodynamic parameters. In order to understand how the whole system thermodynamics work it is important to characterise the individual contribution of each of the systems components. While the change in conformational entropy of the protein can be determined using QENS complementary methods are necessary in order to characterise all components. This paper will describe the challenges that can occur when combining the different methods, as well as how they can be overcome.

Structural basis of DegP protease temperature-dependent activation

Protein quality control is an essential cellular function mainly executed by a vast array of different proteases and molecular chaperones. One of the bacterial high temperature requirement A (HtrA) protein family members, the homo-oligomeric DegP protease, plays a crucial role in the Escherichia coli protein quality control machinery by removing unfolded proteins or preventing their aggregation and chaperoning them to their final folded state within the periplasm. DegP contains two regulatory PDZ domains, which play key roles in substrate recognition and in the transformation of DegP between inactive hexameric and proteolytic active cage-like structures. Here, we analyze the interaction and dynamics of the DegP PDZ domains underlying this transformation by high-resolution NMR spectroscopy complemented with biochemical cleavage assays. We identify an interdomain molecular lock, which controls the interactions between the two PDZ domains, regulated by fine-tuned temperature-dependent protein dynamics, and which is potentially conserved in proteins harboring tandem PDZ domains.

Design of Swollen Lipidic Cubic Phase to Increase Transcutaneous Penetration of Biomacromolecules

Lipidic cubic phase (LCP) is a self-assembled system composed of lipids with interpenetrated aqueous channels, and its potential in drug delivery systems has been investigated. Although LCP was shown to improve transcutaneous penetration of hydrophilic molecules of up to 1203 Da so far, the transcutaneous delivery of larger molecules such as proteins has not been achieved. This is likely because proteins are usually larger than the aqueous channels of LCP (∼37.2 Å in diameter), which limits the molecular diffusion in LCP. In this report, we overcome this issue by adding N-octyl-β-d-glucopyranoside to glyceryl monooleate-water-based LCP to give swollen LCP (SLCP), which has larger aqueous channel diameters (∼65.6 Å). First, we systemically evaluated the effect of swelling on drug diffusion in LCP/SLCP. The release kinetics of various peptides and proteins whose sizes ranged from 9.14 to 55.28 Å in diameter were evaluated, and the diffusion coefficients (D) were calculated by the Fickian diffusion model. As expected, all peptides and proteins diffused faster in SLCP than in LCP. A more detailed analysis revealed a negative linear relationship between log D and the ratio of the radius of gyration of the proteins to the aqueous channel radius, indicating that swelling of a cubic nanostructure is an effective strategy to enhance D. Next, the skin penetration of proteins encapsulated in LCP and SLCP was evaluated. The skin penetration of ovalbumin (42.9 kDa), for example, was enhanced by SLCP but not by LCP, and a positive correlation between D and the amount of skin penetration was found. Collectively, this study provides an effective measure for designing LCP systems that enhance transcutaneous penetration of biomacromolecules.

Microfluidics: A Novel Approach for Dehydration Protein Droplets

The equation of state of colloids plays an important role in the modelling and comprehension of industrial processes, defining the working conditions of processes such as drying, filtration, and mixing. The determination of the equation is based on the solvent equilibration, by dialysis, between the colloidal suspension and a reservoir with a known osmotic pressure. In this paper, we propose a novel microfluidic approach to determine the equation of state of a lysozyme solution. Monodispersed droplets of lysozyme were generated in the bulk of a continuous 1-decanol phase using a flow-focusing microfluidic geometry. In this multiphasic system and in the working operation conditions, the droplets can be considered to act as a permeable membrane system. A water mass transfer flow occurs by molecule continuous diffusion in the surrounding 1-decanol phase until a thermodynamic equilibrium is reached in a few seconds to minutes, in contrast with the standard osmotic pressure measurements. By changing the water saturation of the continuous phase, the equation of state of lysozyme in solution was determined through the relation of the osmotic pressure between protein molecules and the volume fraction of protein inside the droplets. The obtained equation shows good agreement with other standard approaches reported in the literature.

Multiscale Approach for Computing Gated Ligand Binding from Molecular Dynamics and Brownian Dynamics Simulations

We develop an approach to characterize the effects of gating by a multiconformation protein consisting of macrostate conformations that are either accessible or inaccessible to ligand binding. We first construct a Markov state model of the apo-protein from atomistic molecular dynamics simulations from which we identify macrostates and their conformations, compute their relative macrostate populations and interchange kinetics, and structurally characterize them in terms of ligand accessibility. We insert the calculated first-order rate constants for conformational transitions into a multistate gating theory from which we derive a gating factor γ that quantifies the degree of conformational gating. Applied to HIV-1 protease, our approach yields a kinetic network of three accessible (semi-open, open, and wide-open) and two inaccessible (closed and a newly identified, "parted") macrostate conformations. The parted conformation sterically partitions the active site, suggesting a possible role in product release. We find that the binding kinetics of drugs and drug-like inhibitors to HIV-1 protease falls in the slow gating regime. However, because γ = 0.75, conformational gating only modestly slows ligand binding. Brownian dynamics simulations of the diffusional association of eight inhibitors to the protease─having a wide range of experimental association constants (∼10⁴-10¹⁰ M^-1 s^-1)─yields gated rate constants in the range of ∼0.5-5.7 × 10⁸ M^-1 s^-1. This indicates that, whereas the association rate of some inhibitors could be described by the model, for many inhibitors either subsequent conformational transitions or alternate binding mechanisms may be rate-limiting. For systems known to be modulated by conformational gating, the approach could be scaled computationally efficiently to screen association kinetics for a large number of ligands.

Transcription-dependent confined diffusion of enzymes within subcellular spaces of the bacterial cytoplasm

BACKGROUND

Knowledge on the localization and mobility of enzymes inside bacterial cells is scarce, but important for understanding spatial regulation of metabolism. The four central enzymes (Rib enzymes) of the riboflavin (RF) biosynthesis pathway in the Gram positive model bacterium Bacillus subtilis have been studied extensively in vitro, especially the heavy RF synthase, a large protein complex with a capsid structure formed by RibH and an encapsulated RibE homotrimer, which mediates substrate-channeling. However, little is known about the behavior and mobility of these enzymes in vivo.

RESULTS

We have investigated the localization and diffusion of the Rib enzymes in the cytoplasm of B. subtilis. By characterizing the diffusion of Rib enzymes in live cells using single particle tracking (SPT) we provide evidence for confined diffusion at the cell poles and otherwise Brownian motion. A majority of RibH particles showed clear nucleoid occlusion and a high degree of confined motion, which is largely abolished after treatment with Rifampicin, revealing that confinement is dependent on active transcription. Contrarily, RibE is mostly diffusive within the cell, showing only 14% encapsulation by RibH nanocompartments. By localizing different diffusive populations within single cells, we find that fast diffusion occurs mostly across the nucleoids located in the cell centers, while the slower, confined subdiffusion occurs at the crowded cell poles.

CONCLUSIONS

Our results provide evidence for locally different motion of active enzymes within the bacterial cytoplasm, setting up metabolic compartmentalization mostly at the poles of cells.

A molecular dynamics study of protein denaturation induced by sulfonate-based surfactants

Microsecond timescale explicit-solvent atomistic simulations were carried out to investigate how anionic surfactants modulate protein structure and dynamics. We found that lysozyme undergoes near-complete denaturation at the high concentration (> 0.1 M) of sodium pentadecyl sulfonate (SPDS), while only partial denaturation occurs at the concentration slightly below 0.1 M. In large part, protein denaturation is structurally manifested by disappearance of helical segments and loss of tertiary interactions. The computational prediction of the extent of burial of cysteine residues was experimentally validated by measuring the accessibility of the respective sulfhydryl groups. Overall, our work indicates an interesting synergy between electrostatic and hydrophobic contributions to lysozyme's denaturation process by anionic surfactants. In fact, first disulfide bridges and hydrogen bonds from protein surface to SPDS head groups loosen the protein globule followed by fuller denaturation via insertion of the surfactant's hydrophobic tails into the protein core.

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.

NMR structure of the Vibrio vulnificus ribosomal protein S1 domains D3 and D4 provides insights into molecular recognition of single-stranded RNAs

The ribosomal S1 protein (rS1) is indispensable for translation initiation in Gram-negative bacteria. rS1 is a multidomain protein that acts as an RNA chaperone and ensures that mRNAs can bind the ribosome in a single-stranded conformation, which could be related to fast recognition. Although many ribosome structures were solved in recent years, a high-resolution structure of a two-domain mRNA-binding competent rS1 construct is not yet available. Here, we present the NMR solution structure of the minimal mRNA-binding fragment of Vibrio Vulnificus rS1 containing the domains D3 and D4. Both domains are homologues and adapt an oligonucleotide-binding fold (OB fold) motif. NMR titration experiments reveal that recognition of miscellaneous mRNAs occurs via a continuous interaction surface to one side of these structurally linked domains. Using a novel paramagnetic relaxation enhancement (PRE) approach and exploring different spin-labeling positions within RNA, we were able to track the location and determine the orientation of the RNA in the rS1-D34 bound form. Our investigations show that paramagnetically labeled RNAs, spiked into unmodified RNA, can be used as a molecular ruler to provide structural information on protein-RNA complexes. The dynamic interaction occurs on a defined binding groove spanning both domains with identical β2-β3-β5 interfaces. Evidently, the 3'-ends of the cis-acting RNAs are positioned in the direction of the N-terminus of the rS1 protein, thus towards the 30S binding site and adopt a conformation required for translation initiation.

Thermophoretic Random Walks and Enhancement of Diffusion

The contribution of the stochastic thermodiffusion to the diffusion enhancement is studied. The thermodiffusion of particles suspended in a liquid may hold place when the spontaneous endo- or exothermal nanoscale events similar to elementary acts of enzymatic reactions occur as the random series in the space and time. In these events, the energy can be emitted or absorbed at nanoscale during few to hundreds of picoseconds. It may cause local spontaneous temperature spikes spreading quickly in the space and decaying with time. The random local temperature spikes create local transient temperature gradients, where thermodiffusion of the molecules and particles holds place as well as the change in the physical properties of the suspending medium due to heating. These thermodiffusion random walks may appear as the enhanced usual Stokes-Einstein diffusion when the energy absorption/generation is high enough. The evaluated relative contribution of the mentioned effect to the molecular mobility is shown to be in agreement with experimental data on enzymatic reactions from the literature.

Small Angle X-ray Scattering, Molecular Modeling, and Chemometric Studies from a Thrombin-Like (Lmr-47) Enzyme of Lachesis m. rhombeata Venom

INTRODUCTION

Snakebite envenomation is considered a neglected tropical disease, and SVTLEs critical elements are involved in serious coagulopathies that occur on envenoming. Although some enzymes of this group have been structurally investigated, it is essential to characterize other proteins to better understand their unique properties such as the Lachesis muta rhombeata 47 kDa (Lmr-47) venom serine protease.

METHODS

The structure of Lmr-47 was studied in solution, using SAXS, DLS, CD, and in silico by homology modeling. Molecular docking experiments simulated 21 competitive inhibitors.

RESULTS

At pH 8.0, Lmr-47 has an Rg of 34.5 ± 0.6 Å, Dmax of 130 Å, and SR of 50 Å, according to DLS data. Kratky plot analysis indicates a rigid shape at pH 8.0. Conversely, the pH variation does not change the center of mass's intrinsic fluorescence, possibly indicating the absence of fluorescent amino acids in the regions affected by pH variation. CD experiments show a substantially random coiled secondary structure not affected by pH. The low-resolution model of Lmr-47 presented a prolate elongated shape at pH 8.0. Using the 3D structure obtained by molecular modeling, docking experiments identified five good and three suitable competitive inhibitors.

CONCLUSION

Together, our work provided insights into the structure of the Lmr-47 and identified inhibitors that may enhance our understanding of thrombin-like family proteins.