Umbrella Sampling Method Research Articles

The Solvation of Na+ Ions by Ethoxylate Moieties Enhances Adsorption of Sulfonate Surfactants at the Air-Water Interface

HypothesisExperiments show pronounced synergy in the reduction of surface tension when the nonionic surfactant Poly(oxy-1,2-ethanediyl), .alpha.-tris(1-phenylethyl)phenyl-.omega.-hydroxy- (Ethoxylated tristyrylphenol, EOT) is mixed with the anionic surfactant Sodium 4-dodecylbenzenesulfonate (NaDDBS). We hypothesize that the synergism is due to counterion (cation) effects. This would be unusual as one of the surfactants is nonionic. To test this hypothesis, the molecular mechanisms responsible need to be probed using experiments and simulations. ApproachThe interfacial properties of mixtures comprising EOT and NaDDBS are investigated using equilibrium molecular dynamics (MD) simulations. Free energy calculation using thermodynamic integration and umbrella sampling methods are employed to analyze the molecular interactions at surface and reveal the role of counterion solvation on the results observed. Simulation snapshots and trajectories are interrogated to confirm the findings. FindingsSimulation results indicate that the ethoxylate moieties solvate Na+ ions, forming long-lasting cation-EOT complexes. Free energy calculations suggest that these complexes are more stable at the interface than in the bulk, likely because of changes in the dielectric properties of water. The cation-EOT complexes, in turn, cause a stronger affinity between the interface and NaDDBS when EOT is present. Similar studies conducted for mixtures of EOT and cationic surfactant Dodecylammonium chloride (DAC) do not show evidence of Cl− ions solvation via the ethoxylate moieties, while the DAC headgroup was found to form hydrogen bonds with the EOT headgroup. This suggests that the mechanisms observed are somewhat ion specific.

Probing the Molecular Basis of Aminoacyl-Adenylate Affinity With Mycobacterium tuberculosis Leucyl-tRNA Synthetase Employing Molecular Dynamics, Umbrella Sampling Simulations and Site-Directed Mutagenesis.

Leucyl-tRNA synthetase (LeuRS) is clinically validated molecular target for antibiotic development. Recently, we have reported several classes of small-molecular inhibitors targeting aminoacyl-adenylate binding site of Mycobacterium tuberculosis LeuRS with antibacterial activity. In this work, we performed in silico site-directed mutagenesis of M. tuberculosis LeuRS synthetic site in order to identify the most critical amino acid residues for the interaction with substrate and prove binding modes of inhibitors. We carried out 20-ns molecular dynamics (MD) simulations and used umbrella sampling (US) method for the calculation of the binding free energy (ΔGb) of leucyl-adenylate with wild-type and mutated forms of LeuRS. According to molecular modeling results, it was found that His89, Tyr93, and Glu660 are essential amino acid residues both for aminoacyl-adenylate affinity and hydrogen bond formation. We have selected His89 for experimental site-directed mutagenesis since according to our previous molecular docking results this amino acid residue was predicted to be important for inhibitor interaction in adenine-binding region. We obtained recombinant mutant M. tuberculosis LeuRS H89A. Using aminoacylation assay we have found that the mutation of His89 to Ala in the active site of M. tuberculosis LeuRS results in significant decrease of inhibitory activity for compounds belonging to three different chemical classes-3-phenyl-5-(1-phenyl-1H-[1,2,3]triazol-4-yl)-[1,2,4]oxadiazoles, N-benzylidene-N'-thiazol-2-yl-hydrazines, and 1-oxo-1H-isothiochromene-3-carboxylic acid (4-phenyl-thiazol-2-yl)-amide derivatives. Therefore, the interaction with His89 should be taken into account during further M. tuberculosis LeuRS inhibitors development and optimization.

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an SN2 Reaction in Water.

