Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36Linked to original sources
PubMed2004

Investigation of glucose binding sites on insulin.

Possible insulin binding sites for D-glucose have been investigated theoretically by docking and molecular dynamics (MD) simulations. Two different docking programs for small molecules were used; Multiple Copy Simultaneous Search (MCSS) and Solvation Energy for Exhaustive Docking (SEED) programs. The configurations resulting from the MCSS search were evaluated with a scoring function developed to estimate the binding free energy. SEED calculations were performed using various values for the dielectric constant of the solute. It is found that scores emphasizing non-polar interactions gave a preferential binding site in agreement with that inferred from recent fluorescence and NMR NOESY experiments. The calculated binding affinity of -1.4 to -3.5 kcal/mol is within the measured range of -2.0 +/- 0.5 kcal/mol. The validity of the binding site is suggested by the dynamical stability of the bound glucose when examined with MD simulations with explicit solvent. Alternative binding sites were found in the simulations and their relative stabilities were estimated. The motions of the bound glucose during molecular dynamics simulations are correlated with the motions of the insulin side chains that are in contact with it and with larger scale insulin motions. These results raise the question of whether glucose binding to insulin could play a role in its activity. The results establish the complementarity of molecular dynamics simulations and normal mode analyses with the search for binding sites proposed with small molecule docking programs.

Binding Sites↗
PubMed1999

Structure, force, and energy of a double-stranded DNA oligonucleotide under tensile loads.

The end-to-end stretching of a duplex DNA oligonucleotide has been studied using potential of mean force (PMF) calculations based on molecular dynamics (MD) simulations and atomic force microscopy (AFM) experiments. Near quantitative agreement between the calculations and experiments was obtained for both the extension length and forces associated with strand separation. The PMF calculations show that the oligonucleotide extends without a significant energetic barrier from a length shorter than A-DNA to a length 2.4 times the contour length of B-DNA at which the barrier to strand separation is encountered. Calculated forces associated with the barrier are 0.09 +/- 0.03 nN, based on assumptions concerning tip and thermal-activated barrier crossing contributions to the forces. Direct AFM measurements show the oligonucleotide strands separating at 2.6 +/- 0.8 contour lengths with a force of 0.13 +/- 0.05 nN. Analysis of the energies from the MD simulations during extension reveals compensation between increases in the DNA-self energy and decreases in the DNA-solvent interaction energy, allowing for the barrierless extension of DNA beyond the canonical B form. The barrier to strand separation occurs when unfavorable DNA interstrand repulsion cannot be compensated for by favorable DNA-solvent interactions. The present combination of single molecule theoretical and experimental approaches produces a comprehensive picture of the free energy surface of biological macromolecular structural transitions.

DNA↗
PubMed2006

Self-assembling cyclic peptides: molecular dynamics studies of dimers in polar and nonpolar solvents.

The self-assembly of cyclic D,L-alpha-peptides into hollow nanotubes is a crucial mechanistic step in their application as antibacterial and drug-delivery agents. To understand this process, molecular dynamics (MD) simulations were performed on dimers of cyclic peptides formed from cyclo [(-L-Trp-D-N-MeLeu-)4-]2 and cyclo [(-L-Trp-D-Leu-)4-]2 subunits in nonpolar (nonane) and polar (water) solvent. The dimers were observed to be stable only in nonpolar solvent over the full 10 ns length of the MD trajectory. The behavior of the dimers in different solvents is rationalized in terms of the intersubunit hydrogen bonding, hydrogen bonding with the solvent, and planarity of the rings. It is shown that the phi and psi dihedral angles of a single uncapped ring in nonane lie in the beta-sheet region of the Ramachandran plot, and the ring stays in a flat conformation. Steered MD (SMD) simulations based on Jarzynski's equality were performed to obtain the potential of mean force as a function of the distance between the two rings of the capped dimer in nonane. It is also shown that a single peptide subunit prefers to reside close to the nonane/water interface rather than in bulk solvent because of the amphiphilic character of the peptide ring. The present MD results build the foundation for using MD simulations to study the mechanism of the formation of cyclic peptide nanotubes in lipid bilayers.

Dimerization↗
PubMed2006

First-principles theoretical analysis of silyl radical diffusion on silicon surfaces.

