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.

Influence of pH and Salts on Partial Molar Volume of Lysozyme and Bovine Serum Albumin in Aqueous Solutions

The partial molar volume of lysozyme and bovine serum albumin in aqueous solutions at different pH values and in aqueous solutions containing sodium chloride, ammonium chloride, sodium sulfate, or ammonium sulfate at different concentrations at pH 7.0 was investigated experimentally at 298.15 K and 1 bar. It was found that the influence of the pH value and the salts on the partial molar volume of the proteins is small, but trends were measurable. Furthermore, the partial molar volume of lysozyme in pure water at different pH values and in aqueous solutions with different sodium chloride concentrations at pH 7.0 was predicted by molecular simulations. The predictions are in good agreement with the experimental data.

Molecular Explanation for the Abnormal Flux of Material into a Hot Spot in Ester Monolayers

Langmuir monolayers of certain surfactants show a negative derivative of the surface pressure with respect to temperature. In these monolayers, a local temperature gradient leads to local yielding of the solid phase to a kinetically flowing liquid, so that the material flows toward the hotter regions that act as sinks. The accumulation of material leads to the formation of nonequilibrium multilamellar bubbles of different sizes. Here we investigate the molecular factors leading to such a peculiar behavior. First, we identify the required structural molecular moieties, and second we vary the composition of the subphase in order to analyze its influence. We conclude that esters appear to be unique in two key aspects: they form monolayers whose compression isotherms shift to lower areas as the temperature increases, and thus collapse into a hot spot; and they bind weakly to the aqueous subphase, i.e., water does not attach to the monolayer at the molecular level, but only supports it. Molecular simulations for a selected system confirm and help explain the observed behavior: surfactant molecules form a weak hydrogen bonding network, which is disrupted upon heating, and also the molecular tilting changes with temperature, leading to changes in the film density.

MD simulations and multivariate studies for modeling the antileishmanial activity of peptides

Leishmaniasis, a protozoan-caused disease, requires alternative treatments with minimized side-effects and less prone to resistance development. Antimicrobial peptides represent a possible choice to be developed. We report on the prospection of structural parameters of 23 helical antimicrobial and leishmanicidal peptides as a tool for modeling and predicting the activity of new peptides. This investigation is based on molecular dynamic simulations (MD) in mimetic membrane environment, as most of these peptides share the feature of interacting with phospholipid bilayers. To overcome the lack of experimental data on peptides' structures, we started simulations from designed 100% α-helices. This procedure was validated through comparisons with NMR data and the determination of the structure of Decoralin-amide. From physicochemical features and MD results, descriptors were raised and statistically related to the minimum inhibitory concentration against Leishmania by the multivariate data analysis technique. This statistical procedure confirmed five descriptors combined by different loadings in five principal components. The leishmanicidal activity depends on peptides' charge, backbone solvation, volume, and solvent-accessible surface area. The generated model possesses good predictability (q²  = 0.715, r²  = 0.898) and is indicative for the most and the least active peptides. This is a novel theoretical path for structure-activity studies combining computational methods that identify and prioritize the promising peptide candidates.

In silico characterization of protein partial molecular volumes and hydration shells

In this paper we present a computational approach, based on NVT molecular dynamics trajectories, that allows the direct evaluation of the protein partial molecular volume. The results obtained for five different globular proteins demonstrate the accuracy of this computational procedure in reproducing protein partial molecular volumes, providing quantitative characterization of the hydration shell in terms of the protein excluded volume, hydration shell ellipsoidal volume and related solvent density. Remarkably, our data indicate for the hydration shell a ≈10% solvent density increase with respect to the liquid water bulk density, in excellent agreement with the available experimental data.

Exploring volume, compressibility and hydration changes of folded proteins upon compression

Understanding the physical basis of the structure, stability and function of proteins in solution, including extreme environmental conditions, requires knowledge of their temperature and pressure dependent volumetric properties. One physical-chemical property of proteins that is still little understood is their partial molar volume and its dependence on temperature and pressure. We used molecular dynamics simulations of aqueous solutions of a typical monomeric folded protein, staphylococcal nuclease (SNase), to study and analyze the pressure dependence of the apparent volume, Vapp, and its components by the Voronoi-Delaunay method. We show that the strong decrease of Vapp with pressure (βT = 0.95 × 10(-5) bar(-1), in very good agreement with the experimental value) is essentially due to the compression of the molecular volume, VM, ultimately, of its internal voids, V. Changes of the intrinsic volume (defined as the Voronoi volume of the molecule), the contribution of the solvent to the apparent volume, and of the contribution of the boundary voids between the protein and the solvent have also been studied and quantified in detail. The pressure dependences of the volumetric characteristics obtained are compared with the temperature dependent behavior of these quantities and with corresponding results for a natively unfolded polypeptide.

