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 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.     

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.     

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.     

Computer modeling of substrate binding to lipases from Rhizomucor miehei, Humicola lanuginosa, and Candida rugosa.

The substrate-binding sites of the triacyl glyceride lipases from Rhizomucor miehei, Humicola lanuginosa, and Candida rugosa were studied by means of computer modeling methods. The space around the active site was mapped by different probes. These calculations suggested 2 separate regions within the binding site. One region showed high affinity for aliphatic groups, whereas the other region was hydrophilic. The aliphatic site should be a binding cavity for fatty acid chains. Water molecules are required for the hydrolysis of the acyl enzyme, but are probably not readily accessible in the hydrophobic interface, in which lipases are acting. Therefore, the hydrophilic site should be important for the hydrolytic activity of the enzyme. Lipases from R. miehei and H. lanuginosa are excellent catalysts for enantioselective resolutions of many secondary alcohols. We used molecular mechanics and dynamics calculations of enzyme-substrate transition-state complexes, which provided information about molecular interactions important for the enantioselectivities of these reactions.     

Computer analysis of protein-protein interaction

An automatic procedure which generates possible modes of protein-protein association is developed and applied to the bovine pancreatic trypsin inhibitor-trypsin complex as a test case. Using a simplified model in which each residue is replaced by one interaction center, all possible modes of interaction between the inhibitor and the active center of the enzyme are generated systematically. The non-bonded interactions between the molecules and the protein surface area buried in the generated interfaces are evaluated and used as criteria for selecting stable complexes. We show that satisfactory estimates of accessible and buried surface areas can be made using the simplified model.The procedure leads to about nine structures having non-bonded interactions and buried surface areas similar to those of the native complex. This suggests that the major contributions to the free energy of dissociation are taken into account by our selection procedure, though complementarity and specificity are not properly represented in the simplified model. However, it makes it possible to scan a much larger number of configurations than would otherwise be feasible, chiefly through elimination of side-chain detail.     

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.     

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.     

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.     

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.     

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.