We present a tutorial to carry out umbrella-sampling free-energy simulations with a combined quantum mechanical and molecular mechanical (QM/MM) potential, which may also be used in a computational or biophysical chemistry curriculum for first-year graduate and undergraduate students. In this article, we choose the Type II SN2 Menshutkin reaction between ammonia and chloromethane to construct the potential of mean force (PMF) for the reaction in aqueous solution. In this exercise, we wish to accomplish three tasks: (1) an understanding of the concept of PMF and the umbrella-sampling free-energy simulation method, (2) the use of a combined QM/MM potential in molecular dynamics simulation of chemical reactions, and (3) an understanding of solvent effects and intermolecular interactions on chemical reactions through comparison with gas-phase results. Analysis of the simulation results allows students to appreciate the difference between the transition state along a PMF from statistical simulations versus the optimized saddle point in the gas phase, the shift of transition state with respect to solvation, and the significance of electronic polarization. Through this tutorial, students will gain the necessary skills to carry out free energy simulations of chemical reactions in solution and in enzymes both in a classroom and in the laboratory.

Dimerization of Full-Length Aβ-42 Peptide: A Comparison of Different Force Fields and Water Models.

Among the two isoforms of amyloid-β i. e., Aβ-40 and Aβ-42, Aβ-42 is more toxic due to its increased aggregation propensity. The oligomerization pathways of amyloid-β may be investigated by studying its dimerization process at an atomic level. Intrinsically disordered proteins (IDPs) lack well-defined structures and are associated with numerous neurodegenerative disorders. Molecular dynamics simulations of these proteins are often limited by the choice of parameters due to inconsistencies in the empirically developed protein force fields and water models. To evaluate the accuracy of recently developed force fields for IDPs, we study the dimerization of full-length Aβ-42 in aqueous solution with three different combinations of AMBER force field parameters and water models such as ff14SB/TIP3P, ff19SB/OPC, and ff19SB/TIP3P using classical MD and Umbrella Sampling method. This work may be used as a benchmark to compare the performance of different force fields for the simulations of IDPs.

Tunneling Mechanisms of Quinones in Photosynthetic Reaction Center–Light Harvesting 1 Supercomplexes

In photosynthesis, light energy is absorbed and transferred to the reaction center, ultimately leading to the reduction of quinone molecules through the electron transfer chain. The oxidation and reduction of quinones generate an electrochemical potential difference used for adenosine triphosphate synthesis. The trafficking of quinone/quinol molecules between electron transport components has been a long‐standing question. Here, an atomic‐level investigation into the molecular mechanism of quinol dissociation in the photosynthetic reaction center–light‐harvesting complex 1 (RC–LH1) supercomplexes from Rhodopseudomonas palustris, using classical molecular dynamics (MD) simulations combined with self‐random acceleration MD‐MD simulations and umbrella sampling methods, is conducted. Results reveal a significant increase in the mobility of quinone molecules upon reduction within RC–LH1, which is accompanied by conformational modifications in the local protein environment. Quinol molecules have a tendency to escape from RC–LH1 in a tail‐first mode, exhibiting channel selectivity, with distinct preferred dissociation pathways in the closed and open LH1 rings. Furthermore, comparative analysis of free energy profiles indicates that alternations in the protein environment accelerate the dissociation of quinol molecules through the open LH1 ring. In particular, aromatic amino acids form π‐stacking interactions with the quinol headgroup, resembling the key components in a conveyor belt system. This study provides insights into the molecular mechanisms that govern quinone/quinol exchange in bacterial photosynthesis and lays the framework for tuning electron flow and energy conversion to improve metabolic performance.

Entropy-favorable adsorption of polymer-grafted nanoparticles at fluid-fluid interfaces.

