Computer modeling of human movements]

A dynamic model for studying man's movements is proposed. Lagrange equations of the second order are used. Differential equations of the model are presented in the matrix form, and all the coefficients involved are calculated from recurrent formulae. The dynamic model described is easily algorythmized. Differentiating operations can be thus avoided which are realized on electron computers with difficulities.     

Computer modeling of the membrane interaction of FYVE domains

FYVE domains are membrane targeting domains that are found in proteins involved in endosomal trafficking and signal transduction pathways. Most FYVE domains bind specifically to phosphatidylinositol 3-phosphate (PI(3)P), a lipid that resides mainly in endosomal membranes. Though the specific interactions between FYVE domains and the headgroup of PI(3)P have been well characterized, principally through structural studies, the available experimental structures suggest several different models for FYVE/membrane association. Thus, the manner in which FYVE domains adsorb to the membrane surface remains to be elucidated. Towards this end, recent experiments have shown that FYVE domains bind PI(3)P in the context of phospholipid bilayers and that hydrophobic residues on a conserved loop are able to penetrate the membrane interface in a PI(3)P-dependent manner.Here, the finite difference Poisson-Boltzmann (FDPB) method has been used to calculate the energetic interactions of FYVE domains with phospholipid membranes. Based on the computational analysis, it is found that (1) recruitment to membranes is facilitated by non-specific electrostatic interactions that occur between basic residues on the domains and acidic phospholipids in the membrane, (2) the energetic analysis can quantitatively differentiate among the modes of membrane association proposed by the experimentally determined structures, (3) FDPB calculations predict energetically feasible models for the membrane-associated states of FYVE domains, (4) these models are consistent with the observation that conserved hydrophobic residues insert into the membrane interface, and (5) the calculations provide a molecular model for the hydrophobic partitioning: binding of PI(3)P significantly neutralizes positive potential in the region of the hydrophobic residues, which acts as an "electrostatic switch" by reducing the energetic barrier for membrane penetration. Finally, the computational results are extended to FYVE domains of unknown structure through the construction of high quality homology models for human FYVE sequences.     

Computer modeling of the human systemic arterial tree

A model of the human systemic arterial tree has been devised, based on a lumped-parameter-circuit approximate form. This model has been set up and studied on an analog computer. A feature of this simulation is the division of the arterial system into sections whose lengths are inversely proportional (approximately) to their cross-sectional area-or what is termed ‘equal-volume’ modeling.

Great care was exercised in the determination of the model parameters, using expressions for these parameters from a recent paper by Rideout and Dick on fluid flow in distensible tubes, with numerical values based on measurements reported in the medical literature.

The simulated pressure and flow waveforms obtained with the model compare favorably with data recorded from the normal adult human, and exhibit such well-known features as distal delay and peaking of pressure pulses. The aortic input impedance vs. frequency curve checks well against measurements on the human. The model also provides a simple means for determination of cardiac output, cardiac work and cardiac power under various assumed conditions such as variation of heart rate.     

Homology modeling and substrate binding study of human CYP4A11 enzyme

Although both bacterial CYP102 (P450BM3) and mammalian CYP4A isozymes share a common function as fatty acid hydroxylases, distinctly different preferred sites of oxidation are observed with the CYP102 performing the usual non-terminal hydroxylation or epoxidation and the CYP4A enzymes performing the unusual and enigmatic terminal hydroxylation. The origin of this unique product specificity in human CYP4A11 has been explored in this work, focusing on possible differences in the binding site architecture of the two isozymes as the cause. To this end, 3D model structures of the human CYP4A11 enzyme were built and compared to the X-ray structure of CYP102. The substrate-binding channel identified in CYP4A11 was found to have a much more sterically restricted active site than that in CYP102 that could cause limited access of long-chain fatty acid to the ferryl oxygen leading to the preferred omega-hydroxylation. Results of docking of a common substrate, lauric acid, into the binding site of both CYP4A11 and CYP102 and molecular dynamics simulations provided additional support for this hypothesis. Specifically, in the CYP4A11-lauric acid simulations, the omega hydrogens were closest to the ferryl oxygen most of the time. By contrast, in the CYP102-lauric acid complex, the substrate could penetrate further into the active site providing access of the non-terminal (omega-1, omega-2) positions to the ferryl oxygen. These results, taken together, have elucidated the origin of the unusual product specificity of CYP4A11 and illustrated the central role of binding site architecture in subtle modulation of function.     

