Molecular Dynamics of Water at the Protein−Solvent Interface

Water at the nanoscale: From filling or dewetting hydrophobic pores and carbon nanotubes to "sliding" on graphene

In this work, we study the effect of nanoconfinement on the hydration properties of model hydrophobic pores and carbon nanotubes, determining their wetting propensity and the conditions for geometrically induced dehydration. By employing a recently introduced water structural index, we aim at two main goals: (1) to accurately quantify the local hydrophobicity and predict the drying transitions in such systems, and (2) to provide a molecular rationalization of the wetting process. In this sense, we will further discuss the number and strength of the interactions required by the water molecules to promote wetting. In the case of graphene-like surfaces, an explanation for their unexpectedly significant hydrophilicity will also be provided. On the one hand, the structural index will show that the net attraction to the dense carbon network that a water molecule experiences through several simultaneous weak interactions is sufficient to give rise to hydrophilic behavior. On the other hand, we will show that an additional effect is also at play: the hydrating water molecule is retained on the surface by a smooth exchange of such simultaneous weak interactions, as if "sliding" on graphene.

Similarities and Differences in Ligand Binding to Protein and RNA Targets: The Case of Riboflavin

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Design of potential anti-cancer agents as COX-2 inhibitors, using 3D-QSAR modeling, molecular docking, oral bioavailability proprieties, and molecular dynamics simulation

Modeling the structural properties of novel morpholine-bearing 1, 5-diaryl-diazole derivatives as potent COX-2 inhibitor, two proposed models based on CoMFA and CoMSIA were evaluated by external and internal validation methods. Partial least squares analysis produced statistically significant models with Q 2 values of 0.668 and 0.652 for CoMFA and CoMSIA, respectively, and also a significant non-validated correlation coefficient R² with values of 0.882 and 0.878 for CoMFA and CoMSIA, respectively. Both models met the requirements of Golbraikh and Tropsha, which means that both models are consistent with all validation techniques. Analysis of the CoMFA and CoMSIA contribution maps and molecular docking revealed that the R1 substituent has a very significant effect on their biological activity. The most active molecules were evaluated for their thermodynamic stability by performing MD simulations for 100 ns; it was revealed that the designed macromolecular ligand complex with 3LN1 protein exhibits a high degree of structural and conformational stability. Based on these results, we predicted newly designed compounds, which have acceptable oral bioavailability properties and would have high synthetic accessibility.

AI-Driven Design System for Fabrication of Inhalable Nanocatchers for Virus Capture and Neutralization

The global pandemic presents a critical threat to humanity, with no effective rapid-response solutions for early-stage virus dissemination. This study aims to create an AI-driven entry-blocker design system (AIEB) to fabricate inhalable virus-like nanocatchers (VLNCs) fused with entry-blocking peptides (EBPs) to counter pandemic viruses and explore therapeutic applications. This work focuses on developing angiotensin-converting enzyme 2 (ACE2)-mimic domain-fused VLNCs (ACE2@VLNCs) using AIEB and analyzing their interaction with the SARS-CoV-2 receptor binding domain (RBD), demonstrating their potential to hinder SARS-CoV-2 infection. Aerosol-based tests show ACE2@VLNCs persist over 70 min in the air and neutralize pseudoviruses within 30 min, indicating their utility in reducing airborne virus transmission. In vivo results reveal ACE2@VLNCs mitigate over 67% of SARS-CoV-2 infections. Biosafety studies confirm their safety, causing no damage to eyes, skin, lungs, or trachea, and not eliciting significant immune responses. These findings offer crucial insights into pandemic virus prevention and treatment, highlighting the potential of the ACE2@VLNCs system as a promising strategy against future pandemics.

pH-dependent interactions of biologically important metal ions with hen egg white lysozyme based on its hydration properties: Thermodynamic and mechanistic insights

Importance of metal ion selectivity in biomolecules and their key role in proteins are widely explored. However, understanding the thermodynamics of how hydrated metal ions alter the protein hydration and their conformation is also important. In this study, the interaction of some biologically important Ca2+, Mn2+, Co2+, Cu2+, and Zn2+ ions with hen egg white lysozyme at pH 2.1, 3.0, 4.5 and 7.4 has been investigated. Intrinsic fluorescence studies have been employed for metal ion-induced protein conformational changes analysis. Thermostability based on protein hydration has been investigated using differential scanning calorimetry (DSC). Thermodynamic parameters emphasizing on metal ion-protein binding mechanistic insights have been well discussed using isothermal titration calorimetry (ITC). Overall, these experiments have reported that their interactions are pH-dependent and entropically driven. This research also reports the strongly hydrated metal ions as water structure breaker unlike osmolytes based on DSC studies. These experimental results have highlighted higher concentrations of different metal ions effect on the protein hydration and thermostability which might be helpful in understanding their interactions in aqueous solutions.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Temperature-dependent dielectric relaxation and hydrophobicity of aqueous alanine using time domain reflectometry

Physical, chemical and microbiological stability of the materials is affected by the rotational and translational mobility of free and hydrated water. The role of water in areas such as protein hydration and enzyme activity, food technology, lyophilization and polymers hydration is, therefore, important and can be well understood in terms of dielectric relaxation spectroscopy. Concentration and temperature-dependent hydrophobicity of amino acid is reflected in their tendencies to appear in appropriate positions in proteins. Therefore, to gain more insights on the temperature and concentration dependence of hydrophobicity and structural properties of amino acid, dielectric relaxation of aqueous alanine have been studied in the temperature region 303.15 K to 278.15 K. Time domain spectroscopy have been used in the frequency range of 10 MHz to 30 GHz and in the concentration range 0.18708 ≤ c/M ≤ 0.74831. Two relaxation processes namely the low-frequency relaxation (l) and the high-frequency relaxation (h) has been detected for the aqueous alanine. Dielectric parameters such as static dielectric constant (εj), relaxation time (τj) dipole moments (û) and correlation factor (g) have been studied to investigate molecular interaction between alanine and water. The number of water molecules irrotationally bond by the solute molecules (Zib) was also determined to examine the hydrophobicity of alanine which was found more hydrophobic towards low temperatures and concentrations. Thermodynamic parameters calculated are also supported well for the hydrophobic behaviour of alanine towards low temperatures and concentrations.Communicated by Ramaswamy H. Sarma.

