THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Application of terahertz Time-Domain spectroscopy and chemometrics-based whale optimization algorithm in PDE5 inhibitor detection

Combating the illicit use of PDE5 inhibitor drugs is a focal point in forensic science research. In order to achieve rapid identification of such drugs, this study applies terahertz time-domain spectroscopy combined with chemometrics to establish a fast and accurate detection method for PDE5 inhibitors. The optimal detection method is determined by comparing the spectral performance of three optical parameters, namely absorption coefficient, refractive index, and dielectric constant. Linear discriminant models based on different spectral parameters, whale optimization algorithm optimized extreme learning machine models, and whale optimization algorithm optimized random forest models are established. The effectiveness and performance of principal component analysis and competitive adaptive reweighted sampling algorithm for spectral feature data selection are also investigated. The PDE5 inhibitor identification model based on the competitive adaptive reweighted sampling - whale optimization algorithm - random forest (CARS-WOA-RF) model achieves an accuracy of 98.61%, and the identification model for two concentrations of Sildenafil achieves 100% accuracy. The results demonstrate that terahertz time-domain spectroscopy combined with chemometrics can effectively detect various common types of PDE5 inhibitor drugs and different concentrations.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study

Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.

Biomolecular response to hour-long ultralow field microwave radiation: An effective coarse-grained model simulation

Various electronic devices, which we commonly use, radiate microwaves. Such external perturbation influences the functionality of biomolecules. In an ultralow field, the cumulative response of a molecule is expected only over a time scale of hours. To study the structural dynamics of biomolecules over hours, we adopt a simple methodology for constructing the coarse-grained structure of the protein molecule and solve the Langevin equation under different working potentials. In this approach, each amino acid residue of a biomolecule is mapped onto a number of beads, a few for the backbone, and few for the side chain, depending on the complexity of its chemical structure. We choose the force field in such a way that the dynamics of the protein molecule in the presence of ultralow radiation field of microvolt/nm could be followed over the time frame of 2 h. We apply the model to describe a biomolecule, hen egg white lysozyme, and simulate its structural evolution under ultralow strength electromagnetic radiation. The simulation revealed the finer structural details, like the extent of exposure of bioactive residues and the state of the secondary structures of the molecule, further confirmed from spectroscopic measurements [details are available in Phys. Rev. E 97, 052416 (2018)10.1103/PhysRevE.97.052416 and briefly described here]. Though tested for a specific system, the model is quite general. We believe that it harnesses the potential in studying the structural dynamics of any biopolymer under external perturbation over an extended time scale.

Spectroscopic methods for assessing the molecular origins of macroscopic solution properties of highly concentrated liquid protein solutions

In cases of subcutaneous injection of therapeutic monoclonal antibodies, high protein concentrations (>50 mg/ml) are often required. During the development of these high concentration liquid formulations (HCLF), challenges such as aggregation, gelation, opalescence, phase separation, and high solution viscosities are more prone compared to low concentrated protein formulations. These properties can impair manufacturing processes, as well as protein stability and shelf life. To avoid such unfavourable solution properties, a detailed understanding about the nature of these properties and their driving forces are required. However, the fundamental mechanisms that lead to macroscopic solution properties, as above mentioned, are complex and not fully understood, yet. Established analytical methods for assessing the colloidal stability, i.e. the ability of a native protein to remain dispersed in solution, are restricted to dilute conditions and provide parameters such as the second osmotic virial coefficient, B₂₂, and the diffusion interaction coefficient, k_D. These parameters are routinely applied for qualitative estimations and identifications of proteins with challenging solution behaviours, such as high viscosities and aggregation, although the assays are prepared for low protein concentration conditions, typically between 0.1 and 20 mg/ml ("ideal" solution conditions). Quantitative analysis of samples of high protein concentration is difficult and it is hard to obtain information about the driving forces of such solution properties and corresponding protein-protein self-interactions. An advantage of using specific spectroscopic methods is the potential of directly analysing highly concentrated protein solutions at different solution conditions. This allows for collecting/gaining valuable information about the fundamental mechanisms of solution properties of the high protein concentration regime. In addition, the derived parameters might be more predictive as compared to the parameters originating from assays which are optimized for the low protein concentration range. The provided information includes structural data, molecular dynamics at various timescales and protein-solvent interactions, which can be obtained at molecular resolution. Herein, we provide an overview about spectroscopic techniques for analysing the origins of macroscopic solution behaviours in general, with a specific focus on pharmaceutically relevant high protein concentration and formulation conditions.

Equilibrium Ensembles for Insulin Folding from Bias-Exchange Metadynamics

Earliest events in the aggregation process, such as single molecule reconfiguration, are extremely important and the most difficult to characterize in experiments. To this end, we have used well-tempered bias exchange metadynamics simulations to determine the equilibrium ensembles of an insulin molecule under amyloidogenic conditions of low pH and high temperature. A bin-based clustering method that uses statistics accumulated in bias exchange metadynamics trajectories was employed to construct a detailed thermodynamic and kinetic model of insulin folding. The highest lifetime, lowest free-energy ensemble identified consisted of native conformations adopted by a folded insulin monomer in solution, namely, the R-, the R^f-, and the T-states of insulin. The lowest free-energy structure had a root mean square deviation of only 0.15 nm from native x-ray structure. The second longest-lived metastable state was an unfolded, compact monomer with little similarity to the native structure. We have identified three additional long-lived, metastable states from the bin-based model. We then carried out an exhaustive structural characterization of metastable states on the basis of tertiary contact maps and per-residue accessible surface areas. We have also determined the lowest free-energy path between two longest-lived metastable states and confirm earlier findings of non-two-state folding for insulin through a folding intermediate. The ensemble containing the monomeric intermediate retained 58% of native hydrophobic contacts, however, accompanied by a complete loss of native secondary structure. We have discussed the relative importance of nativelike versus nonnative tertiary contacts for the folding transition. We also provide a simple measure to determine the importance of an individual residue for folding transition. Finally, we have compared and contrasted this intermediate with experimental data obtained in spectroscopic, crystallographic, and calorimetric measurements during early stages of insulin aggregation. We have also determined stability of monomeric insulin by incubation at a very low concentration to isolate protein-protein interaction effects.