Homology modeling and substrate binding study of human CYP2C9 enzyme

The CYP2C subfamily of human liver P450 isozymes is of major importance in drug metabolism. The most abundant 2C isozyme, CYP2C9, regioselectively hydroxylates a wide variety of substrates. A major obstacle to understanding this specificity in human CYP2C9 is the absence of a 3D structure. A 3D model of CYP2C9 was built, assessed, and used to characterize explicit enzyme-substrate complexes using methods previously developed in our laboratory. The 3D model was assessed by determining its stability to unconstrained molecular dynamics and by comparison of specific properties with those of known protein structures. The CYP2C9 model was then used to characterize explicit enzyme complexes with three structurally and chemically diverse substrates: (S)-naproxen, phenytoin, and progesterone. Each substrate was found to bind to the enzyme with a favorable interaction energy and to remain in the binding site during unconstrained molecular dynamics. Moreover, the mode of binding of each substrate led to calculated preferred hydroxylation sites consistent with experiment. Binding-site residues identified for the models included Arg 105 and Arg97 as key cationic residues, as well as Asn 202, Asp 293, Pro 101, Leu 102, Gly 296, and Phe 476. Site-specific mutations are proposed for further integrated computational and experimental study.     

Computer modeling of ventricular fibrillation

Electrical activity of a heart in ventricular fibrillation was modeled as a sum of independent pulse streams with various amplitude-frequency and phase characteristics. Results of computer experiments were compared with those of real physiological experiments on rabbits. Identification of the model was carried out by means of the least-squares procedure. The offered technique allows a computer model investigation of internal structure of irregularities of ventricular fibrillation.     

Computer modeling of ventricular fibrillation

Electrical activity of a heart in ventricular fibrillation was modeled as a sum of independent pulse streams with various amplitude-frequency and phase characteristics. Results of computer experiments were compared with those of real physiological experiments on rabbits. Identification of the model was carried out by means of the least-squares procedure. The offered technique allows a computer model investigation of internal structure of irregularities of ventricular fibrillation.     

Three-dimensional modeling of thrombin-fibrinogen interaction

Three-dimensional models of thrombin complexed with large fragments of the fibrinogen Aalpha and Bbeta chains are presented. The models are consistent with the results of recent mutagenesis studies of thrombin and with the information available on naturally occurring fibrinogen mutants. Thrombin recognizes fibrinogen with an extended binding surface, key elements of which are Tyr(76) in exosite I, located about 20 A away from the active site, and the aryl binding site located in close proximity to the catalytic triad. A highly conserved aromatic-Pro-aromatic triplet motif is identified in the primed site region of fibrinogen and other natural substrates of thrombin. The role of this triplet, based on the three-dimensional models, is to correctly orient the substrate for optimal bridge binding to exosite I and the active site. The three-dimensional models suggest a possible pattern of recognition by thrombin that applies generally to other natural substrates.     

Functional modeling of neural-glial interaction

We propose a generalized mathematical model for a small neural-glial ensemble. The model incorporates subunits of the tripartite synapse that includes a presynaptic neuron, the synaptic terminal itself, a postsynaptic neuron, and a glial cell. The glial cell is assumed to be activated via two different pathways: (i) the fast increase of intercellular [K(+)] produced by the spiking activity of the postsynaptic neuron, and (ii) the slow production of a mediator triggered by the synaptic activity. Our model predicts the long-term potentiation of the postsynaptic neuron as well as various [Ca(2+)] transients in response to the activation of different pathways.     

Homology modeling of Mycoplasma pneumoniae enolase and its molecular interaction with human plasminogen