Water molecule ordering on the surface of an intrinsically disordered protein

The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.

Water inside the selectivity filter of a K+ ion channel: structural heterogeneity, picosecond dynamics, and hydrogen-bonding

Water inside biological ion channels regulates the key properties of these proteins such as selectivity, ion conductance, and gating. In this Article we measure the picosecond spectral diffusion of amide I vibrations of an isotope labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100 - 2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D lineshapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent, or non-adjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations were observed on a picosecond timescale. These dynamics are in stark contrast with liquid water that remains highly dynamic even in the presence of ions at high concentrations.

Counteraction Effects of Ammonium-Based Ionic Liquids on Urea-Induced Denaturation of α-Lactalbumin: A Comprehensive Molecular Simulation Study

Ionic liquids (ILs) are known to stabilize protein conformations in aqueous medium. Importantly, ILs can also act as refolding additives in urea-driven denaturation of proteins. However, despite the importance of the problem, detailed microscopic understanding of the counteraction effects of ILs on urea-induced protein denaturation remains elusive. In this work, atomistic molecular dynamics (MD) simulations of the protein α-lactalbumin have been carried out in pure aqueous medium, in 8 M binary urea-water solution and in ternary urea-IL-water solutions containing ammonium-based ethyl ammonium acetate (EAA) as the IL at different concentrations (1-4 M). Attempts have been made to quantify detailed molecular-level understanding of the origin behind the counteraction effects of the IL on urea-induced partial unfolding of the protein. The calculations revealed significant conformational changes of the protein with multiple free energy minima due to its partial unfolding in binary urea-water solution. The counteraction effect of the IL was evident from the enhanced structural rigidity of the protein with propensity to transform into a single native free energy minimum state in ternary urea-IL-water solutions. Such an effect has been found to be associated with preferential direct binding of the IL components with the protein and simultaneous expulsion of urea from the interface, thereby providing additional stabilization of the protein in ternary solutions. Most importantly, modified rearrangement of the hydrogen bond network at the interface due to the formation of stronger protein-cation (PC) and protein-anion (PA) hydrogen bonds by breaking relatively weaker protein-urea (PU) and protein-water (PW) hydrogen bonds has been recognized as the microscopic origin behind the counteraction effects of EAA on urea-induced partial unfolding of the protein.

Variations in the Protein Hydration and Hydrogen-Bond Network of Water Molecules Induced by the Changes in the Secondary Structures of Proteins Studied through Near-Infrared Spectroscopy

This study investigated how the secondary structural changes of proteins in aqueous solutions affect their hydration and the hydrogen-bond network of water molecules using near-infrared (NIR) spectroscopy. The aqueous solutions of three types of proteins, i.e., ovalbumin, β-lactoglobulin, and bovine serum albumin, were denatured by heating, and changes in the NIR bands of water reflecting the states of hydrogen bonds induced via protein secondary structural changes were investigated. On heating, the intermolecular hydrogen bonds between water molecules as well as between water and protein molecules were broken, and protein molecules were no longer strongly bound by the surrounding water molecules. Consequently, the denaturation was observed to proceed depending on the thermodynamic properties of the proteins. When the aqueous solutions of proteins were cooled after denaturation, the hydrogen-bond network was reformed. However, the state of protein hydration was changed owing to the secondary structural changes of proteins, and the variation patterns were different depending on the protein species. These changes in protein hydration may be derived from the differences in the surface charges of proteins. The elucidation of the mechanism of protein hydration and the formation of the hydrogen-bond network of water molecules will afford a comprehensive understanding of the protein functioning and dysfunctioning derived from the structural changes in proteins.

Correlations between defect propensity and dynamical heterogeneities in supercooled water

A salient feature of supercooled liquids consists in the dramatic dynamical slowdown they undergo as temperature decreases while no significant structural change is evident. These systems also present dynamical heterogeneities (DH): certain molecules, spatially arranged in clusters, relax various orders of magnitude faster than the others. However, again, no static quantity (such as structural or energetic measures) shows strong direct correlations with such fast-moving molecules. In turn, the dynamic propensity approach, an indirect measure that quantifies the tendency of the molecules to move in a given structural configuration, has revealed that dynamical constraints, indeed, originate from the initial structure. Nevertheless, this approach is not able to elicit which structural quantity is, in fact, responsible for such a behavior. In an effort to remove dynamics from its definition in favor of a static quantity, an energy-based propensity has also been developed for supercooled water, but it could only find positive correlations between the lowest-energy and the least-mobile molecules, while no correlations could be found for those more relevant mobile molecules involved in the DH clusters responsible for the system's structural relaxation. Thus, in this work, we shall define a defect propensity measure based on a recently introduced structural index that accurately characterizes water structural defects. We shall show that this defect propensity measure provides positive correlations with dynamic propensity, being also able to account for the fast-moving molecules responsible for the structural relaxation. Moreover, time dependent correlations will show that defect propensity represents an appropriate early-time predictor of the long-time dynamical heterogeneity.