The adsorption of polymer-grafted nanoparticles at interfaces is a problem of fundamental interest in physics and soft materials. This adsorption behavior is governed by the interplay between interaction potentials and entropic effects. Here, we use molecular dynamics simulations and umbrella sampling methods to study the adsorption behavior of a Janus-like homopolymer-grafted nanoparticle at fluid-fluid interfaces. By calculating the potential of the mean force as the particle moves from fluid A to the interface, the adsorption energy Ea can be obtained. When two homopolymer chains with types A and B are grafted to the opposite poles of the particle, Ea shows a scaling behavior with respect to chain length N: Ea ∝ N0.598. This is determined by the interactions between polymers and fluids. The enthalpy dominates, and the entropy effects mainly come from the rotational entropy loss of the polymer-grafted nanoparticle at interfaces, which disfavors the stabilization of particles at interfaces. When the grafted polymer number m is large, the adsorption energy exhibits a linear dependence on m. While the enthalpy dominates the behavior, the entropy becomes significant at a larger chain length of N = 15, where the configurational entropy of the polymer chains dominates the entropy of the system. The globule-coil transition occurs when polymers move from poor solvents to good solvents, increasing the configurational entropy and favoring the stabilization of particles at interfaces. Our study provides novel insights into the stabilization mechanism of polymer-grafted nanoparticles at interfaces and reveals the stabilization mechanism favored by the configurational entropy of grafted polymer chains.

Rational Design of Targeted Gold Nanoclusters with High Affinity to Integrin αvβ3 for Combination Cancer Therapy.

The unique attributes of targeted nano-drug delivery systems (TNDDSs) over conventional cancer therapies in suppressing off-target effects make them one of the most promising options for cancer treatment. There is evidence that the density of surface-conjugated ligands is a crucial factor in achieving the desired therapeutic efficacy of TNDDSs, but this is hardly manageable in conventional nanomaterials. In this context, ligand-protected gold nanoclusters (AuNCs) are excellent candidates for developing new TNDDSs with a unique control on their surface functionalities, thus helping to achieve enhanced delivery performance. Here, we study the interactions and binding free energies between ten different functionalized Au144(SR)60 (SR = thiolate ligand) nanoclusters and integrin αvβ3 using molecular dynamics simulations and the umbrella sampling method to obtain the optimal formulations. The AuNCs were functionalized with anticancer drugs (5-fluorouracil or signaling pathways inhibitors, such as capivasertib, linifanib, tanespimycin, and taselisib) and integrin-targeting peptides (RGD4C or QS13), and we identified the optimal mixed ligand layer to enhance their binding affinity to the cancer cell receptor. The results showed that changing the proportions of the same type of ligands on the surface of AuNCs led to differences of up to 38 kcal/mol in computed binding free energies. RGD4C as the targeting peptide resulted in greater affinity for αvβ3, and in most formulations studied, a higher amount of drug than peptide was needed. Polar and charged residues, such as Ser123, Asp150, Tyr178, Arg214, and Asp251 were found to play a significant role in AuNC binding. Our simulations also revealed that Mn2+ cations are crucial for stabilizing the αvβ3-AuNC complex. These findings demonstrate the potential of carefully designing the surface composition of TNDDSs to optimize their target affinity and specificity.

Investigation of the Driving Force of Replacing Adsorbed Hydrocarbons by CO2 in Organic Matter from an Energy Perspective.

Carbon dioxide (CO2) has been widely used to enhance the recovery of adsorbed hydrocarbons from the organic matter (OM) in shale formations. To reveal the driving force of replacing adsorbed hydrocarbons from OM by CO2, we performed molecular dynamics (MD) simulations of the replacement process and calculated the interaction forces between CO2 and hydrocarbons. In addition, based on the umbrella sampling method, steered MD simulations were performed, and the free energy profiles of hydrocarbons were obtained using the weighted histogram analysis method. Results show that the condition of the hydrocarbon replacement by CO2 requires the hydrocarbon to have sufficient kinetic energy or to have a sufficiently large attractive force exerted to ensure that the hydrocarbon escapes the potential well of the OM. The attractive forces exerted on hydrocarbon molecules by CO2 can significantly decrease the energy barrier associated with hydrocarbon movement away from the OM surface. Furthermore, both CO2 and supercritical CO2 can effectively displace adsorbed hydrocarbon gas (methane) on the OM, while supercritical CO2 is required to enhance the recovery of adsorbed hydrocarbon oil (n-dodecane). The results obtained in this study provide guidance for enhancing the recovery of adsorbed hydrocarbons by CO2 in shale formations.

Potential acridinedione derivatives for the development of a heterobifunctional PROTAC for targeted degradation of PCSK9 protein