Vibrational resonance, allostery, and activation in rhodopsin-like G protein-coupled receptors

G protein-coupled receptors are a large family of membrane proteins activated by a variety of structurally diverse ligands making them highly adaptable signaling molecules. Despite recent advances in the structural biology of this protein family, the mechanism by which ligands induce allosteric changes in protein structure and dynamics for its signaling function remains a mystery. Here, we propose the use of terahertz spectroscopy combined with molecular dynamics simulation and protein evolutionary network modeling to address the mechanism of activation by directly probing the concerted fluctuations of retinal ligand and transmembrane helices in rhodopsin. This approach allows us to examine the role of conformational heterogeneity in the selection and stabilization of specific signaling pathways in the photo-activation of the receptor. We demonstrate that ligand-induced shifts in the conformational equilibrium prompt vibrational resonances in the protein structure that link the dynamics of conserved interactions with fluctuations of the active-state ligand. The connection of vibrational modes creates an allosteric association of coupled fluctuations that forms a coherent signaling pathway from the receptor ligand-binding pocket to the G-protein activation region. Our evolutionary analysis of rhodopsin-like GPCRs suggest that specific allosteric sites play a pivotal role in activating structural fluctuations that allosterically modulate functional signals.

A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site

In multi-resolution simulations, different system components are simultaneously modeled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of ligand binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modeled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity by exploring the reproduction of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function. Proteins 2016; 84:1902-1913. © 2016 Wiley Periodicals, Inc.

Role of Ionic Strength in Staphylococcal Cell Aggregation

Cell aggregation plays a key role in biofilm formation and pathogenesis of Staphylococcus species. Although the molecular basis of aggregation in Staphylococci has already been extensively investigated, the influence of environmental factors, such as ionic strength, remains poorly understood. In this paper, we report a new type of cellular aggregation of Staphylococci that depends solely on ionic strength. Seven strains out of 14, all belonging to staphylococcal species, formed large cell clusters within minutes in buffers of ionic strength ranging from 1.5 to 50 mM, whereas isolates belonging to other Gram-positive species did not display this phenotype. Atomic force microscopy (AFM) with chemically functionalized tips provided direct evidence that ionic strength modulates cell surface adhesive properties through changes in cell surface charge. The optimal ionic strength for aggregation was found to be strain dependent, but in all cases, bacterial aggregates formed at an ionic strength of 1.5-50 mM were rapidly dispersed in a solution of higher ionic strength, indicating a reversibility of the cell aggregation process. These findings suggest that some staphylococcal isolates can respond to ionic strength as an external stimulus to trigger rapid cell aggregation in a way that has not yet been reported.

Using THz Spectroscopy, Evolutionary Network Analysis Methods, and MD Simulation to Map the Evolution of Allosteric Communication Pathways in c-Type Lysozymes

It is now widely accepted that protein function is intimately tied with the navigation of energy landscapes. In this framework, a protein sequence is not described by a distinct structure but rather by an ensemble of conformations. And it is through this ensemble that evolution is able to modify a protein's function by altering its landscape. Hence, the evolution of protein functions involves selective pressures that adjust the sampling of the conformational states. In this work, we focus on elucidating the evolutionary pathway that shaped the function of individual proteins that make-up the mammalian c-type lysozyme subfamily. Using both experimental and computational methods, we map out specific intermolecular interactions that direct the sampling of conformational states and accordingly, also underlie shifts in the landscape that are directly connected with the formation of novel protein functions. By contrasting three representative proteins in the family we identify molecular mechanisms that are associated with the selectivity of enhanced antimicrobial properties and consequently, divergent protein function. Namely, we link the extent of localized fluctuations involving the loop separating helices A and B with shifts in the equilibrium of the ensemble of conformational states that mediate interdomain coupling and concurrently moderate substrate binding affinity. This work reveals unique insights into the molecular level mechanisms that promote the progression of interactions that connect the immune response to infection with the nutritional properties of lactation, while also providing a deeper understanding about how evolving energy landscapes may define present-day protein function.

The glassy state of crambin and the THz time scale protein-solvent fluctuations possibly related to protein function

BACKGROUND

THz experiments have been used to characterize the picosecond time scale fluctuations taking place in the model, globular protein crambin.

RESULTS

Using both hydration and temperature as an experimental parameter, we have identified collective fluctuations (<= 200 cm(-1)) in the protein. Observation of the protein dynamics in the THz spectrum from both below and above the glass transition temperature (Tg) has provided unique insight into the microscopic interactions and modes that permit the solvent to effectively couple to the protein thermal fluctuations.

CONCLUSIONS

Our findings suggest that the solvent dynamics on the picosecond time scale not only contribute to protein flexibility but may also delineate the types of fluctuations that are able to form within the protein structure.