We report results from a detailed analysis of the fundamental radical precursor diffusion processes on silicon surfaces and discuss their implications for the surface smoothness of hydrogenated amorphous silicon (a-Si:H) thin films. The analysis is based on a synergistic combination of first-principles density functional theory (DFT) calculations of SiH(3) radical migration on the hydrogen-terminated Si(001)-(2 x 1) surface with molecular-dynamics (MD) simulations of SiH(3) radical precursor migration on surfaces of a-Si:H films. Our DFT calculations yield activation energies for SiH(3) migration that range from 0.18 to 0.89 eV depending on the local electronic environment on the Si(001)-(2 x 1):H surface. In particular, when no substantial surface relaxation (Si-Si bond breaking or formation) accompanies the hopping of the SiH(3) radical the activation barriers are highest, whereas hopping between nearest-neighbor overcoordinated surface Si atoms results in the lowest radical diffusion barrier of 0.18 eV; this low barrier is consistent with the activation barrier for SiH(3) migration through overcoordinated sites on the a-Si:H surface. Specifically, the analysis of the MD simulations of SiH(3) radical migration on a-Si:H surfaces yields an effective diffusion barrier of 0.16 eV, allowing for the rapid migration of the SiH(3) radical prior to its incorporation in surface valleys; rapid migration and subsequent incorporation constitute the two-step mechanism responsible for the smoothness of plasma deposited a-Si:H thin films.

Chemical Phenomena↗
PubMed2006

Collective Langevin dynamics of conformational motions in proteins.

Functionally relevant slow conformational motions of proteins are, at present, in most cases inaccessible to molecular dynamics (MD) simulations. The main reason is that the major part of the computational effort is spend for the accurate description of a huge number of high frequency motions of the protein and the surrounding solvent. The accumulated influence of these fluctuations is crucial for a correct treatment of the conformational dynamics; however, their details can be considered irrelevant for most purposes. To accurately describe long time protein dynamics we here propose a reduced dimension approach, collective Langevin dynamics (CLD), which evolves the dynamics of the system within a small subspace of relevant collective degrees of freedom. The dynamics within the low-dimensional conformational subspace is evolved via a generalized Langevin equation which accounts for memory effects via memory kernels also extracted from short explicit MD simulations. To determine the memory kernel with differing levels of regularization, we propose and evaluate two methods. As a first test, CLD is applied to describe the conformational motion of the peptide neurotensin. A drastic dimension reduction is achieved by considering one single curved conformational coordinate. CLD yielded accurate thermodynamical and dynamical behaviors. In particular, the rate of transitions between two conformational states agreed well with a rate obtained from a 150 ns reference molecular dynamics simulation, despite the fact that the time scale of the transition (approximately 50 ns) was much longer than the 1 ns molecular dynamics simulation from which the memory kernel was extracted.

Biophysics↗
PubMed1999

Structural recognition of an ICAM-1 peptide by its receptor on the surface of T cells: conformational studies of cyclo (1, 12)-Pen-Pro-Arg-Gly-Gly-Ser-Val-Leu-Val-Thr-Gly-Cys-OH.