Proprotein convertase substilisin-like/kexin type 9 (PCSK9) interacts with the low-density lipoprotein (LDL) receptor by protein-protein interaction (PPI). The inhibition of this PPI lowers the LDL cholesterol levels, thereby reducing the risks of cardiovascular diseases. We identified in-house synthesized acridinedione scaffolds as binders of PCSK9 through microsecond timescale conventional and biased molecular dynamics (MD) simulations. The free energy of binding was calculated by the Molecular Mechanics Poisson-Boltzmann Surface Area and umbrella sampling methods for acridinedione molecules and compared with reference molecules (co-crystallized inhibitors). Three acridinedione molecules (DSPD-2, DSPD-4, and DSPD-6) were identified as potential binders of the PCSK9 allosteric binding site from preliminary docking and steered MD simulations. The stability of molecules inside the allosteric pocket was established by the low values of root mean square deviations of backbone C-α atoms. The acridinedione derivatives also showed a shorter distance of hydrogen bonds as compared to the standard molecules. The molecule DSPD-6 showed the lowest binding free energy as compared to other selected and standard molecules. The binding of the selected molecules at the allosteric binding pocket had no effect on structural parameters of LDL receptor-binding site of PCSK9 protein. Finally, a proteolysis targeting chimera (PROTAC) was developed by attaching a proteasome recruiting tag to DSPD-6 (lowest binding affinity) to achieve ligand-induced protein degradation. These results qualify the development of the DSPD-6 molecule as a heterobifunctional PROTAC for targeted degradation of PCSK9 protein. Moreover, our computational strategy provides a framework for the identification of molecules against conventionally undruggable PPI targets.

Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study.

The construction of the stratum corneum (SC) is crucial to the problems of transdermal drug delivery. SC consists of the keratinocyte layers and the lipid matrix surrounding it. Among them, the lipid matrix is the barrier for many exogenous molecules, mainly composed of ceramides (CERs), free fatty acids (FFA), and cholesterol (CHOL). In this work, we developed single-component (CERs, CER-NS, and CER-EOS) and six three-component models, and each model was simulated by using the GROMOS-54A7 force field. Short-period phase (SPP) and long-period phase (LPP) systems were established separately, and area per lipid (APL), thickness, order of carbon chain (SCD), and density distribution were analyzed. The transition of CER-NS and CER-EOS in LPP was observed. The results of hydrogen bonds in the lipid systems indicated that a strong hydrogen-bond network was formed between the skin-lipid bilayers. Umbrella sampling method simulations were performed to calculate the free energy change of ethanol moving into the skin-lipid bilayer. The results revealed that ethanol molecules pulled some water molecules into the membrane when they passed through SPP-1. Our findings provided some insights and models of the stratum corneum that could be used for the subsequent mechanism of macromolecule permeation through membranes in drugs, cosmetics, and so on.

Trapping a non-cognate nucleotide upon initial binding for replication fidelity control in SARS-CoV-2 RNA dependent RNA polymerase.

The RNA dependent RNA polymerase (RdRp) in SARS-CoV-2 is a highly conserved enzyme responsible for viral genome replication/transcription. To understand how the viral RdRp achieves fidelity control during such processes, here we computationally investigate the natural non-cognate vs. cognate nucleotide addition and selectivity during viral RdRp elongation. We focus on the nucleotide substrate initial binding (RdRp active site open) to the prechemical insertion (active site closed) of the RdRp. The current studies were first carried out using microsecond ensemble equilibrium all-atom molecular dynamics (MD) simulations. Due to the slow conformational changes (from open to closed) during nucleotide insertion and selection, enhanced or umbrella sampling methods have been further employed to calculate the free energy profiles of the nucleotide insertion. Our studies find notable stability of noncognate dATP and GTP upon initial binding in the active-site open state. The results indicate that while natural cognate ATP and Remdesivir drug analogue (RDV-TP) are biased toward stabilization in the closed state to facilitate insertion, the natural non-cognate dATP and GTP can be well trapped in off-path initial binding configurations and prevented from insertion so that to be further rejected. The current work thus presents the intrinsic nucleotide selectivity of SARS-CoV-2 RdRp for natural substrate fidelity control, which should be considered in antiviral drug design.