Design of novel anti-cancer drugs targeting TRKs inhibitors based 3D QSAR, molecular docking and molecular dynamics simulation

Tropomyosin receptor kinase (TRK) enzymes are responsible for different types of tumors caused by neurotrophic tyrosine receptor kinase gene fusion and have been identified as an effective target for anticancer therapy. The study of the mechanism between polo-like kinase (PLKs) and pyrazol inhibitors was performed using 3D-QSAR modeling, molecular docking, and MD simulations in order to design high-activity inhibitors. The HQSAR (Q2 = 0.793, R2 = 0.917, R2ext = 0.961), CoMFA (Q2 = 0.582, R2 = 0.722, R2ext = 0.951), CoMSIA/SE (Q2 = 0.603, R2 = 0.801, R2ext = 0.849), and Topomer CoMFA (Q2 = 0.726, R2 = 0.992, R2ext = 0.717) showed good reliability and predictability. All models have been successfully tested by external validation, so all five established models are reliable. The analysis of the different contour maps of different models gives structural information to improve the inhibitory function. Molecular docking results show that the amino acids Met 592, GLU 590, LEU 657, VAL 524, and PHE 589 are the active sites of the tropomyosin receptor TRKs. The results obtained by MD showed that compound 19i could form a more stable complex protein (PDB id: 5KVT). Based on these results, we developed new compounds and their expected inhibitory activities. The results of physicochemical and ADME-Tox properties showed that the four proposed molecules are orally bioavailable, and they are not toxic in the Ames test. Thus, these results would provide modeling information that could help experimental researchers find TRK type I inhibitors more efficiently.Communicated by Ramaswamy H. Sarma.

Elucidating the Sluggish Water Dynamics at the Ice-Binding Surface of the Hyperactive Tenebrio molitor Antifreeze Protein

Quasi-ice-like hydration waters on the ice-binding surface (IBS) of an antifreeze protein (AFP) commonly exhibit sluggish dynamics especially at low temperatures. In this work, we have analyzed molecular dynamics (MD) simulation trajectories at two different temperatures for Tenebrio molitor antifreeze protein (TmAFP) to explore whether the unique quasi-ice-like structuring of hydration water has any impact on making their dynamics slower on the IBS of the protein. Our calculation reveals that, as translational dynamics is coupled with the conformational fluctuations, hydration water on the IBS exhibits sluggish translational motion due to reduced flexibility of the IBS compared to that on the non-ice-binding surface (NIBS) of the protein. Interestingly, it is noticed that rotational motion of hydration water is not coupled with the conformational fluctuations of the surfaces. In that case, structural relaxations of the protein-water (PW) and water-water (WW) hydrogen bonds compete with each other to make the rotational dynamics of hydration water around the IBS either faster or slower with respect to those around the NIBS. At low temperature, the slower structural relaxation of water-water hydrogen bonds dominates and imparts sluggish rotational motion of the hydration water on the IBS of the protein. The slower structural relaxation of water-water hydrogen bonds and hence the retarded rotational dynamics, despite the weak short-lived PW hydrogen bonds on the IBS, is clearly a manifestation of the rigid quasi-ice-like structure of the hydration shell on that surface.

Structure, Thermodynamics, and Dynamics of Thiocyanate Ion Adsorption and Transfer across the Water/Toluene Interface

Molecular dynamics simulations are used to examine in detail the structure, thermodynamics, and dynamics involve in the adsorption and transfer of the thiocyanate ion (SCN^-) across the water/toluene interface. Free energy, hydration structure, and several dynamical properties as a function of the ion location along the interface normal are calculated and contrasted with recent experiments. The free energy profile exhibits a local minimum near the interface corresponding to adsorption free energy relative to bulk water of -6 kJ/mol, in reasonable agreement with experiments. The simulations provide insight into the water surface fluctuations that are coupled to the ion transfer, demonstrating formation of water finger-like structures assisting the transfer process.

Exploring Heterogeneous Dynamical Environment around an Ensemble of Aβ42 Peptide Monomer Conformations

Exploring the conformational properties of amyloid β (Aβ) peptides and the role of solvent (water) in guiding the dynamical environment at their interfaces is crucial for microscopic understanding of Aβ misfolding, which is involved in causing the most common neurodegenerative disorder, i.e., Alzheimer's disease. While numerous studies in the past have emphasized examining the conformational states of Aβ peptides, the role of water has not received much attention. Here, we have performed all-atom molecular dynamics simulations of several full-length Aβ₄₂ peptide monomers with different initial configurations. Our efforts are directed toward probing the origin of the heterogeneous dynamics of water around various segments of the Aβ peptide, identified as the two terminal segments (N-term and C-term) and the two hydrophobic segments (hp1 and hp2), along with the central turn region interconnecting hp1 and hp2. Our results revealed that water hydrating hp1, hp2, and turn (nonterminal segments) and C-term segments exhibit nonuniformly restricted translational as well as rotational motions. The degree of such restriction has been found to be correlated with the hydrogen bond relaxation time scales at the interface. Importantly, it is revealed that the water molecules around hp1 and, to some extent, around hp2, form relatively rigid hydration layers, compared to that around the other segments. Such rigid hydration layers arise due to relatively more solid-like caging motions resulting in relatively lesser hydration entropy. As hp1 and hp2 have been demonstrated to play a central role in Aβ aggregation, we believe that distinct water dynamics in the vicinity of these two segments, as outlined in this study, can provide vital information in understanding the early stages of the onset of the aggregation process of such peptides at higher concentration that can further aid toward advances in AD therapeutics.

Structural, dynamic, and hydration properties of quercetin and its aggregates in solution