The purpose of this study is to elucidate the solution conformation of cyclic peptide 1 (cIBR), cyclo (1, 12)-Pen1-Pro2-Arg3-Gly4-Gly5-Ser6-Val7-Leu8-V al9-Thr10-Gly11-Cys12-OH, using NMR, circular dichroism (CD) and molecular dynamics (MD) simulation experiments. cIBR peptide (1), which is derived from the sequence of intercellular adhesion molecule-1 (ICAM-1, CD54), inhibits homotypic T-cell adhesion in vitro. The peptide hinders T-cell adhesion by inhibiting the leukocyte function-associated antigen-1 (LFA-1, CD11a/CD18) interaction with ICAM-1. Furthermore, Molt-3 T cells bind and internalize this peptide via cell surface receptors such as LFA-1. Peptide internalization by the LFA-1 receptor is one possible mechanism of inhibition of T-cell adhesion. The recognition of the peptide by LFA-1 is due to its sequence and conformation; therefore, this study can provide a better understanding for the conformational requirement of peptide-receptor interactions. The solution structure of 1 was determined using NMR, CD and MD simulation in aqueous solution. NMR showed a major and a minor conformer due to the presence of cis/trans isomerization at the X-Pro peptide bond. Because the contribution of the minor conformer is very small, this work is focused only on the major conformer. In solution, the major conformer shows a trans-configuration at the Pen1-Pro2 peptide bond as determined by HMQC NMR. The major conformer shows possible beta-turns at Pro2-Arg3-Gly4-Gly5, Gly5-Ser6-Val7-Leu8, and Val9-Thr10-Gly11-Cys12. The first beta-turn is supported by the ROE connectivities between the NH of Gly4 and the NH of Gly5. The connectivities between the NH of Ser6 and the NH of Val7, followed by the interaction between the amide protons of Val7 and Leu8, support the presence of the second beta-turn. Furthermore, the presence of a beta-turn at Val9-Thr10-Gly11-Cys12 is supported by the NH-NH connectivities between Thr10 and Gly11 and between Gly11 and Cys12. The propensity to form a type I beta-turn structure is also supported by CD spectral analysis. The cIBR peptide (1) shows structural similarity at residues Pro2 to Val7 with the same sequence in the X-ray structure of D1-domain of ICAM-1. The conformation of Pro2 to Val7 in this peptide may be important for its binding selectivity to the LFA-1 receptor.

Amino Acid Sequence↗
PubMed1999

Structural dynamics of ligand diffusion in the protein matrix: A study on a new myoglobin mutant Y(B10) Q(E7) R(E10).

A triple mutant of sperm whale myoglobin (Mb) [Leu(B10) --> Tyr, His(E7) --> Gln, and Thr(E10) --> Arg, called Mb-YQR], investigated by stopped-flow, laser photolysis, crystallography, and molecular dynamics (MD) simulations, proved to be quite unusual. Rebinding of photodissociated NO, O2, and CO from within the protein (in a "geminate" mode) allows us to reach general conclusions about dynamics and cavities in proteins. The 3D structure of oxy Mb-YQR shows that bound O2 makes two H-bonds with Tyr(B10)29 and Gln(E7)64; on deoxygenation, these two residues move toward the space occupied by O2. The bimolecular rate constant for NO binding is the same as for wild-type, but those for CO and O2 binding are reduced 10-fold. While there is no geminate recombination with O2 and CO, geminate rebinding of NO displays an unusually large and very slow component, which is pretty much abolished in the presence of xenon. These results and MD simulations suggest that the ligand migrates in the protein matrix to a major "secondary site," located beneath Tyr(B10)29 and accessible via the motion of Ile(G8)107; this site is different from the "primary site" identified by others who investigated the photolyzed state of wild-type Mb by crystallography. Our hypothesis may rationalize the O2 binding properties of Mb-YQR, and more generally to propose a mechanism of control of ligand binding and dissociation in hemeproteins based on the dynamics of side chains that may (or may not) allow access to and direct temporary sequestration of the dissociated ligand in a docking site within the protein. This interpretation suggests that very fast (picosecond) fluctuations of amino acid side chains may play a crucial role in controlling O2 delivery to tissue at a rate compatible with physiology.

Animals↗
PubMed2005

Computational studies and peptidomimetic design for the human p53-MDM2 complex.

The interaction between human p53 and MDM2 is a key event in controlling cell growth. Many studies have suggested that a p53 mimic would be sufficient to inhibit MDM2 to reduce cell growth in cancerous tissue. In order to design a potent p53 mimic, molecular dynamics (MD) simulations were used to examine the binding interface and the effect of mutating key residues in the human p53-MDM2 complex. The Generalized Born surface area (GBSA) method was used to estimate free energies of binding, and a computational alanine-scanning approach was used to calculate the relative effects in the free energy of binding for key mutations. Our calculations determine the free energy of binding for a model p53-MDM2 complex to be -7.4 kcal/mol, which is in very good agreement with the experimentally determined values (-6.6--8.8 kcal/mol). The alanine-scanning results are in good agreement with experimental data and calculations by other groups. We have used the information from our studies of human p53-MDM2 to design a beta-peptide mimic of p53. MD simulations of the mimic bound to MDM2 estimate a free energy of binding of -8.8 kcal/mol. We have also applied alanine scanning to the mimic-MDM2 complex and reveal which mutations are most likely to alter the binding affinity, possibly giving rise to escape mutants. The mimic was compared to nutlins, a new class of inhibitors that block the formation of the p53-MDM2 complex. There are interesting similarities between the nutlins and our mimic, and the differences point to ways that both inhibitors may be improved. Finally, an additional hydrophobic pocket is noted in the interior of MDM2. It may be possible to design new inhibitors to take advantage of that pocket.