Infinitely dilute partial molar properties of proteins from computer simulation

A detailed understanding of temperature and pressure effects on an infinitely dilute protein's conformational equilibrium requires knowledge of the corresponding infinitely dilute partial molar properties. Established molecular dynamics methodologies generally have not provided a way to calculate these properties without either a loss of thermodynamic rigor, the introduction of nonunique parameters, or a loss of information about which solute conformations specifically contributed to the output values. Here we implement a simple method that is thermodynamically rigorous and possesses none of the above disadvantages, and we report on the method's feasibility and computational demands. We calculate infinitely dilute partial molar properties for two proteins and attempt to distinguish the thermodynamic differences between a native and a denatured conformation of a designed miniprotein. We conclude that simple ensemble average properties can be calculated with very reasonable amounts of computational power. In contrast, properties corresponding to fluctuating quantities are computationally demanding to calculate precisely, although they can be obtained more easily by following the temperature and/or pressure dependence of the corresponding ensemble averages.

Thermal stability of the hydrogen-bonded water network in the hydration shell of islet amyloid polypeptide

The effect of temperature on the connectivity of hydrogen bonds in the hydration shells of the islet amyloid polypeptides (IAPPs) is studied by means of computer simulations. The hydrogen-bonded network of hydration water homogeneously envelopes a peptide at low temperature and breaks into an ensemble of small clusters upon heating. This thermal break occurs via a percolation transition, which is not found to be sensitive to the chemical modifications of IAPP (IAPP with and without a disulfide bridge, human and rat IAPP). The radius of gyration of IAPP starts to increase when the hydration water network breaks upon heating. The fluctuations of the number of intra-peptide hydrogen bonds show negative correlation with the fraction of molecules in the largest cluster of hydration water. The thermal stability of the network of hydration water is enhanced upon increasing number of intra-peptide hydrogen bonds, which makes the peptide surface more hydrophobic. The thermal stabilities of the hydrogen-bonded water networks in the hydration shells of IAPPs and of several other biomolecules are found to be rather similar: the network breaks between 300 and 330 K, i.e., in the temperature interval where the biological activity of living organisms is maximal.

Comparing the structural properties of human and rat islet amyloid polypeptide by MD computer simulations

Conformational properties of the full-length human and rat islet amyloid polypeptide 1-37 (amyloidogenic hIAPP and non-amyloidogenic rIAPP, respectively) were studied at 310 and 330 K by MD simulations both for the cysteine (reduced IAPP) and cystine (oxidized IAPP) moieties. At all temperatures studied, IAPP does not adopt a well-defined conformation and is essentially random coil in solution, although transient helices appear forming along the peptide between residues 8 and 22, particularly in the reduced form. Above the water percolation transition (at 320 K), the reduced hIAPP moiety presents a considerably diminished helical content remaining unstructured, while the natural cystine moiety reaches a rather compact state, presenting a radius of gyration that is almost 10% smaller and characterized by intrapeptide H-bonds that form many β-bridges in the C-terminal region. This compact conformation presents a short end-to-end distance and seems to form through the formation of β-sheet conformations in the C-terminal region with a minimization of the Y/F distances in a two-step mechanism: the first step taking place when the Y37/F23 distance is ~1.1 nm, and subsequently Y37/F15 reaches its minimum of ~0.86 nm. rIAPP, which does not aggregate, also presents transient helical conformations. A particularly stable helix is located in proximity of the C-terminal region, starting from residues L27 and P28. Our MD simulations show that P28 in rIAPP influences the secondary structure of IAPP by stabilizing the peptide in helical conformations. When this helix is not present, the peptide presents bends or H-bonded turns at P28 that seem to inhibit the formation of the β-bridges seen in hIAPP. Conversely, hIAPP is highly disordered in the C-terminal region, presenting transient isolated β-strand conformations, particularly at higher temperatures and when the natural disulfide bond is present. Such conformational differences found in our simulations could be responsible for the different aggregational propensities of the two different homologues. In fact, the fragment 30-37, which is identical in both homologues, is known to aggregate in vitro, hence the overall sequence must be responsible for the amyloidogenicity of hIAPP. The increased helicity in rIAPP induced by the serine-to-proline variation at residue 28 seems to be a plausible inhibitor of its aggregation.