Quercetin is a flavonoid present in the human diet with multiple health benefits. Quercetin solutions are inhomogeneous even at very low concentrations due to quercetin's tendency to aggregate. We simulate, using molecular dynamics, three systems of quercetin solutions: infinite dilution, 0.22 M, and 0.46 M. The systems at the two highest concentrations represent regions of the quercetin aggregates, in which the concentration of this molecule is unusually high. We study the behavior of this molecule, its aggregates, and the modifications in the surrounding water. In the first three successive layers of quercetin hydration, the density of water and the hydrogen bonds formations between water molecules are smaller than that of bulk. Quercetin has a hydrophilic surface region that preferentially establishes donor hydrogen bonds with water molecules with relative frequencies from 0.12 to 0.46 at infinite dilution. Also, it has two hydrophobic regions above and below the planes of its rings, whose first hydration layers are further out from quercetin (≈0.3 Å) and their water molecules do not establish hydrogen bonds with it. Water density around the hydrophobic regions is smaller than that of the hydrophilic. Quercetin molecules aggregate inπ-stacking configurations, with a distance of ≈0.37 nm between the planes of their rings, and form bonds between their hydroxyl groups. The formation of quercetin aggregates decreases the hydrogen bonds between quercetin and the surrounding water and produces a subdiffusive behavior in water molecules. Quercetin has a subdiffusive behavior even at infinite dilution, which increases with the number of molecules within the aggregates and the time they remain within them.

Dynamic Heterogeneity at the Interface of an Intrinsically Disordered Peptide

It is believed that water around an intrinsically disordered protein or peptide (IDP) in an aqueous environment plays an important role in guiding its conformational properties and aggregation behavior. However, despite its importance, only a handful of studies exploring the correlation between the conformational motions of an IDP and the microscopic properties of water at its surface are reported. Attempts have been made in this work to study the dynamic properties of water present in the vicinity of α-synuclein, an IDP associated with Parkinson's disease (PD). Room temperature molecular dynamics (MD) simulations of eight α-synuclein_1-95 peptides with a wide range of initial conformations have been carried out in aqueous media. The calculations revealed that due to solid-like caging motions, the translational and rotational mobility of water molecules near the surfaces of the peptide repeat unit segments R1 to R7 are significantly restricted. A small degree of dynamic heterogeneity in the hydration environment around the repeat units has been observed with water near the hydrophobic R6 unit exhibiting relatively more restricted diffusivity. The time scales involving the overall structural relaxations of peptide-water and water-water hydrogen bonds near the peptide have been found to be correlated with the time scale of diffusion of the interfacial water molecules. We believe that the relatively more hindered dynamic environment near R6 can help create water-mediated contacts centered around R6 between peptide monomers at a higher concentration, thereby enhancing the early stages of peptide aggregation.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

Interactions of Co, Cu, and non-metal phthalocyanines with external structures of SARS-CoV-2 using docking and molecular dynamics

The new coronavirus, SARS-CoV-2, caused the COVID-19 pandemic, characterized by its high rate of contamination, propagation capacity, and lethality rate. In this work, we approach the use of phthalocyanines as an inhibitor of SARS-CoV-2, as they present several interactive properties of the phthalocyanines (Pc) of Cobalt (CoPc), Copper (CuPc) and without a metal group (NoPc) can interact with SARS-CoV-2, showing potential be used as filtering by adsorption on paints on walls, masks, clothes, and air conditioning filters. Molecular modeling techniques through Molecular Docking and Molecular Dynamics were used, where the target was the external structures of the virus, but specifically the envelope protein, main protease, and Spike glycoprotein proteases. Using the g_MM-GBSA module and with it, the molecular docking studies show that the ligands have interaction characteristics capable of adsorbing the structures. Molecular dynamics provided information on the root-mean-square deviation of the atomic positions provided values between 1 and 2.5. The generalized Born implicit solvation model, Gibbs free energy, and solvent accessible surface area approach were used. Among the results obtained through molecular dynamics, it was noticed that interactions occur since Pc could bind to residues of the active site of macromolecules, demonstrating good interactions; in particular with CoPc. Molecular couplings and free energy showed that S-gly active site residues interacted strongly with phthalocyanines with values ​​of - 182.443 kJ/mol (CoPc), 158.954 kJ/mol (CuPc), and - 129.963 kJ/mol (NoPc). The interactions of Pc's with SARS-CoV-2 may predict some promising candidates for antagonists to the virus, which if confirmed through experimental approaches, may contribute to resolving the global crisis of the COVID-19 pandemic.

Influence of Aqueous Arginine Solution on Regulating Conformational Stability and Hydration Properties of the Secondary Structural Segments of a Protein at Elevated Temperatures: A Molecular Dynamics Study

The effects of aqueous arginine solution on the conformational stability of the secondary structural segments of a globular protein, ubiquitin, and the structure and dynamics of the surrounding water and arginine were examined by performing atomistic molecular dynamics (MD) simulations. Attempts have been made to identify the osmolytic efficacy of arginine solution, and its influence in guiding the hydration properties of the protein at an elevated temperature of 450 K. The similar properties of the protein in pure water at elevated temperatures were computed and compared. Replica exchange MD simulation was performed to explore the arginine solution's sensitivity in stabilizing the protein conformations for a wide range of temperatures (300-450 K). It was observed that although all the helices and strands of the protein undergo unfolding at elevated temperature in pure water, they exhibited native-like conformational dynamics in the presence of arginine at both ambient and elevated temperatures. We find that the higher free energy barrier between the folded native and unfolded states of the protein primarily arises from the structural transformation of α-helix, relative to the strands. Our study revealed that the water structure around the secondary segments depends on the nature of amino acid compositions of the helices and strands. The reorientation of water dipoles around the helices and strands was found hindered due to the presence of arginine in the solution; such hindrance reduces the possibility of exchange of hydrogen bonds that formed between the secondary segments of protein and water (PW), and as a result, PW hydrogen bonds take longer time to relax than in pure water. On the other hand, the origin of slow relaxation of protein-arginine (PA) hydrogen bonds was identified to be due to the presence of different types of protein-bound arginine molecules, where arginine interacts with the secondary structural segments of the protein through multiple/bifurcated hydrogen bonds. These protein-bound arginine formed different kinds of bridged PA hydrogen bonds between amino acid residues of the same secondary segments or among multiple bonds and helped protein to conserve its native folded form firmly.