Binding Sites↗
PubMed2005

Insights into the recognition and association of transmembrane alpha-helices. The free energy of alpha-helix dimerization in glycophorin A.

The free energy of alpha-helix dimerization of the transmembrane (TM) region of glycophorin A was estimated from a 125-ns molecular dynamics (MD) simulation in a membrane mimetic. The free energy profile was obtained by allowing the TM helical segments to diffuse reversibly along the reaction pathway. Partition of the potential of mean force into free energy components illuminates the critical steps of alpha-helix recognition and association. At large separations, the TM segments are pushed together by the solvent, allowing initial, but not necessarily native, interhelical interactions to occur. This early recognition stage precedes the formation of native contacts, which is accompanied by a tilt of the helices, characteristic of the dimeric structure. This step is primarily driven by the van der Waals helix-helix interactions. Free energy perturbation calculations of the L75A and I76A point mutations reveal a disruption in helix-helix association due to a loss of favorable dispersion interactions. Additional MD simulations of the native TM dimer and of a single alpha-helix confirm that, prior to association, individual alpha-helices are independently stable, in agreement with the "two-stage" model of integral membrane protein folding.

Amino Acid Sequence↗
PubMed2005

Theoretical study on interaction of hydrogen with single-walled boron nitride nanotubes. II. Collision, storage, and adsorption.

Collision and adsorption of hydrogen with high incident kinetic energies on a single-walled boron nitride (BN) nanotube have been investigated. Molecular-dynamics (MD) simulations indicate that at incident energies below 14 eV hydrogen bounces off the BN nanotube wall. On the other hand, at incident energies between 14 and 22 eV each hydrogen molecule is dissociated at the exterior wall to form two hydrogen atoms, but only one of them goes through the wall. However, at the incident energies between 23 and 26 eV all of the hydrogen atoms dissociated at the exterior wall are found to be capable of going inside the nanotube and then to recombine to form hydrogen molecules inside the nanotube. Consequently, it is determined that hydrogen should have the incident energy >22 eV to go inside the nanotube. On the other hand, we find that the collisions using the incident energies >26 eV could result in damaging the nanotube structures. In addition our MD simulations find that hydrogen atoms dissociated at the wall cannot bind to either boron or nitrogen atoms in the interior wall of the nanotube.

Journal Article↗
PubMed2006

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding.

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.

Amino Acid Sequence↗
PubMed2006

Molecular dynamics study of the inclusion of cholesterol into cyclodextrins.

The interaction of three cyclodextrins (CDs), viz. beta-CD, heptakis (2,6-di-O-methyl)-beta-CD (DM-beta-CD), and 2-hydroxypropyl-beta-CD (HP-beta-CD), with cholesterol was investigated using molecular dynamics (MD) simulations. The free energy along the reaction pathway delineating the inclusion of cholesterol into each CD was computed using the adaptive biasing force method. The association constant and the corresponding association free energy were derived by integrating the potential of mean force (PMF) over a representative ordering parameter. The results show that the free energy profiles possess two local minima corresponding to roughly equally probable binding modes. Among the three CDs, DM-beta-CD exhibits the highest propensity to associate with cholesterol. Ranking for binding cholesterol, viz. DM-beta-CD > HP-beta-CD > beta-CD, agrees nicely with experiment. Partitioning of the PMF into free energy components illuminates that entering of cholesterol into the CD cavity is driven mainly by electrostatic interactions, whereas deeper inclusion results from van der Waals forces and solvation effects. Additional MD simulations were performed to investigate the structural stability of the host-guest complexes near the free energy minima. The present results demonstrate that association of cholesterol and CDs follows two possible binding modes. Although the latter are thermodynamically favorable for all CDs, one of the two inclusion complexes appears to be preferred kinetically in the case of DM-beta-CD.

Cholesterol↗
PubMed1995

Conformation of the circular dumbbell d : structure determination and molecular dynamics.