Alpha (α)-enolase (e), a glycolytic enzyme, has an alternative role as a surface receptor of several bacteria mediating plasminogen (pg) binding. It is also recognized as a virulence factor of some pathogenic bacteria facilitating plasminogen activation and host cell invasion. A mycoplasmal α-enolase is also a plasminogen binding protein. Molecular interactions of enolase from Mycoplasma pneumoniae with host plasminogen would be useful for exploring the pathogen-host interaction. In an attempt to identify plasminogen binding sites of M. pneumoniae enolase, homology modeling and docking studies were conducted to obtain modeled structures of the M. pneumoniae enolase-plasminogen complex. The refined model was validated further by standard methods. Molecular docking revealed hydrogen bonding of eLys70-pgTyr50, eAsn165-pgThr66, eAla168-pgGlu21, eAsp17-pgLys70, and eAsn213-pgPro68/pgAsn69. Substantial decreases in accessible surface area (ASA) were observed and in concurrence with hydrogen bond pattern. These findings provide a detailed prediction of key residues that interact at the protein-protein interface. Our theoretical prediction is consistent with known biochemical data. The predicted interaction complex can be of great assistance in understanding structural insights, which is necessary to pathogen and host-component interaction. The ability of M. pneumoniae enolase to bind plasminogen may be indicative of an important role in invasion of this pathogen to host.     

Computer modeling of the effective search

A computer model of vector-Brownian processes is proposed, which can be applied to some biophysical problems including the problem of search efficiency. The results obtained suggest that the changes in the ratio of two components (the vector and the Brownian ones) in any process essentially improve the efficiency of search.     

Computer modeling of muscle phosphofructokinase kinetics

The kinetics of the phosphofructokinase reaction were studied by computer modeling. A general random order, two-state allosteric model, of which the Monod--Wyman--Changeux model is a limiting case, was found to most accurately reproduce the experimental observations of Pettigrew & Frieden (1979 a,b). A simplified model with Hill coefficients was found to fit almost as well. In these models substrates bind preferentially to and stabilize the enzyme in the R state, and ATPH3-, the inhibitory species, binds preferentially to and stabilizes the enzyme in the T state. Enzymatic activity is regulated by conversion from the R to the T state, which is effected by protonation, especially of the uncomplexed enzyme, but the experimental data are inadequate for accurate estimation of the pKa of the enzyme. Random order binding of substrates is an important cause of sigmoidal kinetics. Additional experiments that would aid in the discrimination among rival models are described.     

Computer modeling of small heat-shock metalloprotease of the human malaria parasite Plasmodium vivax.

We present here computer generated model of N-terminal fragment, amino acids (aa) 36-245, of a Plasmodium vivax heat shock metalloprotease called PVHSP28, whose gene was cloned and characterised earlier. The fragment showed homology with HSPs from many organisms, including Escherichia coli and Haemophilus influenzae. PVHSP28 had the signature sequence 'HEXXH' and 'EXXXD' of Zinc metalloproteases. Being the first malarial HSP possessing metalloprotease activity, PVHSP28 is an ideal target for the design of new anti-malarial drugs. However, except for a small region (aa 62-132) which had 24.6% sequence similarity with 1TAQ (a DNA polymerase), it did not show sequence similarity with any published structures in protein data bank. Hence it could not be modelled using any automated modeling programs. We modelled 36-245 aa of PVHSP28 using predicted secondary structure as well as experimentally determined and predicted properties of the protein on the basis of its amino acid sequence, using various Internet tools and in-house package MODEL. The model was energy minimised using Sander's module of AMBER 5.0, working on a Silicon Graphics machine, with all atom force field.     

In vitro modeling of cell-scaffold interaction

A simple method of controlling the efficiency of surface ligand-cell receptor interaction has been developed in the course of modeling the specific adhesion of cells on a support with their subsequent proliferation and bone tissue formation, using affinity chromatography on macroporous monolithic sorbents. The biospecific peptide GRGDSP played the role of an active ligand on the support, whereas cells were simulated by polymeric (polystyrene) microparticles with the peptide EDYPVDIYYLMDLSYSMKDD immobilized on their surface. The latter peptide is part of the active site of the integrin molecule responsible for binding the RGD sequence. Thus, the monolithic ultrashort column (CIM® disk) represented a simplified model of the support (structural scaffold) possessing biospecific properties. The parameters of the interaction of affinity partners were quantitatively estimated by frontal analysis involving the construction of adsorption isotherms, followed by their linearization and mathematical processing. The data obtained indicate a high specificity of biological pairing, which is supported by the results of cell culture experiments.     