Site-Specific Water Dynamics in the First Hydration Layer of an Anti-Freeze Glyco-Protein: A Simulation Study

Antifreeze glycoproteins (AFGPs) inhibit ice recrystallization by a mechanism remaining largely elusive. Dynamics of AFGPs’ hydration water and its involvement in the antifreeze activity, for instance, have not been identified...

Prolonged Association between Water Molecules under Hydrophobic Nanoconfinement

We present an investigation of the dynamics of water confined among rigid carbon rods and between parallel graphene sheets with molecular dynamics simulations. Diffusion coefficients, activation energy of diffusion, and residence-time correlation functions as a function of confinement geometry reveal a retardation of water dynamics under hydrophobic confinement compared to bulk water. In fact, water under various confinements possesses longer associations with its neighbors and exhibits diffusion dynamics characteristic of a lower temperature. Analysis of the residence-time correlation functions reveals long and short residence times, which we relate to the diffusion coefficient and activation energy of diffusion, respectively. Additional investigations reveal how the level of confining surface hydrophobicity affects water dynamics, further broadening our understanding of water diffusion inside diverse media. Overall, this study sheds light on the physical origin of retarded water dynamics under hydrophobic confinement and the close relationship between residence times and diffusion behavior.

Digital Microfluidic Thermal Control Chip-Based Multichannel Immunosensor for Noninvasively Detecting Acute Myocardial Infarction

Rapid and automated detection of acute myocardial infarction (AMI) at its developing stage is very important due to its high mortality rate. To quantitatively diagnose AMI, Myo, CK-MB, and cTnI are chosen as three biomarkers, which are usually detected through an immunosorbent assay, such as the enzyme-linked immunosorbent assay. However, the approach poses many drawbacks, such as long detection time, the cumbersome process, the need for professionals, and the difficulty of realizing automatic operation. Here, a multichannel digital microfluidic (DMF) thermal control chip integrated with a sandwich-based immunoassay strategy is proposed for the automated, rapid, and sensitive detection of AMI biomarkers. A miniaturized temperature control module is integrated on the back of the DMF chip, meeting the temperature requirement for the immunoassay. With this DMF thermal control chip, sample and reagent consumption are reduced to several microliters, significantly alleviating reagent consumption and sample dependence, and the automated and multichannel detection of biomarkers can be achieved. In this work, the simultaneously noninvasive detection of the human serum sample containing the three biomarkers of AMI is also achieved within 30 min, which improves the diagnostic accuracy of AMI. Due to the features of automation and miniaturization, the multichannel immunosensor can be used in community hospitals to increase the speed of diagnosis of patients with various acute diseases.

Structure and dynamics of nanoconfined water and aqueous solutions

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Cyclosporin Structure and Permeability: From A to Z and Beyond

Cyclosporins are natural or synthetic undecapeptides with a wide range of actual and potential pharmaceutical applications. Several members of the cyclosporin compound family have remarkably high passive membrane permeabilities that are not well-described by simple structural metrics. Here we review experimental studies of cyclosporin structure and permeability, including cyclosporin-metal complexes. We also discuss models for the conformation-dependent permeability of cyclosporins and similar compounds. Finally, we identify current knowledge gaps in the literature and provide recommendations regarding future avenues of exploration.

Effect of charged mutation on aggregation of a pentapeptide: Insights from molecular dynamics simulations

Aggregation of therapeutic monoclonal antibodies (mAbs) can negatively affect their chemistry, manufacturing, and control attributes and lead to undesirable immune responses in patients. Therefore, optimization of lead mAb drug candidates during discovery stages to mitigate aggregation is increasingly becoming an integral part of their developability assessments. The disruption of short sequence motifs called aggregation prone regions (APRs) found in amino acid sequences of mAb candidates can potentially mitigate their aggregation. In this work, we have performed molecular dynamics simulations to study the aggregation of an APR (VLVIY) found in λ light chains of human antibodies and its single point mutant KLVIY. Eighteen different multicopy peptide simulation systems of "VLVIY" and "KLVIY" were constructed by varying their concentrations, temperatures, termini capping, and flanking gate-keeper regions. Within 20 ns of the simulation, peptide "VLVIY" formed an aggregate of 100 peptides at ~0.1 M concentration with a 60% reduction in solvent accessible surface area (SASA). Furthermore, analysis of the SASA change, peptide cluster distribution, and water residence time demonstrated how Val ➔ Lys mutation resists aggregation and improves solubility. Presence of Lys slows down aggregation kinetics via charge-charge repulsions and by raising the kinetic barrier to formation of large oligomers. However, the effect of the Val ➔ Lys mutation is dependent on sequence and structural contexts around the APR. This mutation also alters the solvation shell around the peptide by favoring solute-solvent interactions, thereby increasing its solubility. This work has provided a detailed mechanistic explanation of how APR disruption can mitigate aggregation in biotherapeutics and improve their developability.

Solvent Relaxation NMR: A Tool for Real-Time Monitoring Water Dynamics in Protein Aggregation Landscape