The circular DNA decamer 5'-d -3' was studied in solution by means of NMR spectroscopy and molecular dynamics in H2O. At a temperature of 269 K, a 50/50 mixture of two dumbbell structures (denoted L2L2 and L2L4) is present. The L2L2 form contains three Watson-Crick C-G base pairs and two two-residue loops is opposite parts of the molecule. On raising the temperature from 269 K to 314 K, the L2L4 conformer becomes increasingly dominant (95% at 314 K). This conformer has a partially disrupted G(anti)-C(syn) closing base pair in the 5'-GTTC-3' loop with only one remaining (solvent-accessible) hydrogen bond between NH alpha of the cytosine dC(1) and O6 of the guanine dG(8). The opposite 5'-CTTG-3' loop remains stable. The two conformers occur in slow equilibrium (rate constant 2-20 s-1). Structure determination of the L2L2 and L2L4 forms was performed with the aid of a full relaxation matrix approach (IRMA) in combination with restrained MD. Torsional information was obtained from coupling constants. Coupling constant analysis (3JHH, 3JHP, 3JCP) gave detailed information about the local geometry around backbone torsion angles beta, gamma, delta, and epsilon, revealing a relatively high flexibility of the 5'-GTTC-3' loop. The values of the coupling constants are virtually temperature-independent. 'Weakly constrained' molecular dynamics in solvent was used to sample the conformational space of the dumbbell. The relaxation matrices from the MD simulation were averaged over to predict dynamic NOE volumes. In order to account for the 1:1 conformational mixture of L2L2 and L2L4 present at 271 K, we also included S2 factors and averaging of the -averaged relaxation matrices. On matrix averaging, the agreement of NOE volumes with experiment improved significantly for protons located in the thermodynamically less stable 5'-GTTC-3' loop. The difference in stability of the 5'-CTTG-3' and 5'-GTTC-3' loops is mainly caused by differences in the number of potential hydrogen bonds in the minor groove and differences in stacking overlap of the base pairs closing the minihairpin loops. The syn conformation for dC(1), favored at high temperature, is stabilized by solvation in the major groove. However, the conformational properties of the dC(1) base, as deduced from R-factor analysis and MD simulations, include a large flexibility about torsion angle chi.

DNA, Circular↗
PubMed2002

Coupling interactions of distal residues enhance dihydrofolate reductase catalysis: mutational effects on hydride transfer rates.

Recently, the participation in catalysis of residues spatially removed from the enzyme's active site has received considerable attention. The influence of the distal Gly-121 residue on the chemical step of hydride transfer in dihydrofolate reductase (DHFR) catalysis had been demonstrated previously [Cameron, C. E., and Benkovic, S. J. (1997) Biochemistry 36, 15792-15800]. In our continuing effort to identify catalytically important residues that are distal from the active site, we used sequence conservation information, kinetic data on site-directed mutants, dynamic motion information from NMR methods, and correlated motions from MD simulations to identify a subset of residues. Among them, the region spanning positions 41-45 is distal to the active site and was chosen as the focus for the mutagenesis and kinetic studies reported here. Specifically, the highly conserved Met-42 was selected for site-directed mutagenesis. While the reaction kinetics for the M42F mutant enzyme did not deviate from wild-type behavior, a 41-fold reduction in the forward hydride transfer rate was found for the M42W mutant. Given the established role of Gly-121 in the hydride transfer process, double mutant enzymes involving positions 42 and 121 were constructed and characterized. These double mutant enzymes generally showed little changes in substrate and cofactor binding but synergistic decreases in forward hydride transfer rates, while the decreases in reverse rates were additive. Along with supporting information from mixed quantum/classical MD simulations [Agarwal, P. K., Billeter, S. R., Rajagopalan, P. T., Benkovic, S. J., and Hammes-Schiffer, S. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 2794-2799], the data suggest that a coupled dynamic motion of these distal residues enhances crossing of the chemical reaction barrier and imply an expanded nonstatic role for the protein fold in catalysis.

Bacterial Proteins↗
PubMed2006

Differential membrane packing of stereoisomers of bis(monoacylglycero)phosphate.