Insights from in silico study of receptor energetics of SARS-CoV-2 variants.

The emergence of new variants of the novel coronavirus SARS-CoV-2 with increased infectivity, superior virulence, high transmissibility, and unmatched immune escape has demonstrated the adaptability and evolutionary fitness of the virus. The subject of relative order of the binding affinity of SARS-CoV-2 variants with the human ACE2 (hACE2) receptor is hotly debated and its resolution has implications for drug design and development. In this work, we have investigated the energetics of the binding of receptor binding domain (RBD) of SARS-CoV-2 variants of concern (VOCs): Beta (B.1.351), Delta (B.1.617.2), Omicron (B.1.1.529), variant of interest (VOI): Kappa (B.1.617.1), and Delta Plus (B.1.617.2.1) variant with the human ACE2 receptor by using the umbrella sampling (US) method. Our work indicates that Delta and Delta Plus variants have greater values of the US binding free energy than Wild-type (WT), whereas Beta, Kappa, and Omicron variants have lower values. Further analysis of hydrogen bonding, salt bridges, non-bonded interaction energy, and contact surface area at the RBD-hACE2 interface establish Delta as the variant with the highest binding affinity among these variants.

Chemistry of the interaction between Imidazole derivatives as corrosion inhibitors molecules and copper/brass/zinc surfaces: A DFT, reactive and classical molecular force fields study

In the present research, we explored the corrosion inhibition mechanisms on metals and alloys by imidazole, methylimidazole, and benzimidazole applying density functional theory (DFT), reactive (ReaxFF), and classical force fields. Our investigation focused on the interaction of these molecules with different types of bare metal surfaces. The results of our studies have been compiled to compare the protection efficiency of the three organic molecules for copper, zinc, and copper-zinc alloys (Cu63%-Zn37%) exposed to the same corrosion conditions in an acid environment. The calculated interaction energy quantities are strongly correlated with the experimental corrosion inhibition efficiencies BIM > MIM > IM. The imidazole derivative studied, significantly better protect copper than zinc, which is attributed to the higher bond strength of Cu–N (formation of chemical bonds) than of Zn–N, confirms the increase in the protection effectiveness against corrosion of Zn < Cu63%-Zn37% < Cu. This present study demonstrates the enormous potential of ReaxFF calculations to accurately simulate the adsorption of inhibitor molecules onto metal surfaces at a fraction of the computational cost of DFT, as well as the performance of the umbrella sampling method to assess the adsorption free energies.

Adsorption of zwitterionic oligomer-grafted silica nanoparticles on rock surface in enhanced oil recovery studied by molecular dynamics simulations

The application of surface-functionalized nanoparticles (NPs) in enhanced oil recovery (EOR) has been rapidly increasing. Adsorption of NPs on the rock surface is the key step. In this study, we systematically investigated the adsorption and desorption mechanisms of various brush-like zwitterionic oligomers functionalized silica NPs on the hydrophilic rock surface in an aqueous solution using molecular dynamics simulations. Our simulations show that the adsorption is primarily influenced by the oligomer chains, while the neutral silica NP core contributes very weakly only when it has direct contact with the rock surface at the grafting density of 5 or below. The interactions from SO3⁻ groups located at the tail of zwitterionic side chains account for more than 72% of the non-bonded interaction energy Enb between the NPs and the rock. The adsorption becomes stronger for NPs with shorter (4 monomer units) oligomer chains but slightly weaker for NPs with longer (9 monomer units) oligomer chains when increasing the grafting density, which is well explained by the change in the number of adsorbed SO3⁻ groups. Desorption simulations with the umbrella sampling method indicate that the free energy of desorption decreases with increasing the grafting density. The silica NPs grafted with longer oligomer chains (9 monomer units) require higher free energies of desorption than the ones with shorter oligomer chains. The adsorption of NPs on the rock surface enhances surface hydrophilicity, and the effect of longer oligomer chains is more profound than that of short oligomer chains. This research provides valuable insights into the adsorption and desorption mechanisms of surface-modified silica NPs on the rock surface at the molecular level and thus gives hints to the design of nanofluids.