Solvent dynamics strongly induce the fibrillation of an amyloidogenic system. Probing the solvation mechanism is crucial as it enables us to predict different proteins' functionalities, such as the aggregation propensity, structural flexibility, and toxicity. This work shows that a straightforward NMR method in conjunction with phenomenological models gives a global and qualitative picture of water dynamics at different concentrations and temperatures. Here, we study amyloid system Aβ40 and its fragment AV20 (A21-V40) and G37L (mutation at Gly37 → Leu of AV20), having different aggregation and toxic properties. The independent validation of this method is elucidated using all-atom classical MD simulation. These two state-of-the-art techniques are pivotal in linking the effect of solvent environment in the near hydration-shell to their aggregation nature. The time-dependent modulation in solvent dynamics probed with the NMR solvent relaxation method can be further adopted to gain insight into amyloidogenesis and link with their toxicity profiles.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Structural aspects of an energy-based water classification index and the structure-dynamics link in glassy relaxation

An energy-based structural indicator for water, [Formula: see text], has been recently introduced by our group. In turn, in this work we aim at: (1) demonstrating that [Formula: see text] is indeed able to correctly classify water molecules between locally structured tetrahedral (T) and locally distorted (D) ones, circumventing the usual problem of certain previous indicators of overestimating the distorted state; (2) correlating [Formula: see text] with dynamic propensity, a measure of the molecular mobility tendency, in order to seek for the existence of a connection between structure and dynamics within the supercooled regime. More specifically, in the first part of this work we will show that [Formula: see text] accurately discriminates between merely thermally deformed local molecular arrangements and truly distorted molecules (defects). This fact will be made evident not only from radial distribution function results but also from the dynamic propensity distributions of the different kinds of molecules. In turn, we shall devote the second part of this work to finding correlations between T and D molecules with low- and high-dynamic-propensity molecules, respectively, thus revealing the existence of a link between local structure and dynamics, while also making evident the dominant role of the D molecules (defects) in the structural relaxation. Moreover, the availability of a proper molecular classification technique will enable us to study the timescale of such influence of structure on dynamics by defining a modified dynamic propensity measure and by applying it to the structured and unstructured water molecular states.

Small Glycols Discover Cryptic Pockets on Proteins for Fragment-Based Approaches

Cryptic pockets are visible in ligand-bound protein structures but are occluded in unbound structures. Utilizing these pockets in fragment-based drug-design provides an attractive option for proteins not tractable by classical binding sites. However, owing to their hidden nature, they are difficult to identify. Here, we show that small glycols find cryptic pockets on a diverse set of proteins. Initial crystallography experiments serendipitously revealed the ability of ethylene glycol, a small glycol, to identify a cryptic pocket on the W6A mutant of the RBSX protein (RBSX-W6A). Explicit-solvent molecular dynamics (MD) simulations of RBSX-W6A with the exposed state of the cryptic pocket (ethylene glycol removed) revealed closure of the pocket reiterating that the exposed state of cryptic pockets in general are unstable in the absence of ligands. Also, no change in the pocket was observed for simulations of RBSX-W6A with the occluded state of the cryptic pocket, suggesting that water molecules are not able to open the cryptic pocket. "Cryptic-pocket finding" potential of small glycols was then supported and generalized through additional crystallography experiments, explicit-cosolvent MD simulations, and protein data set construction and analysis. The cryptic pocket on RBSX-W6A was found again upon repeating the crystallography experiments with another small glycol, propylene glycol. Use of ethylene glycol as a probe molecule in cosolvent MD simulations led to the enhanced sampling of the exposed state of experimentally observed cryptic sites on a test set of two proteins (Niemann-Pick C2, Interleukin-2). Further, analyses of protein structures with validated cryptic sites showed that ethylene glycol molecules bind to sites on proteins (Bcl-xL, G-actin, myosin II, and glutamate receptor 2), which become apparent upon binding of biologically relevant ligands. Our study thus suggests potential application of the small glycols in experimental and computational fragment-based approaches to identify cryptic pockets in apparently undruggable and/or difficult targets, making these proteins amenable to drug-design strategies.

Unraveling the origin of interactions of hydroxychloroquine with the receptor-binding domain of SARS-CoV-2 in aqueous medium

Interactions of hydroxychloroquin (HCQ) with the receptor binding domain (RBD) of SARS-CoV-2 were studied from atomistic simulation and ONIOM techniques. The key-residues of RBD responsible for the human transmission are recognized to be blocked in a heterogeneous manner with the favorable formation of key-residue:HCQ (1:1) complex. Such heterogeneity in binding was identified to be governed by the differential life-time of the hydrogen bonded water network anchoring HCQ and the key-residues. The intermolecular proton transfer facilitates the most favorable Lys417:HCQ complexation. The study demonstrates that off-target bindings of HCQ need to be minimized to efficiently prevent the transmission of SARS-CoV-2.

Time-Resolved Fluorescence and Essential Dynamics Study on the Structural Heterogeneity of p53DBD Bound to the Anticancer p28 Peptide

Time-resolved fluorescence emission was combined with molecular dynamics (MD) simulations to investigate the DNA-binding domain (DBD) of the tumor suppressor p53 alone and its complex with the anticancer peptide p28 (DBD/p28). The fluorescence emission decay of the lone Trp residue, from both DBD and DBD/p28, was well-described by a stretched exponential function. Such a behavior was ascribed to heterogeneity in the Trp relaxation behavior, likely due to the coexistence of different conformational states. The increase of the stretching parameter, on passing from DBD to DBD/p28, indicates a reduced heterogeneity in the Trp146 environment for DBD/p28. Moreover, the effects of p28 on the global dynamics of DBD were analyzed by the essential dynamics method on 30 ns long MD trajectories of both DBD and DBD/p28. We found the establishment of wide-amplitude anharmonic modes throughout the DBD molecule, with a particularly high amplitude being detected in the DNA-binding region. These modes are significantly reduced when DBD is bound to p28, consistently with a structure stabilization. In summary, the results indicate that p28 binding has a strong effect on both the local and global heterogeneity of DBD, thus providing some hints to the understanding of its anticancer activity.