Bis(monoacylglycero)phosphate (BMP) reveals an unusual sn-1,sn-1' stereoconfiguration of glycerophosphate. We synthesized sn-(3-myristoyl-2-hydroxy)glycerol-1-phospho-sn-1'-(3'-myristoyl-2'-hydroxy)glycerol (1,1'-DMBMP) and characterized the thermotropic phase behavior and membrane structure, in comparison with those of the corresponding sn-3:sn-1' stereoisomer (3,1'-DMBMP), by means of differential scanning calorimetry (DSC), small- and wide-angle X-ray scattering (SAXS and WAXS, respectively), pressure-area (pi-A) isotherms, epifluorescence microscopy of monolayers, and molecular dynamics (MD) simulations. In DSC, these lipids exhibited weakly energetic broad peaks with an onset temperature of 9 degrees C for 1,1'-DMBMP and 18 degrees C for 3,1'-DMBMP. In addition, a highly cooperative, strongly energetic transition peak was observed at approximately 40 degrees C for 1,1'-DMBMP and approximately 42 degrees C for 3,1'-DMBMP. These results are supported by the observation that 1,1'-DMBMP exhibited a larger phase transition pressure (pi(c)) than 3,1'-DMBMP. Small- and wide-angle X-ray scattering measurements identified these small and large energetic transitions as a quasi-crystalline (L(c1))-quasi-crystalline with different tilt angle (L(c2)) phase transition and an L(c2)-L(alpha) main phase transition, respectively. X-ray measurements also revealed that these DMBMPs undergo an unbinding at the main phase transition temperature. The MD simulations estimated stronger hydrogen bonding formation in the 3,1'-DMBMP membrane than in 1,1'-DMBMP, supporting the experimental data.

Calorimetry, Differential Scanning↗
PubMed2026

Sanguinarine as a multi-target therapeutic candidate for laryngeal cancer: insights from network pharmacology, molecular dynamics and in vitro validation.

OBJECTIVE: To identify the core targets and elucidate the potential molecular mechanisms of sanguinarine (SA) against laryngeal squamous cell carcinoma (LSCC), and to validate its antitumor effects in vitro. METHODS: Potential targets of SA were predicted using SwissTargetPrediction, TargetNet, and SuperPred and intersected with LSCC-related targets obtained from the GeneCards, OMIM, and DISEASES databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. A protein-protein interaction (PPI) network was constructed using the STRING database (combined score > 0.900), and topological parameters including degree centrality (DC), betweenness centrality (BC), closeness centrality (CC), eigenvector centrality (EC), and local average connectivity (LAC) were calculated in Cytoscape to identify core genes based on median thresholds. Molecular docking and 100-ns molecular dynamics (MD) simulations were conducted for epidermal growth factor receptor (EGFR), Phosphatidylinositide-3-kinase catalytic subunit alpha (PIK3CA), phosphatidylinositol-4,5-biphosphate 3-kinase catalytic subunit β (PIK3CB), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit delta (PIK3CD), and Non-Receptor Tyrosine Kinase (SRC). The effects of SA on LSCC were evaluated using CCK-8, colony formation, Transwell migration, and wound-healing assays in TU177 cells and TU212. RESULTS: A total of 213 common targets were identified, which were significantly enriched in PI3K-Akt signaling, EGFR tyrosine kinase inhibitor resistance, and adhesion- and migration-related pathways. The PPI network comprised 259 nodes and 259 edges, from which five core genes-PIK3CA, PIK3CB, PIK3CD, EGFR, and SRC-were identified. Molecular docking revealed strong binding affinities between SA and the PI3K family proteins (- 9.79 to - 10.96 kcal/mol), as well as EGFR (- 8.58 kcal/mol) and SRC (- 6.77 kcal/mol). MD simulations indicated greater stability of SA complexes with EGFR and PI3K family members compared with SRC. In vitro assays demonstrated that SA significantly inhibited TU177 cell and TU212 cell proliferation, colony formation, and migration. CONCLUSION: SA may exert anti-laryngeal cancer effects through synergistic multi-target inhibition centered on the EGFR/SRC/PI3K signaling axis, highlighting its potential as a promising therapeutic candidate for LSCC.

Humans↗
PubMed2001

A unified picture of protein hydration: prediction of hydrodynamic properties from known structures.

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.

Nuclear Magnetic Resonance, Biomolecular↗
PubMed1999

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5.

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.

Animals↗