Particle Deformability Enables Control of Interactions between Membrane-Anchored Nanoparticles.

Nanoparticles adsorbed on a membrane can induce deformations of the membrane that give rise to effective interactions between the particles. Previous studies have focused primarily on rigid nanoparticles with fixed shapes. However, DNA origami technology has enabled the creation of deformable nanostructures with controllable shapes and mechanical properties, presenting new opportunities to modulate interactions between particles adsorbed on deformable surfaces. Here we use coarse-grained molecular dynamics simulations to investigate deformable, hinge-like nanostructures anchored to lipid membranes via cholesterol anchors. We characterize deformations of the particles and membrane as a function of the hinge stiffness. Flexible particles adopt open configurations to conform to a flat membrane, whereas stiffer particles induce deformations of the membrane. We further show that particles spontaneously aggregate and that cooperative effects lead to changes in their shape when they are close together. Using umbrella sampling methods, we quantify the effective interaction between two particles and show that stiffer hinge-like particles experience stronger and longer-ranged attraction. Our results demonstrate that interactions between deformable, membrane-anchored nanoparticles can be controlled by modifying mechanical properties of the particles, suggesting new ways to modulate the self-assembly of particles on deformable surfaces.

Molecular Dynamics Simulation of Polyacrylamide Adsorption on Calcite.

In poorly consolidated carbonate rock reservoirs, solids production risk, which can lead to increased environmental waste, can be mitigated by injecting formation-strengthening chemicals. Classical atomistic molecular dynamics (MD) simulation is employed to model the interaction of polyacrylamide-based polymer additives with a calcite structure, which is the main component of carbonate formations. Amongst the possible calcite crystal planes employed as surrogates of reservoir rocks, the (1 0 4) plane is shown to be the most suitable surrogate for assessing the interactions with chemicals due to its stability and more realistic representation of carbonate structure. The molecular conformation and binding energies of pure polyacrylamide (PAM), hydrolysed polyacrylamide in neutral form (HPAM), hydrolysed polyacrylamide with 33% charge density (HPAM 33%) and sulfonated polyacrylamide with 33% charge density (SPAM 33%) are assessed to determine the adsorption characteristics onto calcite surfaces. An adsorption-free energy analysis, using an enhanced umbrella sampling method, is applied to evaluate the chemical adsorption performance. The interaction energy analysis shows that the polyacrylamide-based polymers display favourable interactions with the calcite structure. This is attributed to the electrostatic attraction between the amide and carboxyl functional groups with the calcite. Simulations confirm that HPAM33% has a lower free energy than other polymers, presumably due to the presence of the acrylate monomer in ionised form. The superior chemical adsorption performance of HPAM33% agrees with Atomic Force Microscopy experiments reported herein.

Fluctuation-Dominated Ligand Binding in Molten Globule Protein.

A molten globule (MG) state is an intermediate state of protein observed during the unfolding of the native structure. In MG states, milk protein α-lactalbumin (aLA) binds to oleic acid (OLA). This MG-aLA-OLA complex, popularly known as XAMLET, performs cytotoxic activities against cancer cell lines. However, the microscopic understanding of ligand recognition ability in the MG state of the protein has not yet been explored. Motivated by this, we explore the binding of bovine aLA with OLA using all-atom molecular dynamics (MD) simulations. We find the binding mode between MG-aLA and OLA using the conformational thermodynamics method. We also estimate the binding free energy using the umbrella sampling (US) method for both the MG state and the neutral state. We find that the binding free energy obtained from US is comparable with earlier experimental results. We characterize the dihedral fluctuations as the ligand is liberated from the active site of the protein using steered MD. The low energy fluctuations occur near the ligand binding site, which eventually transfer toward the Ca2+-binding site as the ligand is taken away from the protein.