First-passage fingerprints of water diffusion near glutamine surfaces

The extent to which biological interfaces affect the dynamics of water plays a key role in the exchange of matter and chemical interactions that are essential for life. The density and the mobility of water molecules depend on their proximity to biological interfaces and can play an important role in processes such as protein folding and aggregation. In this work, we study the dynamics of water near glutamine surfaces-a system of interest in studies of neurodegenerative diseases. Combining molecular-dynamics simulations and stochastic modelling, we study how the mean first-passage time and related statistics of water molecules escaping subnanometer-sized regions vary from the interface to the bulk. Our analysis reveals a dynamical complexity that reflects underlying chemical and geometrical properties of the glutamine surfaces. From the first-passage time statistics of water molecules, we infer their space-dependent diffusion coefficient in directions normal to the surfaces. Interestingly, our results suggest that the mobility of water varies over a longer length scale than the chemical potential associated with the water-protein interactions. The synergy of molecular dynamics and first-passage techniques opens the possibility for extracting space-dependent diffusion coefficients in more complex, inhomogeneous environments that are commonplace in living matter.

Anomalously Slow Conformational Change Dynamics of Polar Groups Anchored to Hydrophobic Surfaces in Aqueous Media

Water molecules within a thin hydration layer, spontaneously generated on hydrophobic protein surfaces, are reported to form a poorly dynamic network structure. However, how such a water network affects the conformational change dynamics of polar groups has never been explored, although such polar groups play a critical role in protein-protein and protein-ligand interactions. In the present work, we utilized as model protein surfaces a series of self-assembled monolayers (SAMs) appended with polar (Fmoc) or ionic (FITC) fluorescent head groups that were tethered via a 1.5-nm-long flexible oligoether chain to a hydrophobic silicon wafer surface, which was densely covered with paraffinic chains. We found that, not only in deionized water but also in aqueous buffer, these oligoether-appended head groups at ambient temperatures both displayed an anomalously slow conformational change, which required ∼10 h to reach a thermodynamically equilibrated state. We suppose that these behaviors reflect the poorly dynamic and low-permittivity natures of the thin hydration layer.

Protein-Water and Water-Water Long-Time Relaxations in Protein Hydration Water upon Cooling-A Close Look through Density Correlation Functions

We report results on the translational dynamics of the hydration water of the lysozyme protein upon cooling obtained by means of molecular dynamics simulations. The self van Hove functions and the mean square displacements of hydration water show two different temperature activated relaxation mechanisms, determining two dynamic regimes where transient trapping of the molecules is followed by hopping phenomena to allow to the structural relaxations. The two caging and hopping regimes are different in their nature. The low-temperature hopping regime has a time scale of tenths of nanoseconds and a length scale on the order of 2–3 water shells. This is connected to the nearest-neighbours cage effect and restricted to the supercooling, it is absent at high temperature and it is the mechanism to escape from the cage also present in bulk water. The second hopping regime is active at high temperatures, on the nanoseconds time scale and over distances of nanometers. This regime is connected to water displacements driven by the protein motion and it is observed very clearly at high temperatures and for temperatures higher than the protein dynamical transition. Below this temperature, the suppression of protein fluctuations largely increases the time-scale of the protein-related hopping phenomena at least over 100 ns. These protein-related hopping phenomena permit the detection of translational motions of hydration water molecules longly persistent in the hydration shell of the protein.

Molecular Dynamics Studies on the Effect of Surface Roughness and Surface Tension on the Thermodynamics and Dynamics of Hydronium Ion Transfer Across the Liquid/Liquid Interface

Molecular dynamics simulations are used to examine the effect of surface roughness and surface tension on the transfer of the classical hydronium ion (H₃O⁺) across the water/1,2-dichloroethane interface. Free energy of transfer, hydration structure, and dynamics as a function of the ion location along the interface normal are calculated with six different values of a control parameter whose variation modifies the surface tension without impacting the bulk properties of the two solvents. Transfer of the classical hydronium ion across the water/1,2-dichloroethan interface involves the cotransfer of three hydration shell water molecules, independent of the surface tension. However, as the interaction between the two liquids weakens, a rise in interfacial tension and decrease in intrinsic water fingering and capillary fluctuations result in fewer water molecules cotransported with the ion in the second shell and a reduction in the length of the finger that the ion is attached to, consistent with the reduced size of the second hydration shell. First shell water residence time and lateral ion diffusion constants vary with the surface tension in a way that is consistent with the abovementioned structural insight.

Insights into the Sensitivity of Arginine Concentration to Preserve the Folded Form of Insulin Monomer under Thermal Stress

