The influence of correlated protein-water volume fluctuations on the apparent compressibility of proteins determined by ultrasonic velocimetry

Interaction of the new inhibitor paxlovid (PF-07321332) and ivermectin with the monomer of the main protease SARS-CoV-2: A volumetric study based on molecular dynamics, elastic networks, classical thermodynamics and SPT

The COVID-19 pandemic has accelerated the study of drugs, most notably ivermectin and more recently Paxlovid (PF-07321332) which is in phase III clinical trials with experimental data showing covalent binding to the viral protease Mpro. Theoretical developments of catalytic site-directed docking support thermodynamically feasible non-covalent binding to Mpro. Here we show that Paxlovid binds non-covalently at regions other than the catalytic sites with energies stronger than reported and at the same binding site as the ivermectin B1a homologue, all through theoretical methodologies, including blind docking. We volumetrically characterize the non-covalent interaction of the ivermectin homologues (avermectins B1a and B1b) and Paxlovid with the mMpro monomer, through molecular dynamics and scaled particle theory (SPT). Using the fluctuation-dissipation theorem (FDT), we estimated the electric dipole moment fluctuations at the surface of each of complex involved in this study, with similar trends to that observed in the interaction volume. Using fluctuations of the intrinsic volume and the number of flexible fragments of proteins using anisotropic and Gaussian elastic networks (ANM+GNM) suggests the complexes with ivermectin are more dynamic and flexible than the unbound monomer. In contrast, the binding of Paxlovid to mMpro shows that the mMpro-PF complex is the least structurally dynamic of all the species measured in this investigation. The results support a differential molecular mechanism of the ivermectin and PF homologues in the mMpro monomer. Finally, the results showed that Paxlovid despite beingbound in different sites through covalent or non-covalent forms behaves similarly in terms of its structural flexibility and volumetric behaviour.

Characterization of the ergometric properties of commercial bioactive dairy peptides

The thermodynamic properties of bioactive peptides provide insights into their functional behavior and their biological efficacy. We conducted precise analyses of the density, the ultrasonic velocity and the relative attenuation of serial dilutions of three commercial dairy peptides prepared by enzymatic methods. From these we determined the partial specific volume and the partial specific adiabatic compressibility coefficient for the peptides. At concentrations greater than ~2.5 ​mg ​mL⁻¹, the apparent values for specific volume and adiabatic compressibility were constant, differing between the three peptides at ±3% for specific volume and ±70% for compressibility. Both specific volume and adiabatic compressibility were highly dependent on concentration, indicating the importance of precise low concentration measurements to obtain correct values for these thermodynamic parameters. From these parameters it was apparent that restructuring of water molecules around the peptides (and their associated counterions) led to compact solutes that were also incompressible. These thermodynamic analyses are critical for understanding how the properties and the beneficial effects of bioactive peptides are influenced by their chemical environment.

•

Dissolved dairy peptide properties distinguishable from ergometric analyses.

•

Specific volume and adiabatic compressibility evaluate bioactive peptide hydration.

•

Commercial bioactive dairy peptides are compact and incompressible.

•

Compactness and incompressibility of peptide affected by hydrogen-bonding amino acids.

•

Solution concentration affects values of measured thermodynamic parameters.

Pressure-Dependent Conformation and Fluctuation in Folded Protein Molecules

Hydrostatic pressure leads to nonuniform compression of proteins. The structural change is on average only about 0.1 Å kbar(-1), and is therefore within the range of fluctuations at ambient pressure. The largest changes are around cavities and buried water molecules. Sheets distort much more than helices. Hydrogen bonds compress about 0.012 Å kbar(-1), although there is a wide range, including some hydrogen bonds that lengthen. In the presence of ligands and inhibitors, structural changes are smaller. Pressure has little effect on rapid fluctuations, but with larger scale slower motions, pressure increases the population of excited states (if they have smaller overall volume), and slows the fluctuations. In barnase, pressure is shown to be a useful way to characterise fluctuations on the timescale of microseconds, and helps to show that fluctuations in barnase are hierarchical, with the faster fluctuations providing a platform for the slower ones. The excited states populated at high pressure are probably functionally important.

Pressure-dependent structure changes in barnase on ligand binding reveal intermediate rate fluctuations

In this work we measured 1H NMR chemical shifts for the ribonuclease barnase at pressures from 3 MPa to 200 MPa, both free and bound to d(CGAC). Shift changes with pressure were used as restraints to determine the change in structure with pressure. Free barnase is compressed by approximately 0.7%. The largest changes are on the ligand-binding face close to Lys-27, which is the recognition site for the cleaved phosphate bond. This part of the protein also contains the buried water molecules. In the presence of d(CGAC), the compressibility is reduced by approximately 70% and the region of structural change is altered: the ligand-binding face is now almost incompressible, whereas changes occur at the opposite face. Because compressibility is proportional to mean square volume fluctuation, we conclude that in free barnase, volume fluctuation is largest close to the active site, but when the inhibitor is bound, the fluctuations become much smaller and are located mainly on the opposite face. The timescale of the fluctuations is nanoseconds to microseconds, consistent with the degree of ordering required for the fluctuations, which are intermediate between rapid uncorrelated side-chain dynamics and slow conformational transitions. The high-pressure technique is therefore useful for characterizing motions on this relatively inaccessible timescale.