Molecular modeling of human BAD and its interaction with PKAc or PP1c

To build up the structure of human BAD (Bcl-2 antagonist of cell death), subsequently combined with PKAc or PP1c (protein phosphatase 1), to investigate the interaction relationship between BAD and its kinase/PTPese at the molecular level. Additionally, it is concerned with the search for all optimal positions and orientations of a set of amino acid residues of BAD, while its binding sites include N-termini (Glu19, Ala27, and Ser34-Lys35), BH3-located helical domain (Arg98-Lys126), and C-termini (Trp154-Ser163 and Ser167-Gln168). The related sites of PKAc are mainly assembled in C-terminal α/β-domain of PKAc, which comprises the KTL motif (47-49), Glu203 residue, a helical region (Asp241-Arg256), and the span from 328 to 333; while the interaction sites with BAD converge at C-terminal β-domain of PP1c, which includes the DEK motif (166-168), the stretch from 179 to 197 including a helix (Glu184-Arg188), Glu230-Asp242 segment containing Val232-His237 helix, and Glu287-Leu289 loop. In conclusion, analysis of the complex between BAD and PKAc or PP1c provides a novel viewpoint on the structural origins of molecular recognition. And the complex models suggest that BH3 domain of BAD interact with PKAc or PP1c by electrostatic, van der Waals contacts, hydrogen bond and salt bridge. This is helpful for our development and research of some new drugs, especially mimetic BH3 peptides and inspires scientists with BAD complex and molecular mechanism of its integrating glycolysis and apoptosis.     

Computer modeling of parametric structures of water

The algorithm for the arrangement of hydrogen atoms in twist-hexacycle parametric structures of bound water is developed. The calculation of energetic properties is carried out using the TIP3P and Poltev-Malenkov potentials. Optimization of energy for these structures is fulfilled.     

Molecular interaction of human serum albumin with paracetamol: Spectroscopic and molecular modeling studies

The interaction between paracetamol and human serum albumin (HSA) under physiological conditions has been investigated by fluorescence, circular dichroism (CD) and docking. Fluorescence data revealed that the fluorescence quenching of HSA by paracetamol was the result of the formed complex of HSA–paracetamol, and the binding constant (K_a) and binding number obtained is 1.3 × 10⁴ at 298 K and 2, respectively for the primary binding site. Circular dichorism spectra showed the induced conformational changes in HSA by the binding of paracetamol. Moreover, protein–ligand docking study indicated that paracetamols (two paracetamols bind to HSA) bind to residues located in the subdomain IIIA.     

Computer modeling 16 S ribosomal RNA

A three-dimensional structure for 16 S RNA has been produced with a computer protocol that is not dependent on human intervention. This protocol improves upon traditional modeling techniques by using distance geometry to fold the molecule in an objective and reproducible fashion. The method is based on the secondary structure of RNA and treats the molecule as a set of double-stranded helices that are linked by flexible single-strands of variable length. Data derived from chemical cross-linking studies of 16 S RNA and tertiary phylogenetic relationships provide the constraints used to fold the molecule into a compact three-dimensional form. Possibly subjective evaluation of the input data are transformed into verifiable quantitative parameters. Relationships based on general locations within the 30 S subunit or on protein-RNA interactions have been specifically excluded. The resolution of the model exceeds that of electron micrographs and approaches that obtained in preliminary X-ray crystal structures. The model size of 245 x 190 x 140 A is compatible with that of the 30 S subunit as determined by electron microscopy. The volume of the model is 1.87 x 10(6) A which is similar to that of the small subunit in a preliminary X-ray crystal structure. The radius of gyration of the model structure of 76 A is intermediate to that seen for partially denatured and fully folded 16 S RNA. Computer graphics are used to display the results in a manner that maximizes the opportunities for human visual interpretation of the models. A format for displaying the structures has been developed that will make it possible for researchers who have not devoted themselves to ribosomal modeling to comprehend and make use of the information that the models embody. On this basis the computer-generated models are compared with models developed by other researchers and with structural data not included in the folding parameter data set.