Arginine, although popularly known as aggregation suppressor additive, has been found to quench proteins' structure and function by destabilizing their conformations. Driven by such controversial evidence, in this work we performed a series of atomistic molecular dynamics simulations of insulin monomer, a biologically active hormone protein, in arginine solution of varying concentrations (0.5, 1, and 2 M) at ambient and elevated temperature (400 K) to explore the arginine concentration driven structure-based stability of the protein. Our study reveals that the flexibility of the protein's structure is dependent on the arginine concentration, and among all the used solutions, 2 M arginine, a "neutral crowder" that mimics the cellular environment, can preserve the native folded form of the protein at ambient temperature in an excellent manner. Further, while the protein unfolds at 400 K in pure water, this solution worked satisfactorily to preserve the protein's folded conformation more firmly than the other solutions. The replica-exchange MD of insulin in 2 M arginine solution further supports the fact. In this aspect an important issue in molecular pharmacology is to identify and recognize the physical origin of the stability of a protein, i.e, in this case, how arginine directs the conformational flexibility of the protein and preserves its native folded form. We identified that the exclusion of arginine from the protein surface increases the local structuration of water around the protein, thereby preserving its "biological water" layer, and makes the protein more hydrated at 2 M concentration as compared to the other arginine solutions. Additionally, our microscopic investigation on the interactions of the protein-solvation layer revealed that the structural heterogeneity of the protein surface, arising from the differential physicochemical nature of the amino acid residues, controls the favorable formation of sluggish water-arginine mixed solvation layer at higher arginine concentration that helps the protein to maintain its structural rigidity. Importantly, apart from the protein-solvent hydrogen-bonding interactions, the anion-pi interactions, established between the carboxyl group of arginine and the aromatic amino acid residues of insulin, were recognized to facilitate the protein to maintain its native folded form at the experimental temperatures.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

Capturing the Effects of Explicit Waters in Implicit Electrostatics Modeling: Qualitative Justification of Gaussian-Based Dielectric Models in DelPhi

Our group has implemented a smooth Gaussian-based dielectric function in DelPhi (J. Chem. Theory Comput. 2013, 9 (4), 2126-2136) which models the solute as an object with inhomogeneous dielectric permittivity and provides a smooth transition of dielectric permittivity from surface-bound water to bulk solvent. Although it is well-understood that the protein hydrophobic core is less polarizable than the hydrophilic protein surface, less attention is paid to the polarizability of water molecules inside the solute and on its surface. Here, we apply explicit water simulations to study the behavior of water molecules buried inside a protein and on the surface of that protein and contrast it with the behavior of the bulk water. We selected a protein that is experimentally shown to have five cavities, most of which are occupied by water molecules. We demonstrate through molecular dynamics (MD) simulations that the behavior of water in the cavity is drastically different from that in the bulk. These observations were made by comparing the mean residence times, dipole orientation relaxation times, and average dipole moment fluctuations. We also show that the bulk region has a nonuniform distribution of these tempo-spatial properties. From the perspective of continuum electrostatics, we argue that the dielectric "constant" in water-filled cavities of proteins and the space close to the molecular surface should differ from that assigned to the bulk water. This provides support for the Gaussian-based smooth dielectric model for solving electrostatics in the Poisson-Boltzmann equation framework. Furthermore, we demonstrate that using a well-parametrized Gaussian-based model with a single energy-minimized configuration of a protein can also reproduce its ensemble-averaged polar solvation energy. Thus, we argue that the Gaussian-based smooth dielectric model not only captures accurate physics but also provides an efficient way of computing ensemble-averaged quantities.

Rings, Hexagons, Hetals, and Dipolar Moment Sink-Sources: The Fanciful Behavior of Water around Cyclodextrin Complexes

The basket-like geometry of cyclodextrins (CDs), with a cavity able to host hydrophobic groups, makes these molecules well suited for a large number of fundamental and industrial applications. Most of the established CD-based applications rely on trial and error studies, often ignoring key information at the atomic level that could be employed to design new products and to optimize their use. Computational simulations are well suited to fill this gap, especially in the case of CD systems due to their low number of degrees of freedom compared with typical macromolecular systems. Thus, the design and validation of solid and efficient methods to simulate and analyze CD-based systems is key to contribute to this field. The behavior of supramolecular complexes critically depends on the media where they are embedded, so the detailed characterization of the solvent is required to fully understand these systems. In the present work, we use the inclusion complex formed by two α-CDs and one sodium dodecyl sulfate molecule to test eight different parameterizations of the GROMOS and AMBER force fields, including several methods aimed to increase the conformational sampling in computational molecular dynamics simulation trajectories. The system proved to be extremely sensitive to the employed force field, as well as to the presence of a water/air interface. In agreement with previous experiments and in contrast to the results obtained with AMBER, the analysis of the simulations using GROMOS showed a quick adsorption of the complex to the interface as well as an extremely exotic behavior of the water molecules surrounding the structure both in the bulk aqueous solution and at the water surface. The chirality of the CD molecule seems to play an important role in this behavior. All together, these results are expected to be useful to better understand the behavior of CD-based supramolecular complexes such as adsorption or aggregation driving forces, as well as to introduce new methods able to speed up general MD simulations.

Anomalous proton conduction behavior across a nanoporous two-dimensional conjugated aromatic polymer membrane

We investigate aqueous proton penetration behavior across a newly synthesized nanoporous two-dimensional conjugated aromatic polymer (2D-CAP) membrane using extensive ReaxFF reactive molecular dynamics simulations. We found that the proton penetration energy barrier across 2D-CAP is twice as high as that of graphtetrayne, even though 2D-CAP exhibits a larger pore size. Detailed analysis indicates that the anomalous high proton conduction energy barrier of 2D-CAP originates from its unique atomic nanopore structure. The hydrogen atoms at the periphery of the 2D-CAP nanopores can form a stable local hydrogen bond network with water molecules inside or surrounding the nanopores. The mobility of water molecules involved in this local hydrogen bond network will be significantly lowered, and the proton transportation process across the nanopores will thus be impeded. Our results show that the proton penetration behavior across nanoporous 2D materials is influenced not only by the pore size, but also by the decorated atoms or functional groups at the pore edges. Hydrogen atoms at the periphery of nanopores with certain geometry can form a stable local hydrogen bond network with neighboring water molecules, further hampering the proton conductivity.