Combined Antibodies Evusheld against the SARS-CoV-2 Omicron Variants BA.1.1 and BA.5: Immune Escape Mechanism from Molecular Simulation.

The Omicron lineage of SARS-CoV-2, which was first reported in November 2021, has spread globally and become dominant, splitting into several sublineages. Experiments have shown that Omicron lineage has escaped or reduced the activity of existing monoclonal antibodies, but the origin of escape mechanism caused by mutation is still unknown. This work uses molecular dynamics and umbrella sampling methods to reveal the escape mechanism of BA.1.1 to monoclonal antibody (mAb) Tixagevimab (AZD1061) and BA.5 to mAb Cilgavimab (AZD8895), both mAbs were combined to form antibody cocktail, Evusheld (AZD7442). The binding free energy of BA.1.1-AZD1061 and BA.5-AZD8895 has been severely reduced due to multiple-site mutated Omicron variants. Our results show that the two Omicron variants, which introduce a substantial number of positively charged residues, can weaken the electrostatic attraction between the receptor binding domain (RBD) and AZD7442, thus leading to a decrease in affinity. Additionally, using umbrella sampling along dissociation pathway, we found that the two Omicron variants severely impaired the interaction between the RBD of SARS-CoV-2's spike glycoprotein (S protein) and complementary determining regions (CDRs) of mAbs, especially in CDR3H. Although mAbs AZD8895 and AZD1061 are knocked out by BA.5 and BA.1.1, respectively, our results confirm that the antibody cocktail AZD7442 retains activity against BA.1.1 and BA.5 because another antibody is still on guard. The study provides theoretical insights for mAbs interacting with BA.1.1 and BA.5 from both energetic and dynamic perspectives, and we hope this will help in developing new monoclonals and combinations to protect those unable to mount adequate vaccine responses.

Diversity and composition of island insects of the Bejaia region (Algeria).

The study was conducted on four islands located west of Bejaia (north-eastern Algeria) to collect baseline information about diversity and composition of insects. Data were collected during summer and spring of 2015 using traps, sweep nets and the Japanese umbrella sampling methods. A total of 1691 individuals belonging to 112 species were collected. Totals of 72, 56, 41 and 34 species were recorded in Pisans Island, El Euch Island, Sahel Islet and Ail Islet, respectively. Of the 1,691 individuals recorded, 39.85 %, 30.51 %, 18.03 % and 11.59 % were recorded in Pisans Island, El Euch Island, Sahel islet and Ail islet, respectively. Shannon-weaver index ranged from 2.96 to 3.34 and Equitability between 0.77 and 0.87. The composition of the insect community differed between sampled islands (NMDS, Stress =. 0.013). Eleven (11) insect species composed over 50 % of dissimilarity between islands. Two species (Cassida viridis (L.) and Pheidole pallidula (Nyl.) accounted for more than 5% of the diversity on all islands. Sorensen's similarity coefficient showed that Pisans and El Euch islands have the highest similarity (45.37 %), while the lowest similarity was recorded between Pisans Island and Sahel Islet (15.53 %).

Understanding the sources of error in MBAR through asymptotic analysis.

Many sampling strategies commonly used in molecular dynamics, such as umbrella sampling and alchemical free energy methods, involve sampling from multiple states. The Multistate Bennett Acceptance Ratio (MBAR) formalism is a widely used way of recombining the resulting data. However, the error of the MBAR estimator is not well-understood: previous error analyses of MBAR assumed independent samples. In this work, we derive a central limit theorem for MBAR estimates in the presence of correlated data, further justifying the use of MBAR in practical applications. Moreover, our central limit theorem yields an estimate of the error that can be decomposed into contributions from the individual Markov chains used to sample the states. This gives additional insight into how sampling in each state affects the overall error. We demonstrate our error estimator on an umbrella sampling calculation of the free energy of isomerization of the alanine dipeptide and an alchemical calculation of the hydration free energy of methane. Our numerical results demonstrate that the time required for the Markov chain to decorrelate in individual states can contribute considerably to the total MBAR error, highlighting the importance of accurately addressing the effect of sample correlation.